Kurt A. Kistler

Three-state conical intersections in cytosine and pyrimidinone bases

Three-state conical intersections have been located and characterized for cytosine and its analog 5-methyl-2-pyrimidinone using multireference configuration-interaction ab initio methods. The potential energy surfaces for each base contain three different three-state intersections: two different S(0)-S(1)-S(2) intersections (gs/pi pi(*)/n(N)pi(*) and gs/pi pi(*)/n(O)pi(*)) and an S(1)-S(2)-S(3) intersection (pi pi(*)/n(N)pi(*)/n(O)pi(*)). Two-state seam paths from these intersections are shown to be connected to previously reported two-state conical intersections. Nonadiabatic coupling terms have been calculated, and the effects of the proximal third state on these quantities are detailed. In particular, it is shown that when one of these loops incorporates more than one seam point, there is a profound and predictable effect on the phase of the nonadiabatic coupling terms, and as such provides a diagnostic for the presence and location of additional seams. In addition, it is shown that each of the three three-state conical intersections located on cytosine and 5-methyl-2-pyrimidinone is qualitatively similar between the two bases in terms of energies and character, implying that, like with the stationary points and two-state conical intersections previously reported for these two bases, there is an underlying pattern of energy surfaces for 2-pyrimidinone bases, in general, and this pattern also includes three-state conical intersections.

Absorption, circular dichroism, and photoluminescence in perylene diimide bichromophores: polarization-dependent H- and J-aggregate behavior.

Using a single-mode Holstein Hamiltonian with through-space excitonic couplings evaluated quantum mechanically, the absorption, circular dichroism, and photoluminescence spectral line shapes of a chiral perylene diimide dimer complex were accurately reproduced. In general, a dimer consisting of two chromophores related through a C(2) rotation is neither a J- nor an H-aggregate because oscillator strength is divided between the top and bottom of the exciton band. The division gives rise to the two Davydov components per vibronic band in the absorption spectrum. Nevertheless, it is shown that the vibronic structure of the absorption component polarized in the same direction as the lower (upper) Davydov component is identical to what one would obtain from an ideal J- (H-) aggregate. Emission generally contains both polarization components, but the component polarized in the same direction as the lower (upper) Davydov component behaves similarly to the emission from an ideal J- (H-) aggregate. The basic photophysical behavior also applies to molecular crystals containing two molecules per unit cell in which the interactions between inequivalent molecules dominate over interactions between equivalent molecules.

Solvatochromic Shifts of Uracil and Cytosine Using a Combined Multireference Configuration Interaction/Molecular Dynamics Approach and the Fragment Molecular Orbital Method†

A recently developed combined quantum mechanics/molecular mechanics (QM/MM) approach has been applied to the calculation of solvatochromic shifts of the excited states of the pyrimidine nucleobases uracil and cytosine in aqueous solution. In this procedure the quantum mechanical solute is described using a multireference configuration interaction method while molecular dynamics simulations are used to obtain the structure of the solvent around the solute. The fragment molecular orbital multiconfiguration self-consistent field (FMO-MCSCF) method of Fedorov and Kitaura is also used and compared with the QM/MM results. The two methods give similar results. The solvatochromic shifts in uracil are found to be +0.41 (+0.44) eV for the S(1) excited state and -0.05 (-0.19) eV for the S(2) state at the QM/MM (FMO-MCSCF) level. Solvatochromic shifts in cytosine are calculated to be +0.25 (+0.19), +0.56 (+0.62), and +0.83 (+0.83) eV for the S(1), S(2), and S(3) states, respectively, at the QM/MM (FMO-MCSCF) level.

Excited-State Energies and Electronic Couplings of DNA Base Dimers

