The Azulene S1 State Decays via a Conical Intersection: A CASSCF Study with MMVB Dynamics

Iminothioindoxyl Donors with Exceptionally High Cross Section for Protein Vibrational Energy Transfer

Various protein functions are related to vibrational energy transfer (VET) as an important mechanism. The underlying transfer pathways can be experimentally followed by ultrafast Vis-pump/IR-probe spectroscopy with a donor-sensor pair of non-canonical amino acids (ncAAs) incorporated in a protein. However, so far only one donor ncAA, azulenylalanine (AzAla), exists, which suffers from a comparably low Vis extinction coefficient. Here, we introduce two novel donor ncAAs based on an iminothioindoxyl (ITI) chromophore. The dimethylamino-ITI (DMA-ITI) and julolidine-ITI (J-ITI) moieties overcome the limitation of AzAla with a 50 times higher Vis extinction coefficient. While ITI moieties are known for ultrafast photoswitching, DMA-ITI and J-ITI exclusively form a hot ground state on the sub-ps timescale instead, which is essential for their usage as vibrational energy donor. In VET measurements of donor-sensor dipeptides we investigate the performance of the new donors. We observe 20 times larger signals compared to the established AzAla donor, which opens unprecedented possibilities for the study of VET in proteins.

A framework for multiexcitonic logic

Exciton science sits at the intersection of chemical, optical and spin-based implementations of information processing, but using excitons to conduct logical operations remains relatively unexplored. Excitons encoding information could be read optically (photoexcitation-photoemission) or electrically (charge recombination-separation), travel through materials via exciton energy transfer, and interact with one another in stimuli-responsive molecular excitonic devices. Excitonic logic offers the potential to mediate electrical, optical and chemical information. Additionally, high-spin triplet and quintet (multi)excitons offer access to well defined spin states of relevance to magnetic field effects, classical spintronics and spin-based quantum information science. In this Roadmap, we propose a framework for developing excitonic computing based on singlet fission (SF) and triplet-triplet annihilation (TTA). Various molecular components capable of modulating SF/TTA for logical operations are suggested, including molecular photo-switching and multi-colour photoexcitation. We then outline a pathway for constructing excitonic logic devices, considering aspects of circuit assembly, logical operation synchronization, and exciton transport and amplification. Promising future directions and challenges are identified, and the potential for realizing excitonic computing in the near future is discussed.

Excited-State (Anti)Aromaticity Explains Why Azulene Disobeys Kasha's Rule

Fluorescence exclusively occurs from the lowest excited state of a given multiplicity according to Kasha's rule. However, this rule is not obeyed by a handful of anti-Kasha fluorophores whose underlying mechanism is still understood merely on a phenomenological basis. This lack of understanding prevents the rational design and property-tuning of anti-Kasha fluorophores. Here, we propose a model explaining the photophysical properties of an archetypal anti-Kasha fluorophore, azulene, based on its ground- and excited-state (anti)aromaticity. We derived our model from a detailed analysis of the electronic structure of the ground singlet, first excited triplet, and quintet states and of the first and second excited singlet states using the perturbational molecular orbital theory and quantum-chemical aromaticity indices. Our model reveals that the anti-Kasha properties of azulene and its derivatives result from (i) the contrasting (anti)aromaticity of its first and second singlet excited states (S1 and S2, respectively) and (ii) an easily accessible antiaromaticity relief pathway of the S1 state. This explanation of the fundamental cause of anti-Kasha behavior may pave the way for new classes of anti-Kasha fluorophores and materials with long-lived, high-energy excited states.

Quantifying the Relaxation Dynamics of Higher Electronic Excited States in Perylene

