Wei Ji1, Hui Yuan2,3, Bin Xue4, Sarah Guerin5, Hui Li6, Lei Zhang7, Yanqing Liu8, Linda J W Shimon9, Mingsu Si8, Yi Cao4, Wei Wang4, Damien Thompson5, Kaiyong Cai1, Rusen Yang3, Ehud Gazit2. 1. Key Laboratory of Biorheological Science and Technology, Ministry of Education, The National "111" Project for Biomechanics and Tissue Repair Engineering, College of Bioengineering, Chongqing University, Chongqing, 400044, P. R. China. 2. School of Molecular Cell Biology and Biotechnology, George S. Wise Faculty of Life Sciences, Tel Aviv University, Tel Aviv, 6997801, Israel. 3. School of Advanced Materials and Nanotechnology, Xidian University, Xi'an, 710126, China. 4. National Laboratory of Solid State Microstructure, Department of Physics, Nanjing University, Nanjing, 210093, Jiangsu, China. 5. Department of Physics, Bernal Institute, University of Limerick, Limerick, V94 T9PX, Ireland. 6. Science and Technology on Combustion and Explosion Laboratory, Xi'an Modern Chemistry Research Institute, Xi'an, 710065, China. 7. CAEP Software Center for High Performance Numerical Simulation, Beijing, 100088, China. 8. Key Laboratory for Magnetism and Magnetic Materials of the Ministry of Education, Lanzhou University, Lanzhou, 730000, China. 9. Department of Chemical Research Support, Weizmann Institute of Science, Rehovot, 7610001, Israel.
Abstract
The physical characteristics of supramolecular assemblies composed of small building blocks are dictated by molecular packing patterns in the solid-state. Yet, the structure-property correlation is still not fully understood. Herein, we report the unexpected cofacial to herringbone stacking transformation of a small aromatic bipyridine through co-assembly with acetylated glutamic acid. The unique solid-state structural transformation results in enhanced physical properties of the supramolecular organizations. The co-assembly methodology was further expanded to obtain diverse molecular packings by different bipyridine and acetylated amino acid derivatives. This study presents a feasible co-assembly approach to achieve the solid-state stacking transformation of supramolecular organization and opens up new opportunities to further explore the relationship between molecular arrangement and properties of supramolecular assemblies by crystal engineering.
The physical characteristics of supramolecular assemblies composed of small building blocks are dictated by molecular packing patterns in the solid-state. Yet, the structure-property correlation is still not fully understood. Herein, we report the unexpected cofacial to herringbone stacking transformation of a small aromatic bipyridine through co-assembly with acetylated glutamic acid. The unique solid-state structural transformation results in enhanced physical properties of the supramolecular organizations. The co-assembly methodology was further expanded to obtain diverse molecular packings by different bipyridine and acetylated amino acid derivatives. This study presents a feasible co-assembly approach to achieve the solid-state stacking transformation of supramolecular organization and opens up new opportunities to further explore the relationship between molecular arrangement and properties of supramolecular assemblies by crystal engineering.
Hierarchical supramolecular self‐assembly with long‐range ordered spatial arrangement is a ubiquitous process in nature mediated by multiple non‐covalent intermolecular forces.
In particular, the organization of simple building blocks into ordered supramolecular structures by a bottom‐up approach has gained increasing interest due to the unique properties of the assemblies, along with the advantages of easy preparation and low cost.
Typically, the molecular packing motifs, which dictate the physical properties of supramolecular assemblies, can be divided into four types, namely cofacial, brickwork, herringbone, and slipped stacks.
For example, while cofacial structure with large surface overlap facilitates efficient charge transport and attenuates the transport anisotropy, the herringbone structure has been revealed as an excellent molecular arrangement for very high charge mobilities in semiconductor materials.[
,
] Yet, the structure–property correlations between the molecular packing and physical characteristics are still not fully understood as polymorphism is not easy to avoid in supramolecular assemblies. Pure, solid‐state single crystalline phases are quite challenging to prepare,
and this hinders the successful design and control of properties in the desired materials. To overcome this challenge, supramolecular chemists have devoted much effort to chemically modify the molecular structures via small changes that can directly and drastically affect the molecular packing by intermolecular interactions.
Different solid‐state packing modes were obtained and determined by single crystal structures, precisely establishing the relationships between molecular arrangement and properties of the crystals.
However, this approach usually comprises the disadvantages of intricate organic synthesis, static molecular structures, and non‐adjustable physical properties, which hinders the development of high‐performance organic materials.Recently, supramolecular co‐assembly of multiple component systems inspired by biological materials has emerged as a promising approach to further expand the conformational space of supramolecular assemblies in terms of structural and functional complexity.[
,
] This methodology paves the way to provide an attractive alternative to time‐consuming and inaccessible chemical synthesis and was demonstrated to entail several advantages, such as the control of nanostructure physical dimensions,
adjustable mechanical strength,
the fabrication of non‐canonical complex topologies,
switchable fluorescent emissions,
and the regulation of the helicity and orientation of nanoarchitectures.
In particular, the crystallization of supramolecular co‐assembly in the solid state has been demonstrated to be an ideal way to precisely decipher the structure‐properties relationships and prepare multifunctional and high performance optoelectronic materials.
Although the early stage of co‐assembly from small organic molecules was studied recently,
the basic interaction mechanism and principles underlying the co‐crystallization process are still unclear. Furthermore, not every potential pair of donor and acceptor can co‐assemble and crystallize as the criteria for co‐crystal formation are quite restrictive.
The effective tuning of the physical properties of the co‐crystals by designing and controlling the transformation of molecular stacking modes remains a challenge.Herein, we report an unexpected transformation from cofacial to herringbone stacking of a small aromatic bipyridine (BPY) by co‐assembly with a non‐aromatic acetylated glutamic acid (AcE) (Figure 1). The co‐assembly between BPY and AcE was thoroughly characterised by fluorescence microscopy, powder X‐ray diffraction (PXRD), circular dichroism (CD), Raman spectroscopy, and nuclear magnetic resonance (NMR). The solid‐state stacking modes of supramolecular assemblies were determined by single‐crystal X‐ray diffraction. In the single crystal state, BPY could form intermolecular hydrogen bonding with water molecules resulting in cofacial molecular arrangements. However, by adding AcE as the conformer, BPY preferably formed intermolecular hydrogen bonding with AcE rather than with solvent molecules, leading to herringbone molecular arrangements. The unique stacking transformation resulted in modulation of the physical properties of the co‐crystal, including morphology, fluorescence emission, supramolecular chirality, thermal stability, mechanical property, band structure and piezoelectricity. Furthermore, diverse molecular packing modes were obtained through the co‐assembly between different bipyridine and acetylated amino acid derivatives. This study directly demonstrates the co‐assembly induced supramolecular stacking transformation in the solid state by crystal engineering and promotes the precise understanding of the relationship between molecular packing and properties of supramolecular assemblies.
Figure 1
Schematic presentation of co‐assembly induced stacking transformation from cofacial to herringbone packing in molecular crystals, leading to enhanced physical properties of the crystals.
Schematic presentation of co‐assembly induced stacking transformation from cofacial to herringbone packing in molecular crystals, leading to enhanced physical properties of the crystals.
Results and Discussion
Based on the affinity of intermolecular hydrogen bonding between the carboxylic acid group and the pyridine ring,
a molar ratio of BPY/AcE (1 : 1) was chosen for the co‐assemblies, denoting the equal number of two pyridine and two carboxylic acid groups in the chemical structures of BPY and AcE, respectively. BPY self‐assembled to form long needle‐shaped crystals, which could emit blue fluorescence under UV light excitation due to the aromatic nature of BPY (Figure 2a). Needle‐shaped crystals were also observed in the self‐assembly of AcE, that could not emit fluorescence (Figure 2b). However, bulk irregular needles were observed for the BPY/AcE mixture (Figure 2c), and red‐shifted emission from blue to green was found for the BPY/AcE mixture, suggesting the different packing mode of BPY in the new crystalline state. Compared to the pure BPY and AcE crystals, the appearance of new peaks (11.24°, 28.36°, 30.30°) and the absence of several original peaks (10.47°, 12.41°, 12.74°, 13.72°, 21.62°, 23.61°) was observed in the PXRD of BPY/AcE, indicating a new molecular arrangement (Figure 2d). Furthermore, the supramolecular chirality of the crystal samples was evaluated by CD (Figure 2e). A positive Cotton effect was observed at 226 nm for pure AcE crystals while pure BPY crystals did not produce a CD signal. In contrast, a much higher intensity of ellipticity and different Cotton effects (two positive peaks at 208 and 252 nm, two negative peaks at 288 and 342 nm) were detected in the co‐assembly state, reflecting the enhanced supramolecular chirality of the co‐crystals with high asymmetric supramolecular organizations.
Figure 2
a)–c) Microscopy images of a) BPY, b) AcE, and c) BPY/AcE crystals in bright field (DIC, differential interference contrast, on the left) and in different fluorescence channels, as indicated. d) PXRD analysis of the BPY, AcE, and BPY/AcE crystals. e) CD spectra of the BPY, AcE, and BPY/AcE crystals. f) Raman spectra of the BPY, AcE, and BPY/AcE crystals. g) 1H NMR spectra of BPY, AcE, and BPY/AcE. h) Chemical shift of 1H NMR of BPY/AcE, compared to the single component BPY and AcE.
a)–c) Microscopy images of a) BPY, b) AcE, and c) BPY/AcE crystals in bright field (DIC, differential interference contrast, on the left) and in different fluorescence channels, as indicated. d) PXRD analysis of the BPY, AcE, and BPY/AcE crystals. e) CD spectra of the BPY, AcE, and BPY/AcE crystals. f) Raman spectra of the BPY, AcE, and BPY/AcE crystals. g) 1H NMR spectra of BPY, AcE, and BPY/AcE. h) Chemical shift of 1H NMR of BPY/AcE, compared to the single component BPY and AcE.Raman spectrum and 1H NMR spectra were examined to elucidate the intermolecular hydrogen bonding between the carboxylic acid group of AcE and the pyridine ring of BPY during co‐assembly. The peak at 1693 cm−1 ascribed to the C=O stretching vibration of the carboxylic acid in AcE was absent due to the formation of hydrogen bonding between carboxylic acid and pyridine in the BPY/AcE co‐assembly (Figure 2f). The labels of hydrogen atoms in the chemical structures of BPY (Ha, Hb) and AcE (Hc, Hd, He, Hf) are shown in Figure S1 for the analysis of 1H NMR experiments. An obvious downfield shift was observed for the hydrogen atoms of BPY after co‐assembly with AcE. In contrast, an upfield shift was found for the hydrogen atoms of AcE (Figure 2g, h). These results confirmed the formation of a BPY/AcE complex with strong intermolecular hydrogen bonding between carboxylic acid and pyridine.To better understand the molecular packing modes in the pure crystals and co‐crystal at the atomic level, diffraction quality crystals were grown from methanol and their structures were thoroughly analysed by single crystal X‐ray diffraction.
Detailed crystal data collection and refinement parameters of the BPY crystal and BPY/AcE co‐crystal are summarized in Table S1. BPY crystallized in the monoclinic space group P21 with two BPY molecules and four water molecules in the asymmetric unit (Figure 3a, b, Figure S3). Intermolecular hydrogen bonding was found between the N atom in BPY and H atom in water molecules (Figure S4a). First, each N atom interacted with one water molecule (N1⋅⋅⋅H2A−O2) through hydrogen bond with a N⋅⋅⋅H distance of 2.096 Å, which further connected the next BPY molecule via another water molecule (O2−H2B⋅⋅⋅O1) with a O⋅⋅⋅H distance of 1.894 Å. The other connected intermolecular N2⋅⋅⋅H1B−O1 hydrogen bond length was 2.128 Å, forming a 1D molecular chain. Second, the water‐mediated single molecular chains of BPY further formed 3D crystal networks by cofacial π–π stacking interactions between pyridine rings with distances of 3.705 and 3.712 Å and hydrogen bonds between water molecules with O⋅⋅⋅H distance of 1.929 Å (Figure S4b). Thus, pure BPY molecules could form classic cofacial crystal stacking through the formation of intermolecular hydrogen bonding with solvent molecules (Figure 3c). AcE crystallized in the orthorhombic space group P212121 with one AcE molecule in the asymmetric unit (Figure 3d and S2).
The AcE molecule was directly hydrogen bonded to four adjacent molecules by a total of six hydrogen bonds (Figure 3e and Figure S5). Two cyclic dimer hydrogen bonds were formed by two N−H6⋅⋅⋅O25 and O15−H5⋅⋅⋅O21 contacts with bond lengths of 1.799 (O25⋅⋅⋅H6) Å and 2.105 Å (O21⋅⋅⋅H5), respectively. Furthermore, an additional tight hydrogen bond O11‐H1⋅⋅⋅O5 with a shorter distance of 1.528 Å (O5⋅⋅⋅H1) facilitated the 3D crystal network formation. A layer‐by‐layer arrangement was observed in the higher order crystal packing of AcE (Figure 3f).
Figure 3
a) Chemical structures of BPY and water molecules in the crystals. b) Intermolecular hydrogen bonding formation between BPY and water molecules. c) Higher order of cofacial crystal stacking between BPY and water molecules. d) The AcE molecule asymmetric building unit. e) The formation of hydrogen bonds between adjacent AcE molecules. f) Supramolecular layer‐by‐layer arrangement in the higher order crystal packing of AcE. The two molecular packing layers are denoted in yellow and blue. g) BPY and AcE molecules combine via h) intermolecular hydrogen bonding between the BPY pyridine ring and AcE carboxylic acid groups. i) Higher order of herringbone crystal stacking between BPY and AcE molecules.
a) Chemical structures of BPY and water molecules in the crystals. b) Intermolecular hydrogen bonding formation between BPY and water molecules. c) Higher order of cofacial crystal stacking between BPY and water molecules. d) The AcE molecule asymmetric building unit. e) The formation of hydrogen bonds between adjacent AcE molecules. f) Supramolecular layer‐by‐layer arrangement in the higher order crystal packing of AcE. The two molecular packing layers are denoted in yellow and blue. g) BPY and AcE molecules combine via h) intermolecular hydrogen bonding between the BPY pyridine ring and AcE carboxylic acid groups. i) Higher order of herringbone crystal stacking between BPY and AcE molecules.Compared to the individual crystal structures, co‐crystallization of BPY and AcE resulted in the formation of a new crystal containing one molecule of BPY and one molecule of AcE in the asymmetric unit with a completely distinct molecular arrangement, which crystallized in the monoclinic space group P21 (Figure 3g, Figure S6). As indicated by Raman and 1H NMR experiments, the hydrogen bonds between the carboxylic acid group of AcE and the pyridine ring of BPY were directly observed in the crystal structures (Figure 3h). BPY formed a strong trimer with two molecules of AcE stabilized through two intermolecular hydrogen bonds (O2−H2⋅⋅⋅N2 and O4−H4⋅⋅⋅N3, Figure S7a) with distances of 1.814 Å (N2⋅⋅⋅H2) and 1.817 Å (N3⋅⋅⋅H4), respectively. Compared to the hydrogen bond between BPY and water molecules in the single crystal (N⋅⋅⋅H distance: 2.096 Å and 2.128 Å), the shorter hydrogen bonds between BPY and AcE indicates strong pairing in the co‐crystal. The trimer of BPY with two AcE molecules showed an angle of ≈58° directed by the two intermolecular hydrogen bonds, which propagated to form a classic herringbone stacking mode in the higher order packing of the co‐crystal structure (Figure 3i). No face‐to‐face π–π stacking between aromatic rings was found due to the long centre‐to‐centre distance (5.747 Å) (Figure S7b). Thus, the direct observation of solid‐state cofacial to herringbone stacking transformation was validated in the crystal structure of BPY through co‐assembly with AcE. We further performed density functional theory (DFT) calculations to visualise the valence charge density distribution and Hirshfeld surface of the crystals using the HASEM code. The result and analysis were shown in Figure S8.The thermal stability of the single‐ and co‐crystals was evaluated by thermogravimetric analysis (TGA) and differential scanning calorimetry (DSC) analyses (Figure S9a–c). The TGA results revealed that the BPY and AcE single crystals remained stable until the temperature reached 107.5°C and 191.5 °C, respectively. After co‐crystallization, the BPY/AcE crystal degraded only after the temperature reached 155.6 °C. Endothermic peaks were observed in the DSC curve: BPY (70.3 °C, 112.1 °C, 188.0 °C), AcE (192.8 °C, 235.5 °C), and BPY/AcE (171.7 °C, 223.1 °C). The morphological changes and thermal events of the crystals were further investigated using hot‐stage optical microscopy (HSOM) to visually examine thermal transitions (Figure S9d–f). Melting of the three crystals was observed at a decomposition temperature range of approximately 111.2–112 °C, 193–200.4 °C, and 162–172 °C for BPY, AcE, and BPY/AcE, respectively. Thus, the thermal stability order is [AcE>BPY/AcE>BPY], suggesting that co‐assembly with the more thermally‐stable conformer of AcE, producing a new herringbone packing, could improve the thermal stability relative to pure BPY with a cofacial packing.Investigating the mechanical properties of molecular crystals has long been a subject of considerable interest to chemists due to their versatile applications in pharmaceutical technology, biophysics, polymer science, and optoelectronic devices.
To examine the effect of molecular packing modes on the mechanical properties of the crystals, atomic force microscopy (AFM)‐based nanoindentation experiments were performed by applying the AFM cantilever tip to the surface of the crystalline needles. Typically, the surface of the crystals was scanned using the AFM tapping mode to locate the smaller size flat area (Figure S10). Then the cantilever was forced to press the flat area and retracted at a constant speed of 30 μm s−1 (Figure 4a), during which the force‐displacement traces were recorded. The Young's modulus of the crystals at different positions could be determined by fitting the retracting traces using the Hertz model (Figure 4b and Figure S11). Both Young's modulus and point‐stiffness were evaluated to study the mechanical rigidity. The topographic Young's modulus maps showed the Young's moduli for the pure crystals of BPY (Figure 4c), AcE (Figure 4d), and the BPY/AcE co‐crystal (Figure 4e). The histogram distribution of Young's modulus was further analysed, showing the statistical values of Young's modulus along the thickness direction to be 15.3±4.9 GPa (Figure 4f), 16.6±2.6 GPa (Figure 4g), and 83.4±10.7 GPa (Figure 4h) for BPY, AcE, and BPY/AcE, respectively. Moreover, the value of point stiffness could also be calculated according to the force‐displacement traces (Figure S12). The statistical point stiffness values of single‐ and co‐crystals were 188±27 N m−1 (Figure 4i), 200±36 N m−1 (Figure 4j), and 546±61 N m−1 (Figure 4k) for BPY, AcE, and BPY/AcE, respectively. Thus, the rigidity order is [BPY/AcE>AcE>BPY], indicating an enhanced mechanical stiffness by co‐crystallization due to the tightly‐packed supramolecular organization of herringbone stacking through the combination of intermolecular interactions.
Figure 4
a) Schematic diagram of AFM indentation experiments performed on crystal samples. b) Typical force‐displacement traces on the BPY/AcE co‐crystal. c)–e) Topographic modulus maps of c) BPY, d) AcE, and e) BPY/AcE. Scale bar is 100 nm. f)–h) The statistical Young's modulus distributions of f) BPY, g) AcE, and h) BPY/AcE. i)–k) The statistical point stiffness distributions of i) BPY, j) AcE, and k) BPY/AcE.
a) Schematic diagram of AFM indentation experiments performed on crystal samples. b) Typical force‐displacement traces on the BPY/AcE co‐crystal. c)–e) Topographic modulus maps of c) BPY, d) AcE, and e) BPY/AcE. Scale bar is 100 nm. f)–h) The statistical Young's modulus distributions of f) BPY, g) AcE, and h) BPY/AcE. i)–k) The statistical point stiffness distributions of i) BPY, j) AcE, and k) BPY/AcE.We further explored the electric band structures and density of states (DOS) properties of the three crystals by first‐principle calculations. The corresponding 3D packing structures in different directions are shown in Figure S13. The first Brillouin zones of the primitive cells of the three crystals are displayed in Figure S14. The calculated electronic band structures and the corresponding DOS of the three crystals as shown in Figure 5a and Figure S15 demonstrate the indirect band gap of the BPY/AcE co‐crystal to be approximately 2.96 eV, in the semiconductor range. In contrast, both single crystals of BPY and AcE showed larger indirect band gap values of approximately 4.79 eV and 4.94 eV, respectively. The plotted orbital charge density of conduction band minimum (CBM) and valence band maximum (VBM) surfaces (Figures S16–S18) show that the VBM and CBM of BPY/AcE were both not located at high symmetry points. The charge densities of CBM were only located in the structures of the BPY molecules, while in the VBM, the charge densities were distributed in the structure of the AcE molecules (Figure 5b). The intermolecular binding of BPY with AcE decreased its band gap relative to the individual components, suggesting potential future application for the BPY/AcE co‐crystal as an indirect band gap semiconductor.
Figure 5
a) Band structures and density of states of the BPY/AcE co‐crystal. b) Computed orbital charge densities of CBM and VBM bands of the BPY/AcE co‐crystal along the a axis. c)–e) Optimized crystal structures of c) monoclinic BPY, d) orthorhombic AcE, and e) monoclinic BPY/AcE. f) Predicted absolute piezoelectric strain constants, d of all three crystals. g) Comparison of the maximum piezoelectric coefficients of selected different biological (brown) and non‐biological (blue) materials, as well as BPY, AcE, and BPY/AcE crystals (pink). 1: BPY, 2: AcE, 3: BPY/AcE, 4: lysozyme, 5: gelatin nanofiber, 6: chitosan, 7: asparagine, 8: aluminium nitrate, 9: poly (lactic acid), 10: GaN nanowires, 11: PVDF. The relevant references are shown in Table S6. h) Schematic presentation of the preparation and sandwich structure of the piezoelectric nanogenerators. i) Piezoelectric energy conversion process of the sandwich‐structured generator by press and release using molecular crystals as active components. j), k) Open‐circuit voltage of j) AcE and k) BPY/AcE crystal‐based generators obtained by applying a force of 31 N. l) Linear fitting of the open‐circuit voltage of AcE and BPY/AcE crystal‐based nanogenerator. m), n) Short‐circuit current of m) AcE and n) BPY/AcE crystal‐based generator obtained by applying a force of 31 N. o) Linear fitting of the short‐circuit current of AcE and BPY/AcE crystal‐based nanogenerator.
a) Band structures and density of states of the BPY/AcE co‐crystal. b) Computed orbital charge densities of CBM and VBM bands of the BPY/AcE co‐crystal along the a axis. c)–e) Optimized crystal structures of c) monoclinic BPY, d) orthorhombic AcE, and e) monoclinic BPY/AcE. f) Predicted absolute piezoelectric strain constants, d of all three crystals. g) Comparison of the maximum piezoelectric coefficients of selected different biological (brown) and non‐biological (blue) materials, as well as BPY, AcE, and BPY/AcE crystals (pink). 1: BPY, 2: AcE, 3: BPY/AcE, 4: lysozyme, 5: gelatin nanofiber, 6: chitosan, 7: asparagine, 8: aluminium nitrate, 9: poly (lactic acid), 10: GaN nanowires, 11: PVDF. The relevant references are shown in Table S6. h) Schematic presentation of the preparation and sandwich structure of the piezoelectric nanogenerators. i) Piezoelectric energy conversion process of the sandwich‐structured generator by press and release using molecular crystals as active components. j), k) Open‐circuit voltage of j) AcE and k) BPY/AcE crystal‐based generators obtained by applying a force of 31 N. l) Linear fitting of the open‐circuit voltage of AcE and BPY/AcE crystal‐based nanogenerator. m), n) Short‐circuit current of m) AcE and n) BPY/AcE crystal‐based generator obtained by applying a force of 31 N. o) Linear fitting of the short‐circuit current of AcE and BPY/AcE crystal‐based nanogenerator.The non‐centrosymmetric structure of the three crystals suggests piezoelectric responses owing to the supramolecular internal polarization. We predicted the piezoelectric properties of the crystals by DFT calculations. The optimised structures of the three crystals used for prediction are shown in Figure 5c–e. The orthorhombic symmetry of the AcE crystal allows for three shear piezoelectric tensor components only, compared to eight for the monoclinic BPY and BPY/AcE crystals (Tables S3–S5). A slight increase in the piezoelectric charge constant e23 was observed in the co‐crystals with a value of 0.03 C m−2. The maximum piezoelectric strain constants (d
max) did not increase with co‐crystallisation, showing values in‐between those of BPY and AcE crystals (Figure 5f). The highly flexible AcE single crystal showed the highest piezoelectric response, with a large predicted d
25 value of 42.9 pC N−1. The AcE crystal has favourable packing of its carboxyl groups that facilitates strong intermolecular interactions in all directions, and thus large ionic displacement under an applied shear force. The maximum response of the co‐crystal is significant, with d25 of 22.4 pC N−1, and is slightly higher than the highest maximum response in BPY single crystals. The d
max values of the three crystals were comparable to different biological and non‐biological materials (Figure 5g and Table S6), including biological materials such as lysozyme (6.5 pC N−1), gelatin nanofiber (20 pC N−1), chitosan (8 pC N−1), and asparagine (12 pC N−1), and non‐biological materials such as aluminium nitrate (5 pC N−1), poly (lactic acid) (10 pC N−1), GaN nanowires (12.8 pC N−1), and PVDF (22 pC N−1). The relevant references are shown in Table S6. Furthermore, all three crystals showed significant piezoelectric voltage constants, indicating their potential for energy harvesting applications. The calculated g
max values were 1151, 1670, and 784 mV m/N for BPY, AcE and BPY/AcE, respectively, higher than that of lead zirconium titanate (PZT) and many of its lead‐free counterparts (g
max=250 mV m N−1), the piezoelectric materials most widely exploited for sensing and actuation.
A summary of the knowledge on the structural and chemical features that modulate piezoelectric performance in molecular co‐crystals was shown in Supporting Information (Page S20).We further studied the piezoelectric properties of the three crystals powders by designing and fabricating sandwich power generators. The schematic of a fabricated crystal material‐based nanogenerator and the energy conversion process are depicted in Figure 5h and i. The applied force‐dependent open‐circuit voltages and short‐circuit currents of crystal film‐based nanogenerators were both tested in the pressing force range from 0 N to 31 N. By applying a pressing force of 31 N, the values of open‐circuit voltage reached up to 0.42 V, 0.73 V, and 0.47 V for BPY, AcE and BPY/AcE, respectively, confirming the DFT‐predictions of significant voltage output from the co‐crystal and highest piezoelectric response of AcE single crystals (Figure 5j–k and Figures S19–S21). Particularly, the applied force‐dependent open‐circuit voltages from 4 N to 31 N showed a good linear fit (Figures 5l, S19f, S20f, and S21f), indicating the stable piezoelectric response of the three crystals for the voltage signals. These values exceed the value achieved by other piezoelectric nanogenerators based on inorganic materials (e.g. ZnO, monolayer MoS2) and biological materials (e.g. M13 bacteriophage).
Moreover, the open‐circuit voltage produced by the power nanogenerator maintained at a constant level at different pressing and releasing points, suggesting useful dielectric property and negligible current leakage during the energy conversion process. By applying a pressing force of 31 N, the values of short‐circuit currents were found to be 0.57 nA, 17.40 nA, and 27.86 nA for BPY, AcE and BPY/AcE, respectively (Figure 5m, n and Figures S22–S24). Furthermore, a good linear fit was also observed for the applied force‐dependent short‐circuit currents (Figure 5o, Figures S22f, S23f, and S24f), indicating the piezoelectric response of these molecular crystals for the current output. The BPY/AcE co‐crystal generated the highest short‐circuit current output due to its lower resistance, which is in good agreement with the calculated band structures and DOS profiles of the crystals. These results revealed the potential application of the fabricated nanogenerators based on organic molecular crystals in energy harvesting with tunable piezoelectricity.The criteria for co‐crystal formation are quite restrictive since not every pair of donor and acceptor can co‐assemble and crystallize, thus, we further changed the chemical structures of BPY and AcE into different bipyridine and acetylated amino acid derivatives and investigated whether more co‐crystal structures can be obtained by using this co‐assembly methodology (Figure S25). From these experiments, nine diffraction‐quality co‐crystals were grown from methanol and their structures were thoroughly analyzed, including APY/AcA, BPY/AcD, BPY/AcN, BPA/AcA, BPA/AcD, 2,2′‐BPE/AcE, BPE/AcA,
BPE/AcD, and BPE/AcE, demonstrating the generality of this co‐assembly methodology. Detailed crystal data collection and refinement parameters of these co‐crystals are summarized in Supporting Information Tables S7 and S8. The hydrogen bonds and higher order of crystal stacking between bipyridine and acetylated amino acid molecules were analysed and are shown in Figure S26 and Figure 6. Intermolecular hydrogen bonds between the pyridine ring of bipyridines and the carboxylic acid group of acetylated amino acids was observed in all the crystal structures as the driving force for co‐assembly. Different crystal packing modes were generated by the co‐crystallization of the two types of molecules denoted by different colours (Figure 6). From this detailed investigation of the molecular platform, we found that bipyridine derivatives like BPY, BPA, BPE, showed a good tendency to form co‐crystals with acetylated amino acids (AcA, AcD, AcE). Interestingly, strong AcA dimers stabilized through two face‐to‐face intermolecular hydrogen bonds were observed in all AcA‐based co‐crystals (APY/AcA, BPA/AcA, BPE/AcA), which acted as the bridge to cage the bipyridine molecules and further produced similar supramolecular packing modes (Figure S26a, S26d, S26 g and Figure 6a, d, g). Compared to other acetylated amino acids, the molecular structures of AcD and AcE have two carboxylic acid groups with different geometry, increasing the possibility of molecular co‐assembly between the pyridine ring and the carboxylic acid to form higher‐order supramolecular organizations (Figure S26b, S26e, S26f, S26h, S26i and Figure 6b, e, f, h, i).
Figure 6
Higher order of different crystal stacking between bipyridine and acetylated amino acid molecules viewed along their a axes: a) APY/AcA, b) BPY/AcD, c) BPY/AcN, d) BPA/AcA, e) BPA/AcD, f) 2,2′‐BPE/AcE, g) BPE/AcA, h) BPE/AcD, and i) BPE/AcE. The different colours indicate different molecules to more clearly depict the different molecular packing.
Higher order of different crystal stacking between bipyridine and acetylated amino acid molecules viewed along their a axes: a) APY/AcA, b) BPY/AcD, c) BPY/AcN, d) BPA/AcA, e) BPA/AcD, f) 2,2′‐BPE/AcE, g) BPE/AcA, h) BPE/AcD, and i) BPE/AcE. The different colours indicate different molecules to more clearly depict the different molecular packing.
Conclusion
In summary, we present the co‐assembly induced solid‐state cofacial to herringbone stacking transformation, as validated by X‐ray crystallography, using very simple molecular building blocks. The effect of supramolecular packing on the physical properties of the resulting crystals is thoroughly investigated by various experimental methods and DFT calculations. Compared to BPY hydrate crystals with cofacial stacking, co‐crystallization with AcE molecules to form a new herringbone stacking resulted in modified physicochemical properties of the crystals, such as red‐shifted fluorescence emission, enhanced supramolecular chirality, improved thermal stability, increased mechanical strength and reduced band gap. Moreover, high values of force‐dependent open‐circuit voltage and short‐circuit current were generated by fabricating sandwich piezoelectric devices using AcE and BPY/AcE crystals, respectively, making them promising candidates for energy harvesting. Moreover, the co‐assembly methodology was expanded to obtain diverse molecular packings using different bipyridine and acetylated amino acid derivatives. This study not only demonstrates a feasible co‐assembly strategy to regulate the solid‐state packing transformation of minimal molecular building blocks into supramolecular organizations with enhanced physical properties, but also provides new insights to precisely establish the structure–property correlations in multiple‐component crystals and further inspire the rational molecular design of high‐performance organic optoelectronic materials for next‐generation sustainable and energy‐efficient technologies.
Experimental Section
Experimental procedures and supplementary data were described in the Supporting Information.
Conflict of interest
The authors declare no conflict of interest.As a service to our authors and readers, this journal provides supporting information supplied by the authors. Such materials are peer reviewed and may be re‐organized for online delivery, but are not copy‐edited or typeset. Technical support issues arising from supporting information (other than missing files) should be addressed to the authors.Supporting InformationClick here for additional data file.Supporting InformationClick here for additional data file.Supporting InformationClick here for additional data file.Supporting InformationClick here for additional data file.Supporting InformationClick here for additional data file.Supporting InformationClick here for additional data file.Supporting InformationClick here for additional data file.Supporting InformationClick here for additional data file.Supporting InformationClick here for additional data file.Supporting InformationClick here for additional data file.Supporting InformationClick here for additional data file.
Authors: Ayala Lampel; Scott A McPhee; Hang-Ah Park; Gary G Scott; Sunita Humagain; Doeke R Hekstra; Barney Yoo; Pim W J M Frederix; Tai-De Li; Rinat R Abzalimov; Steven G Greenbaum; Tell Tuttle; Chunhua Hu; Christopher J Bettinger; Rein V Ulijn Journal: Science Date: 2017-06-09 Impact factor: 47.728