Efficient Photochemical Merocyanine-to-Spiropyran Ring Closure Mechanism through an Extended Conical Intersection Seam. A Model CASSCF/CASPT2 Study

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.

Asymmetric Dressing of WSe2 with (Macro)molecular Switches: Fabrication of Quaternary-Responsive Transistors

The forthcoming saturation of Moore's law has led to a strong demand for integrating analogue functionalities within semiconductor-based devices. As a step toward this goal, we fabricate quaternary-responsive WSe₂-based field-effect transistors (FETs) whose output current can be remotely and reversibly controlled by light, heat, and electric field. A photochromic silane-terminated spiropyran (SP) is chemisorbed on SiO₂ forming a self-assembled monolayer (SAM) that can switch from the SP to the merocyanine (MC) form in response to UV illumination and switch back by either heat or visible illumination. Such a SAM is incorporated at the dielectric-semiconductor interface in WSe₂-based FETs. Upon UV irradiation, a drastic decrease in the output current of 82% is observed and ascribed to the zwitterionic MC isomer acting as charge scattering site. To provide an additional functionality, the WSe₂ top surface is coated with a ferroelectric co-polymer layer based on poly(vinylidene fluoride-co-trifluoroethylene). Because of its switchable inherent electrical polarization, it can promote either the accumulation or depletion of charge carriers in the WSe₂ channel, thereby inducing a current modulation with 99% efficiency. Thanks to the efficient tuning induced by the two components and their synergistic effects, the device polarity could be modulated from n-type to p-type. Such a control over the carrier concentration and device polarity is key to develop 2D advanced electronics. Moreover, the integration strategy of multiple stimuli-responsive elements into a single FET allows us to greatly enrich its functionality, thereby promoting the development for More-than-Moore technology.

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.

A generalized vibronic-coupling Hamiltonian for molecules without symmetry: Application to the photoisomerization of benzopyran

We present a model for the lowest two potential energy surfaces (PESs) that describe the photoinduced ring-opening reaction of benzopyran taken as a model compound to study the photochromic ring-opening reaction of indolinobenzospiropyran and its evolution toward its open-chain analog. The PESs are expressed in terms of three effective rectilinear coordinates. One corresponds to the direction between the equilibrium geometry in the electronic ground state, referred to as the Franck-Condon geometry, and the minimum of conical intersection (CI), while the other two span the two-dimensional branching space at the CI. The model correctly reproduces the topography of the PESs. The ab initio calculations are performed with the extended multiconfiguration quasidegenerate perturbation theory at second order method. We demonstrate that accounting for electron dynamic correlation drastically changes the global energy landscape since some zwitterionic states become strongly stabilized. Quantum dynamics calculations using this PES model produce an absorption spectrum that matches the experimental one to a good accuracy.

Directional and regioselective hole injection of spiropyran photoswitches intercalated into A/T-duplex DNA

The hole electron transfer of UV excited spiropyran intercalated in dsDNA is directional, asymmetric and regioselective, as shown by quantitative multiscale computations.

Discovery of conical intersection mediated photochemistry with growing string methods

Conical intersections (CIs) are important features of photochemistry that determine yields and selectivity. Traditional CI optimizers require significant human effort and chemical intuition, which typically restricts searching to only a small region of the CI space. Herein, a systematic approach utilizing the growing string method is introduced to locate multiple CIs. Unintuitive MECI are found using driving coordinates that can be generated using a combinatorial search, and subsequent optimization allows reaction pathways, transition states, products, and seam-space pathways to be located. These capabilities are demonstrated by application to two prototypical photoisomerization reactions and the dimerization of butadiene. In total, many reaction pathways were uncovered, including the elusive stilbene hula-twist mechanism, and a previously unidentified product in butadiene dimerization. Overall, these results suggest that growing string methods provide a predictive strategy for exploring photochemistry.

Ultrafast Spintronics: Dynamics of the Photoisomerization-Induced Spin-Charge Excited-State (PISCES) Mechanism in Spirooxazine-Based Photomagnetic Materials

The optical control of spin state is of interest in the development of spintronic materials for data processing and storage technologies. Photomagnetic effects at the single-molecule level have recently been observed in the thin film state at 300 K in photochromic cobalt dioxolenes. Visible light excitation leads to ring-closure of a photochromic spirooxazine bound to a cobalt dioxolene, which leads to generation of a high magnetization state. Formation of the photomagnetic state occurs through a photoisomerization-induced spin-charge excited-state process and is dictated by the spirooxazine ligand dynamics. Here, we report a mechanistic investigation by ultrafast spectroscopy in the UV-vis region of the photochemical ring-closing process in the parent spirooxazine, azahomoadamantylphenanthroline spirooxazine, and the photomagnetic spirooxazine cobalt-dioxolene complex. The cobalt appears to stabilize a photomerocycanine transient intermediate, presumably the TCC isomer, formed along the ground-state potential energy surface (PES). Structural changes associated with the TCC isomer induces formation of the high-spin Co(II) form, suggesting that magnetization dymanics can occur along the excited-state PES, leading to ultrafast switching on the ps time scale. We demonstrate the full ring closure of the spiro-oxazine ligand is not required to switch magnetization states which can be induced with a higher yielding isomerization reaction. The ability of this system to undergo optically induced spin state switching on the ps time scale in the solid state makes it a promising canididate for resistive nonvolatile memory technologies.