The singlet excited electronic states of two pi-stacked thymine molecules and their splittings due to electronic coupling have been investigated with a variety of computational methods. Focus has been given on the effect of intermolecular distance on these energies and couplings. Single-reference methods, CIS, CIS(2), EOM-CCSD, TDDFT, and the multireference method CASSCF, have been used, and their performance has been compared. It is found that the excited-state energies are very sensitive to the applied method but the couplings are not as sensitive. Inclusion of diffuse functions in the basis set also affects the excitation energies significantly but not the couplings. TDDFT is inadequate in describing the states and their coupling, while CIS(2) gives results very similar to EOM-CCSD. Excited states of cytosine and adenine pi-stacked dimers were also obtained and compared with those of thymine dimers to gain a more general picture of excited states in pi-stacked DNA base dimers. The coupling is very sensitive to the relative position and orientation of the bases, indicating great variation in the degree of delocalization of the excited states between stacked bases in natural DNA as it fluctuates.

Distinguishing Between Relaxation Pathways by Combining Dissociative Ionization Pump Probe Spectroscopy and Ab Initio Calculations: A Case Study of Cytosine

We present a general method for tracking molecular relaxation along different pathways from an excited state down to the ground state. We follow the excited state dynamics of cytosine pumped near the S(0)-S(1) resonance using ultrafast laser pulses in the deep ultraviolet and probed with strong field near infrared pulses which ionize and dissociate the molecules. The fragment ions are detected via time of flight mass spectroscopy as a function of pump probe delay and probe pulse intensity. Our measurements reveal that different molecular fragments show different timescales, indicating that there are multiple relaxation pathways down to the ground state. We interpret our measurements with the help of ab initio electronic structure calculations of both the neutral molecule and the molecular cation for different conformations en route to relaxation back down to the ground state. Our measurements and calculations show passage through two seams of conical intersections between ground and excited states and demonstrate the ability of dissociative ionization pump probe measurements in conjunction with ab initio electronic structure calculations to track molecular relaxation through multiple pathways.

A benchmark of excitonic couplings derived from atomic transition charges.

In this report we benchmark Coulombic excitonic couplings between various pairs of chromophores calculated using transition charges localized on the atoms of each monomer chromophore, as derived from a Mulliken population analysis of the monomeric transition densities. The systems studied are dimers of 1-methylthymine, 1-methylcytosine, 2-amino-9-methylpurine, all-trans-1,3,5-hexatriene, all-trans-1,3,5,7-octatetraene, trans-stilbene, naphthalene, perylenediimide, and dithia-anthracenophane. Transition densities are taken from different single-reference electronic structure excited state methods: time-dependent density functional theory (TDDFT), configuration-interaction singles (CIS), and semiempirical methods based on intermediate neglect of differential overlap. Comparisons of these results with full ab initio calculations of the electronic couplings using a supersystem are made, as are comparisons with experimental data. Results show that the transition charges do a good job of reproducing the supersystem couplings for dimers with moderate to long-range interchromophore separation. It is also found that CIS supermolecular couplings tend to overestimate the couplings, and often the transition charges approach may be better, due to fortuitous cancellation of errors.

The Fluorescence Mechanism of 5-Methyl-2-Pyrimidinone: An Ab Initio Study of a Fluorescent Pyrimidine Analog†

The photophysically important potential energy surfaces of the fluorescent pyrimidine analog 5‐methyl‐2‐pyrimidinone have been explored using multireference configuration‐interaction ab initio methods at three levels of dynamical correlation, all of which support a fluorescence mechanism. At vertical excitation S1 (dark, nNπ*) and S2 (bright, ππ*) are almost degenerate at 4.4 eV, with S3 (dark, nOπ*) at 5.1 eV. The excited system can follow the S1–S2 seam of conical intersections, accessible from the Franck–Condon region, to its minimum and then evolve from this conical intersection on the S1 (ππ*) surface to a global minimum. At lower levels of correlation, the S1 surface shows two minima separated by a barrier of up to 0.18 eV. The secondary minimum found at the lower levels of correlation becomes the global minimum with higher correlation. The S1 population at this minimum can be trapped from accessing the lowest energy S0–S1 (ππ*/gs) conical intersection by an energy gap at least 0.3–0.4 eV higher than the S1 minimum. The calculated emission energy from this minimum is 2.80 eV. Gradient pathways connecting important S1 geometries are presented, as well as other excited state conical intersections.

