Photodynamics of Schiff Base Salicylideneaniline: Trajectory Surface-Hopping Simulations

Simulating the Excited-State Dynamics of Polaritons with Ab Initio Multiple Spawning

Over the past decade, there has been a growth of interest in polaritonic chemistry, where the formation of hybrid light-matter states (polaritons) can alter the course of photochemical reactions. These hybrid states are created by strong coupling between molecules and photons in resonant optical cavities and can even occur in the absence of light when the molecule is strongly coupled with the electromagnetic fluctuations of the vacuum field. We present a first-principles model to simulate nonadiabatic dynamics of such polaritonic states inside optical cavities by leveraging graphical processing units (GPUs). Our first implementation of this model is specialized for a single molecule coupled to a single-photon mode confined inside the optical cavity but with any number of excited states computed using complete active space configuration interaction (CASCI) and a Jaynes-Cummings-type Hamiltonian. Using this model, we have simulated the excited-state dynamics of a single salicylideneaniline (SA) molecule strongly coupled to a cavity photon with the ab initio multiple spawning (AIMS) method. We demonstrate how the branching ratios of the photodeactivation pathways for this molecule can be manipulated by coupling to the cavity. We also show how one can stop the photoreaction from happening inside of an optical cavity. Finally, we also investigate cavity-based control of the ordering of two excited states (one optically bright and the other optically dark) inside a cavity for a set of molecules, where the dark and bright states are close in energy.

Prediction of Photochromism of Salicylideneaniline Crystals Using a Data Mining Approach

Salicylideneanilines (SAs) are photochromic compounds that undergo enol-keto photoisomerization in the solid state. Research over the past 60 years has revealed empirically that SAs with steric and planar conformations tend to be photochromic and nonphotochromic, respectively. However, increasing counterexamples in the recent literature raise questions about the nature of the relationship between structure and photochromism in SA crystals and whether the photochromism of SA crystals is predictable. This study is the first to construct a data set on SA crystals and conduct a comprehensive analysis to investigate the relationship between molecular and crystal structures and photochromism. A data mining approach revealed that the dihedral angle is the most dominant structural parameter for photochromism, followed by the Hirshfeld surface volume. SAs with neutral bulky hydrocarbon groups, such as the tert-butyl group, tend to be photochromic because such SAs have steric conformation and a loosely packed structure. In contrast, SAs with fluorine, pyridine, and pyrazine are less likely to be photochromic due to their planar conformation and densely packed structures. The photochromism of the SA crystals in our data set was predicted with high accuracy (>85%) using machine learning. The results of this study provide a useful reference for designing SA crystals with desired photochromic properties.

Finding Excited-State Minimum Energy Crossing Points on a Budget: Non-Self-Consistent Tight-Binding Methods

The automated exploration and identification of minimum energy conical intersections (MECIs) is a valuable computational strategy for the study of photochemical processes. Due to the immense computational effort involved in calculating non-adiabatic derivative coupling vectors, simplifications have been introduced focusing instead on minimum energy crossing points (MECPs), where promising attempts were made with semiempirical quantum mechanical methods. A simplified treatment for describing crossing points between almost arbitrary diabatic states based on a non-self-consistent extended tight-binding method, GFN0-xTB, is presented. Involving only a single diagonalization of the Hamiltonian, the method can provide energies and gradients for multiple electronic states, which can be used in a derivative coupling-vector-free scheme to calculate MECPs. By comparison with high-lying MECIs of benchmark systems, it is demonstrated that the identified geometries are good starting points for further MECI refinement with ab initio methods.

Computational Approaches for Organic Semiconductors: From Chemical and Physical Understanding to Predicting New Materials

