Keegan
McGehee
ab,
Koichiro
Saito
b,
Dennis
Kwaria
b,
Hiroyuki
Minamikawa
c and
Yasuo
Norikane
*ab
aGraduate School of Pure and Applied Science, University of Tsukuba, Tsukuba, Ibaraki 305-8571, Japan
bResearch Institute for Advanced Electronics and Photonics, National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Ibaraki 305-8565, Japan. E-mail: y-norikane@aist.go.jp
cInterdisciplinary Research Center for Catalytic Chemistry, National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Ibaraki 305-8565, Japan
First published on 8th February 2024
Here we present a study on the solid state properties of trans tetra-ortho-bromo azobenzene (4Br-Azo). Two distinct crystal polymorphs were identified: the α-phase and β-phase. Notably, only the β-phase exhibited an extraordinary photosalient effect (jumping/breaking) upon exposure to a wide range of visible light. Powder X-ray diffraction and Raman spectroscopy revealed that the β-phase is metastable and can transition to the α-phase when subjected to specific stimuli like heat and light. Furthermore, single crystal X-ray diffraction and density functional theory calculations highlighted the significance of a highly strained conformer in the β-phase, showing that the metastability of the phase potentially arises from relieving this strain. This metastability leads to a light induced phase transition, which appears to be the cause of the photosalient effect in these crystals. Interestingly the polymorphism at the core of 4Br-Azo's dynamic behavior is based on different arrangements of halogen based intermolecular interactions. It is possible that continued study on combining visible light capturing chromophores with halogen interaction-based polymorphism will lead to the discovery of even more visible light controlled dynamic crystal materials.
Polymorphism of molecular crystals represents an important area for diversification of dynamic crystal responses. This is because both the physicochemical properties of the constituent molecules, and the way they are able react due to their crystal environment contribute to photomechanical responses. For this reason researchers have shown very different photomechanical responses can be achieved in different polymorphs of the same material.20–22 This makes continuing to study the diverse in bulk photoresponse between polymorphs an important goal for crystal engineering research.
One unique class of motions exhibited by some dynamic crystals are so called salient effects, where crystals jump when stimulated. Salient effects have been reported in response to heat (thermosalient)23 and light irradiation (photosalient).24 Several underlying molecular mechanisms like bond rearrangement,24 photochemical reaction,25–28 and polymorphic phase transitions29,30 have been found to trigger a salient response. Some materials even display salient versus nonsalient polymorphs.20 One thing all these events have in common is they result in a sudden anisotropic shift in the unit cell in those crystals. This is thought to cause significant strain within the crystal that can manifest in a rapid jump or hop motion. Cracking or complete breaking of the crystals is also commonly observed in this process.
Azobenzene derivatives represent one of the most well studied classes of dynamic crystal materials.15,17,31,32 Owing to the photochromic nature of azobenzene they have been found to exhibit motions such as crawling33,34 and bending17,31,35 from light irradiation. Activating these motions often requires UV irradiation to trigger the well-known trans to cis isomerization of azobenzene. Though with the rise in azobenzenes that can be isomerized purely with visible light, researchers have reported photomechanical azobenzene crystals that do not require UV irradiation.17,36 Purely visible light responsive systems are ideal for the widest range of uses because UV light can cause unwanted chemical reactions and damage biological tissue.37 One of the most popular molecular design strategies for creating visible light responsive azobenzenes is substitution of halogen atoms on the ortho positions of the aromatic rings.36,38,39 Though this molecular design strategy has been applied to various solution state systems, there are few reports on the solid-state properties of ortho halogenated azobenzenes.40–42 Thus, further investigation of these materials is expected to be valuable.
Here we present a study on the solid-state properties of tetra-ortho-bromo azobenzene (4Br-Azo, Fig. 1a). The synthesis of this compound has been reported previously,43,44 but, to the best of our knowledge, this is the first solid-state study of the material. During preliminary experiments two distinct type of crystals were observed. They have since been identified as two crystal polymorphs of trans 4Br-Azo, which will be referred to as α-phase and β-phase from now on. Although most recrystallization conditions yielded mixed crystals of the polymorphs, they were successfully crystallized individually by optimizing the solvent and concentration (see ESI†). Due to their striking differences in colour and shape (Fig. 1b) they immediately drew interest for further study. Upon exposing the α-phase to visible light irradiation (filtered Hg lamp), no change was detected. However, within seconds of exposure to the same irradiation the β-phase crystals showed a photosalient effect (Fig. 1c and Video S1, ESI†). This effect was reproducibly observed using three different wavelengths of visible light. To the best of our knowledge this relative independence on incident wavelength is unique among reported photosalient materials.24–28,45–48 Understanding the underlying molecular structure differences between the polymorphs and the reason for the photosalient effect in the β-phase crystals then became an important question. Through characterization with X-ray diffraction (XRD), Raman spectroscopy, and computational chemistry methods we have determined the critical factors to 4Br-Azo's photosalient effect and proposed its molecular mechanism. The findings of this study highlight the fascinating potential in combining halogen interaction-based polymorphism with visible light chromophores.
![]() | ||
Fig. 1 (a) Molecular structure of 4Br-Azo with discussed dihedral angles labelled. (b) Microscope photographs of the two isolated polymorphs of trans 4Br-Azo (left) alongside their respective unit cells viewed normal to the [010] face (right). Scale bars indicate 1 mm. (c) Representative microscope photos of each polymorph before (left) and after (right) strong visible light irradiation (Video S1, ESI†). Dashed circle indicates the irradiated area. Scale bars indicate 200 μm. |
The clearest way to discuss the effect of these irradiation conditions is the photosalient response time, which we define as the time after beginning irradiation when at least one piece of the crystal jumped away from its starting position (Fig. 2a, b and Video S2, ESI†). At least one piece typically travelled on the order of millimeters. Within the intensity range of approximately 600–3900 mW cm−2 it was found that a photosalient response could be reliably produced in less than 30 s and that higher light intensities accelerated the response (Fig. 2c). The deviation in response time also lowers significantly with increased intensity. On the other hand, wavelength supplied appeared to have no clear relationship with the viability or response time of 4Br-Azo's photosalient effect. One possible exception is the measured faster response with green versus blue light at 600 mW cm−2 and 700 mW cm−2 respectively. However, between experimental error in the intensity estimations, uncontrollable variations in sample morphology, and the large standard deviation in photosalient response time at these lower intensities, it is not reasonable to conclude that there is a statistically significant difference in the average response time at those two conditions.
This relative independence of visual light response is well rationalized by the wide visual light absorption band measured for the β-phase (Fig. S1, ESI†). It was also noted that irradiation with 200 mW cm−2 light could induce minor to major cracking of the crystals, suggesting that thinner crystals may exhibit photosalience at these intensities. Only the lowest intensity used (<100 mW cm−2) showed no visible change to the crystals. The observed thickness dependence suggests that this photosalient effect is primarily due to a surface reaction on the β-phase which produces strain that becomes explosive when the bulk crystal is not thick enough to absorb it. The intensity dependence should be directly related to the rate of the photoreaction taking place and follows a similar trend to the earliest report of the photosalient effect.24 As previously noted this wide wavelength response is unique though.
To discuss the geometry differences in these conformations we will highlight the ring twisting dihedral angles (ω) as detailed in Fig. 1a. The aromatic ring twisting dihedral angles of 52.0° in the α-phase conformation are completely symmetrical and are relatively large compared to the nearly planar conformation of unsubstituted azobenzene.49 They are also slightly larger than those reported for tetra-ortho-chloro azobenzene,38 which is a logical trend given the increase in size of the ortho substituent. While the βlow conformation is not completely symmetrical in the same way, the difference in ω values is still small, with ω1 = 62.0° and ω2 = 63.2°. Such small asymmetry in the ring twisting is also common in crystals of azobenzene derivatives.38,40,50
On the other hand, the βhigh conformer stands out among previously reported azobenzene crystals and the other 4Br-Azo conformations. Firstly, to properly model the observed disorder of the azo bond and the aromatic ring closest to the βlow conformer, it was necessary to split the involved atoms into two parts (Fig. 4c). Part 1 (colored by element in Fig. 4c) has 0.55 occupancy and part 2 (magenta atoms in Fig. 4c) has 0.45. This treatment of the βhigh conformer was found to give good quality refinement statistics for the β-phase structure. The need to describe the βhigh conformer in this unique way offers the first indication that it is key to the dynamic behavior of the β-phase, because it may indicate that there is more flexibility of the molecules in these crystal sites compared to the α and βlow sites. Additionally, the twisting angles of the aromatic rings are significantly different from those of the other conformers. When considering part 1 of the βhigh conformation the angles are ω1 = 29.0 and ω2 = 39.8, with ω1 = 59.7 and ω2 = 68.4 using part 2.
In addition to the differences in the molecular conformations present in the two crystal phases there are notable differences in their packing. The α-phase has a tightly ordered structure consisting of end-to-end chains of molecules along the a-axis which have alternating ring stacking orientations (Fig. 3d and e). Sets of these chains align nearly in parallel to create sheets roughly in the (010) plane. Interestingly, these sheets are slightly interdigitated which leads to an alternating two-sheet pattern along the b-axis (Fig. 3c and d). This is significantly different than the β-phase, where the two-fold rotation axis in the [010] direction is lost due to the presence of two conformers. Instead, there is a the two-fold screw axis in the [010] direction, which creates a more complex packing network. Viewed normal to the [010] face (Fig. 4d) this leads to the appearance of sheets made of alternating columns of βhigh and βlow conformers but looking to the c-axis projection (Fig. 4e) it is clear that the interdigitation from the screw axis disrupts any clear sheet-like description of the structure. However, while there is not a sheet analogue within the β-phase structure there is still a chain with similar orientation to the one seen in the α-phase. As can be seen most clearly in Fig. 4a, it has a βhigh–βlow–βlow–βhigh pattern and follows a twisted distortion from the straight chains of molecules seen in the α-phase. Visualizing the differences in these two polymorphs as a collection of linear chains organized into slightly offset sheets (α-phase) versus a collection of twisted chains with alternating twisting direction (β-phase) is helpful for interpreting their relative behaviors. Particularly for the case of a potential polymorphic phase transition as we will discuss in detail later.
Interestingly the only difference in the types of close contact (<sum of VdW radii-0.01 Å) intermolecular interactions present in the two polymorphs is that Br–Br interactions are present in the β-phase, but not the α-phase (Fig. S13 and S14, ESI†). Otherwise, both contain Br–H, Br–C, and C–C (π–π) interactions with no clear differences in magnitude among them. This makes identifying the most dominating interactions in the coordination of each conformer a point of significant interest. To do this Hirshfeld surface analysis51 was employed for each conformer. Hirshfeld surfaces of each conformer are shown in Fig. 5a–c, which have been used to identify the dominating intermolecular interactions by mapping the dnorm property52 on them. The areas of red coloration on these surfaces highlight the regions of most dominant close contact interactions for each conformer. By filtering the regions of the surface by atom type we can identify which types of atom-to-atom contacts these red regions correspond to. This highlights that the α-phase contains some dominant π–π stacking interactions (C5–C9 (3.340 Å)) along with C–Br interactions (C5–Br2 (3.498 Å), C6–Br2(3.514 Å)). These contacts are all shorter than the sum of the vdW radii of the participating atoms by 0.060, 0.052, and 0.036 respectively (Table S4, ESI†). On the other hand, the β-phase is dominated only by Br based intermolecular interactions (Br2–C3 (3.524 Å), Br6–C17(3.219 Å), Br1–Br2(3.446 Å), Br2–Br2(3.618 Å), Br3–H5(2.954 Å), Br5–H9(3.011 Å)). They are shorter than the vdW radii sums by 0.026, 0.331, 0.254, 0.082, 0.096, and 0.039 respectively (Table S5, ESI†). Effectively this allows one to summarize that the relative effects of π–π stacking versus Br contacts lead to the differences in packing for each polymorph.
The Br–Br interactions within the β-phase serve as an important case study for its unique behavior. Based on the commonly accepted geometry-based definitions53 the three unique interactions have been classified as type I or type II (Table S5, ESI†). Specifically, the Br2–Br2 contact is type I (dispersion interaction) while the Br1–Br2 interaction is type II (true halogen bond). There is also a second type II interaction, Br3–Br4 (3.64 Å) which is filtered out from the Hirshfeld surface analysis due to its large length. That contact is only 0.060 less than the sum of Br vdW radii, which is much smaller than the 0.254 difference for the Br1–Br2 contact. This indicates a much stronger interaction in the Br1–Br2 case, which is further supported by interaction angles of 89.9° and 174.6°. These are much closer to the ideal 90°/180° pair than 105.5° and 157.6° as seen for the Br3–Br4 contact. All of these interactions are between two different βlow molecules which appears to provide a backbone of βlow molecules in the unit cell of that phase. The stabilization from these strong interhalogen interactions and the coordination of the βhigh sites that is further made possible by the Br–H interaction between the βlow and βhigh molecules must be key to the relatively high metastability of this phase. Which ultimately leads to its photosalient behavior.
The DSC observations were further confirmed to be related to a monotropic β-to-α transition by Raman and powder (P) XRD measurements (Fig. 6b and c). For the Raman measurements, it was found that the molecular conformation differences in the different polymorphs produce a distinct vibrational signature for each phase in the molecular fingerprint region. Most notably are peaks in the region of coupled azo bond and aromatic ring stretching vibrations (∼1100–1600 cm−1).54 In the α-phase there is a single strong sharp peak at 1498 cm−1. The corresponding area for the β-phase spectrum features two overlapped peaks centered at 1485 cm−1 and 1508 cm−1. It is likely that this difference arises from the presence of two conformers in the β-phase. Also, there is a clear difference in relative intensity of the peaks near 1100 cm−1 for the two polymorphs. For the α-phase the peak is centered at 1070 cm−1 and is noticably weaker than the 1498 cm−1 peak. Meanwhile the 1067 cm−1 centered peak in the β-phase spectrum has a similiar relative intensity to its pair in the 1500 cm−1 region.
We also note that the largest peak in the region of CCC/CCBr bending (∼200–400 cm−1) has comparable relative intensity to the azo/aromatic stretching region in the β-phase. On the other hand, the α-phase features no peaks there with similar intensity to its higher wavenumber regions. By considering the combination of these features it was simple to quickly identify the change of an originally β-phase crystal to the α-phase after heating to 130 °C for 1 minute (Fig. 6b). This spectral shift was further confirmed to be related to a periodic structure change by looking at the PXRD patterns of the two polymorphs before and after the same heating procedure (Fig. 6c). Though it can be challenging to confidently identify polymorphs by their powder patterns due to potential overlapping of peaks it was possible to identify at least two key peaks for each polymorph. For the β-phase we use the (50–) and (102) peaks. Characteristic peaks for the α-phase were chosen as (310) and (130). These peaks were chosen by comparing the relative intensities of the predicted powder patterns (Fig. 6c and Fig. S9, S10, ESI†) with the measured patterns. Using these characteristic peaks, it is clear to see that the β-phase peaks disappeared, and the α-phase peaks emerged after heating the powder.
Given the fact that polymorphic phase transitions have been reported as a potential driving force of some thermosalient effects in the literature,30,55 the observation of this transition in 4Br-Azo did raise the question of whether it might exhibit a thermosalient effect too. However, even under a wide variety of heating conditions no displacement of the β-phase crystals was observed by heating. Only a distinct color change and minor morphology change were observed (Fig. 6a).
Even though no thermosalient effect was observed for β-phase 4Br-Azo the question of whether the metastability of β-phase was related to the photosalient effect remained open. To test for the possibility of a light induced phase transition a thin layer of β-phase powder was placed between two layers of colorless polyimide film and irradiated with a 465 nm LED. In between 5 min periods of irradiation the Raman spectrum of the irradiated region was measured. After about 10 min of irradiation it was found that the Raman spectrum started to take on distinctly α-phase character (Fig. 7b). No additional changes in the spectrum could be detected after 15 min, so it was considered that the sample was at a steady state. The final sample was also measured with XRD for further characterization (Fig. 7c). The PXRD pattern shows characteristic diffraction peaks from the α-phase as expected, but some characteristic β-phase peaks also remain present. Since the α-phase Raman peaks showed higher intensity, it is reasonable that the mixed powder only showed an α-phase signature at a certain point. It is possible that this light induced phase transition occurs due to radiationless decay through the crystals after excitation, effectively creating regions of “hot” molecules, similiar to Förster's description of class I photoreactions.56 This explanation could offer a consistant explanation for both heat and light induced transitions even though the ensemble temperature of the powder did not exceed 45 °C during light irradiation.
In all, the observation of the light induced β to α transition highlights an important aspect of the metastability of the β-phase and presents the light induced transition as a potential driving force for the photosalient effect. Considering the previously discussed packing of the two polymorphs, the β-to-α transition can basically be described as the twisted chains of the β-phase being “pulled” straight. In the heat induced case this would be relatively concerted, so only observing slight expansion and damage (Fig. 6a) is reasonable. On the other hand, the more scattered and gradual light induced transition could create significant stress within single crystals producing the photosalient effect. Even though one would expect the amount of bulk transition on the seconds timescale needed for the photosalient effect to be small previous studies have shown that relatively little photoproduct is required to produce the photosalient effect.26
Conformer | ω 1(°) | ω 2(°) | ϕ(°) | E (kJ mol−1)a |
---|---|---|---|---|
a Single point calculation with B3LYP-GD3//cc-pvtz (relative to trans-opt). βhigh value is the weighted average of calculations considering both parts. b Geometry from B3LYP-GD3//6-311(d) optimization. c Geometry extracted from structures solved by SCXRD. | ||||
trans-optb | 51.7 | 51.7 | 175.4 | 0 |
α | 52.0 | 52.0 | 178.0 | 18.8 |
β low | 62.0 | 63.3 | 176.6 | 9.02 |
β high | Part 1:29.0 | Part 1:39.8 | Part 1:174.6 | 45.9 |
Part 2:59.7 | Part 2:68.4 | Part 2:174.1 |
While molecular geometries offer direct comparison with XRD measurements, the energy values from DFT calculations offer further insight to the 4Br-Azo crystal systems. Single point energy calculations (B3LYP-GD3//cc-pvtz) were performed for each calculated and measured 4Br-Azo geometry. These energy calculations clearly show the extent of the combined distortions on the βhigh molecules. For clarity it is useful to discuss these energies relative to the lowest value, which is trans-opt. The α and βlow geometries have relatively low energies at 18.8 kJ mol−1 and 9.02 kJ mol−1 respectively. In contrast, the βhigh geometry's energy is significantly higher at 45.9 kJ mol−1. Between the geometry distortions and the energy difference, one can clearly see that alleviating the strain on the βhigh molecule must be the driving force behind the observed metastability of the β-phase. The monotropic behavior of this transition is also well rationalized by the geometry and energy considerations because expecting half the molecules within an α-phase crystal to twist back to the βhigh conformation by cooling does not seem reasonable. Instead, the β-phase must only form when distorted molecules get trapped during crystallization to form a metastable solid. This makes quantifying the intermolecular interactions holding the β-phase together a high point of interest. Particularly since one may assume the β-phase to be much less stable than observed if only these single molecule properties are considered. Indeed, one can imagine that this could be a key factor to the photosalient effect in β-phase 4Br-Azo.
To investigate the relative binding energy in the crystal environments of 4Br-Azo, CrystalExplorer2160 (CE) was used. This was done using the interaction energies calculated by CE to compute the lattice energy of each phase through the literature method.61 The calculated lattice energies are −119.563 kJ mol−1 for α-phase and −124.840 kJ mol−1 for β-phase. Firstly, the relatively strong lattice energy of the β-phase is logical given the fact that it is stable until close to the melting point of the compound despite the expected instability of the βhigh conformation. Furthermore, transition from the more tightly bound β-phase to the slightly less bound α-phase is consistent with the endothermic nature of the phase transition peak observed by DSC. The magnitude of the difference between calculated lattice energies is also consistent with the DSC peak integration (∼5.3 kJ mol−1vs. ∼1.2 kJ mol−1). Furthermore, this small difference is consistent with commonly observed differences in energy of organic polymorphs.62
Considering these computational results alongside the experimental evidence a clear picture of what could be driving the photosalient effect in β-phase 4Br-Azo emerged. If one were to imagine the strained βhigh conformation as analogous to a loaded spring trapped by the strong intermolecular interactions present in that crystal structure, then it is reasonable to expect that releasing this spring-like tension could lead to the observed cracking and jumping behavior. Interestingly this does not manifest in a thermosalient effect, as previously discussed. Instead the mixed phase end point of the light-induced phase transition suggests that inducing a phase mixture in a single crystal might be the key factor to the photosalient effect.
Footnote |
† Electronic supplementary information (ESI) available: Detailed conditions for experiments/calculations, full crystallographic details, and additional characterization. CCDC 2295252–2295254. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d3cp04691e |
This journal is © the Owner Societies 2024 |