Revealing Deactivation Pathways Hidden in Time-Resolved Photoelectron Spectra.

Time-resolved photoelectron spectroscopy is commonly employed with the intention to monitor electronic excited-state dynamics occurring in a neutral molecule. With the help of theory, we show that when excited-state processes occur on similar time scales the different relaxation pathways are completely obscured in the total photoionization signal recorded in the experiment. Using non-adiabatic molecular dynamics and Dyson norms, we calculate the photoionization signal of cytosine and disentangle the transient contributions originating from the different deactivation pathways of its tautomers. In the simulations, the total signal from the relevant keto and enol tautomers can be decomposed into contributions either from the neutral electronic state populations or from the distinct mechanistic pathways across the multiple potential surfaces. The lifetimes corresponding to these contributions cannot be extracted from the experiment, thereby illustrating that new experimental setups are necessary to unravel the intricate non-adiabatic pathways occurring in polyatomic molecules after irradiation by light.

Strong Field Molecular Ionization in the Impulsive Limit: Freezing Vibrations with Short Pulses.

We study strong-field molecular ionization as a function of pulse duration. Experimental measurements of the photoelectron yield for a number of molecules reveal competition between different ionization continua (cationic states) which depends strongly on pulse duration. Surprisingly, in the limit of short pulse duration, we find that a single ionic continuum dominates the yield, whereas multiple continua are produced for longer pulses. Using calculations which take vibrational dynamics into account, we interpret our results in terms of nuclear motion and nonadiabatic dynamics during the ionization process.

Photoelectron spectra of 2-thiouracil, 4-thiouracil, and 2,4-dithiouracil.

Ground- and excited-state UV photoelectron spectra of thiouracils (2-thiouracil, 4-thiouracil, and 2,4-dithiouracil) have been simulated using multireference configuration interaction calculations and Dyson norms as a measure for the photoionization intensity. Except for a constant shift, the calculated spectrum of 2-thiouracil agrees very well with experiment, while no experimental spectra are available for the two other compounds. For all three molecules, the photoelectron spectra show distinct bands due to ionization of the sulphur and oxygen lone pairs and the pyrimidine π system. The excited-state photoelectron spectra of 2-thiouracil show bands at much lower energies than in the ground state spectrum, allowing to monitor the excited-state population in time-resolved UV photoelectron spectroscopy experiments. However, the results also reveal that single-photon ionization probe schemes alone will not allow monitoring all photodynamic processes existing in 2-thiouracil. Especially, due to overlapping bands of singlet and triplet states the clear observation of intersystem crossing will be hampered.

Efficient and Flexible Computation of Many-Electron Wave Function Overlaps.

A new algorithm for the computation of the overlap between many-electron wave functions is described. This algorithm allows for the extensive use of recurring intermediates and thus provides high computational efficiency. Because of the general formalism employed, overlaps can be computed for varying wave function types, molecular orbitals, basis sets, and molecular geometries. This paves the way for efficiently computing nonadiabatic interaction terms for dynamics simulations. In addition, other application areas can be envisaged, such as the comparison of wave functions constructed at different levels of theory. Aside from explaining the algorithm and evaluating the performance, a detailed analysis of the numerical stability of wave function overlaps is carried out, and strategies for overcoming potential severe pitfalls due to displaced atoms and truncated wave functions are presented.

Nonadiabatic photodynamics of a retinal model in polar and nonpolar environment.

The nonadiabatic photodynamics of the all-trans-2,4-pentadiene-iminium cation (protonated Schiff base 3, PSB3) and the all-trans-3-methyl-2,4-pentadiene-iminium cation (MePSB3) were investigated in the gas phase and in polar (aqueous) and nonpolar (n-hexane) solutions by means of surface hopping using a multireference configuration-interaction (MRCI) quantum mechanical/molecular mechanics (QM/MM) level. Spectra, lifetimes for radiationless deactivation to the ground state, and structural and electronic parameters are compared. A strong influence of the polar solvent on the location of the crossing seam, in particular in the bond length alternation (BLA) coordinate, is found. Additionally, inclusion of the polar solvent changes the orientation of the intersection cone from sloped in the gas phase to peaked, thus enhancing considerably its efficiency for deactivation of the molecular system to the ground state. These factors cause, especially for MePSB3, a substantial decrease in the lifetime of the excited state despite the steric inhibition by the solvent.

Quantum dynamics of ultrafast charge transfer at an oligothiophene-fullerene heterojunction.