Gating logical operations through high-lying electronic excited states presents opportunities for developing ultrafast, subnanometer computational devices. A lack of molecular systems with sufficiently long-lived higher excited states has hindered practical realization of such devices, but recent studies have reported intriguing photophysics from high-lying excited states of perylene. In this work, we use femtosecond spectroscopy supported by quantum chemical calculations to identify and quantify the relaxation dynamics of monomeric perylene's higher electronic excited states. The 21B2u state is accessed through single-photon absorption at 250 nm, while the optically dark 21Ag state is excited via the 11B3u state. Population of either state results in subpicosecond relaxation to the 11B3u state, and we quantify 21Ag and 21B2u state lifetimes of 340 and 530 fs, respectively. These lifetimes are significantly longer than the singlet fission time constant from the perylene 21B2u state, suggesting that the higher electronic states of perylene may be useful for gating logical operations.

Synthesis and Properties of Carbo- and Heterocyclic Benz[a]azulenes

A new and convenient synthesis of aryl-substituted naphtho[2,1-a]azulenes by the combination of Suzuki-Miyaura, Sonogashira, and cycloisomerization reactions is reported. The methodology was applied to the synthesis of hitherto unknown azuleno[1,2-h]quinolines, cyclohepta[1,2]indeno[4,5-b]thiophenes, and cyclohepta[1,2]indeno[4,5-c]thiophenes. The impact of different fused-heterocyclic rings on the photophysical and electrochemical properties of these azulene derivatives was studied by experimental and theoretical methods and hence provides a rationale for the preparation of novel azulene derivatives with improved properties for application as organic materials.

Photophysics of uracil: an explicit time-dependent generating function-based method combining both nonadiabatic and spin-orbit coupling effects

We present a composite framework for calculating the rates of non-radiative deactivation processes, namely internal conversion (IC) and intersystem crossing (ISC), on an equal footing by explicitly computing the non-adiabatic coupling (NAC) and spin-orbit coupling (SOC) constants, respectively. The stationary-state approach uses a time-dependent generating function based on Fermi's golden rule. We validate the applicability of the framework by computing the rate of IC for azulene, obtaining comparable rates to experimental and previous theoretical results. Next, we investigate the photophysics associated with the complex photodynamics of the uracil molecule. Interestingly, our simulated rates corroborate experimental observations. Detailed analyses using Duschinsky rotation matrices, displacement vectors and NAC matrix elements are presented to interpret the findings alongside testing the suitability of the approach for such molecular systems. The suitability of the Fermi's golden rule based method is explained qualitatively in terms of single-mode potential energy surfaces.

From a one-mode to a multi-mode understanding of conical intersection mediated ultrafast organic photochemical reactions

This review discusses how ultrafast organic photochemical reactions are controlled by conical intersections, highlighting that decay to the ground-state at multiple points of the intersection space results in their multi-mode character.

Rules for Magnetic Exchange in Azulene-Bridged Biradicals: Quo Vadis?

Electronic coupling through organic bridges facilitates magnetic exchange interactions and controls electron transfer and single-molecule device electron transport. Electronic coupling through alternant π-systems (e.g., benzene) is better understood than the corresponding coupling through nonalternant π-systems (e.g., azulene). Herein, we examine the structure, spectroscopy, and magnetic exchange coupling in two biradicals (1,3-SQ₂Az and 1,3-SQ-Az-NN; SQ = the zinc(II) complex of spin-1/2 semiquinone radical anion, NN = spin-1/2 nitronylnitroxide; Az = azulene) that possess nonalternant azulene π-system bridges. The SQ radical spin density in both molecules is delocalized into the Az π-system, while the NN spin is effectively localized onto the five-atom ONCNO π-system of NN radical. The spin distributions and interactions are probed by EPR spectroscopy and magnetic susceptibility measurements. We find that J = +38 cm^-1 for 1,3-SQ₂Az and J = +9 cm^-1 for 1,3-SQ-Az-NN (H=-2JS^SQ·S^SQorNN). Our results highlight the differences in exchange coupling mediated by azulene compared to exchange coupling mediated by alternant π-systems.

Azulene as an ingredient for visible-light- and stimuli-responsive photoswitches

