Sebastiano
Romi
a,
Samuele
Fanetti
*ab,
Frederico
Alabarse
c,
Antonio M.
Mio
d and
Roberto
Bini
*abe
aLENS, European Laboratory for Non-linear Spectroscopy, Via N. Carrara 1, I-50019, Sesto Fiorentino, Firenze, Italy. E-mail: fanetti@lens.unifi.it; roberto.bini@unifi.it; Fax: +390554572489; Tel: +390554572436
bICCOM-CNR, Istituto di Chimica dei Composti OrganoMetallici, Via Madonna del Piano 10, I-50019, Sesto Fiorentino, Firenze, Italy
cELETTRA, Elettra Sincrotrone Trieste S.C.p.A, in AREA Science Park 34149 Basovizza, Trieste, Italy
dIMM-CNR, Istituto per la Microelettronica e Microsistemi, VIII Strada 5 - Zona Industriale 95121, Catania, Italy
eDipartimento di Chimica “Ugo Schiff”, Università di Firenze, Via della Lastruccia 3, I-50019, Sesto Fiorentino, Italy
First published on 13th April 2021
Carbon nanothreads are likely the most attracting new materials produced under high pressure conditions. Their synthesis is achieved by compressing crystals of different small aromatic molecules, while also exploiting the applied anisotropic stress to favor nontopochemical paths. The threads are nanometric hollow structures of saturated carbon atoms, reminiscent of the starting aromatic molecule, gathered in micron sized bundles. The examples collected so far suggest that their formation can be a general phenomenon, thus enabling the design of functionalities and properties by suitably choosing the starting monomer on the basis of its chemical properties and crystal arrangement. The presence of heteroatoms or unsaturation within the thread is appealing for improving the processability and tuning the electronic properties. Suitable simple chromophores can fulfill these requirements and their controlled insertion along the thread would represent a considerable step forward in tailoring the optical and electronic properties of these mechanically extraordinary materials. Here, we report the synthesis and extensive characterization of double core nanothreads linked by azo groups. This is achieved by compressing azobenzene in a diamond anvil cell, the archetype of a wide class of dyes, and represents a fundamental step in the realization of nanothreads with tailored photochemical and photophysical properties.
The nanothread formation from benzene crystals has been fully rationalized identifying the crucial non-topochemical reaction step as due to the uniaxial stress which can be realized in the DAC avoiding the use of any hydrostatic compression medium.15,16 When the direction of the applied stress coincides or has an important component along the crystal direction where the rings interact in a parallel configuration, the formation of C–C interring bonds along the stack becomes possible. This reaction path is sharply evidenced by the alignment of the resulting nanothreads along the applied stress direction.10,16 The threads arrange parallel in a pseudohexagonal close-packed 2D lattice (a = 6.5 Å) characterized by a correlation length of several tens of nm, but definitely lower, ≤2.5 nm, along the threads.17,18 A consequence of the mechanism described above leading to the nanothread formation is that the reaction cannot be quantitative, competing with the topochemical one which leads to the formation of amorphous hydrogenated carbon (a-C:H).19,20
After the discovery of nanothread formation from benzene, computational and experimental studies revealed the possibility of synthesising nanothreads from a large number of aromatics and heteroaromatics. Ab initio DFT calculations and classical MD simulations suggested the formation of fully saturated carbon nanothreads from polycyclic aromatic hydrocarbons21 and from substituted benzenes (–CH3, –NH2, –OH, and –F).22 The presence of such groups is extremely attracting in view of exploiting them for functionalization or for tuning the electronic properties without reducing the mechanical properties which are not substantially altered upon functionalization with respect to benzene derived nanothreads. The presence of unsaturation can also play an important role in tuning the electronic properties by changing the gap from the UV to the near infrared, but obviously decreasing both strength and stiffness with respect to fully saturated nanothreads.23 Experimentally, nanothreads have been synthesized from substituted aromatics and heteroaromatics like pyridine,24,25 aniline26 and thiophene,27 and recently also from arene-perfluoroarene co-crystals28,29 and highly strained fully saturated molecules like cubane.30 In spite of the general similarity regarding the mechanism of the nanothread formation, consistent differences in the synthesis conditions, in the thread packing and in the reaction yield can characterize the different systems and are ascribable to the crystal arrangement of the starting monomer. Temperature is an important parameter in these reactions because through its change we can modulate the amplitude of the molecular motions in the crystal. With increasing the temperature above 400 K, the reaction threshold in pyridine has been lowered to 14 GPa improving the quality and yield (≥30%) of the synthesized nanothreads and also modifying the packing arrangement.25 The importance of temperature in the nanothread formation is even more pronounced in the case of aniline26 and cubane,30 where temperatures in excess of 500 K are required for inducing the reaction in spite of pressure conditions around or exceeding 30 GPa. In the cubane case, an excess of energy is required to release the internal strain of the molecule,30 whereas in aniline the reaction threshold experimentally observed, 550 K and 33 GPa,26 nicely fits with a topochemical mechanism where the reactivity is driven by a complex cooperation between lattice motions19,31 and H-bonding32,33 that causes a strongly anisotropic crystal response to compression and favors inter stack bonding while preserving the NH2 groups. With further increasing the temperature and reducing the pressure, the increased amplitude of both thermal oscillation and reaction cavity35 favors the topochemical path, which leads to the formation of amorphous 2D or 3D carbon nitrides both in the case of aniline34 as well as pyridine.25
The presence of heteroatoms or substituents within the threads can be exploited for tuning the electronic properties or to act as potential linkers and active sites for doping. Aromatic or heteroaromatic molecules containing chemically stable or well protected functional groups represent appealing systems for the realization of diamond-like nanothreads where the optical, thermal, electrical and conductivity properties can be suitably adjusted by choosing appropriate substituents. From this point of view, the azobenzene family is extremely attracting in view of the efficient tuning of their absorption, emission and photochemical properties by changing the substituents in the phenyl rings.36 This is the reason of their widespread use as dyes even though the most important characteristic is the efficient photo-isomerization, which makes them attracting building block structures for realizing molecular devices and functional materials,37 molecular photoswitches38–40 and even graphene composites.41 Azobenzene, the simplest member of this family, is characterized, in its trans form, by an excellent chemical stability of the crystal. The structural data, available only at ambient pressure, indicate a monoclinic P21/a structure with two molecules in the asymmetric unit cell (see ESI1†).42 Both molecules lie at inversion centers with one of them showing dynamic orientational disorder.43 Following the behavior observed on cooling, the disorder is expected to reduce with increasing the pressure as indeed supported by Raman experiments under hydrostatic conditions.44 The shortest distance between the centers of mass of equivalent molecules, an arrangement suggestive of a potentially selective Diels–Alder [4 + 2] polymerization, is realized along the b axis and it is 5.756 Å at ambient pressure. The formation of double nanothreads, in which the NN chromophore is preserved and protected, let us envisage the possibility to synthesize diamond-like bundles of nanothreads, where the azo groups, eventually present in low concentration by suitably choosing the starting mixed crystal and the ring substituents, can act as efficient light absorbers giving rise to a nanomaterial that combines the mechanical properties of carbon nanothreads with the tunable optical properties of azocompounds to realize functional ultrahard coatings or ultrafast heat sinks.
In this manuscript we report the formation of double nanothreads by compressing trans-azobenzene crystals. Different pressure and temperature reaction paths have been explored to identify the best conditions for producing high quality bidimensional crystalline packing of double nanothreads linked by unaltered azo groups.
IR absorption spectra have been recorded at ambient temperature as a function of pressure up to 31.5 GPa, and at constant pressure (10 and 15 GPa) as a function of temperature. In all these cases, we detected a chemical transformation identifying the reaction threshold conditions and monitoring the kinetics of the transformation. In Fig. 1, we report some representative spectra collected in a room temperature compression and an isobaric heating experiment up the reaction threshold conditions.
As it can be seen from the spectra collected during room temperature compression, the system is chemically stable at least up to 20 GPa. Careful analysis of the pressure evolution of the spectral components and of the relative parameters (peak frequency, FWHM and band area) shows some changes between 10 and 11 GPa (ESI2†), which are potential hints of a structural change as suggested by Raman spectra.45 However, the fact that most of the peaks are saturated and overlapped does not allow for a clear detection of a phase change requiring, in order to solve this issue, a direct structural characterization by X-ray diffraction.
The determination of the reaction's onset is not trivial because of the broadness of the product bands at high pressure, which makes its identification by spectral changes really challenging. To overcome this problem, we increased the pressure in small steps (0.1–0.2 GPa) subtracting from the last registered spectrum the previous one. In this way the reaction occurrence is readily evident especially from the reduction of the monomer bands, which appear as sharp negative peaks in the difference spectrum. An example is reported in Fig. 2 together with the crystal instability boundary accordingly identified. The higher temperature necessary to induce the reactivity at lower pressure is characteristic of the vast majority of molecular systems and is related to the larger amplitude of the thermal motion required, on decreasing the pressure, to balance the increase of intermolecular distances.19,31,34
Once the reaction onset was identified, we measured the kinetics of the process by acquiring IR absorption spectra as a function of time under constant pressure and temperature conditions. In Fig. 3, we report the kinetics measured at ambient temperature and two different pressures, 25.0 and 31.5 GPa. Although the reaction onset was identified at 300 K around 20 GPa, the transformation rate was found to be extremely slow at this pressure; therefore the kinetic studies were performed at higher pressures, conditions realized by quickly (2–4 minutes) increasing the pressure once the threshold was reached. The kinetics were analyzed by using the Avrami model developed to describe the crystal growth from the melt48–50 and then extended to the study of diffusion controlled solid–state reactions.51 The fraction of the azobenzene monomer R reacted as a function of time is described by the following relationship:
(1) |
Fig. 3 Kinetic evolution of the reactions performed at ambient temperature and 25.0 (lower) and 31.5 GPa (upper). The fraction of the monomer reacted, determined by the absorption of selected azobenzene bands (see the text), is plotted as a function of time. Eqn (1) is employed to reproduce these data; the fitting parameters are reported in both cases in the inset. The data collected at 31.5 GPa present a clear discontinuity after 270–300 minutes and have been reproduced (black line) by using the sum of two kinetic processes (red and blue lines). |
The main result of the kinetics' analysis regards the n values, all smaller than 1, thereby indicating a one-dimensional growth of the product as is expected for the formation of nanothreads. The need for using two different kinetic laws for the higher pressure reaction can be easily explained as due to the activation of competing polymerization paths along different directions because of the decreased intermolecular contacts at higher pressure. An analogous behaviour has been observed in the pressure induced polymerization of the ethylene crystal, where the anisotropic reduction of the lattice parameters with pressure activates a competing reaction path leading to gauche defected polymers.52 This additional reaction path, characterized by a rate constant one order of magnitude smaller, is activated when the faster step of the reaction has already progressed. The main (faster) kinetic contribution well agrees with the one measured at 25.0 GPa, exhibiting the same n value and a two times higher rate constant as expected according to the higher pressure value.
Further insight into the reaction mechanism is provided by the analysis of the molecularity of the process. The two sets of kinetic data just presented provide, as described in the ESI4,† a molecularity value consistent with a process, whose rate limiting step involves two units (molecules or oligomers) and therefore with a linear accretion of the product.
Clear evidence of the nanothread formation is the characteristic, generally hexagonal, spotty XRD pattern which attests to the 2D ordered lattice in which the nanothreads are arranged.16,24,25 In fact, since the threads develop along the applied stress direction, the same direction of the probing X-ray beam, we have from this technique direct evidence of the 2D arrangement of the threads. In Fig. 5, we have reported the 2D images recorded with synchrotron light of the samples recovered from some of the reactions performed together with the integrated diffraction patterns as a function of the d spacing of the samples synthesized at ambient temperature. Diffraction patterns have been acquired in air with the samples contained in the gasket without the diamonds. The 2D images of the samples produced at ambient temperature exhibit the characteristic pattern of the pseudohexagonal lattice although showing some differences between them, but especially with respect to pyridine24,25 and benzene,10,16 where a single diffraction peak was detected. Here the strongest and sharper (FWHM ∼0.4 Å) peak is observed at 5.6 Å, as in pyridine and benzene, but another peak is observed at a d value (11.2–11.5 Å depending on the sample) slightly larger than the double, and having a FWHM of about 0.9–1.0 Å. In the case of the sample produced at the lower pressure (21.7 GPa), we have a clear doubling of the hexagonal pattern likely due to the presence of slightly differently packed 2D crystallites. It is worth remembering that the X-ray beam employed for this analysis, 80 μm in diameter, probes almost the entire sample, thus suggesting a remarkable order all over the sample. The doubling in the 2D image corresponds to a visible splitting of the main peak (blue trace in Fig. 5) with two maxima at 5.34 and 5.54 Å. This splitting is not observed in the products from the reactions performed at 25.0 and 31.5 GPa, where only a peak at 5.60 Å is observed. The peak present at larger d values is broader, exhibits a considerable dispersion (0.3–0.4 Å) and, in contrast to the lower d one, derives from a spotty powder-like diffraction ring, thus suggesting some kind of disorder in the product along this direction. The fourth 2D image reported in Fig. 5(d) is relative to the recovered sample from the experiment performed at 530 K and 10 GPa. In this case, as also for the other high temperature experiment (15 GPa and 350 K), no diffraction peaks can be identified therefore suggesting the formation of an amorphous material.
The thread formation is confirmed by the high resolution transmission electron microscopy (HRTEM) analysis of the recovered materials. In Fig. 6, we report selected images collected at 60 keV primary electron beam energy on the samples recovered from the ambient temperature reactions performed at 25.0 and 31.5 GPa. The fragments analyzed have been obtained by breaking with a needle the sample contained in the gasket and selecting the thinner ones which are likely obtained from the peripheral areas of the recovered product. In these fragments are encountered both areas dominated by the presence of amorphous hydrogenated carbon (see ESI6†) and large portions showing perfectly parallel nanothreads extending for tens of nanometers. The occurrence of two products is in line with the competing mechanisms characterizing this kind of reaction. In some regions, these threads are bent (see the upper left panel of Fig. 6) producing curved bundles. A similar occurrence is observed in the micrographs obtained by Nobrega et al. on the nanothreads recovered by the compression of aniline26 and could likely be ascribed to the polycrystalline nature of the starting samples, which nevertheless allows the propagation of the thread through differently oriented adjacent crystals. The distance between the threads is very reproducible varying between 3.3 and 3.5 Å, a distance remarkably lower, about 30%, than those reported for benzene,10,16 pyridine24,25 and aniline,26 thus evidencing a denser and more compact material.
The XRD and HRTEM data define the constrains for the identification of the 2D crystalline arrangement of the synthesized nanothreads. XRD data unambiguously indicate a pseudo-hexagonal 2D arrangement with an interlayer distance of 5.6 Å corresponding to a lattice parameter a = 6.46 Å. However, the observation of the higher d peak provides additional constrains to identify the molecular arrangement in the 2D unit cell, which must also be reconciled with the small separation between the nanothreads, which develop as usual in the direction perpendicular to this plane, identified by HRTEM (3.3–3.5 Å). A monoclinic cell (a = 6.46 Å; c = 13.05 Å; β = 120°) as that reported in Fig. 7, with the molecule's dimensions as determined by X-ray diffraction,54 well accounts for all the experimental observations exhibiting diffraction peaks at 5.6 Å (100) and 11.2 Å (001) (see Fig. 5). The distance between the threads, obtained by piling up the phenyl rings along the direction perpendicular to the 2D lattice plane, depends on the viewing direction. In Fig. 7, we have highlighted how planes separated by 3.3 Å can be obtained rotating the threads by 90° around axes lying in the ac plane. This side view is also the most probable in the TEM experiment because the thin sample portions prepared for transmission experiments derive from a reproducible cleavage of the sample occurring along the thread direction due to the large Young's modulus of the nanothreads.12 Therefore, XRD experiments performed on the bulk sample contained in the gasket and TEM analysis performed on fragments obtained by breaking the sample provide complementary information related to the nearly orthogonal observation directions.
Our experimental data unambiguously indicate the obtainment of nanothreads, where the tubular structures connected by untouched azo groups are formed through the bonding among phenyl rings arranged in stacks along the azobenzene b crystal axis. In spite of this evidence, we do not have a clue of the nanothread structure. In the attempt of bridging such a gap, we considered the four polymeric structures recognized so far as the most favorable:15,55 tube (3,0), polymer I, polytwistane and zipper. Among these four structures, we can immediately rule out the polytwistane since the formation of this 1D helical structure along the stack direction would require a relative rotation of the neighbouring molecules, thereby becoming incompatible with the formation of double core nanothreads and the preservation of the azo groups. Also the formation of the zipper polymer is unlikely because a relative rotation of neighbouring molecules is required, thus increasing the activation energy of the reaction, and in spite of the possibility of realizing a double core polymer (see Fig. 8), a considerable strain characterizes the thread, thus making the internal energy considerably larger than that of the two remaining structures. According to these preliminary considerations, we have computed, using density functional theory (DFT), the energetics of isolated fully sp3 nanothread fragments composed of six molecular units having the tube (3, 0) and polymer I structures. In Fig. 8, we show the optimized structures of these nanothreads and their Gibbs energies after subtracting the electronic contribution. The detailed results of the calculation are reported in ESI7.† The energy of the polymer I structure is about 2.4 eV lower than that of the tube (3, 0) but it must be considered that in the latter polymorph 6 more hydrogen atoms are required to saturate the terminations thus increasing of 18 units the number of vibrational modes contributing to the zero point energy. By suitably subdivide these vibrations between C–H stretching (6) and bending (12) modes and assuming an average frequency for these modes of 2900 and 1200 cm−1, respectively, an additional contribution to the zero point energy of 1.97 eV occurs. Therefore, although the polymer 1 structure still presents an energy lower by approximately 0.5 eV, the small difference does not allow a certain selection of the thread structure. For this purpose, we have computed the IR spectra of these two nanothread fragments in order to compare them with those recorded on the recovered materials. After the structural optimization, no imaginary vibrational frequencies were obtained, thus indicating the achievement of the minimum of the potential surface for the resulting geometries. The comparison between computed and experimental spectra is shown in Fig. 8. The similarities between the computed polymer 1 spectrum and the experimental ones are impressive also in view of the reduced chain length considered. All the following characteristic features are nicely reproduced: the doublet below 800 cm−1, the reduced intensity between 900 and 1000 cm−1, the absorption density between 1000 and 1400 cm−1 and the other characteristic doublet between 1500 and 1600 cm−1. For the latter feature, the discrepancy in the intensity of the higher frequency component can be due to the presence of unsaturation in the recovered materials. The strong band around 1450 cm−1 present in the computed spectrum of the tube (3, 0) deserves a comment: it is ascribable to the –CH2 scissoring modes given the large number of such groups present in the simulated tube terminations. A schematic drawing of the possible reaction path leading to the polymer 1 thread, assuming an isotropic compression of the ambient pressure crystal structure and considering the monomer-thread structural relationships, is reported in the ESI8.† However, accurate XRD studies of the pressure evolution of the azobenzene crystal structure up to the reaction threshold pressure are mandatory to elucidate the reaction mechanism.
Angle-dispersive synchrotron X-ray diffraction measurements were performed on the recovered products at the Xpress beamline at the Elettra Italian synchrotron with a beam energy of 25.0 KeV, a focal spot diameter of about 80 μm (FWHM) and an image plate detector (MAR345). In-house X-ray diffraction measurements were performed with a custom made, laboratory diffractometer, equipped with a focused Xenocs-GeniX Mo Small Spot microsource with a wavelength of 0.7107 Å and a beam diameter on the focal plane of 150 μm. A PI-SCX 4300 CCD was employed as the detector.
Fourier transform infrared absorption spectra were recorded using a Bruker-IFS 120 HR spectrometer suitably modified for experiments in a diamond anvil cell, with an instrumental resolution set to 1 cm−1.57 The sample pressure was measured by the ruby fluorescence method using few mW of a 532 nm laser line from a doubled Nd:YAG laser source.
Transmission electron microscopy (TEM) investigation was performed with a JEOL ARM200F Cs-corrected microscope, equipped with a cold-field emission gun having an energy spread of 0.3 eV and operating at 60 keV. In conventional TEM (CTEM) mode, micrographs have been acquired in bright field (BF), while in STEM mode imaging has also been performed in Z-contrast mode using a high-angle annular dark field (HAADF) detector, with a probe size of 1.0 Å. The recovered material was mechanically removed from the gasket with a needle and placed directly onto the surface of a standard TEM lacey carbon film, covering a 300-mesh copper grid.
The structural optimization and the infrared spectrum calculation of the nanothread fragment composed of six molecular units have been performed using density functional theory approximation (DFT) with the Orca program,58 with the REVPBE density functional, the def2-TZVPP basis set, and the def2-TZVPP/C auxiliary basis set. The infrared frequencies are obtained by numerical Hessian computation.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d0sc06968j |
This journal is © The Royal Society of Chemistry 2021 |