Semiempirical Hamiltonian for simulation of azobenzene photochemistry

The Role of Double Excitations in Exciton Dynamics of Multiazobenzenes: Trisazobenzenophane as a Test Case

Molecular exciton dynamics underlie energy and charge transfer processes in organic multichromophoric systems. A particularly interesting class of the latter is multiphotochromic systems made of molecules capable of photochemical transformations. Exciton dynamics in assemblies of photoswitches have been recently investigated using either the molecular exciton model or supermolecular configuration interaction (CI) singles, both approaches being based on a semiempirical Hamiltonian and combined with surface hopping molecular dynamics. Here, we study how inclusion of double excitations in nonadiabatic dynamics simulations affects exciton dynamics of multiazobenzenes, using trisazobenzenophane as an example. We find that both CI singles and CI singles and doubles yield virtually the same time scale of dynamical exciton localization, ∼50 fs for the studied multiazobenzene. However, inclusion of double excitations considerably affects the excited state lifetimes and isomerization quantum yields.

Photoisomerization of Spiropyran in Solution: A Surface Hopping Investigation

We performed a computational investigation of the photoisomerization of an unsubstituted spiropyran in solution using surface hopping molecular dynamics simulations in a semiempirical framework. To bring out the solvent effects on the excited state dynamics, we have run simulations in three different environments: chloroform, methanol, and ethylene glycol. We show that the interaction with a moderately polar solvent such as chloroform has little effect on the dynamics when compared to vacuum results previously obtained by our group. At variance, the interaction with the protic solvents considered considerably affects the reaction mechanism, the quantum yield, and the excited state lifetimes.

Topology of Conical Intersection Seams and the Geometric Phase

In this paper, we demonstrate two topological properties of crossing seams, that is, the sets of points in the N-dimensional space of nuclear coordinates where two electronic eigenstates are degenerate. We shall examine the typical case of states of the same spin with accidental degeneracies, whereby the crossing seam is of dimension N - 2. The first property we demonstrate is that a crossing seam has no boundary, therefore, it must either extend asymptotically to infinite values of one or more coordinates or wrap on itself. The second property is that two (or more) crossing seams can intersect each other but in such a way that neither of them ends at the intersection. When N = 3, the crossing seam is a line in a 3D space; this is so in triatomic molecules but also in reduced dimensionality treatments of larger polyatomics. The above-mentioned rules then mean that the crossing seam is a line of infinite length or a closed loop and can split into three branches but not in two. The example of the first two excited 1A' states of H2Cl+ illustrates these rules and shows their usefulness for computational search and characterization of crossing seams.

Visible Light Induced Exciton Dynamics and trans-to-cis Isomerization in Azobenzene Aggregates: Insights from Surface Hopping/Semiempirical Configuration Interaction Molecular Dynamics Simulations

Assemblies of photochromic molecules feature exciton states, which govern photochemical and photophysical processes in multichromophoric systems. Understanding the photoinduced dynamics of the assemblies requires nonadiabatic treatment involving multiple exciton states and numerous nuclear degrees of freedom, thus posing a challenge for simulations. In this work, we address this challenge for aggregates of azobenzene, a prototypical molecular switch, performing on-the-fly surface hopping calculations combined with semiempirical configuration interaction electronic structure and augmented with transition density matrix analysis to characterize exciton evolution. Specifically, we consider excitation of azobenzene tetramers in the nπ* absorption band located in the visible (blue) part of the electromagnetic spectrum, thus extending our recent work on dynamics after ππ* excitation corresponding to the ultraviolet region [Titov, J. Phys. Chem. C2023, 127, 13678-13688]. We find that the nπ* excitons, which are initially strongly localized by ground-state conformational disorder, undergo further (very strong) localization during short-time photodynamics. This excited-state localization process is extremely ultrafast, occurring within the first 10 fs of photodynamics. We observe virtually no exciton transfer of the localized excitons in the nπ* manifold. However, the transfer may occur via secondary pathways involving ππ* states or the ground state. Moreover, we find that the nπ* quantum yields of the trans-to-cis isomerization are reduced in the aggregated state.

Investigating the Photodynamics of trans-Azobenzene with Coupled Trajectories