While a complete understanding of organic semiconductor (OSC) design principles remains elusive, computational methods─ranging from techniques based in classical and quantum mechanics to more recent data-enabled models─can complement experimental observations and provide deep physicochemical insights into OSC structure-processing-property relationships, offering new capabilities for in silico OSC discovery and design. In this Review, we trace the evolution of these computational methods and their application to OSCs, beginning with early quantum-chemical methods to investigate resonance in benzene and building to recent machine-learning (ML) techniques and their application to ever more sophisticated OSC scientific and engineering challenges. Along the way, we highlight the limitations of the methods and how sophisticated physical and mathematical frameworks have been created to overcome those limitations. We illustrate applications of these methods to a range of specific challenges in OSCs derived from π-conjugated polymers and molecules, including predicting charge-carrier transport, modeling chain conformations and bulk morphology, estimating thermomechanical properties, and describing phonons and thermal transport, to name a few. Through these examples, we demonstrate how advances in computational methods accelerate the deployment of OSCsin wide-ranging technologies, such as organic photovoltaics (OPVs), organic light-emitting diodes (OLEDs), organic thermoelectrics, organic batteries, and organic (bio)sensors. We conclude by providing an outlook for the future development of computational techniques to discover and assess the properties of high-performing OSCs with greater accuracy.

Excited-State Intramolecular Proton Transfer in Salicylidene-α-Hydroxy Carboxylate Derivatives: Direct Detection of the Triplet Excited State of the cis-Keto Tautomer

Excited-state intramolecular hydrogen transfer on the triplet surface of salicylideneaniline derivatives has received much less attention than the corresponding ultrafast process on the singlet surface. To enhance the understanding of this triplet reactivity, the photochemical properties of a series of salicylidene-α-hydroxy acid salts with different substituents on the phenol moiety (1-3) were characterized. UV/vis absorption and phosphorescence measurements in ethanol revealed that 1-3 exist as both enol and keto tautomers, with the enol form being predominant. Irradiation of 1 at 310 nm in ethanol glass (77 K) yielded an absorption band with a λ_max at ∼405 nm, which was assigned to the trans-keto tautomer (trans-1K). In contrast, laser flash photolysis of 1-3 in methanol or acetonitrile resulted in a transient absorption with λ_max at 440-460 nm. This transient, which decayed on the microsecond timescale and was significantly shorter lived in methanol than in acetonitrile, was assigned to the triplet excited state (T₁) of the cis-keto tautomer (cis-1K-3K) and residual absorption of trans-1K-3K by comparison with TD-DFT calculations. The assignment of the T₁ of cis-1K was further supported by quenching studies with anthracene and 2,5-dimethyl-2,4-hexadiene. Laser flash photolysis of 1 in the temperature range of 173-293 K gave an activation barrier of 6.7 kcal/mol for the decay of the T₁ of cis-1K. In contrast, the calculated activation barrier for cis-1K to undergo a 1,5-H atom shift to reform 1 was smaller, indicating that intersystem crossing of the T₁ of cis-1K is the rate-determining step in the regeneration of 1.

The vertical excitation energies and a lifetime of the two lowest singlet excited states of the conjugated polyenes from C2 to C22: Ab initio, DFT, and semiclassical MNDO-MD simulations

Electronic excited states in the series of polyene molecules were explored. Optimal ground-state geometry was used for the evaluation of vertical excitation energies. Results of a chosen set of functionals were compared to post-HF methods (EOM-CCSD, NEVPT2, CASPT2, and MRCI). In addition, the semiempirical OM2/MNDO method using MRCISD computational level was confronted with the above-mentioned techniques. Despite the fact that the first excited state has a significant double-excitation character some functionals were able to qualitatively determine the correct state order (where the lowest excited state has a A g - character). The most successful functionals in transition energies predictions were PBE, TPSS and BLYP in Tamm-Dancoff approach (TDA), which had the smallest root-mean-square deviation (RMSD) scoring towards the experimental values. Regarding RMSD scoring, the OM2/MNDO method performed fairly well, too. Besides absorption spectra, lifetimes of the first two excited states were estimated based on a stochastic approach exploring a swarm of OM2/MNDO hopping dynamics using the Tully fewest switch algorithm for each molecule. The longest lifetime of the first excited state (S₁ ) was found for decapentaene (about 5 ps). Further elongation of the conjugated chain caused a mild decrease of this value to ca 1.5 ps for docosaundecaene.

ESIPT in the pyrrol pyridine molecule: mechanism, timescale and yield revealed using dynamics simulations

