Oliver Weingart

Conical intersection dynamics of the primary photoisomerization event in vision

Ever since the conversion of the 11-cis retinal chromophore to its all-trans form in rhodopsin was identified as the primary photochemical event in vision, experimentalists and theoreticians have tried to unravel the molecular details of this process. The high quantum yield of 0.65 (ref. 2), the production of the primary ground-state rhodopsin photoproduct within a mere 200 fs (refs 3–7), and the storage of considerable energy in the first stable bathorhodopsin intermediate all suggest an unusually fast and efficient photoactivated one-way reaction. Rhodopsins unique reactivity is generally attributed to a conical intersection between the potential energy surfaces of the ground and excited electronic states enabling the efficient and ultrafast conversion of photon energy into chemical energy. But obtaining direct experimental evidence for the involvement of a conical intersection is challenging: the energy gap between the electronic states of the reacting molecule changes significantly over an ultrashort timescale, which calls for observational methods that combine high temporal resolution with a broad spectral observation window. Here we show that ultrafast optical spectroscopy with sub-20-fs time resolution and spectral coverage from the visible to the near-infrared allows us to follow the dynamics leading to the conical intersection in rhodopsin isomerization. We track coherent wave-packet motion from the photoexcited Franck–Condon region to the photoproduct by monitoring the loss of reactant emission and the subsequent appearance of photoproduct absorption, and find excellent agreement between the experimental observations and molecular dynamics calculations that involve a true electronic state crossing. Taken together, these findings constitute the most compelling evidence to date for the existence and importance of conical intersections in visual photochemistry.

MOLCAS 8: New Capabilities for Multiconfigurational Quantum Chemical Calculations across the Periodic Table

In this report, we summarize and describe the recent unique updates and additions to the Molcas quantum chemistry program suite as contained in release version 8. These updates include natural and spin orbitals for studies of magnetic properties, local and linear scaling methods for the Douglas–Kroll–Hess transformation, the generalized active space concept in MCSCF methods, a combination of multiconfigurational wave functions with density functional theory in the MC‐PDFT method, additional methods for computation of magnetic properties, methods for diabatization, analytical gradients of state average complete active space SCF in association with density fitting, methods for constrained fragment optimization, large‐scale parallel multireference configuration interaction including analytic gradients via the interface to the Columbus package, and approximations of the CASPT2 method to be used for computations of large systems. In addition, the report includes the description of a computational machinery for nonlinear optical spectroscopy through an interface to the QM/MM package Cobramm. Further, a module to run molecular dynamics simulations is added, two surface hopping algorithms are included to enable nonadiabatic calculations, and the DQ method for diabatization is added. Finally, we report on the subject of improvements with respects to alternative file options and parallelization.

Cooperating Dinitrogen and Phenyl Rotations in trans-Azobenzene Photoisomerization

Semiempirical OM2/MRCI surface-hopping simulations have been performed to study the trans-to-cis photoisomerization of azobenzene upon excitation to the S1 state. The decay dynamics to the ground state shows an oscillatory pattern that can be attributed to an out-of-plane rotation of the N2 moiety. The reaction is thus initially driven by N2 rotation which triggers phenyl rotations around the C-N bonds. The cis isomer is produced most effectively when the phenyl rings rotate in phase. Mode-specific excitations cause variations in the computed decay times and product yields.

Electron promotion and electronic friction in atomic collision cascades

We present a computer simulation model for the space- and time-resolved calculation of electronic excitation energy densities in atomic collision cascades. The model treats electronic friction as well as electron promotion as a source term of electronic energy that is carried away from the original point of excitation according to a nonlinear diffusion equation. While the frictional source is treated within the Lindhard model of electronic stopping, electron promotion is described using diabatic correlation curves derived from ab initio molecular orbital energy level calculations in combination with the Landau–Zener curve crossing model. Results calculated for two selected collision cascades show that the electron promotion mechanism may contribute significantly to the excitation energy density in the cascade volume, giving rise to distinct peaks of the local electron temperature at the surface. This contribution is essentially restricted to the first 100 fs after the projectile impact and may therefore be of significance for either external or internal kinetic electron emission. At later times, where the bombardment-induced particle kinetics lead to the sputter ejection of material from the surface, the excitation is shown to be primarily governed by electronic friction. This finding is important in light of excitation and ionization probabilities of sputtered particles.

Interfacial States in Donor–Acceptor Organic Heterojunctions: Computational Insights into Thiophene-Oligomer/Fullerene Junctions

