Nonadiabatic Decay Dynamics of a Benzylidene Malononitrile

Advanced theoretical design of light-driven molecular rotary motors: enhancing thermal helix inversion and visible-light activation

In this study, we have advanced the field of light-driven molecular rotary motors (LDMRMs) by achieving two pivotal goals: lowering the thermal helix inversion (THI) barrier and extending the absorption wavelength into the visible spectrum. This study involves the structural reengineering of a second-generation visible LDMRM, resulting in the synthesis of a novel class, specifically, 2-((2S)-5-methoxy-2-methyl-2,3-dihydro-1H-cyclopenta[a]naphthalen-1-yl)-3-oxo-2,3-dihydro-1H-dibenzo[e,g]indole-6,9-dicarbonitrile. This redesigned motor stands out with its two photoisomerization stages and two thermal helix inversions, featuring exceptionally low THI barriers (4.00 and 2.05 kcal mol-1 at the OM2/MRCI level for the EM → EP and ZM → ZP processes, respectively). Moreover, it displays absorption wavelengths in the visible light range (482.98 and 465.76 nm for the EP and ZP isomers, respectively, at the TD-PBE0-D3/6-31G(d,p) level), surpassing its predecessors in efficiency, as indicated by the narrow HOMO-LUMO energy gap. Ultrafast photoisomerization kinetics (approximately 0.8-1.6 ps) and high quantum yields (around 0.3-0.6) were observed through trajectory surface hopping simulations. Additionally, the simulated time-resolved fluorescence emission spectrum indicates a significantly reduced "dark state" duration (0.09-0.26 ps) in these newly designed LDMRMs compared to the original ones, marking a substantial leap forward in the design and efficiency of LDMRMs.

On-the-fly simulation of time-resolved fluorescence spectra and anisotropy

We combine on-the-fly trajectory surface hopping simulations and the doorway-window representation of nonlinear optical response functions to create an efficient protocol for the evaluation of time- and frequency-resolved fluorescence (TFRF) spectra and anisotropies of the realistic polyatomic systems. This approach gives the effective description of the proper (e.g., experimental) pulse envelopes, laser field polarizations, and the proper orientational averaging of TFRF signals directly from the well-established on-the-fly nonadiabatic dynamic simulations without extra computational cost. To discuss the implementation details of the developed protocol, we chose cis-azobenzene as a prototype to simulate the time evolution of the TFRF spectra governed by its nonadiabatic dynamics. The results show that the TFRF is determined by the interplay of several key factors, i.e., decays of excited-state populations, evolution of the transition dipole moments along with the dynamic propagation, and scaling factor of the TFRF signals associated with the cube of emission frequency. This work not only provides an efficient and effective approach to simulate the TFRF and anisotropies of realistic polyatomic systems but also discusses the important relationship between the TFRF signals and the underlining nonadiabatic dynamics.

Effect of Temperature on Photoisomerization Dynamics of a Newly Designed Two-Stroke Light-Driven Molecular Rotary Motor

The working mechanism of conventional light-driven molecular rotary motors, especially Feringa-type motors, contains two photoisomerization steps and two thermal helix inversion steps. Due to the existence of a thermal helix inversion step, both the ability to work at lower temperatures and the rotation speed are limited. In this work, a two-stroke light-driven molecular rotary motor, 2-(1,5-dimethyl-4,5-dihydrocyclopenta[b]pyrrol-6(1H)-ylidene)-1,2-dihydro-3H-pyrrol-3-one (DDPY), is proposed, which is capable of performing unidirectional and repetitive rotation by only two photoisomerization (EP→ZP and ZP→EP) steps. With trajectory surface-hopping simulation at the semi-empirical OM2/MRCI level, the EP→ZP and ZP→EP nonadiabatic dynamics of DDPY were systematically studied at different temperatures. Both EP→ZP and ZP→EP photoisomerizations are on an ultrafast timescale (ca. 200–300 fs). The decay mode of EP→ZP photoisomerization is approximately bi-exponential, while that of ZP→EP photoisomerization is found to be periodic. For EP and ZP isomers of DDPY, after the S0→S1 excitation, the dynamical processes of nonadiabatic decay are both followed by twisting about the central C=C double bond and the pyramidalization of the C atom at the stator-axle linkage. The effect of temperature on the nonadiabatic dynamics of EP→ZP and ZP→EP photoisomerizations of DDPY has been systematically investigated. The average lifetimes of the S1 excited state and quantum yields for both EP→ZP and ZP→EP photoisomerization are almost temperature-independent, while the corresponding unidirectionality of rotation is significantly increased (e.g., 74% for EP→ZP and 72% for ZP→EP at 300 K vs 100% for EP→ZP and 94% for ZP→EP at 50 K) with lowering the temperature.

Design and Nonadiabatic Photoisomerization Dynamics Study of a Three-Stroke Light-Driven Molecular Rotary Motor