Excited State Intramolecular Proton Transfer in pyrrol pyridine is theoretically investigated using non-adiabatic dynamics simulations. The photochemical process is completely characterised: the reaction time, the total yield and the accessibility of the conical intersection are evaluated. Finally, new mechanistic interpretation are extracted: the proton transfer reaction in this molecule is shown to be driven by two complementary mechanisms.

Wavelength-dependent photochemistry of a salicylimine derivative studied with cryogenic and ultrafast spectroscopy approaches

Salicylimines are versatile compounds in which an excited-state intramolecular proton transfer and torsional motions may set in upon photoexcitation. Here, we study N-(α-phenylethyl)salicylimine (PESA) to elucidate how the photochemical reaction pathways depend on the excitation wavelength and to what extent the relative photoproduct distribution can be steered towards a desired species. DFT structure and potential energy calculations disclose that the most stable ground-state conformer is an enol species and that the photodynamics may proceed differently depending on the excited state that is reached. With matrix isolation infrared spectroscopy, the predominance of the enol conformer of PESA is confirmed. Illumination of the cryogenic sample with different wavelengths shifts the ratio of enol and keto products, and by sequential irradiation a selective re- and depopulation is possible. Femtosecond transient absorption spectroscopy further reveals that also at room temperature, the outcome of the photoreaction depends on excitation wavelength, and in combination with the calculations, it can be rationalized that the decisive step occurs within the first hundred femtoseconds. Since the ultrafast dynamics mostly match those of similar salicylimines, our findings might also apply to those systems and provide additional insight into their reported sensitivity on excitation energy.

Interplay of conformational relaxation and hydrogen bond dynamics in the excited states of fluorescent Schiff base anions

Time resolved fluorescence spectroscopic investigation of four Schiff base anions has established that their excited state dynamics is governed by several solvent properties: polarity, viscosity and hydrogen bond donating ability. With viscous protic solvents like glycerol, fluorescence lifetimes of anions have been found to be markedly longer than those in ethanol, implying that conformational relaxation of molecules plays a key role in their nonradiative relaxation. Surprisingly, the lifetimes in less viscous aprotic solvents, like acetonitrile, are found to be even longer. The only plausible rationalization of this observation is in the light of hydrogen bond-assisted nonradiative phenomena that are operative in protic solvents. This contention draws support from a time evolution of the emission in the red end of the spectrum in low to moderately hydrogen bond donating protic solvents, with regard to an absence of such a rise time in aprotic solvents and strongly hydrogen bond donating solvents, viz., 2,2,2-trifluoroethanol. Rudimentary quantum chemical calculations provide a preliminary idea about the nature of excited state hydrogen bond redistribution involved in the process.

Fast Screening of Minimum Energy Crossing Points with Semiempirical Tight-Binding Methods

The investigation of photochemical processes is a highly active field in computational chemistry. One research direction is the automated exploration and identification of minimum energy conical intersection (MECI) geometries. However, due to the immense technical effort required to calculate nonadiabatic potential energy landscapes, the routine application of such computational protocols is severely limited. In this study, we will discuss the prospect of combining adiabatic potential energy surfaces from semiempirical quantum mechanical calculations with specialized confinement potential and metadynamics simulations to identify S₀/T₁ minimum energy crossing point (MECP) geometries. It is shown that MECPs calculated at the GFN2-xTB level can provide suitable approximations to high-level S₀/S₁ab initio conical intersection geometries at a fraction of the computational cost. Reference MECIs of benzene are studied to illustrate the basic concept. An example application of the presented protocol is demonstrated for a set of photoswitch molecules.

Exploring the Photophysics of a Zn2+ Fluorescence Sensor and Its Sensing Mechanism