The azulene molecule features a unique combination of optical, luminescence, and stimuli-responsive properties. This makes the azulene motif a promising functional group to be introduced in photoswitches. Recent challenges in the field of photochromic compounds require the development of new approaches to molecules that are switched by visible light (400-760 nm), are proton responsive and have advanced luminescent properties. Merging azulene with photoswitches opens prospects for fulfilling these requirements. Herein, we highlight recent results on the application of this hydrocarbon motif in various photochromic systems, such as stilbenes, diarylethenes, and azobenzenes.

A systematic model study quantifying how conical intersection topography modulates photochemical reactions

Despite their important role in photochemistry and expected presence in most polyatomic molecules, conical intersections have been thoroughly characterized in a comparatively small number of systems. Conical intersections can confer molecular photoreactivity or photostability, often with remarkable efficacy, due to their unique structure: at a conical intersection, the adiabatic potential energy surfaces of two or more electronic states are degenerate, enabling ultrafast decay from an excited state without radiative emission, known as nonadiabatic transfer. Furthermore, the precise conical intersection topography determines fundamental properties of photochemical processes, including excited-state decay rate, efficacy, and molecular products that are formed. However, these relationships have yet to be defined comprehensively. In this article, we use an adaptable computational model to investigate a variety of conical intersection topographies, simulate resulting nonadiabatic dynamics, and calculate key photochemical observables. We varied the vibrational mode frequencies to modify conical intersection topography systematically in four primary classes of conical intersections and quantified the resulting rate, total yield, and product yield of nonadiabatic decay. The results reveal that higher vibrational mode frequencies reduce nonadiabatic transfer, but increase the transfer rate and resulting photoproduct formation. These trends can inform progress toward experimental control of photochemical reactions or tuning of molecules' photochemical properties based on conical intersections and their topography.

Cyclization from Higher Excited States of Diarylethenes Having a Substituted Azulene Ring

The cyclization reaction of diarylethenes having an azulene ring occurs only via higher excited states. Novel diarylethenes having an azulene ring with a strong donor or acceptor were synthesized and examined in these reactions. A derivative having an electron-donating 1,3-benzodithiol-2-ylidenemethyl group at the 1-position of the azulene ring showed photochromism, whereas neither a derivative having a π-conjugated electron-donating group at the 3-position of the azulene ring nor derivatives having a π-conjugated electron-withdrawing group at the 1- or 3-position of the azulene ring showed any photochromism. The photoreactivities of these compounds were explained by calculating forces and bond orders on the excited states using density functional theory (DFT) and time-dependent (TD)-DFT.

Semiclassical Approach to Photophysics Beyond Kasha's Rule and Vibronic Spectroscopy Beyond the Condon Approximation. The Case of Azulene

Azulene is a prototypical molecule with an anomalous fluorescence from the second excited electronic state, thus violating Kasha's rule, and with an emission spectrum that cannot be understood within the Condon approximation. To better understand the photophysics and spectroscopy of azulene and other nonconventional molecules, we developed a systematic, general, and efficient computational approach combining the semiclassical dynamics of nuclei with ab initio electronic structure. First, to analyze the nonadiabatic effects, we complement the standard population dynamics by a rigorous measure of adiabaticity, estimated with the multiple-surface dephasing representation. Second, we propose a new semiclassical method for simulating non-Condon spectra, which combines the extended thawed Gaussian approximation with the efficient single-Hessian approach. S₁ ← S₀ and S₂ ← S₀ absorption and S₂ → S₀ emission spectra of azulene, recorded in a new set of experiments, agree very well with our calculations. We find that accuracy of the evaluated spectra requires the treatment of anharmonicity, Herzberg-Teller, and mode-mixing effects.

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.

Azulene-Derived Fluorescent Probe for Bioimaging: Detection of Reactive Oxygen and Nitrogen Species by Two-Photon Microscopy

