Nonadiabaticab initiomolecular dynamics study of photoisomerization inN-salicilydenemethylfurylamine (SMFA)

Spectroscopic Identification of Hydrogen Bond Vibrations and Quasi-Isostructural Polymorphism in N-Salicylideneaniline

The ortho-hydroxy aryl Schiff base 2-[(E)-(phenylimino)methyl]phenol and its deutero-derivative have been studied by the inelastic incoherent neutron scattering (IINS), infrared (IR) and Raman experimental methods, as well as by Density Functional Theory (DFT) and Density-Functional Perturbation Theory (DFPT) simulations. The assignments of vibrational modes within the 3500–50 cm⁻¹ spectral region made it possible to state that the strong hydrogen bond in the studied compound can be classified as the so-called quasi-aromatic bond. The isotopic substitution supplemented by the results of DFT calculations allowed us to identify vibrational bands associated with all five major hydrogen bond vibrations. Quasi-isostructural polymorphism of 2-[(E)-(phenylimino)methyl]phenol (SA) and 2-[(E)-(phenyl-D₅-imino)methyl]phenol (SA-C₆D₅) has been studied by powder X-ray diffraction in the 20–320 K temperature range.

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.

Position of the Benzene Ring Substituent Regulates the Excited-State Deactivation Process of the Benzyluracil Systems

A systematic theoretical study of the regulating effect of the substituent position on the photoinduced deactivation process of the benzyluracil systems has been performed based on the high-level static electronic structure calculations and on-the-fly full-dimensional excited-state dynamics simulations. Similarities and differences coexist for the two systems by comparative studies on the photoinduced deactivation process of the 5-benzyluracil (5-BU) and 6-benzyluracil (6-BU) systems. They both obey an S₂ → S₁ → S₀ two-step decay pattern, and the decay coordinates of the S₂ → S₁ and S₁ → S₀ processes are mainly driven by the elongation of the bridging bond and the out-of-plane ring deformation motion, respectively. However, the puckering motion occurring at the C2 atom in the uracil fragment dominates the decay pathway of the 5-BU system. On the contrary, the puckering motion at the C5 atom in the benzene fragment mainly drives the decay coordinate of the 6-BU system. Therefore, the substituent position could play significant roles in the deactivation process of the benzyluracil systems. Moreover, the S₁ → S₀ decay process of the 6-BU system consists of five pathways, possessing a more complex deactivation picture than the 5-BU system. The fitted time scale of the puckering motion is compatible with the experimentally observed lifetimes. This work provides a fundamental understanding of the photophysical and photochemical properties of the benzyluracil systems and can give rational suggestions to further design or regulate the bionic molecular systems.

Theoretical studies of the ultrafast deactivation mechanism of 8-oxo-guanine on the S1 and S2 electronic states in gas phase

The 8-oxo-deoxyguanosine is the most abundant specie of the DNA oxidative damage. Despite the deleterious effects such as gene mutation it may cause, the 8-oxodG was also reported to have beneficial effect such as repairing the nearby cyclobutane pyrimidine dimer (CPD) after photoexcitation. Due to its strong biological relevance, the photoinduced excited state dynamics behavior of 8-oxo-deoxyguanosine is of particular interest. In this work, a theoretical investigation by combination of complete active space self-consistent field (CASSCF) ab initio calculations and on-the-fly nonadiabatic dynamics simulations are implemented to provide intrinsic deactivation mechanism of its free base 8-oxoguanine after being excited to the S₁ and S₂ states. Two minimum energy conical intersections (MECIs) characterized by the C3-puckered motion with attractive chiral character are located, which contribute appreciably to the S₁ state deactivation process. When the system is being excited to the S₂ state directly, a S₂ → S₁ → S₀ two-step decay pattern is proposed. A nearly planar S₂/S₁ intersection plays a significant role in the S₂ → S₁ decay process. The subsequent S₁ state relaxation process is also dominated by the C3-puckered deformation motion. One decay time is estimated to be 704 fs, which compares well with the experimental observation of 0.9 ± 0.1 ps in solvents. Particular illustration is the fact that the MECIs configurations we located bear an exceptional resemblance with previous reported thymine, cytosine and guanine, suggesting that the current work could lend support for better understanding of the non-natural nucleobases and derivatives.

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.

The photoinduced isomerization mechanism of the 2-(1-(methylimino)methyl)-6-chlorophenol (SMAC): Nonadiabatic surface hopping dynamics simulations

The photophysical properties of the Schiff base family are crucial for their applications such as molecular switches and molecular memories. However, it was found that the photophysical behavior is not uniform for all Schiff base molecules, which shows a significant substituent dependent property. In this article, we studied the photoisomerization mechanism of one Schiff base chlorosubstituted derivative 2-(1-(methylimino)methyl)-6-chlorophenol by employing geometrical optimization, energy profiles scanning, and on-the-fly dynamical simulations. Three types of minimum energy conical intersections were located on the S₁/S₀ crossing seam, with two characterized by twisting motion of the C=N bond and one featured with the excited state intramolecular proton transfer process and then twisting motion around the C=C bond [excited-state intramolecular proton transfer process (ESIPT)-then-twisting]. By a combination of the dynamics simulation results with the energy profiles scanned along with the ESIPT coordinate, it was found that the photophysical property of the targeted molecule is different from that of most Schiff base members, which prefer to decay by a twisting motion around the C=N bridge bond rather than the ESIPT-then-twisting channel. The minor ESIPT channel is probably governed by a tunneling mechanism. The proposed deactivation mechanism can provide a reasonable explanation for the observations in the experiment and would provide fundamental indications for further design of new and efficient photochromic products.

Nonadiabatic dynamics simulation of photoisomerization mechanism of the second stablest isomer of N-salicilydenemethylfurylamine

The photoisomerization processes of the second stablest isomer in the aromatic Schiff base, N-salicilydenemethylfurylamine, in the gas phase have been studied by static electronic structure calculations and surface-hopping dynamics simulations based on the Zhu-Nakamura theory. Various stable structures are obtained in the optimization because of different orientations of methyl-furyl part with respect to the salicylaldimine part and different orientations of hydroxy group with respect to the benzene ring. Upon photoexcitation into the first excited state, bond isomerization in the salicylaldimine part is completely suppressed until the strong excited-state hydrogen bond is broken. The decay pathway involves two excited-state minima, one in cis-enol form and the other in cis-keto form. After the excited-state proton transfer, twists of bonds lead to a conical intersection between the ground and excited states. After internal conversion around a conical intersection, the molecule is stabilized in cis- or trans-keto form. If the reverse hydrogen transfer process occurs in the ground state, the molecule will finally end up in the cis-enol region. The cis-keto and trans-keto isomers are observed as photoproducts. According to our full-dimensional nonadiabatic dynamics simulations, we find the excited-state intramolecular proton transfer and torsions of three single bonds in the chain to be responsible for photoisomerization of the second stablest isomer of N-salicilydenemethylfurylamine.

Performance of TDDFT with and without spin-flip in trajectory surface hopping dynamics: cis–trans azobenzene photoisomerization

Quantum yields simulated by LR-TDDFT and SF-TDDFT methods for azobenzene photoisomerizations.