Theoretical study of photochromic compounds. 1. Bond length alternation and absorption spectra for the open and closed forms of 29 diarylethene derivatives

QM/MM Study of the Catalytic Mechanism and Substrate Specificity of the Aromatic Substrate C-Methyltransferase Fur6

In the field of medical chemistry and other organic chemistry, introducing a methyl group into a designed position has been difficult to achieve. However, owing to the vigorous developments in the field of enzymology, methyltransferases are considered potential tools for addressing this problem. Within the methyltransferase family, Fur6 catalyzes the methylation of C3 of 1,2,4,5,7-pentahydroxynaphthalene (PHN) using S-adenosyl-l-methionine (SAM) as the methyl donor. Here, we report the catalytic mechanism and substrate specificity of Fur6 based on computational studies. Our molecular dynamics (MD) simulation studies reveal the reactive form of PHN and its interactions with the enzyme. Our hybrid quantum mechanics/molecular mechanics (QM/MM) calculations suggest the reaction pathway of the methyl transfer step in which the energy barrier is 8.6 kcal mol-1. Our free-energy calculations with a polarizable continuum model (PCM) indicate that the final deprotonation step of the methylated intermediate occurs after it is ejected into the water solvent from the active center pocket of Fur6. Additionally, our studies on the protonation states, the highest occupied molecular orbital (HOMOs), and the energy barriers of the methylation reaction for the analogs of PHN demonstrate the mechanism of the specificity to PHN. Our study provides valuable insights into Fur6 chemistry, contributing to a deeper understanding of molecular mechanisms and offering an opportunity to engineer the enzyme to achieve high yields of the desired product(s).

Unsymmetrical benzothiazole-based dithienylethene photoswitches

Herein, we investigate the structure-property relationships in a new series of benzothiazole based unsymmetrical hexafluorocyclopentene dithienylethenes (DTEs) and compare the results with the known facts for symmetric diarylethenes (DAEs). We reveal high photocyclization efficiency resulting from a significant shift of ground state equilibrium to the antiparallel conformation and a barrierless excited state pathway to conical intersection, which remains unperturbed even in polar solvents for most of the prepared DTEs. Furthermore, we uncover that the rate of back thermal cycloreversion correlates clearly more with the central C-C bond-length in the transition state than with the central C-C bond-length in the ground state of the cyclic form. Finally, our detailed vibrational spectral analysis of studied DTEs points out significant changes in Raman and infrared spectra during photoswitching cycles which pave the way for a non-destructive readout of stored information.

Excited state dynamics of normal dithienylethene molecules either isolated or deposited on an argon cluster

Real-time dynamics of the electronically excited open-ring isomer of 1,2-bis(2-methylbenzo[b]thiophen-3-yl)perfluorocyclopentene (BTF6) and 1,2-bis(2,4-dimethyl-5-phenyl-3-thienyl)perfluorocyclopentene (PTF6) molecules was investigated using a set-up that associates a molecular beam, femtosecond lasers and velocity map imaging. The molecules were either free in the gas phase or bound to an argon cluster. DFT and TDDFT calculations were performed on BTF6. The calculated vertical excitation energies indicate an excitation by the pump laser towards a superposition of S₅ and S₆ states. The free molecule dynamics was found to follow a three wavepacket model. One describes the parallel conformer (P) of these molecules. It is unreactive with respect to the ring closure reaction which is responsible for the photochromic property of these molecules. It has no observable decay at the experiment time scale (up to 350 ps). The other two wavepackets describe the reactive antiparallel conformer (AP). They are formed by an early splitting of the wavepacket that was launched initially by the pump laser. They can be considered as generated by excitation of different, essentially uncoupled, deformation modes. They subsequently evolve along independent pathways. One is directed ballistically towards a conical intersection (CI) and decays through the CI to a potential energy surface where it can no longer be detected. The other fraction of the wavepacket decays also towards undetected states but in this case the driving mechanism is a non-adiabatic electronic relaxation within a potential well of the energy surfaces where it was launched. When BTF6 and PTF6 molecules are bound to an argon cluster, the same three wavepacket model applies. The vibronic relaxation timespan is enhanced by a factor 5 and a larger fraction of AP conformers follows this pathway. In contrast, the time constant associated with the ballistic movement is enhanced by only a factor of 2.