Two-photon fluorescence microscopy has become an indispensable technique for cellular imaging. Whereas most two-photon fluorescent probes rely on well-known fluorophores, here we report a new fluorophore for bioimaging, namely azulene. A chemodosimeter, comprising a boronate ester receptor motif conjugated to an appropriately substituted azulene, is shown to be an effective two-photon fluorescent probe for reactive oxygen species, showing good cell penetration, high selectivity for peroxynitrite, no cytotoxicity, and excellent photostability.

Interactions of benzene, naphthalene, and azulene with alkali-metal and alkaline-earth-metal atoms for ultracold studies

We consider collisional properties of polyatomic aromatic hydrocarbon molecules immersed into ultracold atomic gases and investigate intermolecular interactions of exemplary benzene, naphthalene, and azulene with alkali-metal (Li, Na, K, Rb, and Cs) and alkaline-earth-metal (Mg, Ca, Sr, and Ba) atoms. We apply the state-of-the-art ab initio techniques to compute the potential energy surfaces (PESs). We use the coupled cluster method restricted to single, double, and noniterative triple excitations to reproduce the correlation energy and the small-core energy-consistent pseudopotentials to model the scalar relativistic effects in heavier metal atoms. We also report the leading long-range isotropic and anisotropic dispersion and induction interaction coefficients. The PESs are characterized in detail, and the nature of intermolecular interactions is analyzed and benchmarked using symmetry-adapted perturbation theory. The full three-dimensional PESs are provided for the selected systems within the atom-bond pairwise additive representation and can be employed in scattering calculations. The present study of the electronic structure is the first step toward the evaluation of prospects for sympathetic cooling of polyatomic aromatic molecules with ultracold atoms. We suggest azulene, an isomer of naphthalene which possesses a significant permanent electric dipole moment and optical transitions in the visible range, as a promising candidate for electric field manipulation and buffer-gas or sympathetic cooling.

Ultrafast Vibrational Cooling Inside of a Molecular Container

Vibrational cooling of azulene encapsulated in a hemicarcerand molecular container was studied by pump-probe spectroscopy. Within 1.5 ps of excitation of azulene to the S₁ state, rapid internal conversion through a conical intersection leads to the formation of a vibrationally hot (∼1080 K) ground state, the subsequent cooling of which can be monitored by tracking the evolution of the red-shifted hot band at the edge of the ground-state absorption. It was found that the cooling of the hot S₀ state of azulene in the host-guest complex (hemicarceplex) is 2-4 times faster than that in common organic solvents. Such large acceleration points to a high density of matching vibrational modes and efficient mechanical coupling between the guest and the host. The experimental observations were fully corroborated by the results of molecular dynamics simulations.

An efficient approximate algorithm for nonadiabatic molecular dynamics

We propose a modification to the nonadiabatic surface hopping calculation method formulated in a paper by Yu et al. [Phys. Chem. Chem. Phys. 16, 25883 (2014)], which is a multidimensional extension of the Zhu-Nakamura theory with a practical diabatic gradient estimation algorithm. In our modification, their diabatic gradient estimation algorithm, which is based on a simple interpolation of the adiabatic potential energy surfaces, is replaced by an algorithm using the numerical derivatives of the adiabatic gradients. We then apply the algorithm to several models of nonadiabatic dynamics, both analytic and ab initio models, to numerically demonstrate that our method indeed widens the applicability and robustness of their method. We also discuss the validity and limitations of our new nonadiabatic surface hopping method while considering in mind potential applications to excited-state dynamics of biomolecules or unconventional nonadiabatic dynamics such as radiation decay processes in ultraintense X-ray fields.

Excited-state dynamics of a molecular dyad with two orthogonally-oriented fluorophores