Synchronised photoreversion of spirooxazine ring opening in thin crystals to uncover ultrafast dynamics

A ‘recover before destroy’ approach to minimise photoproduct build-up in solid state enables ultrafast studies of chemical reactions.

Rationalization and Design of Enhanced Photoinduced Cycloreversion in Photochromic Dimethyldihydropyrenes by Theoretical Calculations

This study presents a computational investigation of the initial step of the dimethyldihydropyrene (DHP) to cyclophanediene (CPD) photoinduced ring-opening reaction using time-dependent density functional theory (TD-DFT). In particular, the photochemical path corresponding to the formation of the CPD precursor (CPD*) on the zwitterionic state is scrutinized. The TD-DFT approach was first validated on the parent compound against accurate ab initio calculations. It confirms that CPD* formation is efficiently quenched in this system by an easily accessible S2/S1 conical intersection located in the vicinity of the CPD* minimum and leading to a locally excited state minimum responsible for DHP luminescence. Increased ring-opening quantum yields were observed in benzo[e]-fused-DHP (DHP-1), isobutenyl-DHP (DHP-2), and naphthoyl-DHP (DHP-3). The calculations show that CPD* formation is much more favorable in these systems, either due to an inversion of electronic states in DHP-1, suppressing the formation of the locally excited state, or due to efficient stabilization of CPD* on the S1 potential energy surface in DHP-2 and DHP-3. Both effects can be combined in a rationally designed benzo[e]-fused-naphthoyl-DHP (DHP-4) for which we anticipate an unprecedented efficiency.

Unravelling the S → O linkage photoisomerization mechanisms in cis- and trans-[Ru(bpy)2(DMSO)2](2+) using density functional theory

A mechanistic study of the intramolecular S → O linkage photoisomerization in the cis and trans isomers of [Ru(bpy)2(DMSO)2](2+) was performed using density functional theory. This study reveals that for the cis isomer the linkage photoisomerization of the two DMSO ligands occurs sequentially in the lowest triplet excited state and can either be achieved by a one-photon or by a two-photon mechanism. A mechanistic picture of the S → O photoisomerization of the trans isomer is also proposed. This work especially highlights that both adiabatic and nonadiabatic processes are involved in these mechanisms and that their coexistence is responsible for the rich photophysics and photochemical properties observed experimentally for the studied complexes. The different luminescent behavior experimentally observed at low temperature between the cis and trans isomers is rationalized based on the peculiarity of the topology of the triplet excited-state potential energy surfaces.

One photon yields two isomerizations: large atomic displacements during electronic excited-state dynamics in ruthenium sulfoxide complexes

Photochromic compounds efficiently transduce photonic energy to potential energy for excited-state bond-breaking and bond-forming reactions. A critical feature of this reaction is the nature of the electronic excited-state potential energy surface and how this surface facilitates large nuclear displacements and rearrangements. We have prepared two photochromic ruthenium sulfoxide complexes that feature two isomerization reactions following absorption of a single photon. We show by femtosecond transient absorption spectroscopy that this reaction is complete within a few hundred picoseconds and suggest that isomerization occurs along a conical intersection seam formed by the ground-state and excited-state potential energy surfaces.

Spiropyran to merocyanine conversion: explicit versus implicit solvent modeling

A polarizable continnum model study and explicit solvation in water investigated through GPU-accelerated ab initio molecular dynamics followed by quantum-chemical calculations have been applied to the process of spiropyran to merocyanine isomerization of two spiropyran derivatives. It has been found that interaction with only one or two water molecules is sufficient to stabilize the merocyanine with respect to the spiropyran form. It has been shown that the agreement between energies obtained in implicit and explicit solvent models depends on the structure of the system and possible specific interactions (hydrogen bonds). Both solvent models predict similar effects in absorption spectra.

Computational study of the mechanism of the photochemical and thermal ring-opening/closure reactions and solvent dependence in spirooxazines

The spirooxazine/merocyanine couple constitutes a photochromic system that can change from the colorless spirooxazine to the intensely colored merocyanine by thermal or photochemical activation by a reaction that opens the spiro ring of the oxazine. The mechanisms of the ring-opening/closure reactions that interconnect these two isomers have been elucidated by means of a computational study. First, we have used the CASSCF/CASPT2 method to determine in detail these mechanisms in the gas phase for a small model that contains the photoactive part of the whole system. We have found that the state of spirooxazine excited by the initial absorption changes gradually to a lower excited state that is involved in a conical intersection that connects it with the ground state of merocyanine. The same conical intersection is involved in the backward photochemical reaction. Second, using a larger model that includes all the heteroatoms of the system and using the DFT (B3LYP) method, we have studied the influence of a solvent environment on the thermal equilibrium between the open and the closed species. It has been observed experimentally that the thermal equilibrium between these forms is practically unaltered by polar aprotic solvents, but it can be displaced toward the colored form in mixtures of polar protic and aprotic solvents, even if the first one is found in a very small proportion. To reproduce the experimental environments, we have taken into account the long-range effect of the polar aprotic solvent considering it a polarizable continuum, as done in the PCM method, and the short-range effect of the protic solvent including some explicit water molecules in the cluster studied at the atomic level. The results obtained are in good agreement with the experimental observations and explain the reason for this peculiar behavior.