Photocyclization of Diarylethenes: The Effect of Electron and Proton Acceptors as Additives

The effect of electron and proton acceptors on the photocyclization of diarylethenes has been studied. Without any additives, the deprotonation reaction is predominant, although other processes, including the sigmatropic shift, are not excluded. A deuterium exchange experiment has shown that a strong base (DABCO) facilitates the deprotonation reaction, thereby limiting the sigmatropic shift. In the presence of an oxidizing agent or additional sources of radicals (O₂, I₂, TEMPO), the processes of deprotonation and rearrangement (H-shift) are practically not observed, and the reaction proceeds along a radical pathway with the formation of phenanthrene or its heterocyclic analogue.

"On-The-Fly" Non-Adiabatic Dynamics Simulations on Photoinduced Ring-Closing Reaction of a Nucleoside-Based Diarylethene Photoswitch

Nucleoside-based diarylethenes are emerging as an especial class of photochromic compounds that have potential applications in regulating biological systems using noninvasive light with high spatio-temporal resolution. However, relevant microscopic photochromic mechanisms at atomic level of these novel diarylethenes remain to be explored. Herein, we have employed static electronic structure calculations (MS-CASPT2//M06-2X, MS-CASPT2//SA-CASSCF) in combination with non-adiabatic dynamics simulations to explore the related photoinduced ring-closing reaction of a typical nucleoside-based diarylethene photoswitch, namely, PS-IV. Upon excitation with UV light, the open form PS-IV can be excited to a spectroscopically bright S₁ state. After that, the molecule relaxes to the conical intersection region within 150 fs according to the barrierless relaxed scan of the C1–C6 bond, which is followed by an immediate deactivation to the ground state. The conical intersection structure is very similar to the ground state transition state structure which connects the open and closed forms of PS-IV, and therefore plays a crucial role in the photochromism of PS-IV. Besides, after analyzing the hopping structures, we conclude that the ring closing reaction cannot complete in the S₁ state alone since all the C1–C6 distances of the hopping structures are larger than 2.00 Å. Once hopping to the ground state, the molecules either return to the original open form of PS-IV or produce the closed form of PS-IV within 100 fs, and the ring closing quantum yield is estimated to be 56%. Our present work not only elucidates the ultrafast photoinduced pericyclic reaction of the nucleoside-based diarylethene PS-IV, but can also be helpful for the future design of novel nucleoside-based diarylethenes with better performance.

Improving photosensitivity without changing thermal reactivity in photochromic diarylbenzenes based on accurate prediction by DFT calculations

1,2-Diarylbenzenes (DABs) have been developed as a new family of fast T-type photochromic switches. However, the molecular design strategy for DABs with desired optical and thermal properties is not established. In this work, we explored the best functional in quantum chemical calculations to predict the properties of DABs. Furthermore, we newly designed and synthesized DABs based on the calculation using the best functional, resulting in the improvement of the photosensitivity in the UV-A region (i.e. a shift of absorption to lower energies and an increase in the absorption coefficient) without changing the thermal back-reaction rate.

Reversible Shifting of a Chemical Equilibrium by Light: The Case of Keto-Enol Tautomerism of a β-Ketoester

Manipulating the equilibrium between a ketone and an enol by light opens up ample opportunities in material chemistry and photopharmacology. By incorporating β-ketoester into the ethene bridge of a photoactive diarylethene, we achieved reversible light-induced tautomerization to give thermally stable enol. In a pristine state, the tautomeric equilibrium is almost completely shifted toward the ketone. Photocyclization of diarylethene results in a new equilibrium containing a significant fraction of the enol tautomer.

Energy- and conformer-dependent excited-state relaxation of an E/Z photoswitchable thienyl-ethene

Bis(bithienyl)-1,2-dicyanoethene (4TCE) is a photoswitch that operates via reversible E/Z photoisomerization following absorption of visible light. cis-to-trans photoisomerization of 4TCE requires excitation below 470 nm, is relatively inefficient (quantum yield < 5%) and occurs via the lowest-lying triplet. We present excitation-wavelength dependent (565-420 nm) transient absorption (TA) studies to probe the photophysics of cis-to-trans isomerization to identify sources of switching inefficiency. TA data reveals contributions from more than one switch conformer and relaxation cascades between multiple states. Fast (∼4 ps) and slow (∼40 ps) components of spectral dynamics observed at low excitation energies (>470 nm) are readily attributed to deactivation of two conformers; this assignment is supported by computed thermal populations and absorption strengths of two molecular geometries (P_A and P_B) characterized by roughly parallel dipoles for the thiophenes on opposite sides of the ethene bond. Only the P_B conformer is found to contribute to triplet population and the switching of cis-4TCE: high-energy excitation (<470 nm) of P_B involves direct excitation to S₂, relaxation from which prepares an ISC-active S₁ geometry (ISC QY 0.4-0.67, k_ISC∼ 1.6-2.6 × 10^-9 s^-1) that is the gateway to triplet population and isomerization. We ascribe low cis-to-trans isomerization yield to excitation of the nonreactive P_A conformer (75-85% loss) as well as loses along the P_B S₂→ S₁→ T₁ cascade (10-20% loss). In contrast, electrocyclization is inhibited by the electronic character of the excited states, as well as a non-existent thermal population of a reactive "antiparallel" ring conformation.

Is energy transfer limiting multiphotochromism? answers from ab initio quantifications

Dithienylethenes (DTEs) can be assembled to form supramolecular multiphotochromic systems that are highly functional molecular architectures of potential interest for building complex optoelectronic devices. Yet even simple DTE dimers relying on an organic linker may suffer from a partial photoactivity, i.e., only one of the two switches does isomerise. Experimentally, this limited photochromism has been attributed to an excited state energy transfer (EET) between the two DTEs of the multimer; this EET taking place instead of the desired photoinduced cyclisation of the DTE. However, no clear evidences of this phenomenon have been provided so far. In this work we propose the first rationalisation of this potential parasite photoinduced event using a computational approach based on Time-Dependent Density Functional Theory (TD-DFT) for the calculation of the electronic coupling in DTE dimers. Besides quantifying EET in several systems, we dissect the role of through-bond and through-space mechanisms on this process and clarify their dependence on both the nature and length of the bridge separating the two photochromes. The theoretical data obtained in this framework are in full agreement with the experimental outcomes and pave the way toward a molecular design of coupled, yet fully functionals, DTE-based multiswitches.

Theoretical analysis of electrochromism under redox of bis(3-thienyl)/(2-thienyl)hexafluorocyclopentene: effects of charged and substituted systems

Electrochromism with the ring-closing or ring-opening isomerization of substituted and unsubstituted bis(3-thienyl)/(2-thienyl)hexafluorocyclopentene is discussed using the DFT method. In the neutral ground state, bond making and breaking between two reactive C atoms on thienyls are thermodynamically forbidden. Under redox conditions, the gain or loss of electrons can have a significant effect on the frontier molecular orbital distribution of both open- and closed-ring isomers, particularly in reactive sites. Corresponding structural changes show a trend toward isomerization. The reaction energy barrier shows greater reduction for dication than monocation and even becomes barrierless for dianion. During the isomerization in different states, the conjugated system switches distinctively, which is attributed to the special redistribution of molecular orbitals and spin population in each state. In monocation and monoanion, for the involvement of a single electron, isomerization is inclined to proceed sequentially between right and left thienyls, whereas it becomes synchronous in dication. The direction depends on the stabilization achieved by the formation of a global conjugated system and more average spin population on the molecule. The effect of substituents on thienyls is demonstrated in the promotion of the extent of conjugation and the determination of the spin population level on the reactive C atoms. Moreover, according to their electron-donating and withdrawing abilities, they can kinetically support or suppress the electron transfer pattern in the process from isomer to transition state, which leads to the control of reaction efficiency.

Theoretical design of solvatochromism switching by photochromic reactions using donor-acceptor disubstituted diarylethene derivatives with oxidized thiophene rings

We have designed several diarylethene derivatives with oxidized thiophene rings and donor-acceptor substituents, which show the solvatochromism switching by photochromic reactions, using a time-dependent density functional theory (TD-DFT) method using the polarizable continuum model (PCM). It is found that in the UV-vis spectral region examined only the open-ring isomers exhibit the solvatochromism, while the closed-ring isomers do not. The mechanism of the solvatochromism behavior and its switching process are clarified from the viewpoint of the charge-transfer (CT) excitation from the donor to the acceptor substituents. We demonstrate that this CT excitation can be controlled by choosing appropriate pairs of the donor and the acceptor substituents on the basis of the orbital correlation diagram between the diarylethene derivatives and the donor-acceptor substituents, which is constructed from the topologies and the orbital energies of the molecular orbitals primarily contributing to the excitations.

Designing efficient photochromic dithienylethene dyads

Aiming at designing more efficient multiphotochromes, we investigate with the help of ab initio tools the impact of the substitution on a series of dimers constituted of two dithienylethene (DTE) moieties, strongly coupled to each other through an ethynyl linker. The electronic structure and the optical properties of a large panel of compounds, substituted on different positions by various types of electroactive groups, have been compared with the aim of designing a dyad in which the three possible isomers (open-open, closed-open, closed-closed) can be reached. We show that appending the reactive carbons atoms of the DTE core with electroactive groups on one of the two photochromes allows cyclisation to be induced on a specific moiety, which leads to the formation of the desired closed-open isomer. Substituting the lateral positions of the thiophene rings provides further control of the topology of the frontier molecular orbitals, so that the electronic transition inducing the second ring closure stands out in the spectrum of the intermediate isomer.

Dye chemistry with time-dependent density functional theory

In this perspective, we present an overview of the determination of excited-state properties of "real-life" dyes, and notably of their optical absorption and emission spectra, performed during the last decade with time-dependent density functional theory (TD-DFT). We discuss the results obtained with both vertical and adiabatic (vibronic) approximations, choosing relevant examples for several series of dyes. These examples include reproducing absorption wavelengths of numerous families of coloured molecules, understanding the specific band shape of amino-anthraquinones, optimising the properties of dyes used in solar cells, mimicking the fluorescence wavelengths of fluorescent brighteners and BODIPY dyes, studying optically active biomolecules and photo-induced proton transfer, as well as improving the properties of photochromes.

Spectral properties and isomerisation path of retinal in C1C2 channelrhodopsin

Computed torsion profiles along the reactive coordinate in S₁reveal a two-path deactivation mechanism for retinal in C1C2 channelrhodopsin.

Theoretical investigation on a series of novel S,S-dioxide diarylethenes with abnormal photochromic properties and design of new dyads

A novel diarylethene shows an abnormal conjugation system switch in a photochromic reaction, according to which new dyads were designed.

Density functional theory study of electric field effects on the isomerization of a photochromic molecular switch based on 1,2-dithienylethene

Structural and electronic properties of 1,2-bis(5-methyl-[2,2′-bithiophen]-4-yl)cyclopent-1-ene in closed form and open form under various external electric field with strengths, 0, 10 × 10−4, 20 × 10−4, 30 × 10−4, 40 × 10−4, and 50 × 10−4a.u., were studied using the DFT-B3LYP/6-31G* method. As a positive index, structural parameters, length of the photoisomers, and the electronic spatial extents are almost stable at different external electric fields. The UV-Vis electronic spectrum based on time-dependent density functional theory indicated that the HOMO → LUMO transition in the closed form under different electric field strengths is strongly allowed, whereas is very weak in the open form. Electronic response parameters such as the HOMO−LUMO gap, electric dipole moment, and polarizability showed that electric conductivity of the closed form at different field strengths is greater than in the open form. Results of electronic density of states show that at high external electric field, the conductivity of the open form and closed form will be probably equal and switching behavior cannot be observed. Isomerization of the closed form to the open form at different external electric fields can be considered as exothermic and spontaneous.

Two-Photon Absorption Spectrum of a Single Crystal Cyanine-like Dye

The two-photon absorption (2PA) spectrum of an organic single crystal is reported. The crystal is grown by self-nucleation of a subsaturated hot solution of acetonitrile, and is composed of an asymmetrical donor-π-acceptor cyanine-like dye molecule. To our knowledge, this is the first report of the 2PA spectrum of single crystals made from a cyanine-like dye. The linear and nonlinear properties of the single crystalline material are investigated and compared with the molecular properties of a toluene solution of its monomeric form. The maximum polarization-dependent 2PA coefficient of the single crystal is 52 ± 9 cm/GW, which is more than twice as large as that for the inorganic semiconductor CdTe with a similar absorption edge. The optical properties, linear and nonlinear, are strongly dependent upon incident polarization due to anisotropic molecular packing. X-ray diffraction analysis shows π-stacking dimers formation in the crystal, similar to H-aggregates. Quantum chemical calculations demonstrate that this dimerization leads to the splitting of the energy bands and the appearance of new red-shifted 2PA bands when compared to the solution of monomers. This trend is opposite to the blue shift in the linear absorption spectra upon H-aggregation.

Femtosecond dynamics of the ring closing process of diarylethene: a case study of electrocyclic reactions in photochromic single crystals

The cyclization reaction of the photochromic diarylethene derivative 1,2-bis(2,4-dimethyl-5-phenyl-3-thienyl)perfluorocyclopentene was studied in its single crystal phase with femtosecond transient absorption spectroscopy. The transient absorption measurements were performed with a robust acquisition scheme that explicitly exploits the photoreversibility of the molecular system and monitors the reversibility conditions. The crystalline system demonstrated 3 × 10(4) repeatable cycles before significant degradation was observed. Immediately following photoexcitation, the excited state absorption associated with the open-ring conformation undergoes a large spectral shift with a time constant of approximately 200 fs. Following this evolution on the excited state potential energy surface, the ring closure occurs with a time constant of 5.3 ps, which is significantly slower than previously reported measurements for similar derivatives in the solution phase. Time resolved electron diffraction studies were used to further demonstrate the assignment of the transient absorption dynamics to the ring closing reaction. The mechanistic details of the ring closing are discussed in the context of prior computational work along with a vibrational mode analysis using density functional theory to give some insight into the primary motions involved in the ring closing reaction.

Efficient photochromic transformation of a new fluorenyl diarylethene: one- and two-photon absorption spectroscopy

Efficient reversible phototransformation of a new diarylethene-fluorene derivative, 1,2-bis(5-(9,9-didecyl-7-nitro-9H-fluoren-2-yl)-2-methylthiophen-3-yl)cyclopent-1-ene (1), was demonstrated in organic media under low-intensity laser excitation. Linear photophysical characterization of 1 was performed at room temperature in solvents of different polarity and viscosity. Significantly, close to unity quantum yield for the cyclization reaction of 1 was shown in nonpolar solutions. The lifetimes of the excited states of the open (OF) and closed (CF) forms of 1 were measured by a femtosecond transient absorption technique, and corresponding values of ∼0.7 and ∼0.9 ps were shown in dichloromethane (DCM), respectively. Degenerate two-photon absorption (2PA) spectra of the OF and CF of 1 were obtained over a broad spectral range by the open aperture Z-scan method under 1 kHz femtosecond excitation. The values of 2PA cross sections of the OF in DCM (∼50-70 GM) were found to increase up to 1 order of magnitude (∼600 GM) after cyclization to the CF. The nature of cyclization and cylcoreversion processes were investigated by quantum chemistry with employment of DFT-based methods implemented in the Gaussian'09 program. The potential of 1 for application in optical data storage was shown using poly(methyl methacrylate)-doped films and two-photon fluorescence microscopy readout.

DFT/TDDFT investigation of the modulation of photochromic properties in an organoboron-based diarylethene by fluoride ions

The diarylethene derivative 1,2-bis-(5'-dimesitylboryl-2'-methylthieny-3'-yl)-cyclopentene (1) containing dimesitylboryl groups is an interesting photochromic material. The dimesitylboryl groups can bind to F(-), which tunes the optical and electronic properties of the diarylethene compound. Hence, the diarylethene derivative 1 containing dimesitylboryl groups is sensitive to both light and F(-), and its photochromic properties can be tuned by a fluoride ion. Herein, we studied the substituent effect of dimesitylboron groups on the optical properties of both the closed-ring and open-ring isomers of the diarylethene molecule by DFT/TDDFT calculations and found that these methods are reliable for the determination of the lowest singlet excitation energies of diarylethene compounds. The introduction of dimesitylboron groups to the diarylethene compound can elongate its conjugation length and change the excited-state properties from π→π* transition to a charge-transfer state. This explains the modulation of photochromic properties through the introduction of dimesitylboron groups. Furthermore, the photochromic properties can be tuned through the binding of F(-) to a boron center and the excited state of the diarylethene compound is changed from a charge-transfer state to a π→π* transition. Hence, a subtle control of the photochromic spectroscopic properties was realized. In addition, the changes of electronic characteristics by the isomerization reaction of diarylethene compounds were also investigated with theoretical calculations. For the model compound 2 without dimesitylboryl groups, the closed-ring isomer has better hole- and electron-injection abilities, as well as higher charge-transport rates, than the open-ring isomer. The introduction of dimesitylboron groups to diarylethene can dramatically improve the charge-injection and -transport abilities. The closed isomer of compound 1 (1 C) has the best hole- and electron-injection abilities, whereas the charge-transport rates of the open isomer of compound 1 (1 O) are higher than those of 1 C. Importantly, 1 O is an electron-accepting and -transport material. These results show that the diarylethene compound containing dimesitylboryl groups has promising potential to be applied in optoelectronic devices and thus is worth to be further investigated.

Probing photochromic properties by correlation of UV-visible and infra-red absorption spectroscopy: a case study with cis-1,2-dicyano-1,2-bis(2,4,5-trimethyl-3-thienyl)ethene

Quantification of the relative composition of the isomers in a photochromic system at any irradiation time interval is a critical issue in determining absolute quantum yields. For this purpose, we have developed a simple and convenient protocol involving combination of UV-visible and infra-red absorption spectroscopy. Photochromic cyclization reaction of cis-l,2-dicyano-l,2-bis(2,4,5-trimethyl-3-thieny1)ethene (CMTE) is analyzed to demonstrate the efficiency of the proposed methodology. This approach is based on the fact that the two isomers show distinctive infra-red bands. Detailed investigations of the UV-visible and infra-red spectra of the mixture obtained at different irradiation times in CCl(4) supported by quantum chemical computations lead to the unambiguous estimation of molar absorption coefficients of the closed isomer (epsilon(CF) = 4650 L mol(-1) cm(-1) at 512 nm). It facilitates the first determination of absolute quantum yields of this reversible photochromic reaction in CCl(4) by fitting the UV-visible spectral data (Phi(OF-->CF) = 0.41 +/- 0.05 and Phi(CF-->OF) = 0.12 +/- 0.02 at 405 nm and 546 nm, respectively).