The excited-state dynamics of a T-shaped bichromophoric molecule, consisting of two strong fluorophores, diphenyloxazole and diphenylpyrazoline, directly linked in an orthogonal geometry, was investigated. Despite the weak coupling ensured by this geometry and confirmed by the electronic absorption spectra, this dyad exhibits only weak fluorescence in both apolar and polar solvents, with fluorescence lifetimes ranging from 200 ps in CHX to 10 ps in ACN. Ultrafast spectroscopic measurements reveal that the fluorescence quenching in polar solvents is due to the population of a charge-separated state. In non-polar solvents, this process is energetically not feasible, and a quenching due to an efficient intersystem crossing (ISC) to the triplet manifold is proposed, based on quantum-chemical calculations. This process occurs via the spin-orbit charge-transfer (SOCT) ISC mechanism, which is enabled by the charge-transfer character acquired by the S1 state of the dyad upon structural relaxation and by the orthogonal arrangement of the molecular orbitals involved in the transition. The same mechanism is proposed to explain why the recombination of the charge-separated state is faster in medium than in highly polar solvents, as well as to account for the fast decay of the lowest triplet state to the ground state.

Manipulating azobenzene photoisomerization through strong light-molecule coupling

The formation of hybrid light–molecule states (polaritons) offers a new strategy to manipulate the photochemistry of molecules. To fully exploit its potential, one needs to build a toolbox of polaritonic phenomenologies that supplement those of standard photochemistry. By means of a state-of-the-art computational photochemistry approach extended to the strong-coupling regime, here we disclose various mechanisms peculiar of polaritonic chemistry: coherent population oscillations between polaritons, quenching by trapping in dead-end polaritonic states and the alteration of the photochemical reaction pathway and quantum yields. We focus on azobenzene photoisomerization, that encompasses the essential features of complex photochemical reactions such as the presence of conical intersections and reaction coordinates involving multiple internal modes. In the strong coupling regime, a polaritonic conical intersection arises and we characterize its role in the photochemical process. Our chemically detailed simulations provide a framework to rationalize how the strong coupling impacts the photochemistry of realistic molecules.

Manipulation of the photochemistry of molecules is traditionally achieved through synthetic chemical modifications. Here the authors use computational photochemistry to show how to control azobenzene photoisomerization through hybrid light–molecule states (polaritons).

Theoretical Studies of Singlet Fission: Searching for Materials and Exploring Mechanisms

In this Review article, a survey is given for theoretical studies in the subject of singlet fission. Singlet fission converts one singlet exciton to two triplet excitons. With the doubled number of excitons and the longer lifetime of the triplets, singlet fission provides an avenue to improve the photoelectric conversion efficiency in organic photovoltaic devices. It has been a subject of intense research in the past decade. Theoretical studies play an essential role in understanding singlet fission. This article presents a Review of theoretical studies in singlet fission since 2006, the year when the research interest in this subject was reignited. Both electronic structure and dynamics studies are covered. Electronic structure studies provide guidelines for designing singlet fission chromophores and insights into the couplings between single- and multi-excitonic states. The latter provides fundamental knowledge for engineering interchromophore conformations to enhance the fission efficiency. Dynamics studies reveal the importance of vibronic couplings in singlet fission.

Ultrafast isomerization in acetylene dication after carbon K-shell ionization

Ultrafast proton migration and isomerization are key processes for acetylene and its ions. However, the mechanism for ultrafast isomerization of acetylene [HCCH]²⁺ to vinylidene [H₂CC]²⁺ dication remains nebulous. Theoretical studies show a large potential barrier ( > 2 eV) for isomerization on low-lying dicationic states, implying picosecond or longer isomerization timescales. However, a recent experiment at a femtosecond X-ray free-electron laser suggests sub-100 fs isomerization. Here we address this contradiction with a complete theoretical study of the dynamics of acetylene dication produced by Auger decay after X-ray photoionization of the carbon atom K shell. We find no sub-100 fs isomerization, while reproducing the salient features of the time-resolved Coulomb imaging experiment. This work resolves the seeming contradiction between experiment and theory and also calls for careful interpretation of structural information from the widely applied Coulomb momentum imaging method.

