Deepa Rajwara,
Xinfeng Liub,
Zheng Bang Lima,
Sung Ju Choa,
Shi Chenb,
Jens M. H. Thomasc,
Abbie Trewind,
Yeng Ming Lama,
Tze Chien Sum*b and
Andrew C. Grimsdale*a
aSchool of Materials Science and Engineering, Nanyang Technological University, Singapore. E-mail: ACGrimsdale@ntu.edu.sg; Fax: +65-6790 9081; Tel: +65-6790 6728
bDivision of Physics and Applied Physics, School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore. E-mail: tzechien@ntu.edu.sg; Fax: +65-6316 2971; Tel: +65-6795 7981
cInstitute of Integrative Biology, University of Liverpool, Liverpool, L69 7ZB, UK
dDepartment of Chemistry, Lancaster University, Bailrigg, Lancaster, LA1 4YB, UK
First published on 2nd April 2014
The synthesis of linear and trigonal terpyridine bearing molecules (tpys) and their self-assembly to form novel extended self-assembled 2D networks of trigonal tpys with linear tpys through zinc metal ion (Zn2+) co-ordination is reported. The resulting Zn2+ co-ordination networks have been characterized by means of X-ray photoelectron spectroscopy (XPS), small angle X-ray scattering (SAXS), BET, and photophysical methods. The presence of some short range order in these networks has been shown by the SAXS results and these results have been analyzed with the help of molecular modelling studies. These metallo-supramolecular Zn2+ networks revealed the influence of the metal ion on the thermal and optical properties of the synthesized metallo-supramolecular assemblies, similar to the results previously reported for 1D and 3D self-assembled metallo-supramolecular materials. Moreover, these networks have shown high luminescence with a long fluorescence lifetime and good thermal stabilities. The monolayer of such ordered networks can be used as a template to build hierarchical nanostructures. These hierarchical nanostructures could be used as active components in electronic devices and as templates for the formation of hybrid organic–inorganic nanomaterials.
Thus, this strategy of metal–ligand co-ordination represents an alternative method for “bottom-up” technique in nanoelectronics by synthetic approach for designing molecules of required dimensions of few nanometers. For instance, a wide range of 1D, 2D and 3D systems formed by co-ordination bonding have been widely explored in past years.7–12 Hence, these supramolecular co-ordination networks have found a lot of prospective to be utilized for different nanotechnological applications. For example, supramolecular co-ordination networks formed by Fe metal ion co-ordination have been reported with excellent magnetic properties.10 Ru2+ co-ordination complexes with terpyridine ligands have been used in the building of polymer wires on gold.13 Zinc metal ion complexes are very well known for their optoelectronic properties14–18 and are potential candidates for many optoelectronic devices. In another application, non-linear optical properties of zinc metal ion co-ordination networks have been reported recently.19
Terpyridines were proposed as building blocks for supramolecular chemistry decades ago.20,21 Still there remains much scope for development in the supramolecular chemistry of terpyridines and much research work has been performed in this field.11,13,18,22–27 Transition metal-complexes of π-conjugated bis- and tris(terpyridine) have high prospective to be used as new functional materials in electronic and optoelectronic devices. Heavy transition metal ions, such as Pt2+, Ru2+, Os2+, and Ir2+ have been explored thoroughly in the past.20,24,28 In particular, electron-poor square-planar divalent Pt2+ and Pd2+ ions have been widely used in conjunction with electron-rich nitrogen containing moieties in the self-assembly process.25,29 Metal ions with filled electron shells have recently achieved much importance as models for the advancement of new metallopolymers. The low cost and nontoxicity of Zn2+ ions has made them a promising candidate for making highly developed co-ordination networks. The other advantage of using Zn2+ ions for co-ordination is that metal to ligand charge transfer (MLCT) processes do not occur in such systems, because of the full d10 electron configuration of the Zn2+ ions. Hence there is only the possibility of intraligand (IL) charge transfer in these systems. High PL quantum yields and EL performance have been reported30 using such systems. Winter et al. have used Zn2+ as a template to assemble organic building block into polymer chains through coordination to chelating terpyridyl units.18 These linear oligomers with tpy (2,2′:6′,2′′-terpyridine) end-groups self-assembled with Zn2+ to give luminescent polymers. Moving one step forward towards achieving ordered and extended networks, this paper mainly focuses on the synthesis of 2D metallo-supramolecular networks of Zn2+ ions with terpyridine ligands and a comparison of their optoelectronic properties with the 3D networks previously reported by us.12,19 To create metallo-supramolecular networks, monomers with self-assembling end-groups are mandatory. In the previous work in this field, the focus has been on linear and tetrahedral Zn2+ terpyridine networks. By systematic examination on their photophysical properties together with their structural characterization, it has been concluded that by changing conjugated backbone on the terpyridine units it is possible to modify the optoelectronic properties of these networks. Strong emission in the solid state indicates that these materials correspond to the capable candidates with respect to their potential applications in optoelectronic devices such as OLEDs. Hence it was anticipated that the Zn2+ based materials might be made solution processable by the attachment of more solubilizing groups to π-conjugated terpyridines. These complexes could be very promising candidate for solar cell applications. In extension of this work, the present investigations aim the design of zinc metal co-ordinated 2D networks made up of trigonal and linear π-conjugated terpyridines with respect to their possibility in the field of optoelectronic device applications. Various bis(terpyridine) derivatives were used in the synthesis of self-assembled supramolecular 2D networks by Zn2+ metal co-ordination with tpy ligands. Photophysical properties with respect to their device applications along with their thermal, and electrochemical properties, as well as ordering in their self-assembled networks, are discussed, with comparisons made to the similar 3D networks previously made by us.12,19 The chemical structures of basic building blocks for 2D self-assembled networks are shown in Fig. 1.
Synthesis scheme of this trigonal core molecule 5 is outlined in Scheme 1. In the next step towards synthesis of molecule 1, in the presence of Pd(0) catalyst, compound 9 was Suzuki-coupled with triborane 5 to produce compound 1 as an off white solid. The synthetic steps involved in the preparation of key intermediate molecule 9 and new 2D p-type core molecule 1 for 2D self-assembly are outlined in Schemes 2–4. The structure of this novel trigonal-tpy molecule 1 was confirmed by 1H NMR and 13C NMR and its molecular weight was confirmed by high resolution mass spectrometry (HRMS). 2D self-assembled networks of molecule 1 were made through zinc metal ion co-ordination.
Formation of 2D Zn2+ co-ordinated self-assembled network N1 was obtained by refluxing its 2D monomers 1 with zinc acetate (Zn(OAc)2) at the stoichiometric ratio of 1:1.5, in NMP solutions followed by anion exchange using an excess of NH4PF6. Although 2D monomer 1 is soluble in common organic solvents (e.g. ether, DCM, CHCl3, THF, etc.), the 2D networks N1–3 are not only insoluble in those common organic solvents but also in highly polar aprotic solvents (e.g. DMSO, DMF, NMP, acetonitrile etc.) at room temperature. This behaviour is believed to be caused by the highly cross linked or interpenetrated structure of the 2D networks. Therefore, all the unwanted byproducts (including any unreacted starting material 1) were washed away using ether to get the final network N1 in solid form. The formation of new 2D Zn2+ co-ordinated self-assembled supramolecular network N1 from trigonal tpy molecule 1 is shown in Scheme 5.
Using the method previously developed by us for making extended 3D networks,12 extended networks N2 and N3 were formed using the p-type or n-type conjugated tpy linker units 3 and 2 as shown in Scheme 6. 2,5-Bisdodecyloxy-p-phenylene has been employed as the p-type linker (3) in N3. Its being linear and electron rich, it is a good electron donor. In the other extended network N2, 2,1,3-benzothiadiazole (benzo[c][1,2,5]thiadiazole) has been used as the n-type linker (2) as it is known to be a good electron acceptor. The same insoluble behaviour was found for the 2D extended networks N2 and N3. Networks N2 and N3 were obtained by refluxing 2D network N1 with zinc acetate Zn(OAc)2 at the stoichiometric ratio of 1:1.5 in NMP solution to break the network, followed by inserting the linear bis(tpy) monomers 2 and 3 in between the 2D molecules, respectively, to form the 2D extended networks N2 and N3 after anion exchange using an excess of NH4PF6. After forming networks N2 and N3, any residual 2D starting monomer (1) and the linear bis(tpy) monomers (2 and 3) were washed away with ether.
Scheme 6 Synthesis of 2D Zn2+ co-ordinated self-assembled two component extended supramolecular networks N2 and N3. |
The concept of reversible self-assembly was demonstrated in these networks which confirmed our concept for their formation. By adding more metal ions into the 2D network N1, co-ordination complexes can be dissociated again in favour of an open form with each zinc ion only attached to one tpy unit and its other co-ordination sites presumably saturated by solvent molecules. Disassembly of N1 was followed by re-assembly to form extended 2D networks N2 and N3. 2D network N1 was only partially soluble in the NMP solution at 105 °C and the colour of solution was light brown at that moment, after adding 1.5 equivalents of Zn(OAc)2 in NMP, the insoluble solid gradually dissolved in the solution and the colour of solution gradually changed to light yellow before adding the linear bis(tpy) linkers. All these novel 2D Zn2+ co-ordinated self-assembled networks N1–3 were obtained in moderate yields (>50%).
Fig. 2 XPS spectra of 1 and its zinc complexes N1, N2 and N3 on a carbon tape. Wide-scan spectra and zinc 2p region. |
N2 has two very broad weak peaks centred at scattering vector q = 0.38 nm−1 and 2.09 nm−1 respectively. The peak at scattering vector q = 0.38 nm−1 ranges between 0.3 and 0.5 nm−1. This corresponds to a distance of 1.65 nm with a range between 12.6 and 20.9 nm. The peak at scattering vector q = 2.09 nm−1 ranges between 1.00 and 3.10 nm−1. This corresponds to a distance of 3.0 nm with a range between 2.0 and 6.3 nm. N3 shows one broad but clear peak centred at scattering vector q = 2.66 nm−1 and ranges between 2.00 and 3.00 nm−1. This corresponds to a distance of 2.3 nm with a range between 2.1 and 3.0 nm. N1 shows three clear, although relatively broad, peaks compared to N2 and N3, which have fewer peaks that are both weaker and broader. This suggests that N1 has a more defined structure compared to N2 and N3. This reflects the single component self-assembly composition of N1, in contrast to N2 and N3, which are composed of two different basic building blocks and will contain a statistically random combination of each. Due to the insolubility of the networks we could not form films for study by STM or AFM as has been recently done for sheets of a similar network made by Schlüter et al.31 High resolution TEM images were obtained from deposits of N1 (see ESI†) on carbon coated copper grids but were unable to determine any structural data on the networks from them. Even with the use of electron energy loss spectroscopy (EELS) we were unable to see any peaks corresponding to the presence of zinc ions, which we attribute to the low concentration of the zinc in these networks.
Firstly, atomistic models of the building blocks, 1 and 3, and their respective one-component and two-component extended 3D self-assembled networks were generated. Models of 1 and 3 were fully optimised using the COMPASS force field37 in the Forcite module of Materials Studio 5.0 (ref. 38) with charges calculated using the Gasteiger method39 and are shown in Fig. 5. Assessing the optimised geometry of model of 1 shows that the structure deviates from the idealised planar form. There are three angles that can be defined between the centroid of the central phenyl and the central nitrogen atom of two of the tpy ligands for 1. These range between 103° and 138° compared to the ideal of 120°, and is 2 Å out of plane. The geometry of the optimised model for 3 shows a planar rigid molecule with the tpy-N–centroid–tpy-N angle being 177°.
Using these idealised building units, we have rationalised the SAXS data discussed above. For N1, there are four distinct repeat units available. These are, in order of dimension, (i) the distance between the molecular centroid of 1 and a coordinated Zn2+, (ii) the molecular centroid – molecular centroid, and finally, depending upon the topology of the network either (iii) the molecular centroid – molecular centroid across a hexagonal repeat unit, if the network has a 2D graphite-like honeycomb structure, or (iv) a Zn2+ to opposite Zn2+ across a partial hexagonal ring, if the structure has a 3D structure. For N3, the two-component self-assembly results in a number of types of repeat units being available. These are the based on the same repeat units as identified for N1, but for each distinct repeat unit of N1, there are now a number of different combinations available for N3, particularly for repeat unit type (iii) where there will be multiple combinations possible. This variation in the build units results in the likelihood that these types will not be observed in SAXS measurements with only the build units of type (i) now being regular repeat units. These repeat units are identified in Fig. 5(c) and (d) for network N1 and N3 respectively.
A hexagonal and partial hexagon ring for both networks were generated for N1 and N3, shown in Fig. 5(c) and (d) for network N1 and N3 respectively and were fully optimised using the COMPASS force field37 in the Forcite module of Materials Studio 5.0 (ref. 38) with charges calculated using the Gasteiger method.39 For N1, the average distance between repeat units (i) and (ii) is 2.0 nm and 4.0 nm respectively, correlating well to the observed peaks which have corresponding distances centred on 2.5 nm and 5.0 nm. For the idealised hexagon ring the type (iii) distance is 8.0 nm, which does not correlate well with the final peak observed in the SAXS data centred at 12.6 nm. For the partial hexagon ring, the type (iii) distance is 12.0 nm, which correlates well with the final peak observed in the SAXS diffractogram centred at 12.6 nm. This suggests that hexagon rings are not a significant structural component of N1. For N3, the type (i) average distances measure 2.0 nm (same as for N1) and 2.4 nm, which correspond well with the peak centred at a distance of 2.3 nm. No further peaks are observed in the SAXS diffractogram for N3, suggesting that the structural components of type (ii), (iii) or (iv) and are not repeating in dimension. This also suggests that hexagon ring formation is not a significant structural component of N3. The lack of a hexagon ring component in both N1 and N3 further suggests that the building blocks, 1 and 3 are not rigid structurally-directing units but have some conformational flexibility.
To fully explore the conformational freedom of the two building blocks, we have used molecular dynamic simulation at 100 K and 298 K using the Forcite module, the COMPASS force field and charges assigned using the Gasteiger method and with an NPT ensemble. For each simulation, a total of 2 ns simulation time with a 1.0 fs time step was used. The angles and out-of-plane distance are monitored for the final 1 ns of each simulation for 1 and 3. For 1, the three angles range between 64° and 174°. One angle and the out-of-plane distance is monitored throughout the 298 K MD simulation showing an average angle of 111 ± 18° and an average out-of-plane distance of 6.4 ± 1.2 Å. It is noticeable that during the 298 K MD simulation the extended side chains wrap inwards towards the centre of the molecule. For 3, the tpy–centroid–tpy angle ranges between 158° and 180° through out-of plane and in-plane molecular bends with a distance of between 0 Å and 2 Å out of the molecular plane. This demonstrates the high degree of flexibility and conformational freedom available to both 1 and 3 network building blocks.
We then turned our attention to the automated generation of a 3D self-assembled network structure. The 3D self-assembled networks were broken into fragment building units; these are not necessarily chemically representative of the actual monomeric reactants but enable us to automate the building process. For N1, we used two basic building units, Fragments A and B, shown in Fig. S.11(a) and (c) respectively in ESI.† For Fragment A (based on 1), a library of 21 different conformers is available. This includes the idealised optimised structure and an additional 20 structures generated by sampling the two MD simulations. The MD simulations were sampled every 100000 steps and the geometry of the resulting structures were optimised before the tpy end groups were removed. The remaining three unsaturated carbon atoms of the terminal phenyl groups were assigned as bonding atoms. The distance between the bonding atoms and the molecular centroid is approximately 14.0 Å. Fragment B is generated from the Zn2+ and the tpy coordinating ligands with the two PF6− counter ions. The geometry of both the Zn2+-tpy and PF6− were fully optimised with DFT, using NWChem 6.1,40 a def2-tzvp basis set,41 B3LYP functional42 and incorporating the Grimme dispersion correction.43 A hydrogen atom is removed from the two central rings of each tpy group and the remaining unsaturated carbons were allocated as bonding atoms. The bonding atom – bonding atom distance is 9.7 Å. The incorporation of two PF6− groups within Fragment C ensures that these counter ions are accounted for within the network generation process. For N3, constructed from 1 and 3, we included an additional Fragment C (shown in Fig. S.11(b) in ESI†), which is generated from a library of 5 conformers which include the idealised structure and an additional four created by sampling the MD simulation of 3, similarly to 1. The terminal tpy groups were removed and the remaining unsaturated carbons on the terminal phenyl groups are allocated as bonding atoms. The distance between the bonding atoms of Fragment C is 11.5 Å. Representative structures for Fragments A, B, and C are shown in Fig. S.11(c)–(e)† respectively with the bonding atoms highlighted in yellow. For Fragments A and C, Fig. S.11 and S.12† shows an overlay of the library of structures available for each fragment.
A Python code was used to automate the generation of the self-assembled network. The code builds the network following pre-set rules for a series of building steps. The first step is a Seed step where one or more fragments are added to a simulation cell of pre-set dimensions. We assign rules as to which Fragments are permitted to bond as follows; for network N1, Fragment A can bond to Fragment B, but not A to A or B to B, similarly for network N3. Fragment A can bond to B but not A or C and Fragment C can bond to B but not C or A. A Build step bonds a randomly selected additional fragment to a randomly selected existing fragment or partially built network. The fragment is bonded with idealised bond lengths, although the torsional angles of the fragments about the bond are random. Additional bonds between fragments are only permitted if the bond geometry fits within the specified geometrical parameters, the bond distance is within a set range and no non-bonding atoms are closer than 1 Å. This ensures that the initial generation of the network is geometrically sensible.
We generated an idealised cluster for each network N1 and N3. This enables us to assess the presence of the idealised structural components and their respective dimensions within an extended structure. The idealised clusters are generated exclusively from the idealised models of 1 and 3, and do not include any of the structures obtained from the MD simulation. The idealised cluster for N3 was generated following the idealised alternating 1 and 3 building pattern. The clusters were generated by an initial Seed step with one randomly chosen fragment and then a series of Build steps following the bonding rules. Each cluster Fragment is capped with either an uncoordinated tpy or a H atom. After generation, the geometry of the cluster is optimised using the Fast Inertial Relaxation Engine (FIRE) rigid body minimiser within HOOMD-blue.44 The fragment build units are treated as rigid bodies so that only the bonds between fragments are optimised. The Intermolecular forces, bond and angle parameters are taken from the PCFF force field.45 Intermolecular forces are described using a Lennard-Jones potential with additional Zn2+ parameters.46 Charges are calculated using the Gasteiger method and QEq for the Zn2+ coordinated molecule in fragment B.39 The idealised clusters for N1 and N3 are shown in Fig. S.12(a) and (b)† respectively. There are structural features of type (iv), as described above, present in both idealised networks N1 and N3. For N1, the average is 11 nm and range between 10.0 nm and 14.4 nm, in good agreement with the range observed in the SAXs diffractogram for the peak centred at 12.6 nm with a range between 9 nm and 15.7 nm. For N3, the average is 16 nm and range between 19 nm and 14 nm. This does not correlate with any peaks in the diffractogram for N3. There are no structural components of type (iii) due to the automated network generation method. For N1, the type (i) and type (ii) are similar to the values obtained for the hexagon ring and for the partial ring. For the cluster model of idealised N3, the repeating structural components of type (ii) are well defined at 5.0 nm. We would expect to see peaks that correspond to both these distances in the SAXS diffractogram. These peaks are not present, with only peaks corresponding to the type (i) structural components present. This suggests that the regular alternating building unit model for N3 does not represent the real structure. An alternative model for N3 can be generated based on a randomised combination of the building block Fragments A and C.
Finally, we have constructed cluster models that include the structural flexibility observed in the MD simulations and randomised combination for N3. For N1 and N3 we have included the full library of structures for Fragment A and C. For network N3 the choice is random between Fragment A and C with an equal probability for either Fragment. A representative cluster is shown in Fig. 6(a) and (b) for N1 and N3 respectively also showing the underlying net and fragment combination. Five clusters were generated for each network and are shown in Fig. S.13.† The clusters of network N3 are largely planar, as might be expected from the two dimensionality of Fragment C, whereas those of N1 extend in all three dimensions. We see a wide range of different combinations of Fragment A and C leading to a very wide range of structural components. The only repeat units observed in the randomised clusters of N3 are those of type (i). This accounts for the absence of any peaks observed that do not correlate to type (i) distances within the SAXS diffractogram. Each cluster shows an extended structure with no evidence of order. The networks often turn in on themselves, suggesting that network interpenetration is likely. This suggests that an amorphous, efficiently packed solid will be formed with no porosity. This is in agreement with the porosity study discussed above.
Fig. 8 Solid state photoluminescence spectra of core molecules 1, 2, 3 and their 2D self-assembled networks N1–3. |
Moreover, it was expected that the Zn2+ tpy connectivity will increase the electron delocalization on the conjugated and rigid spacers, which can be observed by red shift in the PL spectra of N1–3. The reason behind the red shift in the emission spectrum is the increase in conjugation, which will increase the number of delocalized electrons in the system. Finally, the longer the value of conjugation length, the lower will be the energy gap between closest energy levels and the longer will be the absorption wavelength, and so the emission wavelength too. Femtosecond time-resolved photoluminescence dynamics of 1 and its self-assembled networks N1–3 and comparison with previously reported 3D networks has been performed. Fundamentally, fluorescence is a property of electronically excited states. Thus, fluorescence mechanisms can be better understood by time-resolved fluorescence spectroscopy. To obtain the dynamics of the PL spectra, the fluorescence lifetimes of compound 1 and its Zn2+ self-assembled networks N1–3 were measured. The fluorescence decay behavior in the presence of Zn2+ is shown in Fig. 9, with the exponential fit results. Table 2 shows the fitted lifetimes and the pre-exponential factors. For comparison, we have also calculated the average lifetime. Comparatively, molecule 1 was found to show a shorter lifetime compared to N1, N2 and N3. An increase in the fluorescence lifetime of networks N1–3 was observed after zinc meal ion co-ordination of trigonal and linear terpyridine moiety in 1, 2 and 3. Similar results were observed of networks S1 (similar to N1), S3 (similar to N3) and S4 (similar to N2) made up of a tetragonal tpy molecule (having similar branching units to trigonal tpy molecule 1) previously reported by us (see ESI†).12 The values of fluorescent life times are given in Table S1 in the ESI.† These results further suggest that networks with similar units have similar lifetime values. Enhancement in the lifetime after zinc metal ion co-ordination may be attributed to the increase in conjugation after metal complexation as observed with a red shift in the PL spectra too. These results also suggest a complete modification in the energy levels after zinc metal ion complexation. Networks N2 and S4 were formed by the n-type linker unit 2, which has a lower band gap. Hence the fluorescence lifetime of these networks are also lower compared to the results obtained for N1, S1, N3 and S3 which were formed by the p-type linker units.
Networks | A1 | τ1 (ns) | A2 | τ1 (ns) | τaverage (ns) |
---|---|---|---|---|---|
a An average lifetime, τa, was also calculated using τa = (A1τ12 + A2τ22)/(A1τ1 + A2τ2) for comparison. | |||||
1 | 0.56 | 0.10 ± 0.01 | 0.44 | 0.53 ± 0.02 | 0.45 ± 0.06 |
N1 | 0.63 | 0.17 ± 0.01 | 0.37 | 2.16 ± 0.02 | 1.92 ± 0.07 |
N2 | 0.60 | 0.19 ± 0.01 | 0.40 | 0.84 ± 0.01 | 0.68 ± 0.04 |
N3 | 0.46 | 0.61 ± 0.01 | 0.54 | 2.75 ± 0.02 | 2.41 ± 0.06 |
2D Zn2+ supramolecular networks were compared to the 3D Zn2+ networks in order to understand the influence of the Zn2+ ions on the photophysical properties of the designed structures. The synthesized 2D Zn2+ co-ordination networks revealed marked red shifts in the emission spectra like the 3D Zn2+ networks made up of similar building blocks. It was observed that the presence of the Zn2+ ions in the synthesized model complexes revealed interesting spectroscopic properties, in particular, enhanced fluorescence and fluorescence life time. Thus, from the experimental results it was concluded that the characteristics of the basic building units are considerably influencing the photophysical properties of the Zn2+ co-ordination networks. Such Zn2+ systems could be potentially interesting for the construction of optoelectronic devices.
Footnote |
† Electronic supplementary information (ESI) available: NMR, CV, TGA, DSC, TEM and mass spectral data for new compounds and networks. Structures, PL profiles and PL data for S1, S3 and S4, molecular model for 1, 3, N1, and N3. See DOI: 10.1039/c3ra47377e |
This journal is © The Royal Society of Chemistry 2014 |