In this work, we present the first implementation of coupled-trajectory Tully surface hopping (CT-TSH) suitable for applications to molecular systems. We combine CT-TSH with the semiempirical floating occupation molecular orbital-configuration interaction electronic structure method to investigate the photoisomerization dynamics of trans-azobenzene. Our study shows that CT-TSH can capture correctly decoherence effects in this system, yielding consistent electronic and nuclear dynamics in agreement with (standard) decoherence-corrected TSH. Specifically, CT-TSH is derived from the exact factorization and the electronic coefficients' evolution is directly influenced by the coupling of trajectories, resulting in the improvement of internal consistency if compared to standard TSH.

Theoretical Advances in Polariton Chemistry and Molecular Cavity Quantum Electrodynamics

When molecules are coupled to an optical cavity, new light-matter hybrid states, so-called polaritons, are formed due to quantum light-matter interactions. With the experimental demonstrations of modifying chemical reactivities by forming polaritons under strong light-matter interactions, theorists have been encouraged to develop new methods to simulate these systems and discover new strategies to tune and control reactions. This review summarizes some of these exciting theoretical advances in polariton chemistry, in methods ranging from the fundamental framework to computational techniques and applications spanning from photochemistry to vibrational strong coupling. Even though the theory of quantum light-matter interactions goes back to the midtwentieth century, the gaps in the knowledge of molecular quantum electrodynamics (QED) have only recently been filled. We review recent advances made in resolving gauge ambiguities, the correct form of different QED Hamiltonians under different gauges, and their connections to various quantum optics models. Then, we review recently developed ab initio QED approaches which can accurately describe polariton states in a realistic molecule-cavity hybrid system. We then discuss applications using these method advancements. We review advancements in polariton photochemistry where the cavity is made resonant to electronic transitions to control molecular nonadiabatic excited state dynamics and enable new photochemical reactivities. When the cavity resonance is tuned to the molecular vibrations instead, ground-state chemical reaction modifications have been demonstrated experimentally, though its mechanistic principle remains unclear. We present some recent theoretical progress in resolving this mystery. Finally, we review the recent advances in understanding the collective coupling regime between light and matter, where many molecules can collectively couple to a single cavity mode or many cavity modes. We also lay out the current challenges in theory to explain the observed experimental results. We hope that this review will serve as a useful document for anyone who wants to become familiar with the context of polariton chemistry and molecular cavity QED and thus significantly benefit the entire community.

Dual photoisomerization mechanism of azobenzene embedded in a lipid membrane

The photoisomerization of chromophores embedded in biological environments is of high importance for biomedical applications, but it is still challenging to define the photoisomerization mechanism both experimentally and computationally. We present here a computational study of the azobenzene molecule embedded in a DPPC lipid membrane, and assess the photoisomerization mechanism by means of the quantum mechanics/molecular mechanics surface hopping (QM/MM-SH) method. We observe that while the trans-to-cis isomerization is a slow process governed by a torsional mechanism due to the strong interaction with the environment, the cis-to-trans mechanism is completed in sub-ps time scale and is governed by a pedal-like mechanism in which both weaker interactions with the environment and a different geometry of the potential energy surface play a key role.

Excited state non-adiabatic dynamics of large photoswitchable molecules using a chemically transferable machine learning potential

Light-induced chemical processes are ubiquitous in nature and have widespread technological applications. For example, photoisomerization can allow a drug with a photo-switchable scaffold such as azobenzene to be activated with light. In principle, photoswitches with desired photophysical properties like high isomerization quantum yields can be identified through virtual screening with reactive simulations. In practice, these simulations are rarely used for screening, since they require hundreds of trajectories and expensive quantum chemical methods to account for non-adiabatic excited state effects. Here we introduce a diabatic artificial neural network (DANN), based on diabatic states, to accelerate such simulations for azobenzene derivatives. The network is six orders of magnitude faster than the quantum chemistry method used for training. DANN is transferable to azobenzene molecules outside the training set, predicting quantum yields for unseen species that are correlated with experiment. We use the model to virtually screen 3100 hypothetical molecules, and identify novel species with high predicted quantum yields. The model predictions are confirmed using high-accuracy non-adiabatic dynamics. Our results pave the way for fast and accurate virtual screening of photoactive compounds.