The timescale of isomerization in molecules involving ultrafast migration of constituent atoms is difficult to measure. Here the authors report that sub-100 fs isomerization time on acetylene dication in lower electronic states is not possible and point to misinterpretation of recent experimental results.

A diabatic electronic state system to describe the internal conversion of azulene

A diabatic system of two electronic potential energy surfaces as well as the coupling between them is presented. The system is to be used to study the dynamics of the S₁ → S₀ internal conversion of azulene and is based on single point calculations of the minima of the two surfaces and a dipole-quadrupole (DQ) diabatization. Based on this, a couple of harmonic diabatic surfaces together with a linear coupling surface have been devised. Some preliminary dynamics results are shown.

Towards an understanding of the singlet–triplet splittings in conjugated hydrocarbons: azulene investigated by anion photoelectron spectroscopy and theoretical calculations

Electron density overlaps are correlated with singlet triplet splittings: azulene as a test case.

Seeking small molecules for singlet fission: a heteroatom substitution strategy

We design theoretically small molecule candidates for singlet fission chromophores, aiming to achieve a balance between sufficient diradical character and kinetic persistence. We develop a perturbation strategy based on the captodative effect to introduce diradical character into small π-systems. Specifically, this can be accomplished by replacing pairs of not necessarily adjacent C atoms with isoelectronic and isosteric pairs of B and N atoms. Three rules of thumb emerge from our studies to aid further design: (i) Lewis structures provide insight into likely diradical character; (ii) formal radical centers of the diradical must be well-separated; (iii) stabilization of radical centers by a donor (N) and an acceptor (B) is essential. Following the rules, we propose candidate molecules. Employing reliable multireference calculations for excited states, we identify three likely candidate molecules for SF chromophores. These include a benzene, a napthalene, and an azulene, where four C atoms are replaced by a pair of B and a pair of N atoms.

Modulating structure and properties in organic chromophores: influence of azulene as a building block

The properties of isomeric azulene derivatives, substituted through the 5-membered ring, were examined using a combination of experimentation and theoretical calculations for a series of well-defined electroactive oligomers.

Quantum chemical insights into the aggregation induced emission phenomena: a QM/MM study for pyrazine derivatives

There have been intensive studies on the newly discovered phenomena called aggregation induced emission (AIE), in contrast to the conventional aggregation quenching. Through combined quantum mechanics and molecular mechanics computations, we have investigated the aggregation effects on the excited state decays, both via radiative and nonradiative routes, for pyrazine derivatives 2,3-dicyano-5,6-diphenylpyrazine (DCDPP) and 2,3-dicyanopyrazino phenanthrene (DCPP) in condensed phase. We show that for DCDPP there appear AIE for all the temperature, because the phenyl ring torsional motions in gas phase can efficiently dissipate the electronic excited state energy, and get hindered in aggregate; while for its "locked"-phenyl counterpart, DCPP, theoretical calculation can only give the normal aggregation quenching. These first-principles based findings are consistent with recent experiment. The primary origin of the exotic AIE phenomena is due to the nonradiative decay effects. This is the first time that AIE is understood based on theoretical chemistry calculations for aggregates, which helps to resolve the present disputes over the mechanism.

Communication: A dramatic transition from nonferromagnet to ferromagnet in finite fused-azulene chain

One-dimensional fused-azulene oligomers (n = 2-6) are studied with the effective valence bond as well as density functional theory methods. A nonferromagnetic (closed-shell singlet) to ferromagnetic (triplet) ground state transformation is witnessed with increasing length of oligomers. The computational results are interpreted in terms of spin coupling between the unpaired electrons of two nonbonding molecular orbitals localized, respectively, on the top and bottom chains of the oligomers. The present study provides a theoretical suggestion for understanding the ferromagnetic spin polarizations that has been observed very recently in graphene nanoribbons.