We present an efficient approach for surface hopping-based nonadiabatic dynamics in the condensed phase. For the systems studied, a restricted Kohn-Sham orbital formulation of the delta self-consistent field (ΔSCF) method was used for efficient calculation of excited electronic states. Timedependent density functional theory (DFT) is applied to aid excited-state SCF convergence and provide guess electronic state densities. Aside from that the Landau-Zener procedure simplifies the surface hopping between electronic states. By utilizing the combined Gaussian and plane waves approach with periodic boundary conditions the method is easily applicable to full atomistic DFT simulations of condensedphase systems and was used to study the nonradiative deactivation mechanism of photoexcited diimide in water solution.
The mechanisms of nonradiative deactivation of a phenylalanine residue after near-UV photoexcitation have been investigated in an isolated peptide chain model (N-acetylphenylalaninylamide, NAPA) both experimentally and theoretically. Lifetime measurements at the origin of the first ππ* state of jet-cooled NAPA molecules have shown that (i) among the three most stable conformers of the molecule, the folded conformer NAPA B is ∼50-times shorter lived than the extended major conformer NAPA A and (ii) this lifetime is virtually insensitive to deuteration at the NH(2) and NH sites. Concurrent time-dependent density functional theory (TDDFT) based nonadiabatic dynamics simulations in the full dimensionality, carried out for the NAPA B conformer, provided direct insights on novel classes of ultrafast deactivation mechanisms, proceeding through several conical intersections and leading in fine to the ground state. These mechanisms are found to be triggered either (i) by a stretch of the N(Phe)H bond, which leads to an H-transfer to the ring, or (ii) by specific backbone amide distortions. The potential energy surfaces of the NAPA conformers along these critical pathways have been characterized more accurately using the coupled cluster doubles (CC2) method and shown to exhibit barriers that can be overcome with moderate excess energies. These results analyzed in the light of the experimental findings enabled us to assign the short lifetime of NAPA B conformer to a number of easily accessible exit channels from the initial ππ* surface, most importantly the one involving a transfer of electronic excitation to an nπ* surface, induced by distortions of the backbone peptide bond.
The excitation wavelength dependent photodynamics of pyrrole are investigated by nonadiabatic trajectory-surface-hopping dynamics simulations based on time dependent density functional theory (TDDFT) and the algebraic diagrammatic construction method to the second order (ADC(2)). The ADC(2) results confirm that the N-H bond dissociation occurring upon excitation at the origin of the first excited state, S1(πσ*), is driven by tunnelling [Roberts et al., Faraday Discuss., 2013, 163, 95] as a barrier of ΔE = 1780 cm(-1) traps the population in a quasi-bound minimum. Upon excitation to S1(πσ*) in the wavelength range of 236-240 nm, direct dissociation of the N-H bond takes place with a time constant of 28 fs. The computed time constant is in very good agreement with recently reported measurements. Excitation to the lowest B2 state in the 198-202 nm range returns a time constant for N-H fission of 48 fs at the B3LYP/def2-TZVPD level, in perfect agreement with the experiment [Roberts et al. Faraday Discuss., 2013, 163, 95]. For the same wavelength range the ADC(2)/aug-cc-pVDZ decay constant is more than three times longer than the experimentally reported one. The accuracy of the B3LYP/def2-TZVPD dynamics is checked against reference complete-active-space second-order perturbation theory (CASPT2) calculations and explained in terms of correct topography of the ππ* surface and the lack of mixing between the ππ* and the 3px Rydberg states which occurs in the ADC(2) method.
Irradiation of 2-phenyl-1-naphthol (6) in CH(3) CN/D(2) O (3:1) leads to very efficient incorporation of deuterium at the ortho-positions of the adjacent phenyl ring (overall Φ=0.73±0.07), along with minor incorporation at the naphthalene positions 5 and 8. These finding are explained by excited state intramolecular proton transfer (ESIPT) from the phenolic OH group to the corresponding carbon atoms, the main pathway giving rise to quinone methide (QM) 7, which has been characterized by LFP (τ≈20 ns; 460 nm). The ESIPT reaction paths have been explored with the second order approximate coupled cluster (CC2) method. In nonprotic solvents the ESIPT from the naphthol O-H to the ortho-position of the phenyl ring proceeds in a barrierless manner along the (1) L(a) energy surface via a conical intersection with the S(0) state, delivering 7. In aqueous solvent, clusters with H(2) O are formed wherein proton transfer (PT) to solvent and a H(2) O-mediated relay mechanism gives rise to naphtholates and QMs. The results are compared with 2-phenylphenol (3) that also undergoes barrierless ESIPT giving a QM via a conical intersection. However, due to an unfavorable conformation in the ground state, the quantum efficiency for ESIPT of 3 is significantly lower (Φ for D-exchange=0.041). These results show that ESIPT from phenol to carbon need not be an intrinsically inefficient process.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.