Change in Electronic Structure upon Optical Excitation of 8-Vinyladenosine: An Experimental and Theoretical Study

8-Vinyladenosine (8VA) is an adenosine analog, like 2-aminopurine (2AP), that has a red-shifted absorption and high fluorescence quantum yield. When introduced into double-stranded DNA (dsDNA), its base-pairing and base-stacking properties are similar to those of adenine. Of particular interest, the fluorescence quantum yield of 8VA is sensitive to base stacking, making it a very useful real-time probe of DNA structure. The fundamental photophysics underlying this fluorescence quenching by base stacking is not well understood, and thus exploring the excited state electronic structure of the analog is warranted. In this study, we report on changes in the electronic structure of 8VA upon optical excitation. Stark spectroscopy was performed on 8VA monomer in frozen ethanol glass at 77 K to obtain the direction and degree of charge redistribution in the form of the difference dipole moment, Deltamu(01) = 4.7 +/- 0.3 D, and difference static polarizability, tr(Delta(alpha)01) = 21 +/- 11 A(3), for the S(0)-->S(1) transition. In addition, solvatochromism experiments were performed on 8VA in various solvents and analyzed using Bakhshievs model. High level ab initio methods were employed to calculate transition energies, oscillator strengths, and dipole moments of the ground and excited states of 8VA. The direction of Deltamu(01) was assigned in the molecular frame for the lowest optically accessible state. Our study shows that the angle between ground and excited state dipole moment plays a critical role in understanding the change in electronic structure upon optical excitation. Compared to 2AP, 8VA has a larger difference dipole moment which, with twice the extinction coefficient, suggests that 8VA is superior as a two-photon probe for microscopy studies. To this end, we have measured the ratio of the two-photon fluorescence yields of the two analogs by excitation at the respective monomer absorption maxima. We show that 8VA is indeed a significantly brighter two-photon fluorophore, based on our experimental and computational results.

Photophysical pathways of cytosine in aqueous solution

The effects of aqueous solvation on the photophysical pathways involving the S(1) excited state in cytosine have been studied with a mean-field QM/MM approach. Two main pathways with small barriers were found previously in isolated cytosine, using multireference configuration interaction (MRCI) methods, that facilitate radiationless decay to the ground state. These pathways are examined in solvated cytosine using a mean-field QM/MM combined with MRCI, and it is found that barriers in each direction increase moderately. The barriers in the presence of the solvent are 0.23 eV and 0.31 eV for the two different pathways compared to 0.15 eV and 0.14 eV in the gas phase, indicating that the aqueous environment does not make one of the two directions much more preferable.

2-Aminopurine excited state electronic structure measured by stark spectroscopy.

2-Aminopurine (2AP) is an adenine analogue that has a high fluorescence quantum yield. Its fluorescence yield decreases significantly when the base is incorporated into DNA, making it a very useful real-time probe of DNA structure. However, the basic mechanism underlying 2AP fluorescence quenching by base stacking is not well understood. A critical element in approaching this problem is obtaining an understanding of the electronic structure of the excited state. We have explored the excited state properties of 2AP and 2-amino,9-methylpurine (2A9MP) in frozen solutions using Stark spectroscopy. The experimental data were correlated with high level ab initio (MRCI) calculations of the dipole moments, mu0 and mu1, of the ground and excited states. The magnitude and direction of the dipole moment change, Deltamu01 = mu1 - mu0, of the lowest energy optically allowed transition was determined. While other studies have reported on the magnitude of the dipole moment change, we believe that this is the first report of the direction of Deltamu, a quantity that will be of great value in interpreting absorption spectral changes of the 2AP chromophore. Polarizability changes due to the transition were also obtained.