Matthias Ruckenbauer

Newton-X: a surface-hopping program for nonadiabatic molecular dynamics

The Newton‐X program is a general‐purpose program package for excited‐state molecular dynamics, including nonadiabatic methods. Its modular design allows Newton‐X to be easily linked to any quantum‐chemistry package that can provide excited‐state energy gradients. At the current version, Newton‐X can perform nonadiabatic dynamics using Columbus, Turbomole, Gaussian, and Gamess program packages with multireference configuration interaction, multiconfigurational self‐consistent field, time‐dependent density functional theory, and other methods. Nonadiabatic dynamics simulations with a hybrid combination of methods, such as Quantum‐Mechanics/Molecular‐Mechanics, are also possible. Moreover, Newton‐X can be used for the simulation of absorption and emission spectra. The code is distributed free of charge for noncommercial and nonprofit uses at www.newtonx.org. WIREs Comput Mol Sci 2014, 4:26–33. doi: 10.1002/wcms.1158

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.

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.

The photodynamics of ethylene: A surface-hopping study on structural aspects

Simulations of the photodynamics of ethylene were carried out by employing the semiempirical direct trajectory with surface hopping method in order to assess quantitatively the importance of different regions of the S(2)S(1) and S(1)S(0) crossing seams. The results show that during the first 50 fs after a vertical photoexcitation to the pipi(*) state, the nonadiabatic coupling between the S(1) and the S(2) states produces a recurrence pattern of oscillation of the populations in these states. Within the first 100 fs, the S(1) state population spans a limited region of the configuration space between the initial geometries and the twisted-pyramidalized minimum on the crossing seam (MXS). Depending on the way of counting, about 50% of the S(1)-->S(0) transitions occur in the pyramidalized region of the crossing seam, but not necessarily close to the MXS. The remaining 50% occurs in the H-migration and ethylidene regions. Our analysis shows that the ethylidene region becomes more important in later stages of the dynamics when the flux of trajectories that was not effectively converted to the ground state in the pyramidalized region starts to reach this part of the configuration space. The excited-state nonadiabatic dynamics could be employed to generate suitable initial phase space distributions for the hot-ethylene ground-state kinetic studies.

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.

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.

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.

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.

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.