Donor-acceptor heterojunctions composed of thiophene oligomers and C60 fullerene were investigated with computational methods. Benchmark calculations were performed with time-dependent density functional theory. The effects of varying the density functional, the number of oligomers, the intermolecular distance, the medium polarization, and the chemical functionalization of the monomers were analyzed. The results are presented in terms of diagrams where the electronic states are classified as locally excited states, charge-transfer states, and delocalized states. The effects of each option for computational simulations of realistic heterojunctions employed in photovoltaic devices are evaluated and discussed.

Wavepacket Splitting and Two‐Pathway Deactivation in the Photoexcited Visual Pigment Isorhodopsin

Isorhodopsin is the visual pigment analogue of rhodopsin. It shares the same opsin environment but it embeds 9-cis retinal instead of 11-cis. Its photoisomerization is three times slower and less effective. The mechanistic rationale behind this observation is revealed by combining high-level quantum-mechanical/molecular-mechanical simulations with ultrafast optical spectroscopy with sub-20 fs time resolution and spectral coverage extended to the near-infrared. Whereas in rhodopsin the photoexcited wavepacket has ballistic motion through a single conical intersection seam region between the ground and excited states, in isorhodopsin it branches into two competitive deactivation pathways involving distinct conical intersection funnels. One is rapidly accessed but unreactive. The other is slower, as it features extended steric interactions with the environment, but it is productive as it follows forward bicycle pedal motion.

Nonadiabatic Decay Dynamics of a Benzylidene Malononitrile

The photoinduced nonadiabatic decay dynamics of 2-[4-(dimethylamino)benzylidene]malononitrile (DMN) in the gas phase is investigated at the semiempirical OM2/MRCI level using surface hopping simulations. A lifetime of 1.2 ps is predicted for the S(1) state, in accordance with experimental observation. The dominant reaction coordinate is found to be the twisting around the C7═C8 double bond accompanied by pronounced pyramidalization at the C8 atom. Motion along this coordinate leads to the lowest-energy conical intersection (CI(01α)). Several other S(0)/S(1) conical intersections have also been located by full optimization but play no role in the dynamics. The time-resolved fluorescence spectrum of DMN is simulated by computing emission energies and oscillator strengths along the trajectories. It compares well with the experimental spectrum. The use of different active spaces in the OM2/MRCI calculations yields similar results and thus demonstrates their internal consistency.

Modelling Time-Resolved Two-Dimensional Electronic Spectroscopy of the Primary Photoisomerization Event in Rhodopsin

Time-resolved two-dimensional (2D) electronic spectra (ES) tracking the evolution of the excited state manifolds of the retinal chromophore have been simulated along the photoisomerization pathway in bovine rhodopsin, using a state-of-the-art hybrid QM/MM approach based on multiconfigurational methods. Simulations of broadband 2D spectra provide a useful picture of the overall detectable 2D signals from the near-infrared (NIR) to the near-ultraviolet (UV). Evolution of the stimulated emission (SE) and excited state absorption (ESA) 2D signals indicates that the S1 → SN (with N ≥ 2) ESAs feature a substantial blue-shift only after bond inversion and partial rotation along the cis → trans isomerization angle, while the SE rapidly red-shifts during the photoinduced skeletal relaxation of the polyene chain. Different combinations of pulse frequencies are proposed in order to follow the evolution of specific ESA signals. These include a two-color 2DVis/NIR setup especially suited for tracking the evolution of the S1 → S2 transitions that can be used to discriminate between different photochemical mechanisms of retinal photoisomerization as a function of the environment. The reported results are consistent with the available time-resolved pump–probe experimental data, and may be used for the design of more elaborate transient 2D electronic spectroscopy techniques.

Sampling excited state dynamics: influence of HOOP mode excitations in a retinal model

Zero point energy and classical thermal sampling techniques are compared in semi-classical photodynamics of the pentadienyliminium cation, a minimal retinal model. Using both methods, the effects of vibrational hydrogen-out-of-plane (HOOP) excitations on the photo-reactivity are probed at the ab initio CASSCF level. With 2376 individual trajectories the calculations reveal a clear picture of the relation between the excited state reaction coordinate, surface crossing and product generation. The productivity is strongly coupled with hydrogen torsion and the number of hopping attempts before the molecule finally decays.

Modelling vibrational coherence in the primary rhodopsin photoproduct.

Molecular dynamics simulations of the rhodopsin photoreaction reveal coherent low frequency oscillations in the primary photoproduct (photorhodopsin), with frequencies slightly higher than observed in the experiment. The coherent molecular motions in the batho-precursor can be attributed to the activation of ground state vibrational modes in the hot photo-product, involving out-of-plane deformations of the carbon skeleton. Results are discussed and compared with respect to spectroscopic data and suggested reaction mechanisms.