Solvent dependence of structure, charge distribution, and absorption spectrum in the photochromic merocyanine-spiropyran pair

We have studied the structures and absorption spectra of merocyanine, the photoresponsive isomer of the spiropyran (SP)-merocyanine (MC) pair, in chloroform and in water solvents using a combined hybrid QM/MM Car-Parrinello molecular dynamics (CP-QM/MM) and ZINDO approach. We report remarkable differences in the molecular structure and charge distribution of MC between the two solvents; the molecular structure of MC remains in neutral form in chloroform while it becomes charge-separated, zwitterionic, in water. The dipole moment of MC in water is about 50% larger than in chloroform, while the value for SP in water is in between, suggesting that the solvent is more influential than the conformation itself in deciding the dipole moment for the merocyanine-spiropyran pair. The calculations could reproduce the experimentally reported blue shift in the absorption spectra of MC when going from the nonpolar to the polar solvent, though the actual value of the absorption maximum is overestimated in chloroform solvent. We find that the CP-QM/MM approach is appropriate for structure modeling of solvatochromic and thermochromic molecules as this approach is able to capture the solvent and thermal-induced structural changes within the solute important for an accurate assessment of the properties.

Ring-closure and isomerization capabilities of spiropyran-derived merocyanine isomers

We report the photochemistry of two ring-open isomers, namely TTC and TTT, of a bidirectional photoswitchable spiropyran, 6,8-dinitro-1',3',3'-trimethylspiro[2H-1-benzopyran-2,2'-indoline] (6,8-dinitro BIPS). Both isomers are capable of ring closure after excitation with visible fs laser pulses, as disclosed by pump-wavelength-dependent transient absorption experiments in the visible spectral range. The main isomer TTC has its maximum absorption at 560 nm, whereas the minor isomer TTT is red-shifted (600 nm). The excited-state lifetimes differ strongly (τ ≈ 900 ps for TTT and τ ≈ 95 ps for TTC), nevertheless the quantum efficiencies for ring closure (40% for TTC and 35% for TTT) and isomerization (1-2% for TTC and 1-2% for TTT) are comparable. With regard to the bidirectional photoswitching capabilities, 6,8-dinitro BIPS is the first molecular switch based on a 6π-electrocyclic reaction where both ring-open isomers are capable of ring closure.

Photoreactivity of a push-pull merocyanine in static electric fields: a three-state model of isomerization reactions involving conical intersections

The photochemistry of a prototype push-pull merocyanine is discussed using a simple three-state model. As a derivative of butadiene, the model focuses on two isomerization reactions around the two double bonds of the butadiene backbone. As a molecule substituted by an electron donor and electron acceptor at opposite ends, its structure as well as its photochemistry are expected to be strongly affected by the environment. In polar solvents, a zwitterion transition state for each of the isomerization reactions is stabilized, and its energy is on the same order as that of the biradical one; this leads to the symmetry allowed crossing (S(0)/S(1) conical intersection). It is shown that applying an external electric field or varying the solvent polarity changes the relative energies of the different transition states as well as that of the conical intersection, and thus different photochemical products can be obtained. In particular, the very existence of conical intersections is found to depend on these external parameters. This work provides a theoretical foundation for ideas expressed by Squillacote et al. (J. Am. Chem. Soc. 2004, 126, 1940) concerning the electrostatic control of photochemical reactions.

Opening a Spiropyran Ring by Way of an Exciplex Intermediate

A molecular dyad has been synthesized in which the main chromophore is a 1,4-diethynylated benzene residue terminated with pyrene moieties, this latter unit acting as a single chromophore. A spiropyran group has been condensed to the central phenylene ring so as to position a weak electron donor close to the pyrene unit. Illumination of the pyrene-based chromophore leads to formation of a fluorescent exciplex in polar solvents but pyrene-like fluorescence is observed in nonpolar solvents. The exciplex has a lifetime of a few nanoseconds and undergoes intersystem crossing to the pyrene-like triplet state with low efficiency. Attaching a 4-nitrobenzene group to the open end of the spiropyran unit creates a new route for decay of the exciplex whereby the triplet state of the spiropyran is formed. Nonradiative decay of this latter species results in ring opening to form the corresponding merocyanine species. Rate constants for the various steps have been obtained from time-resolved fluorescence spectroscopy carried out over a modest temperature range. Under visible light illumination, the merocyanine form reverts to the original spiropyran geometry so that the cycle is closed. Energy transfer from the pyrene chromophore to the merocyanine unit leads to an increased rate of ring closure and serves to push the steady-state composition in favor of the spiropyran form.