Misato Suganumaa,
Daichi Kitagawa
*a,
Shota Hamatani
a,
Hikaru Sotome
*b,
Cédric Mittelheisserc,
Michel Sliwa
*cd,
Syoji Ito
b,
Hiroshi Miyasaka
b and
Seiya Kobatake
*a
aDepartment of Chemistry and Bioengineering, Graduate School of Engineering, Osaka Metropolitan University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka 558-8585, Japan. E-mail: kitagawa@omu.ac.jp; kobatake@omu.ac.jp
bDivision of Frontier Materials Science and Center for Promotion of Advanced Interdisciplinary Research, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan. E-mail: sotome@laser.chem.es.osaka-u.ac.jp
cUniversité de Lille, CNRS, UMR 8516, LASIR, Laboratoire de Spectrochimie Infrarouge et Raman, Lille 59000, France. E-mail: michel.sliwa@univ-lille.fr
dLOB, CNRS, INSERM, École Polytechnique, Institut Polytechnique de Paris, Palaiseau 91120, France
First published on 17th January 2025
Photoreactivity in crystals is one of the essential properties for creating photo-functional crystalline materials. This study explores the impact of the dihedral angle in aryl groups on the photocyclization reactivity of inverse-type diarylethenes, both in solution and crystalline phases. By synthesizing various diarylethene derivatives with different dihedral angles, the relationship between structural geometry and photoreactivity is systematically examined. We find that larger dihedral angles between the thiophene and phenyl rings enhance photocyclization reactivity in solution, indicating that destabilized π-conjugation lowers the activation barrier. In fact, ultrafast spectroscopy confirms that the cyclization time constant decreases with larger dihedral angles. In the crystalline phase, X-ray crystallographic analysis shows that all diarylethene derivatives adopt ideally photoreactive anti-parallel conformations, but only crystals with a dihedral angle exceeding approximately 81° exhibit photocyclization. These findings indicate that a certain threshold dihedral angle is essential for photocyclization to occur in crystals. The results of this work provide new insights into the role of molecular geometry in photoreactivity and offer a strategy for designing functional photochromic materials that operate efficiently in the solid state.
With such a background, we focused on the photoreactivity of diarylethene in crystals.42–46 Diarylethene is known as one of the photochromic molecules which reversibly isomerizes between open and closed-ring forms based on a 6π electrocyclic reaction.47 Most of the diarylethene derivatives that have been studied so far have thiophene rings as the aryl groups48–52 and can be classified into two types, normal and inverse, with respect to the orientation of the thiophene ring.53–55 Regarding the photoreactivity of diarylethene in crystals, there is an empirical rule that the normal-type diarylethene derivatives can undergo a reversible photochromic reaction in the crystalline state when the molecules adopt an antiparallel conformation in the crystals and the distance between reactive carbons is less than 4.2 Å, similar to Schmidt's criteria.56,57 In contrast, however, we previously found that most inverse-type diarylethene derivatives cannot undergo photocyclization reactions from the open-ring form to the closed-ring form in crystals even though the molecules satisfy the above reaction conditions in the crystalline state,40,58–61 which are singular cases that do not follow the empirical rule. It has also been confirmed that photoreactivity, i.e. the photocyclization quantum yield, of an inverse-type diarylethene decreases when the viscosity of the solvent is high (from 0.2 in n-hexane to 0.1 in cis-decalin), suggesting that the loss of the photoreactivity in the crystalline state would be due to the property of the molecule itself, rather than being influenced by intermolecular interactions or the size of voids in the crystal.54 It remains an open question as to what controls the photoreactivity of inverse-type diarylethenes in crystals. To gain insight into the answer to this question, we recently focused on the host–guest chemistry of cyclodextrins (CDs), which restricts the molecular geometrical change like in the crystalline phase, and investigated the effect of inclusion into CDs on the photochromic behavior of an inverse-type diarylethene derivative.62 As a result, we found that inclusion of the inverse-type diarylethene into βCD reduced the photocyclization quantum yield from 0.11 to 0.056. A detailed study of the interaction between diarylethene and βCD suggests that the reduced photocyclization reactivity is ascribed to a restriction of the rotational motion of the bond between the thiophene and phenyl rings. Based on this result, we hypothesized that the photoreactivity of inverse-type diarylethenes in crystals depends on the dihedral angle between the thiophene and the phenyl rings.
In this study, we focus on the relationship between the photoreactivity of the inverse-type diarylethene and the dihedral angle. We synthesize various inverse-type diarylethenes (1–6) (Scheme 1) to tune the dihedral angle and examine the photoreactivity in solution and in crystals. We find that the photocyclization quantum yield in solution increases with an increase in the dihedral angle. Ultrafast spectroscopy reveals that the time constant of cyclization tends to decrease with an increase in the dihedral angle, which is interpretable from the viewpoint of the destabilization due to breaking of the π-conjugation. Moreover, we find that the photoreactivity in crystals strongly depends on the dihedral angle and that there is a dihedral angle threshold determining photoreactivity. These results provide new insights into not only the reaction dynamics between hexatriene and cyclohexadiene skeletons but also the molecular design strategy to achieve photoreactivity in the crystalline state.
Open-ring isomer | Closed-ring isomer | Quantum yield | ||||||
---|---|---|---|---|---|---|---|---|
λmax/nm | ε/M−1 cm−1 | φsoln/° | λmax/nm | ε/M−1 cm−1 | φsoln/° | Φo→c | Φc→o | |
a Ref. 63. | ||||||||
1a | 370 | 22![]() |
27 | 438 | 5250 | 31 | 0.17 | 0.48 |
2 | 355 | 19![]() |
41 | 436 | 6820 | 45 | 0.23 | 0.54 |
3 | 350 | 17![]() |
46 | 436 | 6670 | 49 | 0.26 | 0.56 |
4 | 347 | 15![]() |
50 | 435 | 6530 | 51 | 0.30 | 0.60 |
5 | 342 | 14![]() |
72 | 435 | 7440 | 82 | 0.31 | 0.59 |
6 | 340 | 16![]() |
86 | 434 | 7310 | 86 | 0.36 | 0.57 |
It should be noted that the open-ring isomer of diarylethene derivatives is present as reactive and unreactive conformers such as the anti-parallel (AP) and parallel (P) forms, and their relative population affects the apparent cyclization reaction yield. In our previous work, three conformers (AP1, P1, AP2) were demonstrated to be dominant in the ground state of 1o. AP1 is the only reactive conformer, and P1 and AP2 are unreactive.54 We evaluated the relative energies of each conformer using quantum chemical calculation and quantified the relative population of AP1, P1 and AP2 by assuming the thermal equilibrium under the Boltzmann distribution (Fig. S4 and Table S1, ESI†). The relative portion of AP1, P1 and AP2 are respectively 60–65%, 30–35% and ∼5%, which is almost constant in a series of the compounds although 5o shows a slight deviation. This result indicates that the change in Φo→c in 1–6 is less affected by the relative population of the conformers and predominantly originates from the intrinsic cyclization reactivity of the reactive AP form (AP1).
As 6o shows the highest photocyclization quantum yield, we firstly investigated the photocyclization reaction dynamics of 6o using transient absorption spectroscopy with an instrument allowing the recording of transient spectra until 1000 nm as the excited state absorption (ESA) of 1o was reported to extend to the near-infrared region in our precedent work.54 Fig. 2 shows transient absorption spectra of 6o in n-hexane in the femtosecond and later timescales. At the time zero, a positive absorption band was observed at 610 nm concomitantly with a small negative band at around 435 nm. This negative one is ascribable to stimulated emission and is dynamically shifted toward the longer wavelength region within 0.2 ps. It reaches a maximum at 490 nm similar to the steady-state emission maximum shown in Fig. S5 (ESI†), indicating that the molecules undergo structural relaxation on the excited state potential energy surface. The resultant stimulated emission at around 490 nm and ESA at around 635 nm decayed on a picosecond timescale. The transient spectrum at 100 ps is characterized by two positive bands centered at 430 and 610 nm. The first one corresponds to the absorption band of the newly produced closed-ring isomer (strictly speaking, it is a differential spectrum between the open- and closed-ring isomers), and the second one is assigned to the long-lived species of the unreactive conformer such as P1 and AP2. A dip at around 520 nm is attributable to the stimulated emission of these conformers from the spectral agreement with the fluorescence maximum wavelength (Fig. S5, ESI†) and persists from several hundreds of picoseconds to a few nanoseconds. The overall spectral features of 6o are similar to those of 1o with phenyl groups reported in the previous work54 although its timescales become shorter.
As well as for the cyclization dynamics of 1o, increase of transient absorbance corresponding to the formation of the closed-ring isomer was unfortunately elusive due to the spectral overlap with the absorption and stimulated emission of the precursor. Thus, we applied global fitting analysis to the transient absorption data in Fig. 2 to obtain the formation time of the closed-ring isomer, which is essential information for discussion on the topology of the excited state potential. The time evolution of the transient absorption was well reproduced with four components in the 100 ps time window, as shown in Fig. 3a. Similarly to 1o, an oscillatory feature in the first 2 ps is ascribable to coherent vibrations triggered by the pulsed excitation, and detailed analysis with assignment is shown in Fig. S6 and S7 (ESI†). Four decay-associated spectra (DAS) thus yielded are shown in Fig. 3b with the corresponding time constants of 0.12, 2.9, 17 ps and infinity (offset component).
![]() | ||
Fig. 3 (a) Time profiles of transient absorbance of 6o in n-hexane monitored at selected wavelengths. (b) Decay-associated spectra obtained from global fitting analysis. |
The first DAS (0.12 ps) shows negative (<450 nm) and positive components (450–600 nm), which are attributable to the red shift of stimulated emission. A negative signal at 640 nm is also related to the initial geometrical relaxation in the ultrafast timescale. For the second DAS (2.9 ps), a negative component lies at around 440 nm and indicates the formation of a new species and/or the recovery of the stimulated emission. In the present case, this negative feature is well correlated with a positive component at 440 nm in the fourth DAS (offset) showing the absorption of the closed-ring isomer. This spectral agreement shows that the closed-ring isomer is produced with a time constant of 2.9 ps. On the other hand, the third DAS (17 ps) is assignable to the decay of the unreactive conformer in the excited state from the spectral similarity to 1o observed in the previous work.54 The offset component is the superposition of absorption bands due to the closed-ring isomer and long-lived species of the unreactive conformer. In our previous report, time constants were 36 ps (cyclization), 130 and 390 ps for two unreactive (recovery of SE band) conformers. To check if we do not have any extra long time species, we used commercial apparatus (Ultrafast Systems, HELIOS) that allows the measurement of transient absorption until 8 ns (see the ESI†) and found time constant for 6o of 0.11, 2.1, 18 ps and 890 ps, a minor second unreactive conformer (Fig. S9e, ESI†).
We further performed transient absorption measurements of 2o, 3o, 4o and 5o (Fig. S8, ESI†) and extracted the cyclization time constant for investigation of dihedral angle dependence. For 5o, three time constants were needed to model unreactive conformers (7 ps, 26 ps, and 1 ns) and can be rationalized with a tert-butyl group that increases the existence of different conformers. The full datasets of transient absorption data and obtained DAS are described in Fig. S8 and S9 (ESI†). The obtained cyclization timescales of 1o–6o are shown in Table 2, in which the maximum wavelengths of ESA are also shown for comparison. The cyclization time constant decreases with an increase of the dihedral angle reaching a plateau for 5o and 6o. A similar trend can be seen in the dihedral angle dependence of the absorption maximum wavelength of the positive band in 600–700 nm in the transient absorption spectra. These results can be interpreted as follows: approaching the planar geometry with less bulky substituents extends the π-conjugation and leads to stabilization of the energy minimum of the excited state. This is verified by a relaxation time constant of a few ps for 2o and 3o concomitant with the red-shift of the ESA maximum although it is also affected by the excited state of the unreactive conformer. Thus, the activation energy substantially becomes large along the cyclization coordinate toward the energetic barrier, which results in the deceleration of the cyclization process. The comparable time constants in 4o–6o are probably due to sufficiently small activation energy by almost complete breaking of π-conjugation.
Taken together with the above results, we built the cyclization reaction scheme of inverse-type diarylethenes with different dihedral angles. As shown in Fig. 4, photoexcitation first populates the molecule into the Franck–Condon region of the excited state potential energy surface. This excited molecule undergoes slight geometrical relaxation such as the decrease of the dihedral angle, which extends π-conjugation and stabilizes the S1 energy minimum. The calculated molecular geometries theoretically verify that the dihedral angle changes of 1o from the initial 27° to 12° at this minimum.54 From the transient absorption data the timescale of this structural evolution was estimated to be 120 fs. Subsequently, the molecule in the S1 minimum evolves on the excited state potential and overcomes the activation barrier therein. However, this barrier is sufficiently small and the S1 molecule can easily overcome this within a few picoseconds in 6o where the energy stabilization is lower due to shorter π-conjugation. On the other hand, the activation energy of 1o is larger than that of 6o because the S1 minimum of 1o is more stabilized, which leads to the deceleration of the cyclization process (36 ps) and the existence of vibrational relaxation time constant for 2o and 3o. The slower time constant for the cyclization increases the probability of internal conversion, reducing the cyclization quantum yield. Actually, the cyclization reaction of an inverse-type derivative without any aromatic substituents proceeds with a time constant of 0.8 ps according to a previous work64 and this barrierless reaction behavior supports the above consideration. In this way, the dihedral angle of the aromatic groups at the molecular terminals modulates the topology of the excited state potential through the stabilization due to the π-conjugation and plays a crucial role in the evolution of the cyclization process in the solution phase. The slight increase of cyclization quantum yield for 6o in comparison to 5o can be assigned to a lower concentration of unreactive conformer and an increase of internal conversion without the cyclization for 5o due to the tert-butyl group.
![]() | ||
Fig. 4 Schematic representation of the photocyclization reaction dynamics of inverse-type diarylethenes with different dihedral angles. |
Conformation | C–C distance | φcry/°° | Photoreactivity | |
---|---|---|---|---|
a Ref. 40. | ||||
1o-αa | Antiparallel | 3.530/3.441 | 6.42, 23.8/7.18, 8.98 | Inert |
1o-βa | Antiparallel | 3.561 | 10.2, 28.8 | Inert |
2o | Antiparallel | 3.537 | 28.4, 40.79 | Inert |
3o | Antiparallel | 3.536 | 38.09, 41.8 | Inert |
4o | Antiparallel | 3.582 | 47.1, 53.2 | Inert |
5o | Antiparallel | 4.023 | 74.7, 74.7 | Inert |
6o-α | Antiparallel | 3.554 | 87.9, 87.9 | Reactive |
6o-β | Antiparallel | 3.536 | 80.5, 89.8 | Reactive |
7o | Antiparallel | 3.518/3.512 | 85.2, 87.8, 89.3 | Reactive |
To experimentally examine the photochromic reactivity in the crystalline state, the absorption spectral changes of crystals 2o–6o upon UV irradiation were measured. Note that crystals 1o-α and 1o-β are photo-inert as reported previously. The optical microphotographs and the absorption spectra for crystals 2o–6o before and after UV irradiation are shown in Fig. 5b, c and Fig. S11, S12 (ESI†). In the case of crystals 2o–5o, there were no changes in the color of the crystals and the absorption spectra, indicating that the photocyclization reaction did not take place even though they were satisfying the required conditions for the photochromic reaction in the crystalline phase. On the other hand, interestingly, the changes in the color of the crystals and the absorption spectra upon UV irradiation were confirmed for crystals 6o-α and 6o-β as shown in Fig. 5b, c and Fig. S11e, S12e (ESI†). Fig. 5d and Fig. S13 (ESI†) show the polar plots of the polarized absorption spectra at 462 nm on the (001) face for crystal 6o-α and at 455 nm on the (001) face 6o-β, respectively. Note that the face indices were determined by powder X-ray diffraction pattern as shown in Fig. S14 (ESI†). The absorption intensity was dependent on the rotation angle, and there is clear dichroism of the absorption, indicating that the photochromic reaction takes place in the single crystalline phase. Table 3 also summarizes the relationship between the photoreactivity and the φcry. As can be seen, crystals with the φcry less than 74.7° were photo-inert, but crystal 6o-β with the φcry of 80.5° was photo-reactive, indicating that the dihedral angle threshold that determines photoreactivity exists between 74.7° and 80.5°. These results rationalize the previous works that inverse-type diarylethenes bearing hydrogen or methyl groups in place of the aryl ring can undergo photocyclization in solution with a time constant of 0.8 ps and in the crystalline phase, since twisting of the aryl ring is not required.9,64,65
To further verify the effect of the dihedral angle on the photoreactivity in the crystalline state, we synthesized an additional compound 7 shown in Fig. 6a that would have twisted aryl groups. From X-ray crystallographic analysis, it was found that 7o also satisfies the reaction conditions in the crystalline state (the antiparallel conformation in the crystalline state and the distance between the reactive carbons less than 4.2 Å), and the smallest φcry of 7o is 85.2° (Fig. S15 and Table 3, Table S2, ESI†). Upon irradiation with UV light, the crystal color changed from colorless to yellow as shown in Fig. 6b. In addition, the clear dichroism of the absorption was also observed from the polar plots of the polarized absorption spectra shown in Fig. S16 and S17 (ESI†), indicating 7o is photoreactive in the crystalline state. This result follows the relationship between the photoreactivity and the dihedral angle described above, which supports the validity of our explanation of the reaction mechanism of inverse-type diarylethenes.
![]() | ||
Fig. 6 (a) Molecular structure of compound 7 and (b) optical microphotographs of crystal 7o before and after UV irradiation. |
Summing up the results described above, the twisting of the phenyl groups is necessary for inverse-type diarylethenes to overcome the activation barrier in the excited state, resulting in the photocyclization reaction. In particular, it was found that the dihedral angle should be larger than at least 75° to undergo photocyclization as revealed by the investigation of the photoreactivity in the crystalline phase. This characteristic molecular property is the reason why the photoreactivity of inverse-type diarylethene derivatives in crystals does not follow the empirical rule reported previously.
Footnote |
† Electronic supplementary information (ESI) available: Optimized structures obtained by DFT calculations, absorption spectra of 3–6 in n-hexane, the full datasets of transient absorption data and the obtained DASs for 2–5, the detailed data for the investigation of the photoreactivity in the crystalline phase, and the detailed procedures of synthesis (PDF), X-ray crystallographic data (CIF), and cartesian coordinates of the optimized structures (EXCEL). CCDC 2401652–2401657 and 2416015. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d4tc04853a |
This journal is © The Royal Society of Chemistry 2025 |