Surface Hopping Dynamics with the Frenkel Exciton Model in a Semiempirical Framework

We present an implementation of the Frenkel exciton model in the framework of the semiempirical floating occupation molecular orbitals-configuration interaction (FOMO-CI) electronic structure method, aimed at simulating the dynamics of multichromophoric systems, in which excitation energy transfer can occur, by a very efficient approach. The nonadiabatic molecular dynamics is here dealt with by the surface hopping method, but the implementation we proposed is compatible with other dynamical approaches. The exciton coupling is computed either exactly, within the semiempirical approximation considered, or by resorting to transition atomic charges. The validation of our implementation is carried out on the trans-azobenzeno-2S-phane (2S-TTABP), formed by two azobenzene units held together by sulfur bridges, taken as a minimal model of multichromophoric systems, in which both strong and weak exciton couplings are present.

On the Low-Lying Electronically Excited States of Azobenzene Dimers: Transition Density Matrix Analysis

Azobenzene-containing molecules may associate with each other in systems such as self-assembled monolayers or micelles. The interaction between azobenzene units leads to a formation of exciton states in these molecular assemblies. Apart from local excitations of monomers, the electronic transitions to the exciton states may involve charge transfer excitations. Here, we perform quantum chemical calculations and apply transition density matrix analysis to quantify local and charge transfer contributions to the lowest electronic transitions in azobenzene dimers of various arrangements. We find that the transitions to the lowest exciton states of the considered dimers are dominated by local excitations, but charge transfer contributions become sizable for some of the lowest ππ* electronic transitions in stacked and slip-stacked dimers at short intermolecular distances. In addition, we assess different ways to partition the transition density matrix between fragments. In particular, we find that the inclusion of the atomic orbital overlap has a pronounced effect on quantifying charge transfer contributions if a large basis set is used.

Photoisomerization dynamics of spiropyran: A surface-hopping investigation

In the present work, we performed a computational investigation of the photoisomerization of spiro[1,3-dihydroindole-2,2'-chromene] [spiropyran (SP)] to merocyanine. The electronic energies and wavefunctions were obtained from configuration interaction calculations, using the floating occupation molecular orbital method, in a semiempirical framework. The parameters of the semiempirical Hamiltonian were re-optimized to reproduce ab initio literature data for SP. In our dynamics simulations, we considered, besides S₀, the excited states S₁, S₂, and S₃, which are very close in energy in the Franck-Condon region. We obtained a singlet lifetime of 0.67 ps, in line with the experimental results. We found the photoisomerization quantum yield to depend on the electronic state initially populated.

Hybrid theoretical models for molecular nanoplasmonics

The multidisciplinary nature of the research in molecular nanoplasmonics, i.e., the use of plasmonic nanostructures to enhance, control, or suppress properties of molecules interacting with light, led to contributions from different theory communities over the years, with the aim of understanding, interpreting, and predicting the physical and chemical phenomena occurring at molecular- and nano-scale in the presence of light. Multiscale hybrid techniques, using a different level of description for the molecule and the plasmonic nanosystems, permit a reliable representation of the atomistic details and of collective features, such as plasmons, in such complex systems. Here, we focus on a selected set of topics of current interest in molecular plasmonics (control of electronic excitations in light-harvesting systems, polaritonic chemistry, hot-carrier generation, and plasmon-enhanced catalysis). We discuss how their description may benefit from a hybrid modeling approach and what are the main challenges for the application of such models. In doing so, we also provide an introduction to such models and to the selected topics, as well as general discussions on their theoretical descriptions.

Photochemistry in the strong coupling regime: A trajectory surface hopping scheme

The strong coupling regime between confined light and organic molecules turned out to be promising in modifying both the ground state and the excited states properties. Under this peculiar condition, the electronic states of the molecule are mixed with the quantum states of light. The dynamical processes occurring on such hybrid states undergo several modifications accordingly. Hence, the dynamical description of chemical reactivity in polaritonic systems needs to explicitly take into account the photon degrees of freedom and nonadiabatic events. With the aim of describing photochemical polaritonic processes, in the present work, we extend the direct trajectory surface hopping scheme to investigate photochemistry under strong coupling between light and matter.

Full-Dimensional Photodynamics of Bistable Proton Transfer Switches

Excited-state intramolecular proton transfers (ESIPT) are one of the fastest reactions in chemistry (<100 fs) which – among other features like high photostability – makes them an important reaction class for molecular switches. ESIPTs can be coupled with double bond rotation/isomerization, so that molecules can act as “molecular cranes”, facilitating long-range proton transfer. A versatile model system is 7-hydroxy-4-methylquinoline-8-carbaldehyde (HMQCA): it features two proton-accepting sites, two stable ground-state isomers and should allow for easy derivatization. There is also experimental and theoretical reference data available, however, only for static properties, e.g. ground-state IR spectra or potential energy surface scans. In this contribution we show the results of full-dimensional surface-hopping molecular dynamics (MD) of HMQCA after photo-excitation, employing semiempirical quantum mechanics coupled to floating-occupation configuration interaction. The results support the potential of HMQCA as prototype system for directed proton transport by ESIPT.

Nonadiabatic dynamics simulations of singlet fission in 2,5-bis(fluorene-9-ylidene)-2,5-dihydrothiophene crystals

We present simulations of the singlet fission dynamics in 2,5-bis(fluorene-9-ylidene)-2,5-dihydrothiophene (ThBF), a thienoquinoid compound recently investigated experimentally by Kawata et al. The simulation model consisted of two ThBF molecules embedded in their crystal environment. The aim was to understand the singlet fission mechanism, and to predict the excited state lifetimes and the singlet fission quantum yield, hitherto unknown. The simulations were performed by the trajectory surface hopping approach with on-the-fly calculations of the electronic wave functions and energies by the semiempirical FOMO-CI method. We found that the initially photogenerated excitonic bright state decays to the lower dark state with a biexponential behaviour, essentially due to transitions to other close lying states. The dark state in turn decays with a lifetime of about 1 ps to the double triplet 1TT state, which is long-lived, as ascertained by performing a simulation with inclusion of the spin-orbit coupling. The singlet fission quantum yield is predicted to be close to the theoretical maximum of 200%. In view of using this thienoquinoid compound in photovoltaic devices, a major drawback is the low energy of the T1 state at its equilibrium geometry.

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).

