Keigo
Yano
a,
Takumi
Matsuo
*ab and
Shotaro
Hayashi
*ab
aSchool of Engineering Science, Kochi University of Technology, Kami, Kochi 782-8502, Japan. E-mail: hayashi.shotaro@kochi-tech.ac.jp
bFOREST Center, Research Institute, Kochi University of Technology, Kami, Kochi 782-8502, Japan
First published on 24th March 2025
Robust and mechanically deformable organic molecular crystals present a promising alternative to conventional inorganic and polymeric solids as flexible materials. Here, we report the design of elastic crystals based on asymmetric fluorinated donor–acceptor–donor molecular structures involving 2D planes of heteroatom interactions and pitched π-stacking arrangements. 4,7-Bis(2-thienyl)-5,6-difluoro-2,1,3-benzothiadiazole in crystals exhibited a planar structure with minimized torsion angles (within 2.04°) due to intramolecular heteroatom interactions and hydrogen bonds between the fluorinated benzothiadiazole acceptor and thiophene donor groups. The crystal also formed a two-dimensional noncovalent network along the ac-axis, stabilized by F–H hydrogen bonds between the fluorinated benzothiadiazole-thiophene units and CH–π interactions between thiophene moieties. Furthermore, introducing heteroatoms along the molecular short-axis leads to the formation of pitched π-stack structures. These crystals show macroscopic elasticity for applying stress in each (100) and (001) plane, respectively; thus, the material has 2D elasticity. Micro- and nano-indentation tests revealed that the mechanical robustness (elasticity and hardness) of the crystal surpasses that of common polymeric solids. In addition to their mechanical properties, the crystal exhibited yellow-photoluminescence due to charge-transfer interaction (a maximum emission wavelength of 563 nm). This design of donor–acceptor–donor promotes inter- and intra-molecular charge-transfer interactions. The 2D elastic crystals were easily hybridized with polyimide films and exhibited 3D active optical waveguide properties. This molecular-supramolecular strategy for tuning the mechanical and optical properties provides new opportunities for the engineering of crystal functionalities.
Highly planar π-conjugated molecules offer a narrow HOMO–LUMO energy gap due to the delocalization of π-electrons across the main structure. The crystal structures of such molecules exhibit herringbone or brick-like arrangements due to CH–π and π–π interactions between neighboring molecules (Fig. 1).29–34 In our previous work, we designed molecular stacking structures based on π–π interactions driven by a partial repulsive interaction (δ−–δ−) to develop the first flexible, elastic, molecular crystals of π-conjugated molecules.35–38 Pitched π-stacking structures, exemplified by rubrene crystals—one of the most prominent molecular semiconductors—are known for supporting high carrier mobility, a key advantage of π–π interactions. According to the strategy for elastic π-conjugated molecular crystals, 1,4-bis(2-thienyl)-2,3,5,6-tetrafluorobenzene (BTFB) crystals have been designed and synthesized.38 The crystal structure exhibited high planarity due to intramolecular H–F hydrogen bonds and S–F dipole–dipole interactions; thus, pitched π-stacked structures and H–F hydrogen bonds along the short-axis direction are shown. This crystal structure is key to high elasticity in response to mechanical stress on the (001) plane. However, this crystal is limited to one-dimensional elasticity, and its robustness is inferior to that of polymer materials. To address this issue, we focused on the symmetry of the molecular structure. A decrease in molecular symmetry induces disorder in the intramolecular rotation units. Therefore, we replaced the central benzene unit of BTFB, a C2v group, with an asymmetric skeleton, hoping to control the molecular arrangement and mechanical properties by inducing hydrogen–nitrogen competition in intramolecular and intermolecular interactions.
In this study, we developed 2D elastic and yellow-photoluminescent crystals through the following design strategies: (1) asymmetric donor–acceptor–donor (DAD) π-conjugated structures, (2) control over a two-dimensional network through heteroatom interactions, and (3) “pitched π-stack” structures by minimizing CH interactions via the introduction of heteroatoms along the short-axis direction. As an asymmetric skeleton, fluorinated benzothiadiazole was expected to inhibit the intramolecular interaction with thiophene and the successive hydrogen bonds. DAD structures composed of fluorinated benzothiadiazole (A) and thiophene (D) are expected to exhibit a lower energy absorption band due to intramolecular charge-transfer (CT) interactions, as well as facilitate the formation of pitched π-stacked structures via intermolecular interactions. Here, we demonstrate DAD-type elastic crystals by designing 4,7-bis(2-thienyl)-5,6-difluoro-2,1,3-benzothiadiazole (BTFT). Additionally, we investigate the relationship between the molecular structure, crystal structure, and their corresponding mechanical and optical properties by comparing BTFT crystals with those of BTFB and 4,7-bis(2-thienyl)-2,1,3-benzothiadiazole (BTT).39,40
The single-crystal X-ray diffraction data are summarized in Table 1. BTFT crystals were analyzed in the monoclinic system with Z′ = 4, space group P21/c (CCDC dep. no. 2424463†), and the thienyl unit is disordered (Fig. S8, ESI†). Within the molecule, several intermolecular distances—F–H (2.28 Å), N–H (2.40 Å), F–S (2.70 Å), and N–S (2.83 Å)—are shorter than the corresponding van der Waals radii (dFH = rF + rH = 2.67 Å, dNH = rN + rH = 2.75 Å, dFS = rF + rS = 3.27 Å, and dNS = rN + rS = 3.35 Å) (Fig. 3a). The maximum torsion angle between the benzene and thiophene units is 2.04°, indicating a highly planar molecular structure. These results suggest the presence of intramolecular F–H, N–H, F–S, and N–S interactions, which contribute to its overall planarity. The molecules show π–π stacking in an intermolecular plane (center-to-center distance, 4.78 Å) along the b-axis (Fig. 3b). Additionally, intermolecular F–H interactions (2.60 Å) between the fluorinated benzothiadiazole and thiophene units, as well as C–H interactions (2.81 Å) between thiophenes, establish a two-dimensional noncovalent network via CH–F hydrogen bonds along the a-axis and CH–π interactions between thiophenes along the c-axis (Fig. 3c). Thus, by designing an alternating thiophene–fluorinated benzothiadiazole–thiophene structure, it was possible to achieve molecular planarity through intramolecular interactions and to achieve a pitched π-stacking structure through intermolecular interactions. The low symmetry of the fluorinated benzothiadiazole backbone also facilitates anisotropic intermolecular interactions. In the structures of the BTT and BTFB crystals, intramolecular interactions include N–S (2.95 Å) and N–H (2.46 Å) of BTT and F–H (2.23 Å) and F–S (2.74 Å) of BTFB, respectively (Fig. S9 and S10, ESI†). The maximum torsion angles of both BTT and BTFB molecules were 15.94° and 6.48°, respectively, which indicate a highly planar molecular structure. The intermolecular interactions in both BTT and BTFB crystals show both π–π stacking between molecular planes along the b-axis and CH–π contacts between thienyl–thienyl terminals along the c-axis (Fig. S9 and S10, ESI†). On the other hand, in the intermolecular interactions along the a-axis, hydrogen bonds were formed between F–H in the BTFB crystal; however, no atomic short contact was observed in the BTT crystal (Fig. S9, ESI†).
BTFT | BTT | BTFB | |
---|---|---|---|
a R 1 = Σ||Fo| − |Fc||/Σ|Fo|[β] wR2 = {Σ[w(Fo2 − Fc2)2]/Σ[w(Fo2)2]}1/2. | |||
Empirical formula | C14H6F2N2S3 | C14H8N2S3 | C14H6F4S2 |
Formula weight | 336.39 | 300.40 | 314.32 |
Temperature [K] | 103 | 93 | 103 |
Crystal system | monoclinic | monoclinic | monoclinic |
Space group | P21/c | P21/n | P21/n |
a [Å] | 15.7783(9) | 13.6321(4) | 6.3133(3) |
b [Å] | 4.7778(2) | 5.2997(2) | 5.0761(4) |
c [Å] | 17.6686(9) | 18.5171(6) | 18.6666(5) |
α [°] | 90 | 90 | 90 |
β [°] | 104.731(5) | 111.118(4) | 90.463(3) |
γ [°] | 90 | 90 | 90 |
V [Å3] | 1288.18(11) | 1247.94(8) | 598.19(4) |
Z | 4 | 4 | 4 |
ρ calcd [g cm−3] | 1.734 | 1.599 | 1.745 |
Goodness-of-fit on F2 | 1.054 | 1.067 | 1.086 |
R 1 [I ≧ 2σ(I)]a | 0.0674 | 0.0328 | 0.0384 |
wR 2 (for all data) | 0.1482 | 0.0698 | 0.0905 |
The BTFT crystal shows 2D elastic bending flexibility. Thus, we performed bending tests by applying load on both the (100) plane (Fig. 4a) and the (001) plane (Fig. 4b), respectively, resulting in bending deformation under stress and relaxation along both the c-axis and a-axis. On the other hand, BTFB and BTT crystals showed 1D elastic deformation along the c-axis direction (Fig. S11 and S12, ESI†). To investigate the macroscopic deformation of flexible materials, the measurements of stress (σ), strain (ε), mean elastic modulus (Eave), and mean elastic limit (Slim) were performed via indentation tests of BTFT, BTT, and BTFB crystals (Fig. S13, ESI†).41Eave is the resistance to macroscopic elastic deformation, i.e., the material's robustness. Indentation was performed along the c-axis direction of each crystal, resulting in elastic deformation with Eave = 1.168 GPa for BTFT, Eave = 0.757 GPa for BTT, and Eave = 0.306 GPa for the BTFB crystals (Fig. S14 and Table S1, ESI†). The Slim values along the c-axis were 4.29% for BTFT, 1.54% for BTT, and 3.00% for BTFB crystals. The elastic modulus of the PET resin and polystyrene resin measured with the same apparatus was Eave = 0.572 GPa, and Eave = 0.349 GPa, respectively (Fig. S15 and Table S2, ESI†); thus, these results indicate that BTFT and BTT crystals are capable of undergoing more robust mechanical deformation like the polymeric materials and that BTFT crystals can undergo the most significant deformation. On the other hand, indentation tests along the a-axis direction of each crystal showed that the BTFT crystal exhibited macroscopic elastic deformation, with an elastic modulus of Eave = 0.406 GPa (Fig. S16 and Table S3, ESI†). The BTT and BTFB crystals fractured immediately when indentation was applied, estimating elastic moduli of Eave = 0.526 GPa and Eave = 0.374 GPa, respectively. The Slim values along the c-axis were 7.50% for BTFT, 2.40% for BTT, and 3.27% for BTFB crystals. The results of the BTFT and BTT crystals show high Eave in the elastic direction, and the BTFT crystals have two bending directions with good elastic deformation. Thus, it was suggested that reducing the molecular symmetry makes it possible to adjust the intermolecular interactions and molecular arrangement, which will significantly contribute to controlling the mechanical properties of the crystal.
![]() | ||
Fig. 4 Mechanical deformation of the BTFT crystal. (a) Applied stress with loading on the (100) face. (b) Applied stress with loading on the (001) face. |
Furthermore, to investigate the effect of temperature on the mechanical properties of BTFT crystals, the elastic modulus was evaluated under high-temperature conditions (Fig. S13, ESI†). A soldering iron was placed near the fixed BTFT crystals, and macroindentation tests along the c-axis direction were performed at least three times at heater temperatures of 30 °C, 50 °C, 100 °C, 150 °C, 200 °C, and 250 °C (Fig. S17 and Table S4, ESI†). As a result, the elastic modulus of BTFT crystals decreased from 1.168 GPa to 0.433 GPa in the temperature range of 30 °C to 100 °C. Young's modulus remained constant in the temperature range from 100 °C to 250 °C, although there was some variation in the modulus from room temperature to 100 °C. This suggests that the elastic modulus of BTFT crystals decreases with heating. The DSC results indicate that the crystals exhibit endothermic peaks with melting point (Tm) at 216 °C and crystallization exothermic peaks (Tc) at 171 °C, respectively. No phase transition is observed in this temperature range, suggesting that the effect on mechanical properties is due to temperature-dependent lattice deformation.
To investigate the relationship between the mechanical properties and the structure of the BTFT crystal, energy framework calculations (HF/3-21G) were carried out using Crystal Explorer. The simulation results showed that the total pair energy (ETot) of the intermolecular interactions was dominated by π–π interactions (ETot = −61.5 kJ mol−1) (Fig. 5a). The two π-stacked columns were mainly bonded by F–H interactions (ETot = −17.6 kJ mol−1), London dispersion forces (ETot = −13.7 kJ mol−1), and CH–π interactions (ETot = −9.3 kJ mol−1), forming a non-covalent network that extended to the ac plane. In the BTT and BTFB crystals, the π–π interactions were dominant, as in the BTFT crystals (ETot = −48.3 kJ mol−1) (Fig. 5b). In BTT crystals, the two π-stacked rows were mainly bonded by London dispersion forces (ETot = −25.8 kJ mol−1) and CH–π interactions (ETot = −11.0 kJ mol−1, −12.1 kJ mol−1). In the BTTF crystals, the two π-stacked rows were mainly bonded by F–H interactions (ETot = –18.8 kJ mol−1) and CH–π interactions (ETot = –9.2 kJ mol−1), and both the BTT and BTTF crystals formed relatively weak non-covalent networks on the ac-plane.
It has already been reported that the intermolecular interaction forces between orthogonally arranged molecules play an essential role in accessing elastic deformation on two planes.23,42 The BTFT, BTT, and BTFB crystals exhibited elastic deformation in response to mechanical stress on the (001) plane, forming C–π interactions or London dispersion forces along the c-axis corresponding to the loading direction (ETot = −9.2 kJ mol−1 ∼ −14.2 kJ mol−1). Therefore, we focused on the intermolecular interactions along the a-axis perpendicular to the c-axis. In the BTFT and BTT crystals, weak interactions (ETot = −7.4 kJ mol−1, −11.6 kJ mol−1) and strong interactions (ETot = −17.6 kJ mol−1, −25.8 kJ mol−1) were alternately arranged along the a-axis. This is due to the asymmetric structure of the benzothiadiazole skeleton. On the other hand, in the BTFB crystal, strong interactions between CH–F (ETot = −18.8 kJ mol−1) were continuously formed. Comparing the total interaction energy along the a-axis, the BTFT crystal shows −24.6 kJ mol−1, the BTT crystal shows −32.4 kJ mol−1, and the BTFB crystal shows −37.6 kJ mol−1, revealing that the BTFT crystal, which has an asymmetric structure and is mediated by F–H hydrogen bonds, had the lowest value.
From the above results, it can be inferred that the BTFT, BTT, and BTFB crystals formed attractive intermolecular interactions along the c-axis and a-axis. The molecular design significantly affected the interactions along the a-axis, and the total interaction energy of the BTFT crystal was the lowest. These results suggest that adjusting the crystal hardness by attractive intermolecular interactions along the loading direction is essential for accessing the two bending planes. In addition, mechanical stress along the a-axis caused a fracture in the BTT crystal and vertical cracking in the BTFB crystal (Fig. S11 and S12, ESI†). This indicates that the adjustment of strong and weak ETot by CH–F hydrogen bonds and alternating arrangements significantly contributes to the flexibility of the molecular crystal.
To further investigate the mechanical properties of the BTFT, BTT, and BTFB crystals, elastic modulus EIT and the hardness HIT of each crystal were measured using the nanoindentation method (Table 2). The indentation test of each single crystal using a Berkovich indenter was performed on the (001) plane at least 3 times. Fig. 6a and b display the load–displacement curves and indentation marks on the crystal surfaces of each crystal. The measurement results showed that the EIT and HIT of the BTFT crystal were 7.836 GPa and 0.215 GPa, respectively, under a load of 5.0 mN. When compared with the BTT crystal (EIT = 8.811 GPa and HIT = 0.275 GPa) and BTFB crystal (EIT = 2.476 GPa and HIT = 0.101 GPa), both the elastic modulus and hardness are higher in the order of the BTT crystal, the BTFT crystal, and the BTFB crystal. In addition, the BTFB crystals were found to exhibit extremely low values. Nanoindentation is a measurement method that involves the three-dimensional destruction of materials by pressing an indenter into them. Therefore, it can be predicted that the elastic modulus and hardness of molecular crystals are affected by the smallest ETot among the intermolecular interactions that primarily contribute to the formation of the structure. In fact, among the main contributing interactions in the crystals compared in this study, π–π interactions, CH–F interactions, and CH–π interactions, the smallest ETot was ETot = −7.4 kJ mol−1 in the BTFT crystal, ETot = −11.1 kJ mol−1 in the BTT crystal, and ETot = −6.2 kJ mol−1 in the BTFB crystal, suggesting that the intermolecular interaction with the smallest ETot affects the elastic modulus and hardness of the crystal.
Crystals | E IT (GPa) | H IT (GPa) |
---|---|---|
BTFT | 7.836 ± 0.421 | 0.215 ± 0.009 |
BTT | 8.811 ± 0.863 | 0.275 ± 0.037 |
BTFB | 2.476 ± 0.308 | 0.101 ± 0.006 |
Fig. 7 displays the optical properties and the calculated electron distribution of the molecules. The crystals show photoluminescence under UV irradiation (Fig. 7a). The PL spectra of the crystals showed that the maximum wavelengths of photoluminescence (λPL) are 567 nm for the BTFT crystal, 606 nm for the BTT crystal, and 461 nm for the BTFB crystal, respectively. Compared to the BTFB crystal, the BTT and BTFT crystals showed lower energy photoluminescence, indicating charge transfer interaction due to the DA-type molecular structure consisting of benzothiadiazole A and thiophene D. On the other hand, λPL of the BTFB crystal shows higher energy photoluminescence than the BTT crystals despite a similar structure. This is probably due to the cancellation of the transition dipole by the substituted fluorine atoms. The photoluminescence quantum yields ΦPL were 0.27 for BTFT, 0.07 for BTT, and 0.17 for BTFB crystals. Thus, the BTFT crystal exhibits a lower energy photoluminescence than the BTFB crystal despite the highest photoluminescence quantum yield among the compared crystals due to its DA-type molecular structure design and its aggregated state.
In organic optoelectronics, hybridization with flexible substrates is essential.43,44 A BTFT crystal exhibiting 2D elasticity was fixed to a polyimide film tape with excellent flexibility (tensile strength: 0.428 kg mm−1)45 and optical transparency in the near-infrared region to fabricate a crystal–polymer hybrid (Fig. 8a). When the crystal was deformed into an arch shape, it showed a similar arch-shaped deformation without breaking (Video S1, ESI†). The jig was rotated to the left and right in opposite directions, and the tape was slowly twisted. As a result, it was observed that the crystal showed a bending deformation with a strong inclination (Videos S2 and S3, ESI†). In this case, the BTFT crystal was observed before and after the test, and no cracks or fractures were found in the crystal (Fig. 8b). These results reveal that the BTFT crystal exhibits excellent 3D deformation without breaking when subjected to loads from multiple directions (Fig. S18, ESI†). On the one hand, in a similar test performed on a BTT crystal as a comparison, the crystal was cleaved in multiple places, and no sign of 3D deformation was observed (Fig. S19, ESI†). Next, the waveguiding properties of the BTFT crystal were observed by irradiating light. The BTFT crystal-polyimide hybrid film was fixed in a straight and coil-like bend of 360° (Fig. 8c). When the center of the crystal was irradiated with excitation light, the emission at the irradiated site and the crystal end was observed, and the fluorescence spectrum was measured (Fig. 8c–e). Although there has been much research into flexible organic crystal-based optical waveguides,46–51 in the case of crystal–polymer hybrids, it is expected that the optical signal at the end may not be observed due to light diffusion. For the BTFT crystal fixed on the straight film tape, the emission was observed at the irradiated site and terminal end of the crystal (Fig. 8d). The emission spectrum showed a maximum fluorescence wavelength of 553.9 nm at the irradiated site and 568.9 nm at the crystal end (Fig. 8f). The apparent red-shift of the emission band is due to the fluorescence reabsorption within the crystal. In addition, for the crystal bent on the film at 360° like a coil while tilting, the emissions were observed only at the irradiated site and one end of the crystal (Fig. 8e). The emission spectrum showed a maximum fluorescence wavelength of 547.3 nm at the irradiated site and 594.6 nm at the crystal end (Fig. 8g). Therefore, we found that the crystal deformation can be spatially controlled using a polyimide film and that the BTFT crystal can be deformed in two dimensions. In addition, it was shown that BTFT crystals maintain their optical waveguiding properties not only in a linear state but also when deformed in three dimensions.
Footnote |
† Electronic supplementary information (ESI) available. CCDC 2424463. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d5tc00790a |
This journal is © The Royal Society of Chemistry 2025 |