Working cycle of conventional light-driven molecular rotary motors (LDMRMs), especially Feringa-type motors, usually have four steps, two photoisomerization steps, and two thermal helix inversion (THI) steps. THI steps hinder the ability of the motor to operate at lower temperatures and limit the rotation speed of LDMRMs. A three-stroke LDMRM, 2-(2,7-dimethyl-2,3-dihydro-1H-inden-1-ylidene)-1,2-dihydro-3H-pyrrol-3-one (DDIY), is proposed, which is capable of completing an unidirectional rotation by two photoisomerization steps and one thermal helix inversion step at room temperature. On the basis of trajectory surface-hopping simulation at the semi-empirical OM2/MRCI level, the EP→ZP and ZP→EM nonadiabatic photoisomerization dynamics of DDIY were systematically analyzed. Quantum yields of EP→ZP and ZP→EM photoisomerization of DDIY are ca. 34% and 18%, respectively. Both EP→ZP and ZP→EM photoisomerization processes occur on an ultrafast time scale (ca. 100-300 fs). This three-stroke LDMRM may stimulate further research for the development of new families of more efficient LDMRMs.

Intramolecular charge transfer of a push-pull chromophore with restricted internal rotation of an electron donor

Intramolecular charge transfer (ICT) of 4-(dicyanomethylene)-2-methyl-6-[2-(2,3,6,7-tetrahydro-1H,5H-benzo[ij]quinolizin-9-yl)vinyl]-4H-pyran (LD688) in DMSO solution was investigated by femtosecond stimulated Raman spectroscopy (FSRS) with 403 nm excitation. The molecular structure of LD688 is similar to that of a well-known push-pull chromophore, 4-dicyanomethylene-2-methyl-6-(p-dimethylaminostyryl)-4H-pyran (DCM), except that the internal rotation of the electron-donating dimethylamino group is restricted with the introduction of the julolidine moiety. Upon photo-excitation, LD688 shows an ultrafast (1.0 ps) ICT followed by the vibrational relaxation (3-8 ps) in the charge-transfer (CT) state. Two distinct Raman spectra of LD688 in the locally excited (LE) and CT state of the S₁ state were retrieved from FSRS measurements. Based on the time-dependent density functional theory (TDDFT) simulations, a "twisted" julolidine geometry of LD688 was proposed for the ICT state, which was further confirmed in comparison to the spectral changes of several push-pull chromophores with the π-conjugated backbone of stilbene, biphenyl, styrylpyran, styrylpyridinium, and styrene in terms of the skeletal vibrational modes of ν_19b,py, ν_CC,ph, and ν_CN.

Ultrafast Internal Conversion Dynamics through the on-the-Fly Simulation of Transient Absorption Pump-Probe Spectra with Different Electronic Structure Methods

An on-the-fly surface-hopping simulation protocol is developed for the evaluation of transient absorption (TA) pump-probe (PP) signals of molecular systems exhibiting internal conversion to the electronic ground state. We study the nonadiabatic dynamics of azomethane and the associating TA PP spectra at three levels of the electronic-structure theory, OM2/MRCI, SA-CASSCF, and XMS-CASPT2. The impact of these methods on the population dynamics and time-resolved TA PP signals is substantially different. This difference is attributed to the strong non-Condon effects that must be taken into account for the proper understanding and interpretation of time-resolved TA PP signals of nonadiabatic polyatomic systems. This shows that the combination of the dynamical and spectral simulations definitely provides more accurate and detailed information on the microscopic mechanisms of photophysical and photochemical processes. Hence the simulation of time-resolved spectroscopic signals provides another important dimension to examine the accuracy of quantum chemistry methods.

Spectral Fingerprint of Excited-State Energy Transfer in Dendrimers through Polarization-Sensitive Transient-Absorption Pump-Probe Signals: On-the-Fly Nonadiabatic Dynamics Simulations

The time-resolved polarization-sensitive transient-absorption (TA) pump-probe (PP) spectra are simulated using on-the-fly surface-hopping nonadiabatic dynamics and the doorway-window representation of nonlinear spectroscopy. A dendrimer model system composed of two linear phenylene ethynylene units (2-ring and 3-ring) is taken as an example. The ground-state bleach (GSB), stimulated emission (SE), and excited-state absorption (ESA) contributions as well as the total TA PP signals are obtained and carefully analyzed. It is shown that intramolecular excited-state energy transfer from the 2-ring unit to the 3-ring unit can be conveniently identified by employing pump and probe pulses with different polarizations. Our results demonstrate that time-resolved polarization-sensitive TA PP signals provide a powerful tool for the elucidation of excited-state energy-transfer pathways, notably in molecular systems possessing several optically bright nonadiabatically coupled electronic states with different orientations of transition dipole moments.

Ab Initio Surface-Hopping Simulation of Femtosecond Transient-Absorption Pump-Probe Signals of Nonadiabatic Excited-State Dynamics Using the Doorway-Window Representation

An ab initio theoretical framework for the simulation of femtosecond time-resolved transient absorption (TA) pump-probe (PP) spectra with quasi-classical trajectories is presented. The simulations are based on the classical approximation to the doorway-window (DW) representation of third-order four-wave-mixing signals. The DW formula accounts for the finite duration and spectral shape of the pump and probe pulses. In the classical DW formalism, classical trajectories are stochastically sampled from a positive definite doorway distribution, and the signals are evaluated by averaging over a positive definite window distribution. Nonadiabatic excited-state dynamics is described by a stochastic surface-hopping algorithm. The method has been implemented for the pyrazine molecule with the second-order algebraic-diagrammatic construction (ADC(2)) ab initio electronic-structure method. The methodology is illustrated by ab initio simulations of the ground-state bleach, stimulated emission, and excited-state absorption contributions to the TA PP spectrum of gas-phase pyrazine.

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.

Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials

Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.