Simulations of optically switchable molecular machines for particle transport

A promising application for design and deployment of molecular machines is nanoscale transport, driven by artificial cilia. In this contribution, we present several further steps toward this goal, beyond our first-generation artificial cilium (Raeker et al., J. Phys. Chem. A 2012, 116, 11241). Promising new azobenzene-derivatives were tested for use as cilium motors. Using a QM/MM partitioning in on-the-fly photodynamics, excited-state surface-hopping trajectories were calculated for each isomerization direction and each motor version. The methods used were reparametrized semiempirical quantum chemistry together with floating-occupation configuration interaction as the QM part and the OPLSAA-L forcefield as MM part. In addition, we simulated actual particle transport by a single cilium attached to a model surface, with varying attachment strengths and modes, and with transport targets ranging from single atoms to multi-molecule arrangements. Our results provide valuable design guidelines for cilia-driven nanoscale transport and emphasize the need to carefully select the whole setup (not just the cilium itself, but also its surface attachment and the dynamic cilium-target interaction) to achieve true transport. © 2018 Wiley Periodicals, Inc.

Azobenzene as a photoregulator covalently attached to RNA: a quantum mechanics/molecular mechanics-surface hopping dynamics study

Azobenzene covalently attached to RNA undergoes trans-to-cis photo-switching on a time scale of ∼15 picoseconds – 30 times slower than in vacuo.

The photoregulation of nucleic acids by azobenzene photoswitches has recently attracted considerable interest in the context of emerging biotechnological applications. To understand the mechanism of photoinduced isomerisation and conformational control in these complex biological environments, we employ a Quantum Mechanics/Molecular Mechanics (QM/MM) approach in conjunction with nonadiabatic Surface Hopping (SH) dynamics. Two representative RNA–azobenzene complexes are investigated, both of which contain the azobenzene chromophore covalently attached to an RNA double strand via a β-deoxyribose linker. Due to the pronounced constraints of the local RNA environment, it is found that trans-to-cis isomerization is slowed down to a time scale of ∼10–15 picoseconds, in contrast to 500 femtoseconds in vacuo, with a quantum yield reduced by a factor of two. By contrast, cis-to-trans isomerization remains in a sub-picosecond regime. A volume-conserving isomerization mechanism is found, similarly to the pedal-like mechanism previously identified for azobenzene in solution phase. Strikingly, the chiral RNA environment induces opposite right-handed and left-handed helicities of the ground-state cis-azobenzene chromophore in the two RNA–azobenzene complexes, along with an almost completely chirality conserving photochemical pathway for these helical enantiomers.