Sensor X is a turn-on sensor, which is applied in the fluorescence detection of Zn²⁺ ions. Its photophysical process is comprehensively investigated to clarify its weak fluorescence. With the aid of density functional theory (DFT) and time-dependent density functional theory (TDDFT), the potential energy surfaces (PES) of X on both ground and first excited states are studied. Excited-state intramolecular proton transfer (EPT) processes as well as molecule twisting motion are observed, which induces several minima on the excited-state PES. Transition states as well as rate constants for these dynamic processes are obtained to evaluate their occurrences. The twisting motion of the sensor is an ultrafast process, which is initiated by a specific EPT process and leads to a nonemissive twisted intramolecular charge transfer (TICT) state. The fluorescence of the sensor is barely observable because of the easily attainable TICT state on the excited PES. This mechanism is trustworthy and intrinsically different from the previously proposed mechanism. After clarifying the photophysical process of the sensor, the Zn²⁺ sensing mechanism is uncovered. Also, the selectivity against Cd²⁺ and Hg²⁺ is fully discussed.

Quantum mechanics/molecular mechanics studies on the mechanistic photophysics of sunscreen oxybenzone in methanol solution

Herein, we have employed the QM(CASPT2//CASSCF)/MM method to explore the photophysical and photochemical mechanism of oxybenzone (OB) in methanol solution. Based on the optimized minima, conical intersections and crossing points, and minimum-energy reaction paths related to excited-state intramolecular proton transfer (ESIPT) and excited-state decay paths in the ¹ππ*, ¹nπ*, ³ππ*, ³nπ*, and S₀ states, we have identified several feasible excited-state relaxation pathways for the initially populated S₂(¹ππ*) state to decay to the initial enol isomer' S₀ state. The major one is the singlet-mediated and stretch-torsion coupled ESIPT pathway, in which the system first undergoes an essentially barrierless ¹ππ* ESIPT process to generate the ¹ππ* keto species, and finally realizes its ground state recovery through the subsequent carbonyl stretch-torsion facilitating S₁ → S₀ internal conversion (IC) and the reverse ground-state intramolecular proton transfer (GSIPT) process. The minor ones are related to intersystem crossing (ISC) processes. At the S₂(¹ππ*) minimum, an S₂(¹ππ*)/S₁(¹nπ*)/T₂(³nπ*) three-state intersection region helps the S₂ system branch into the T₁ state through a S₂ → S₁ → T₁ or S₂ → T₂ → T₁ process. Once it has reached the T₁ state, the system may relax to the S₀ state via direct ISC or via subsequent nearly barrierless ³ππ* ESIPT to yield the T₁ keto tautomer and ISC. The resultant S₀ keto species significantly undergoes reverse GSIPT and only a small fraction yields the trans-keto form that relaxes back more slowly. However, due to small spin-orbit couplings at T₁/S₀ crossing points, the ISC to S₀ state occurs very slowly. The present work rationalizes not only the ultrafast excited-state decay dynamics of OB but also its phosphorescence emission at low temperature.

The curious case of the photochemistry of 2-hydroxyphenylazo-3,5-dimethylisoxazole: Unravelling the process among tautomer-ization, photoisomerization, and conformational changes

Photoswitching in azo compounds is well established. Typically, the planar trans molecule (native) can undergo isomerization to cis isomer and vice versa in solution by light. However, observing such photochemistry...

Theoretical Studies on the Excited-State Decay Mechanism of Homomenthyl Salicylate in a Gas Phase and an Acetonitrile Solution

Here, we employ the CASPT2//CASSCF and QM(CASPT2//CASSCF)/MM approaches to explore the photochemical mechanism of homomenthyl salicylate (HMS) in vacuum and an acetonitrile solution. The results show that in both cases, the excited-state relaxation mainly involves a spectroscopically "bright" S₁(¹ππ*) state and the lower-lying T₁ and T₂ states. In the major relaxation pathway, the photoexcited S₁ keto system first undergoes an essentially barrierless excited-state intramolecular proton transfer (ESIPT) to generate the S₁ enol minimum, near which a favorable S₁/S₀ conical intersection decays the system to the S₀ state followed by a reverse ground-state intramolecular proton transfer (GSIPT) to repopulate the initial S₀ keto species. In the minor one, an S₁/T₂/T₁ three-state intersection in the keto region makes the T₁ state populated via direct and T₂-mediated intersystem crossing (ISC) processes. In the T₁ state, an ESIPT occurs, which is followed by ISC near a T₁/S₀ crossing point in the enol region to the S₀ state and finally back to the S₀ keto species. In addition, a T₁/S₀ crossing point near the T₁ keto minimum can also help the system decay to the S₀ keto species. However, small spin-orbit couplings between T₁ and S₀ at these T₁/S₀ crossing points make ISC to the S₀ state very slow and make the system trapped in the T₁ state for a while. The present work rationalizes not only the ultrafast excited-state decay dynamics of HMS but also its low quantum yield of phosphorescence at 77 K.

From fluorogens to fluorophores by elucidation and suppression of ultrafast excited state processes of a Schiff base

Strategies have been explored for developing strongly fluorescent species out of a weakly fluorescent Schiff base, 2-(((pyridin-2-ylmethyl)imino)methyl)phenol (salampy). The locally excited enolic state of salampy undergoes an intramolecular proton transfer with a time constant of ca. 200 fs. The emissive cis keto state thus formed decays completely within 50 ps. Its fast decay and miniscule fluorescence quantum yield are attributed to efficient non-radiative channels associated with conformational relaxation. The anionic form, salampy^-, has a significantly longer fluorescence lifetime of 800 ps. Its emissive state evolves in tens of picoseconds, from the locally excited state, by solvent and conformational relaxation. Both the neutral and anionic forms have a fluorescence lifetime of about 6 ns at 77 K, a temperature at which all activated nonradiative channels are blocked. This lifetime is similar to that obtained at room temperature, upon rigidification of the anion by complexation with Zn²⁺. Two such complexes have been studied. The first is binuclear, with acetate bridge between the two Zn²⁺ ions. The second, with ClO₄^- as the counterion, is mononuclear with two salampy ligands ligating the metal ion. Unlike a previous report on a different Schiff base, in which the ligands are π-stacked in its dimeric Zn²⁺ complex, no additional nonradiative deactivation pathway opens up in the Zn complexes of salampy, which are devoid of such stacking. The complex of salampy with Al³⁺ has an even longer fluorescence lifetime of 9 ns, indicating a greater degree of rigidification and consequent suppression of nonradiative processes.

Nonadiabatic Dynamics Simulations on Early-Time Photochemistry of Spirobenzopyran

Photoinduced ring-opening, decay, and isomerization of spirobenzopyran have been explored by the OM2/MRCI nonadiabatic dynamics simulations based on Tully's fewest-switches surface hopping scheme. The efficient S₁ to S₀ internal conversion as observed in experiments is attributed to the existence of two efficient excited-state decay pathways. The first one is related to the C-N dissociation, and the second one is done to the C-O dissociation. The C-O dissociation pathway is dominant, and more than 90% trajectories decay to the S₀ state via the C-O bond-fission related S₁/S₀ conical intersections. Near these regions in the S₀ state, trajectories can either return to spirobenzopyran or proceed to various intermediates including merocyanine via a series of bond rotations. Our nonadiabatic dynamics simulations also demonstrate that the hydrogen-out-of-plane (HOOP) motion is important for efficient and ultrafast excited-state deactivation. On the other hand, we have also found that the replacement of methyl groups by hydrogen atoms in spirobenzopyran can artificially introduce different intramolecular hydrogen transfers leading to hydrogen-transferred intermediates. This finding is important for the community and demonstrates that such a kind of structural truncation, sometimes, could be problematic, leading to incorrect photodynamics. Our present work provides valuable insights into the photodynamics of spirobenzopyran, which could be helpful for the design of spiropyran-based photochromic materials.

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.

Experimental and theoretical investigation of ground state intramolecular proton transfer (GSIPT) in salicylideneaniline Schiff base derivatives in polar protic medium

Ground state intramolecular proton transfer process has been comprehensively investigated in three salicylideneaniline Schiff base derivatives (SB1, SB2, and SB3) using experimental and theoretical methods. It has been confirmed that all the three Schiff base molecules in the ground electronic state exist in the enol form in non-polar and polar aprotic solvents. Keto form is being populated by the polar protic solvent through ground state intramolecular proton transfer (GSIPT) process. Ground state equilibrium between the enol and keto tautomers for SB1 and SB3 is mainly governed by the proton donating ability of the solvent. Ground state equilibria between the enol and keto tautomers of SB2 which is a positional isomer of SB3 is governed by the polarity and proton donating ability of the solvents. Excited state intramolecular proton transfer (ESIPT) process is also evidenced in all the three Schiff base molecules. Theoretical calculations at the B3LYP/cc-pVDZ level in the gas phase and in different solvents using polarisable continuum model (PCM) have failed to establish the GSIPT process. Microsolvation of individual enol and keto conformers has been investigated considering upto three solvent molecules. The energetics of the individual conformers together with the corresponding transition state have been calculated. It has been confirmed that the keto conformer is more stable compared to the enol conformer in microsolvated cluster of three methanol molecules. Lowering of activation energy for the enol to keto tautomerisation in the presence of methanol also supports the experimental observation for GSIPT process. TDDFT/B3LYP/cc-pVDZ single point calculations for microsolvated clusters of enol and keto form of the Schiff base molecules exhibit an excellent agreement with the experimentally obtained absorption spectra. Difference in spectral nature of the Schiff base molecules has been explained using natural bond orbital (NBO) analysis. Quantum theory of atoms in molecules (QTAIM) has also been utilised to understand the GSIPT process in detail.

Photophysical transformations induced by chemical substitution to salicylaldimines

The role of electron-acceptor strength and microenvironment polarity on the photophysical properties of salicylaldimines.

A Spectroscopic Study of Tautomeric Equilibrium of Salicylideneaniline in ZSM-5 Zeolites

Salicylideneaniline (SA) sorbed in cation-exchanged M-ZSM-5 (M = H⁺, Li⁺, Na⁺, K⁺, Rb⁺, Cs⁺ and Zn²⁺) zeolites was studied by spectroscopic techniques assisted by quantum-chemical calculations. The nature of extra-framework cations present in the zeolite void was found to affect the spectral signature of the sorbed SA molecule that points to the shift of tautomeric equilibrium between the enol and keto forms. Small size cations, such as H⁺ and Li⁺, stabilize a cis-keto SA tautomer along with a enol one in the zeolite structure. The calculations indicate that the sorbed cis-keto tautomer may have the dipole large enough to be considered as a zwitterion. New features appearing in the spectra with the increase of the cation size were attributed to the presence of trans-keto SA tautomer, which up to now has been observed only in time-resolved spectroscopic experiments. A strong interaction of the molecule with cations in Zn-ZSM-5 zeolite results in the chelation of enol SA with the divalent Zn²⁺ ions. The results of the study suggest that the tautomeric equilibrium of molecules belonging to the Schiff base family can be tuned by the confinement in the nanoporous materials via a choice of topology of zeolite framework and the nature of extra-framework cations.

GFN2-xTB-An Accurate and Broadly Parametrized Self-Consistent Tight-Binding Quantum Chemical Method with Multipole Electrostatics and Density-Dependent Dispersion Contributions

An extended semiempirical tight-binding model is presented, which is primarily designed for the fast calculation of structures and noncovalent interaction energies for molecular systems with roughly 1000 atoms. The essential novelty in this so-called GFN2-xTB method is the inclusion of anisotropic second order density fluctuation effects via short-range damped interactions of cumulative atomic multipole moments. Without noticeable increase in the computational demands, this results in a less empirical and overall more physically sound method, which does not require any classical halogen or hydrogen bonding corrections and which relies solely on global and element-specific parameters (available up to radon, Z = 86). Moreover, the atomic partial charge dependent D4 London dispersion model is incorporated self-consistently, which can be naturally obtained in a tight-binding picture from second order density fluctuations. Fully analytical and numerically precise gradients (nuclear forces) are implemented. The accuracy of the method is benchmarked for a wide variety of systems and compared with other semiempirical methods. Along with excellent performance for the "target" properties, we also find lower errors for "off-target" properties such as barrier heights and molecular dipole moments. High computational efficiency along with the improved physics compared to its precursor GFN-xTB makes this method well-suited to explore the conformational space of molecular systems. Significant improvements are furthermore observed for various benchmark sets, which are prototypical for biomolecular systems in aqueous solution.

Unravelling the solvent polarity effect on the excited state intramolecular proton transfer mechanism of the 1- and 2-salicylideneanthrylamine. A TD-DFT case study

TD-DFT has been used to investigate the relaxation pathways of 1- and 2-salicylideneanthrylamine in different solvent environments.

A theoretical study of the photodynamics of salicylidene-2-anthrylamine in acetonitrile solution

The ultrafast photoinduced processes of salicylidene-2-anthrylamine (2-AntSA) in acetonitrile solution have been investigated using DFT/TD-DFT static electronic structure calculations and excited state ab initio molecular dynamics simulations. Two different isoenergetic enol ground-state structures with suitable geometry for excited state intramolecular proton transfer (ESIPT) where chosen for the excited-state dynamics. The S1 relaxed potential energy profiles for the excited state intramolecular proton transfer and the N[double bond, length as m-dash]C double bond isomerization reactions predict that both reactions occur over an energy barrier and that they are competitive processes in the deactivation of the Franck-Condon state. The photodynamic simulations show that the ESIPT occurs in the femtosecond time scale for both conformers (77 and 213 fs for IA and IIA, respectively) and that the speed is modulated by the ability of the conformers to evolve toward a planar conformation in the S1 state. The trajectories predict two conical intersections which provide nonradiative relaxation pathways to the S0 state. The first one is located in the twisted enol region, where the proton transfer process is unlikely, and only occurs for the conformer IIA in a time scale ≥600 fs. The second conical intersection is located in the cis-keto region, and represents an effective depopulation channel toward the trans-keto form. All our results are in remarkably good agreement with experiments.

Facile Synthesis of an Organic Solid State Near-Infrared-Emitter with Large Stokes Shift via Excited-State Intramolecular Proton Transfer

An organic solid-state near-infrared (NIR)-emitter (λ_em = 738 nm) exhibiting large Stokes shift (Δλ = 293 nm) through the excited-state intramolecular proton transfer phenomenon has been synthesized and characterized. The present discovery points to the possibility of achieving a new family of solid-state NIR emitters starting from simple aldehyde and amine precursors.

Ultrafast action chemistry in slow motion: atomistic description of the excitation and fluorescence processes in an archetypal fluorescent protein

We report quantum mechanical/molecular mechanical non-adiabatic molecular dynamics simulations on the electronically excited state of green fluorescent protein mutant S65T/H148D. We examine the driving force of the ultrafast (τ < 50 fs) excited-state proton transfer unleashed by absorption in the A band at 415 nm and propose an atomistic description of the two dynamical regimes experimentally observed [Stoner Ma et al., J. Am. Chem. Soc., 2008, 130, 1227]. These regimes are explained in terms of two sets of successive dynamical events: first the proton transfers quickly from the chromophore to the acceptor Asp148. Thereafter, on a slower time scale, there are geometrical changes in the cavity of the chromophore that involve the distance between the chromophore and Asp148, the planarity of the excited-state chromophore, and the distance between the chromophore and Tyr145. We find two different non-radiative relaxation channels that are operative for structures in the reactant region and that can explain the mismatch between the decay of the emission of A* and the rise of the emission of I*, as well as the temperature dependence of the non-radiative decay rate.

Nuclear Ensemble Approach with Importance Sampling

We show that the importance sampling technique can effectively augment the range of problems where the nuclear ensemble approach can be applied. A sampling probability distribution function initially determines the collection of initial conditions for which calculations are performed, as usual. Then, results for a distinct target distribution are computed by introducing compensating importance sampling weights for each sampled point. This mapping between the two probability distributions can be performed whenever they are both explicitly constructed. Perhaps most notably, this procedure allows for the computation of temperature dependent observables. As a test case, we investigated the UV absorption spectra of phenol, which has been shown to have a marked temperature dependence. Application of the proposed technique to a range that covers 500 K provides results that converge to those obtained with conventional sampling. We further show that an overall improved rate of convergence is obtained when sampling is performed at intermediate temperatures. The comparison between calculated and the available measured cross sections is very satisfactory, as the main features of the spectra are correctly reproduced. As a second test case, one of Tully's classical models was revisited, and we show that the computation of dynamical observables also profits from the importance sampling technique. In summary, the strategy developed here can be employed to assess the role of temperature for any property calculated within the nuclear ensemble method, with the same computational cost as doing so for a single temperature.

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.

An Ab initio study on the photophysics of tris(salicylideneaniline)

Sequential excited-state proton transfers result in multiple band fluorescence spectrum.

Photochromic Mechanism of a Bridged Diarylethene: Combined Electronic Structure Calculations and Nonadiabatic Dynamics Simulations

Intramolecularly bridged diarylethenes exhibit improved photocyclization quantum yields because the anti-syn isomerization that originally suppresses photocyclization in classical diarylethenes is blocked. Experimentally, three possible channels have been proposed to interpret experimental observation, but many details of photochromic mechanism remain ambiguous. In this work we have employed a series of electronic structure methods (OM2/MRCI, DFT, TDDFT, RI-CC2, DFT/MRCI, and CASPT2) to comprehensively study excited state properties, photocyclization, and photoreversion dynamics of 1,2-dicyano[2,2]metacyclophan-1-ene. On the basis of optimized stationary points and minimum-energy conical intersections, we have refined experimentally proposed photochromic mechanism. Only an S₁/S₀ minimum-energy conical intersection is located; thus, we can exclude the third channel experimentally proposed. In addition, we find that both photocyclization and photoreversion processes use the same S₁/S₀ conical intersection to decay the S₁ system to the S₀ state, so we can unify the remaining two channels into one. These new insights are verified by our OM2/MRCI nonadiabatic dynamics simulations. The S₁ excited-state lifetimes of photocyclization and photoreversion are estimated to be 349 and 453 fs, respectively, which are close to experimentally measured values: 240 ± 60 and 250 fs in acetonitrile solution. The present study not only interprets experimental observations and refines previously proposed mechanism but also provides new physical insights that are valuable for future experiments.

Tailoring the Schiff base photoswitching – a non-adiabatic molecular dynamics study of substituent effect on excited state proton transfer

Dynamics reveals how to design chemical substitutions to control excited-state proton transfer efficiency.

Excited-State Intramolecular Proton-Transfer Properties of Three Tris(N-Salicylideneaniline)-Based Chromophores with Extended Conjugation

The synthesis and photophysical properties of three tris(N-salicylideneaniline) (TSA) compounds containing 1,3,5-triarylbenzene, -tristyrylbenzene, and -tris(arylethynyl)benzene core units are reported. The TSA compounds underwent efficient excited-state intramolecular proton transfer (ESIPT) in solution and in solid state due to the preformed C=N⋅⋅⋅H-O hydrogen-bonded motifs of the structures. Steady-state fluorescence emission spectra of the TSA molecules revealed dual bands only in DMSO, and large Stokes shifts in other polar aprotic and protic solvents. Femtosecond transient absorption spectroscopic measurements in THF revealed lifetime values in the range of 14-16 ps for the excited-state keto-tautomer. The TSA compounds are also responsive to metal ions (Cu²⁺ and Zn²⁺ ) in DMSO, exhibit enhanced aggregate-induced emission (AIE) properties in DMSO/water mixtures, and are highly luminescent in the solid state.

Intersystem crossing-branched excited-state intramolecular proton transfer for o-nitrophenol: An ab initio on-the-fly nonadiabatic molecular dynamic simulation

The 6SA-CASSCF(10, 10)/6-31G (d, p) quantum chemistry method has been applied to perform on-the-fly trajectory surface hopping simulation with global switching algorithm and to explore excited-state intramolecular proton transfer reactions for the o-nitrophenol molecule within low-lying electronic singlet states (S0 and S1) and triplet states (T1 and T2). The decisive photoisomerization mechanisms of o-nitrophenol upon S1 excitation are found by three intersystem crossings and one conical intersection between two triplet states, in which T1 state plays an essential role. The present simulation shows branch ratios and timescales of three key processes via T1 state, non-hydrogen transfer with ratio 48% and timescale 300 fs, the tunneling hydrogen transfer with ratios 36% and timescale 10 ps, and the direct hydrogen transfer with ratios 13% and timescale 40 fs. The present simulated timescales might be close to low limit of the recent experiment results.