On the Nature of Interplay among Major Flexibility Channels in Molecular Rotors

As a part of our interest in the excited-state dynamics of flexible materials, we have undertaken a theoretical investigation to the photo-induced reactions of 2-[4-(dimethylamino)benzylidene]malononitrile (BMN) by a combination of the density functional theory, its extended time-dependent (TD-DFT) single reference, and ab initio molecular dynamic (MD) simulations. The results showed that double-bond twisting and the neighbor single-bond twisting togetherness in the excited singlet state is the most important nonradiative deactivation channel to the ground state. Double- and single-bond twisting insert clear intersections among the potential energy surfaces of the singlet states (especially S1/S0) leading to fluorescence quenching. Furthermore, effects of molecular dynamic simulations on molecular properties in the femtosecond to picosecond time domain are studied to validate the results. In agreement with the experimental results, the findings conclude the existence of a flexible geometry-dependent single emission band. Such a study may give information on how the molecule could be externally modified/fixed to yield a desired effect, i.e., more fluorescence or more nonradiative decay.

Assessment of approximate computational methods for conical intersections and branching plane vectors in organic molecules

Quantum-chemical computational methods are benchmarked for their ability to describe conical intersections in a series of organic molecules and models of biological chromophores. Reference results for the geometries, relative energies, and branching planes of conical intersections are obtained using ab initio multireference configuration interaction with single and double excitations (MRCISD). They are compared with the results from more approximate methods, namely, the state-interaction state-averaged restricted ensemble-referenced Kohn-Sham method, spin-flip time-dependent density functional theory, and a semiempirical MRCISD approach using an orthogonalization-corrected model. It is demonstrated that these approximate methods reproduce the ab initio reference data very well, with root-mean-square deviations in the optimized geometries of the order of 0.1 Å or less and with reasonable agreement in the computed relative energies. A detailed analysis of the branching plane vectors shows that all currently applied methods yield similar nuclear displacements for escaping the strong non-adiabatic coupling region near the conical intersections. Our comparisons support the use of the tested quantum-chemical methods for modeling the photochemistry of large organic and biological systems.

Photoswitching of salicylidene methylamine: a theoretical photodynamics study

Photoswitching of simple photochromic molecules attracts substantial attention because of its possible role in future photon-driven molecular electronics. Here we model the full photoswitching cycle of a minimal photochromic Schiff base–salicylidene methylamine (SMA). We perform semiempirical nonadiabatic on-the-fly photodynamics simulations at the OM2/MRCI level and thoroughly analyze the structural time evolution and switching efficiency of the system. We also identify and examine in detail the crucial steps in the SMA photochemistry ruled by excited-state intramolecular proton transfer. The results place the investigated model aromatic Schiff base among the promising candidates for novel photoswitching molecular materials. Our study also shows the potential of the semiempirical multireference photodynamics simulations as a tool for early stage molecular photodevice design.

Adiabatic and nonadiabatic contributions to the energy of a system subject to a time-dependent perturbation: complete separation and physical interpretation

When a time-dependent perturbation acts on a quantum system that is initially in the nondegenerate ground state ∣0> of an unperturbed Hamiltonian H(0), the wave function acquires excited-state components ∣k> with coefficients c(k)(t) exp(-iE(k)t/ℏ), where E(k) denotes the energy of the unperturbed state ∣k>. It is well known that each coefficient c(k)(t) separates into an adiabatic term a(k)(t) that reflects the adjustment of the ground state to the perturbation--without actual transitions--and a nonadiabatic term b(k)(t) that yields the probability amplitude for a transition to the excited state. In this work, we prove that the energy at any time t also separates completely into adiabatic and nonadiabatic components, after accounting for the secular and normalization terms that appear in the solution of the time-dependent Schrödinger equation via Dirac's method of variation of constants. This result is derived explicitly through third order in the perturbation. We prove that the cross-terms between the adiabatic and nonadiabatic parts of c(k)(t) vanish, when the energy at time t is determined as an expectation value. The adiabatic term in the energy is identical to the total energy obtained from static perturbation theory, for a system exposed to the instantaneous perturbation λH'(t). The nonadiabatic term is a sum over excited states ∣k> of the transition probability multiplied by the transition energy. By evaluating the probabilities of transition to the excited eigenstates ∣k'(t)> of the instantaneous Hamiltonian H(t), we provide a physically transparent explanation of the result for E(t). To lowest order in the perturbation parameter λ, the probability of finding the system in state ∣k'(t)> is given by λ(2) ∣b(k)(t)∣(2). At third order, the transition probability depends on a second-order transition coefficient, derived in this work. We indicate expected differences between the results for transition probabilities obtained from this work and from Fermi's golden rule.