Photoisomerization of Self-Assembled Monolayers of Azobiphenyls: Simulations Highlight the Role of Packing and Defects

We present surface hopping simulations of the photodynamics of self-assembled monolayers (SAMs) of 4'-(biphenyl-4-ylazo)-biphenyl-4-thiol (ABPT) on Au(111). We show that trans → cis photoisomerization is suppressed because of steric hindrance in a well-ordered SAM. Photoisomerization is instead viable in the presence of defects. Two particularly important defects are the boundaries between domains of trans-ABPT molecules leaning in different directions (a line defect) and single cis molecules embedded in a SAM of trans (a point defect). Our findings explain the cooperative behavior observed during the photoisomerization of a trans-ABPT SAM, leading to large domains of pure cis and trans isomers. The line and point defects are predicted to produce different patterns of cis-ABPT molecules during the early stages of the photoconversion.

Dynamics of Azobenzene Dimer Photoisomerization: Electronic and Steric Effects

While azobenzenes readily photoswitch in solution, their photoisomerization in densely packed self-assembled monolayers (SAMs) can be suppressed. Reasons for this can be steric hindrance and/or electronic quenching, e.g., by exciton coupling. We address these possibilities by means of nonadiabatic molecular dynamics with trajectory surface hopping calculations, investigating the trans → cis isomerization of azobenzene after excitation into the ππ* absorption band. We consider a free monomer, an isolated dimer and a dimer embedded in a SAM-like environment of additional azobenzene molecules, imitating in this way the gradual transition from an unconstrained over an electronically coupled to an electronically coupled and sterically hindered, molecular switch. Our simulations reveal that in comparison to the single molecule the quantum yield of the trans → cis photoisomerization is similar for the isolated dimer, but greatly reduced in the sterically constrained situation. Other implications of dimerization and steric constraints are also discussed.

Surface hopping investigation of benzophenone excited state dynamics

A mechanism of S₁ decay in benzophenone: S₁ → T₁ is the main pathway, although transitions to T₂ and higher triplets play a relevant role.

Theoretical study of the cis–trans isomerization mechanism of a pendant metal-bound azobenzene

Theoretical study of the cis–trans isomerization mechanism of azobenzene substituents in rhenium complexes.

Indandiazocines: Unidirectional molecular switches

We report theoretical investigations on azobenzene-based indandiazocines, novel chiral systems that perform unidirectional cis↔trans isomerizations upon photo-excitation. For three different systems of this kind, we have simulated excited-state surface-hopping trajectories for both isomerization directions, using a configuration-interaction treatment based on system-specifically reparametrized semiempirical AM1 theory. Our results are also compared to experimental and theoretical results for the parent system diazocine. We show that, as intended by design, the trans→cis bending of the azo unit in these indandiazocines can only happen in one of the two possible directions due to sterical constraints, which is a new feature for photoswitches and a necessary prerequisite for directional action at the nanoscale.

The photo-orientation of azobenzene in viscous solutions, simulated by a stochastic model

We report a computational study of the photo-orientation kinetics in a viscous solution of azobenzene in ethylene glycol, under irradiation with linearly polarized light. The development of anisotropy and its interplay with photoisomerization are simulated by a stochastic model. A distinctive feature of the model is that it takes into account the photo-orientation angular distributions, specific for each isomer, obtained by nonadiabatic dynamics simulations at the molecular level. We find that the anisotropy, as measured by optical absorption dichroism, does not necessarily increase monotonously with time. As expected, the photo-orientation turns out to be strongly coupled with photoisomerization, but the latter is not a mandatory ingredient of this phenomenon: we predict that any chromophore undergoing large amplitude geometry relaxation during its excited state dynamics can develop anisotropy under suitable conditions.