Following up on our recent study of ultrafast charge separation at oligothiophene-fullerene interfaces [H. Tamura, I. Burghardt, and M. Tsukada, J. Phys. Chem. C 115, 10205 (2011)], we present here a detailed quantum dynamical perspective on the charge transfer process. To this end, electron-phonon coupling is included non-perturbatively, by an explicit quantum dynamical treatment using the multi-configuration time-dependent Hartree (MCTDH) method. Based upon a distribution of electron-phonon couplings determined from electronic structure studies, a spectral density is constructed and employed to parametrize a linear vibronic coupling Hamiltonian. The diabatic coupling is found to depend noticeably on the inter-fragment distance, whose effect on the dynamics is here investigated. MCTDH calculations of the nonadiabatic transfer dynamics are carried out for the two most relevant electronic states and 60 phonon modes. The electron transfer process is found to be ultrafast and mediated by electronic coherence, resulting in characteristic oscillatory features during a period of about 100 fs.

Strikingly different effects of hydrogen bonding on the photodynamics of individual nucleobases in DNA: comparison of guanine and cytosine.

Ab initio surface hopping dynamics calculations were performed to study the photophysical behavior of cytosine and guanine embedded in DNA using a quantum mechanical/molecular mechanics (QM/MM) approach. It was found that the decay rates of photo excited cytosine and guanine were affected in a completely different way by the hydrogen bonding to the DNA environment. In case of cytosine, the geometrical restrictions exerted by the hydrogen bonds did not influence the relaxation time of cytosine significantly due to the generally small cytosine ring puckering required to access the crossing region between excited and ground state. On the contrary, the presence of hydrogen bonds significantly altered the photodynamics of guanine. The analysis of the dynamics indicates that the major contribution to the lifetime changes comes from the interstrand hydrogen bonds. These bonds considerably restricted the out-of-plane motions of the NH(2) group of guanine which are necessary for the ultrafast decay to the ground state. As a result, only a negligible amount of trajectories decayed into the ground state for guanine embedded in DNA within the simulation time of 0.5 ps, while for comparison, the isolated guanine relaxed to the ground state with a lifetime of about 0.22 ps. These examples show that, in addition to phenomena related to electronic interactions between nucleobases, there also exist relatively simple mechanisms in DNA by which the lifetime of a nucleobase is significantly enhanced as compared to the gas phase.

QM/MM non-adiabatic decay dynamics of formamide in polar and non-polar solvents.

Non-adiabatic on-the-fly dynamics simulations of the photodynamics of formamide in water and n-hexane were performed using a QM/MM approach. It was shown that steric restrictions imposed by the solvent cage do not have an influence on the initial motion which leads to the lowest energy conical intersection seam. The initial deactivation in water is faster than in n-hexane and in the gas phase. However, most of the formamide molecules in water do not reach the ground state. The reason for the deactivation inefficiency in water is traced back to a decrease of close CO···HOH and NH···OH(2) contacts which fall in the range of hydrogen bonds. The energy deposition into H-bond breaking events leaves molecules with less energy for surmounting the CN dissociation barrier. In both solvents, after hopping to the ground state, the solvent cage keeps the HCO and NH(2) fragments or CO and NH(3) products in close proximity. Consequently, the number of trajectories where fast recombination happens is augmented with delayed recombinations that start when the dissociation fragments hit the cage wall and return back. The hot ground state formamide is formed in an internal conversion process identical to the path leading to CN photodissociation. In the case of aqueous formamide, good agreement with experimental results is achieved by combining dynamics simulations starting from the S(1) and the S(2) excited states collecting high and low energy trajectories, respectively.

Azomethane: nonadiabatic photodynamical simulations in solution.

The nonadiabatic deactivation of trans-azomethane starting from the nπ* state has been investigated in gas phase, water, and n-hexane using an on-the-fly surface-hopping method. A quantum mechanical/molecular mechanics (QM/MM) approach was used employing a flexible quantum chemical level for the description of electronically excited states and bond dissociation (generalized valence bond perfect-pairing complete active space). The solvent effect on the lifetime and structural parameters of azomethane was investigated in detail. The calculations show that the nonadiabatic deactivation is characterized by a CNNC torsion, mainly impeded by mechanic interaction with the solvent molecules. The similar characteristics of the dynamics in polar and nonpolar solvent indicate that solvent effects based on electrostatic interactions do not play a major role. Lifetimes increase by about 20 fs for both solvents with respect to the 113 fs found for the gas phase. The present subpicosecond dynamics also nicely show an example of the suppression of C-N dissociation by the solvent cage.

Matrix-controlled photofragmentation of formamide: dynamics simulation in argon by nonadiabatic QM/MM method.

The short-time photodynamics (2 ps) of formamide embedded into an Ar matrix starting from the low-lying singlet excited S(1) (n(0)π*) and S(2) (ππ*) states were explored using a nonadiabatic photodynamics QM/MM approach. The interaction between formamide and the Ar matrix is taken into account at the MM level by means of Lennard-Jones potentials. This is the first example of exploring photodissociation of formamide with full nonadiabatic dynamics in a matrix and it nicely illustrates importance of considering environmental effects on photodissociation behavior of the peptide bond. It is shown that embedding of the formamide molecule in the argon matrix has strong impact on the outcome of the process. This is illustrated by formation of the 1:1 complex between ammonia and CO and prevention of full separation of the NH(2)˙ and HCO˙ subunits in the NH(2)˙ + HCO˙ radical pair. In addition, the argon matrix strongly influences the lifetime of the S(1) state, which increases by 211 fs relative to the gas phase.

