Ultrafasttrans-cisphotoisomerization of the neutral chromophore in green fluorescent proteins: Surface-hopping dynamics simulation

Wavelength-dependent photoisomerization of trans-4,4'-azopyridine: Nonadiabatic dynamics simulation

The trans-cis photoisomerization processes of 4,4'-azopyridine upon S1 and S2 excitations have been investigated by nonadiabatic dynamics simulations based on multi-reference CASSCF calculations. 119 sampling trajectories were simulated starting from the trans form excited to the S1 (S2) state and the cis-isomer quantum yield is evaluated to be (3 ± 2)% ((18 ± 4)%), which is qualitatively in agreement with the recent experimental results in ethanol. We found that rotation around the central N-N bond accompanied by the N-N-C symmetrical bending vibrations is the main mechanism in photoisomerization of the target molecule excited to the S1 and S2 states. Upon S1 excitation, S1-S0 transition occurs earlier along the C-N-N-C torsional coordinate, leading to a low cis-isomer quantum yield. Upon S2 excitation, half of the simulated trajectories are trapped in a potential well on the S2 state, from which the twisted conical intersections are more easily reached in the internal conversion, resulting in a higher cis-isomer quantum yield.

πσ*-Mediated Nonadiabatic Tunneling Dynamics of Thiophenols in S1: The Semiclassical Approaches

The S-H bond tunneling predissociation dynamics of thiophenol and its ortho-substituted derivatives (2-fluorothiophenol, 2-methoxythiophenol, and 2-chlorothiphenol) in S₁ (ππ*) where the H atom tunneling is mediated by the nearby S₂ (πσ*) state (which is repulsive along the S-H bond extension coordinate) have been investigated in a state-specific way using the picosecond time-resolved pump-probe spectroscopy for the jet-cooled molecules. The effects of the specific vibrational mode excitations and the SH/SD substitutions on the S-H(D) bond rupture tunneling dynamics have been interrogated, giving deep insights into the multidimensional aspects of the S₁/S₂ conical intersection, which also shapes the underlying adiabatic tunneling potential energy surfaces (PESs). The semiclassical tunneling rate calculations based on the Wentzel-Kramers-Brillouin (WKB) approximation or Zhu-Nakamura (ZN) theory have been carried out based on the ab initio PESs calculated in the (one, two, or three) reduced dimensions to be compared with the experiment. Though the quantitative experimental results could not be reproduced satisfactorily by the present calculations, the qualitative trends among different molecules in terms of the behavior of the tunneling rate versus the (adiabatic) barrier height or the number of PES dimensions could be rationalized. Most interestingly, the H/D kinetic isotope effect observed in the tunneling rate could be much better explained by the ZN theory compared to the WKB approximation, indicating that the nonadiabatic coupling matrix elements should be invoked for understanding the tunneling dynamics taking place in the proximity of the conical intersection.

Ultrafast Photoisomerization of N-(2-Methoxybenzylidene)aniline: Nonadiabatic Surface-Hopping Study

We investigated the ultrafast photoisomerization of N-(2-methoxybenzylidene)aniline in the gas phase excited into the second singlet (S₂) state by nonadiabatic surface-hopping dynamics calculations. Two trans isomers (1E and 1E') were taken into consideration in our dynamics simulation. Three conical intersections (CIs) were characterized in the optimization. The CI between S₂ and the first singlet (S₁) states presents a nearly planar structure, while the other two CIs (CI_twist-I and CI_twist-II) between S₁ and the ground (S₀) states show nearly perpendicular geometries. After two trans isomers excited to the S₂ state, the torsion of the C-N bond connected the phenyl group and the stretch of the central bridging bond make the molecule reach CI_planar, and the S₂/S₁ hopping occurs. During the S₁-state dynamics, the molecule moves to a S₁/S₀ CI (CI_twist-I or CI_twist-II) by the rotation of the central bridging bond. The cis isomer is obtained through a barrierless pathway in the S₀ state with the torsion of the three bridging bonds. There is a main channel and an alternative one for the photoisomerization process of both trans isomers. CI_twist-I and CI_twist-II act as S₁/S₀ decay funnels in the main isomerization channels of 1E and 1E' isomers, respectively, and the photochemical processes of 1E and 1E' lead to different cis isomers.

Volume-conserving photoisomerization of a nonplanar GFP chromophore derivative: Nonadiabatic dynamics simulation

Nonadiabatic dynamics of a nonplanar green fluorescent protein (GFP) chromophore derivative was examined theoretically by trajectory surface hopping approach at the CASSCF level based on Zhu-Nakamura theory. The geometry optimizations show that there are four ground-state minima and four conical intersections between the ground (S₀) and first excited (S₁) states. Four S₁-state minima were found at a perpendicularly twisted conformation around the imidazolinone-bridged bond. Upon excitation to S₁ state, the main decay pathways of four isomers involve different S₀/S₁ potential energy surface conical intersections. The dominant excited-state relaxation mechanism of this GFP chromophore derivative is the twists of two bridging bonds and the methyl group in the bridge combined with the pyramidalization character of the central carbon atom. Further twists of two bridging bonds and the methyl group occur sequentially in the S₀ state. It is worth to mention that the special volume-conserving motion of this molecule is attributed to twists of two bridging bonds in the same direction during the whole photoisomerization processes. The theoretical investigation presented herein provides important insights into the volume-conserving photoisomerization mechanisms of a nonplanar GFP chromophore derivative. We believe that the present work can benefit the design of the photochromic molecule devices in confined media.