Nonadiabatic excited-state dynamics with hybrid ab initio quantum-mechanical/molecular-mechanical methods: solvation of the pentadieniminium cation in apolar media.

A new implementation of nonadiabatic excited-state dynamics using hybrid methods is presented. The current approach is aimed at the simulation of photoexcited molecules in solution. The chromophore is treated at the ab initio level, and its interaction with the solvent is approximated by point charges within the electrostatic embedding approach and by a Lennard-Jones potential for the nonbonded interactions. Multireference configuration interaction (MRCI) and multiconfiguration self-consistent field (MCSCF) methods can be used. The program implementation has been performed on the basis of the Columbus and Newton-X program systems. For example, the dynamics of penta-2,4-dien-1-iminium (PSB3) and 4-methyl-penta-2,4-dien-1-iminium cations (MePSB3) was investigated in gas phase and in n-hexane solution. The excited-state (S(1)) lifetime and temporal evolution of geometrical parameters were computed. In the case of PSB3 the n-hexane results resemble closely the gas phase data. MePSB3, however, shows a distinct extension of lifetime due to steric hindering of the torsion around the central bond because of solute-solvent interactions.

Does stacking restrain the photodynamics of individual nucleobases?

Nonadiabatic photodynamical simulations of 4-aminopyrimidine (4-APy) used as a model for adenine were performed by embedding it between two stacking methyl-guanine (mGua) molecules to determine the effect of spatial restrictions on the ultrafast photodeactivation mechanism of this nucleobase. A hybrid multiconfigurational ab initio/molecular mechanical approach in combination with surface hopping was used. During the dynamics the formation of a significant fraction of intrastrand hydrogen bonding from 4-APy to mGua above and below is observed. These findings show that this type of hydrogen bond may play an important role for the photodynamics within one DNA strand and that it should be of interest even in irregular segments of double stranded nucleic acids structures. The relaxation mechanism of internal conversion to the ground state is dominated by ring puckering, and an overall elongation of the lifetime of the embedded system by approximately 20% as compared to the isolated 4-APy is computed.

Photodynamics of azomethane: a nonadiabatic surface-hopping study.

The internal conversion and hot ground-state dynamics of trans- and cis-azomethane starting in the S(1) state have been investigated by nonadiabatic ab initio surface hopping dynamics using MCSCF-GVB-CAS and MRCISD methods and by determining energy minima and saddle points, minima on the crossing seam, and minimum energy pathways on the ground and first excited-state surfaces. The lifetimes and photoproducts from the dynamics simulations, geometric properties, excitation energies of selected stationary points and minimum energy pathways between them are reported. Our results favor a statistical model with trans-AZM moving to the ground-state minima before the first CN dissociation takes place. A detailed discussion in comparison to recent experimental and theoretical data is presented.

Nonadiabatic excited-state dynamics of polar pi-systems and related model compounds of biological relevance.

Multireference ab initio dynamics simulations have become available as a tool for the investigation of photochemical processes, mainly for those related to nonadiabatic phenomena taking place in the sub-picosecond time scale. For organic molecules, these phenomena are in many cases deeply dependent on the relaxation of the photoexcited pi-system. We review the latest contributions of our group to this subject and report new results for systems studied previously, grouping them in single pi bonds, chains and aromatic rings. The dynamics of ethylene and substituted ethylenes is discussed mainly in connection to the competition between the two available relaxation paths in the excited states and their relation to the conical intersections in large systems. The trans-cis and the cis-trans dynamics of the pentadieniminium cation is investigated as well. Finally, we discuss the photodynamics of aminopyrimidine starting in the S1 and S2 states and the conclusions, which can be drawn from this for the interpretation of the adenine dynamics.

Multilayer adsorption of water at a rutile TiO2)(110) surface: towards a realistic modeling by molecular dynamics.

We present a model combining ab initio concepts and molecular dynamics simulations for a more realistic treatment of complex adsorption processes. The energy, distance, and orientation of water molecules adsorbed on stoichiometric and reduced rutile TiO(2)(110) surfaces at 140 K are studied via constant temperature molecular dynamics simulations. From ab initio calculations relaxed atomic geometries for the surface and the most probable adsorption sites were derived. The study comprises (i) large two-dimensional surface supercells, providing a realistically low concentration of surface oxygen defects, and (ii) a water coverage sufficiently large to model the onset of the growth of a bulk phase of water on the surface. By our combined approach the influence of both, the metal oxide surface, below, and the bulk water phase, above, on the water molecules forming the interface between the TiO(2) surface and the water bulk layer is taken into account. The good agreement of calculated adsorption energies with experimental temperature programmed desorption spectra demonstrates the validity of our model.