Catherine M.
Aitchison‡
a,
Christopher M.
Kane‡
a,
David P.
McMahon‡
b,
Peter R.
Spackman
bc,
Angeles
Pulido
b,
Xiaoyan
Wang
a,
Liam
Wilbraham
d,
Linjiang
Chen
ac,
Rob
Clowes
a,
Martijn A.
Zwijnenburg
d,
Reiner Sebastian
Sprick
a,
Marc A.
Little
a,
Graeme M.
Day
*b and
Andrew I.
Cooper
*ac
aDepartment of Chemistry and Materials Innovation Factory, University of Liverpool, Liverpool L7 3NY, UK. E-mail: aicooper@liverpool.ac.uk
bComputational Systems Chemistry, School of Chemistry, University of Southampton, Southampton SO17 1BJ, UK. E-mail: g.m.day@soton.ac.uk
cLeverhulme Research Centre for Functional Materials Design, University of Liverpool, Liverpool L7 3NY, UK
dDepartment of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK
First published on 20th March 2020
We show that a hydrogen-bonded framework, TBAP-α, with extended π-stacked pyrene columns has a sacrificial photocatalytic hydrogen production rate of up to 3108 μmol g−1 h−1. This is the highest activity reported for a molecular organic crystal. By comparison, a chemically-identical but amorphous sample of TBAP was 20–200 times less active, depending on the reaction conditions, showing unambiguously that crystal packing in molecular crystals can dictate photocatalytic activity. Crystal structure prediction (CSP) was used to predict the solid-state structure of TBAP and other functionalised, conformationally-flexible pyrene derivatives. Specifically, we show that energy–structure–function (ESF) maps can be used to identify molecules such as TBAP that are likely to form extended π-stacked columns in the solid state. This opens up a methodology for the a priori computational design of molecular organic photocatalysts and other energy-relevant materials, such as organic electronics.
Molecular organic crystals have highly ordered structures that can be prepared in a modular way using solution-processable units. This makes molecular crystals attractive candidates for studying the effect of secondary structure on photocatalytic activity, since it is possible to compare materials that are chemically identical and that differ only in terms of their solid-state packing. By contrast, such structural comparisons are challenging for amorphous polymers and extended organic networks, such as covalent–organic frameworks (COFs), which often have only moderate crystallinity: for example, there are only a handful of single crystal structures reported in the literature for COFs,10,11 and none of those materials have been shown to have photocatalytic activity. On the other hand it is relatively straightforward to grow high-quality single crystals of organic molecules. Until now, however, there are no examples of appreciable photocatalytic hydrogen evolution from molecular organic crystals.
We showed previously that amorphous pyrene-containing polymer networks can produce hydrogen photochemically from water in the presence of a sacrificial hole scavenger.12 Recently, Lotsch et al. reported that a more crystalline material, a 2-D pyrene-based COF,13 shows higher photochemical activity. While this material has in-plane conjugation, out-of-plane conjugation between the close-packed organic layers was also invoked as an important structural feature for high photocatalytic activity. Through-plane conjugation is a relatively common feature of pyrene-containing materials and pyrene itself has been described as “the fruit fly of photochemists”14 because it is known to have interesting packing-related photophysical properties for organic electronics. Here, we chose three pyrene-based molecules; tetraphenylpyrene (TPhP), 1,3,6,8-tetrapyridin-4-yl pyrene (TPyP), and 1,3,6,8-tetra(4′-carboxyphenyl)pyrene (TBAP) (Fig. 1), which all form stable crystalline solids, to investigate the effect of molecular structure and extended crystal packing on photocatalytic hydrogen production rate.
In its native α-form, pyrene has a herringbone-type crystal packing arrangement, and this α-polymorph transforms to the structurally-related β-polymorph below 163 °C.15 By crystallising pyrene from CH2Cl2 at pressures > 0.3 GPa, it is possible to obtain the denser γ-polymorph, which features aligned π-stacked columns of pyrene molecules.16 However, this structure was reported to dissolve in the crystallisation solvent during decompression, making it impossible for us to access this π-stacked polymorph for photochemical experiments. The pyrene derivative, TPhP, crystallises from CH2Cl2/n-hexane,17 to produce a solid that also lacks π–π stacking interactions between pyrene molecules. This is because the phenyl groups frustrate the packing of pyrene cores.
Here, we use TPyP and TBAP to also investigate the effect of pyridyl and benzoic acid groups on the extended packing of pyrene cores. Hydrogen bonding has been shown to have the potential to frustrate dense packing of organic molecules, and to generate electrostatically-stabilised, low-density, hydrogen bonded organic frameworks (HOFs).18–20 Similarly, labile C–H⋯N bonding interactions in molecular crystals, comprising pyridyl functionalised molecules, have been used to stabilise low density crystal packings.21 Hence, TPyP and TBAP both feature functional groups that have the potential to direct low-density, porous crystal packings with extended π-stacks, which is attractive for organic photocatalysts.
When we started this study, there were no reported structures of either TPyP or TBAP in their native uncoordinated states, even though these molecules have been used as struts in metal–organic frameworks (MOFs).22–24 A crystal structure of TBAP was subsequently published while this study was in progress.25
Molecular crystals are not subject to the same intuitive design rules associated with materials such as MOFs26–28 and COFs,29–31 because the crystallisation of organic molecules is governed by an interplay of many weak intra- and intermolecular interactions. It therefore remains a significant challenge to design, from first principles, new organic molecular crystals that feature porosity and/or extended π-conjugation. Porosity is also much less common in organic crystals because of their strong preference for close packing.32 At the atomistic level, discovering a new material involves substantiating a stable minimum on the free energy surface defined by chemical composition and the relative positions of atoms within an extended solid.33 Traditionally, the search for new functional molecular materials has been led by knowledge-guided experiment. However, the role of computation in materials discovery has evolved rapidly in the past few years, from explanatory post-analysis of new materials to genuinely predictive methods that can be applied in advance of, or in tandem with, experiment. Such computationally-led approaches34–36 promise to accelerate experimental searches and to unearth materials that might otherwise have remained undiscovered.37 One area where this applies is the design of porous materials, which have applications in gas storage,38–41 separation,42–44 and organic catalysis,45 as well as photocatalytic water splitting.46–48 For example, we have developed computational methods for crystal structure prediction (CSP) to generate energy–structure–function (ESF) maps,36,49 that have been used to target the discovery of highly porous hydrogen-bonded-frameworks (HOFs)36 and organic semiconductors.50,51 These ESF maps summarize the energetic distribution of stable crystal structures available to a given molecule, along with simulated properties relevant to the desired function.
No prior knowledge is required about the crystal packing of a candidate molecule because these approaches are based on ab initio structure prediction. However, to be reliable, the CSP method must fully explore the configurational space of available crystal packings; this is a high dimensional problem including molecular positions, orientation, and unit cell dimensions, as well as flexible intramolecular degrees of freedom. The associated computational expense has meant that most applications of CSP for the discovery of functional materials have been limited to rigid molecules,36,40,50–52 where a lack of conformational freedom leads to a reduced search space. TPhP, TPyP and TBAP are conformationally flexible and in this study, we have implemented a CSP method that accounts for this flexibility during structure searching, thus expanding the scope of CSP in the de novo design of functional materials. Specifically, we performed CSP calculations with TPhP, TPyP and TBAP, to determine whether these molecules were likely to form porous or π-stacked structures that could increase their photocatalytic activity for hydrogen production from water.
By contrast, the computed crystal energy landscape for TBAP (Fig. 2a) contains several regions of low-energy, low-density predicted structures that fall well below the bulk energy-density trend. These low energy ‘spikes’ are reminiscent of those on the landscapes of the triptycene benzimidazolone molecule, T2,36 and also trimesic acid,37 both of which corresponded to experimentally accessible porous structures. These spikes correspond to isolated, deep basins on the lattice energy surface, separated by a high energy barrier from regions of the lattice energy surface corresponding to dense structures. The most prominent of these spikes has a density of ca. 0.6 g cm−3; the lowest energy structure in this spike, 1 (Fig. 2a) has a lattice energy 57 kJ mol−1 above the dense global minimum on the TBAP landscape.
TBAP structure 1 features 2-dimensional sheets with rhomboid voids held open by acid–acid hydrogen bonds between TBAP molecules (Fig. 3). These sheets are stacked to form infinite pyrene columns (Fig. 3b) and parallel channels that run perpendicular to the hydrogen-bonded sheets. A second spike at a density of ca. 0.75 g cm−3 features similar structures containing one-dimensional channels between π-stacked pyrene columns, but with some collapsed channels (ESI, Fig. S1†). Thus, the CSP results suggested that TBAP had the potential for the construction of porous frameworks, like MOFs22–24 but without the inclusion of metals. During the course of this study, although after these calculations, structure 1 was in fact reported independently by another research group,25 confirming our prediction.
Fig. 4 shows ESF maps that summarize the propensity of TPhP, TPyP, and TBAP to crystallise to form extended π–π stacks with significant overlap between the pyrene cores. In this analysis, we define an extended π–π stack as four or more nearly co-parallel molecules with limits on the separation between centres of mass (see ESI, Section 1.4 for details†). Due to the translational symmetry, this is equivalent to infinite stacks in almost all cases. All three of these molecules have similar energetic distributions of stacked versus unstacked structures. Note that many unstacked structures are hidden in Fig. 4 (see ESI, Fig. S66†). In all three cases, the low-energy edge of the energy-density distribution is dominated by structures containing π-stacked columns. One key difference is the presence of low-energy spikes on the TBAP landscape, discussed above, which are not present for TPhP or TPyP. The structures within these deep basins on the lattice energy surface all show extended π–π stacking; that is, these regions of enhanced stability are explicitly linked to π–π stacking. We would therefore predict a priori that TBAP is a more promising candidate for π-stacked HOFs than either of its two structural analogues, TPhP or TPyP.
Fig. 4 Energy-density distributions of predicted crystal structures. (a) TBAP, (b) TPhP and (c) TPyP. Coloured data points are structures containing extended stacks of coplanar molecules (see ESI for details†) and are coloured by the extent of molecular overlap between stacked molecules (inset, part D). Structures without stacking are shown in grey, but stacked data points are plotted on top of unstacked ones, thus hiding many of the unstacked structures (see ESI, Fig. S66†). Annotations in A–C refer to the observed crystal structures. (d) Histogram of the degree of molecular overlap in predicted crystal structures with extended stacks. Note that the histograms are transparent, so that overlapping regions appear dark green. |
Charge transport is highly sensitive to small changes in molecular packing, which dictates the electronic coupling between molecules. We therefore developed methods to analyse geometric parameters within predicted crystal structures that display extended π-stacked columns. The distribution of predicted structures in Fig. 4 is coloured according to the geometric overlap of neighbouring stacked molecules, which was calculated from the projection of their planes of best fit (Fig. 4d and ESI†). This geometric measure of stacking overlap is used as a proxy for overlap of frontier molecular orbitals and is correlated with efficient charge transport and exciton dissociation in π-conjugated materials.56–58
From a predictive point of view, the distribution of π-overlap on the ESF maps is important. Experimentally accessible structures are expected to be found either near the lattice energy global minimum or, in the case of porous structures that are stabilised by solvent inclusion during crystal growth, along the low energy ‘leading edge’ of the energy-density distribution. All three molecules have predicted structures with nearly perfect π–π overlap (overlap ≈ 1, Fig. 4d), but these occur in high energy regions of their crystal structure landscapes, away from the global energy minimum and leading edge of the ESF maps. The structures in experimentally accessible regions are those with overlaps in the range 0.85–0.95, where TBAP structures show significantly higher overlap than TPhP and TPyP (blue peak in histogram; Fig. 4d). This difference stems from the carboxylic acid groups in TBAP, which direct the molecule into crystal packings that favour π-stacks with strong π–π overlap. This demonstrates a global, predictable preference for TBAP to adopt favourable packing for charge transport compared to TPhP and TPyP, which may enhance charge-transport to the interface with water and the hole scavenger. While the global energy minima structures for both TPhP and TPyP also show π-stacking with a reasonable degree of overlap (0.863 and 0.857, respectively), an experimentally known structure of TPhP17 that sits well above the global energy minimum displays no π-stacking, suggesting that these more ‘optimal’ structures may be kinetically difficult to access, at least for TPhP.
These three crystalline solids, TPhP-α, TPyP-α, TBAP-α, display different π–π stacking modes of pyrene cores, and the TBAP-α structure features a high surface area. These three crystalline solids are therefore good candidates to investigate the effect of crystal structure on photocatalytic activity in molecular crystals.
Material | Hole-scavenger | pH | HERa (μmol h−1 g−1) |
---|---|---|---|
a Catalyst (25 mg) loaded with 1 wt% Pt, from in situ photodeposition of H2PtCl6, suspended in water and scavenger (25 mL), irradiated with a 300 W Xe light source fitted with a λ > 420 nm filter. The HER was determined over five hours. b HER calculated over 20 hours. c Material fully or partially dissolved under these conditions. | |||
TPhP-α | Ascorbic acid 0.1 M | 2.6 | 2b |
TPyP-αc | Ascorbic acid 0.1 M | 2.6 | 18 |
TBAP-α | Ascorbic acid 0.1 M | 2.6 | 1293 |
Amorphous TBAP | Ascorbic acid 0.1 M | 2.6 | 6 |
TPyP-α | Ascorbic acid 0.1 M | 7 | <0.1 |
TBAP-α | Ascorbic acid 0.1 M | 7 | 3108 |
Amorphous TBAP | Ascorbic acid 0.1 M | 7 | 156 |
TPhP-α | Triethylamine 5 vol% | 11.5 | 6 |
TPyP-α | Triethylamine 5 vol% | 11.5 | 40 |
TBAP-αc | Triethylamine 5 vol% | 11.5 | <0.1 |
No hydrogen production is observed in the dark or in the absence of the TBAP-α. When D2O was used as the proton source, D2 production was mostly observed (ESI, Fig. S29†), with a small amount of H2, probably due to H–D exchange with protons in the non-deuterated ascorbic acid.60 Essentially no CO production is observed under these conditions (ESI, Fig. S31†). Taken together, these observations lead us to conclude that the hydrogen production process is indeed photocatalytic.
TPyP-α was much less active and produced hydrogen at a rate of just 18 μmol g−1 h−1, although the basic pyridyl groups in TPyP meant that a significant proportion of the TPyP-α catalyst (>50 wt%) dissolved in the acidic medium during this measurement. Consequently, a direct comparison between the HER activity of TPyP-α and TBAP-α in ascorbic acid is not possible. TPhP-α is stable in ascorbic acid but had an even lower HER of 2 μmol g−1 h−1.
To allow for a direct comparison between TBAP-α and TPyP-α, we tested the materials in a 0.1 M ascorbic acid solution adjusted to pH 7 with NaOH, where neither material dissolves. Under these conditions, TPyP-α produced no measurable amount of hydrogen over a five-hour period. By contrast, TBAP-α was even more active with a HER of 3108 μmol g−1 h−1 (ESI, Fig. S24†). This large change in HER can be rationalised by considering the driving forces of the two half reactions occurring in the system. Density functional theory (DFT)61,62 calculations performed on isolated molecules immersed in water at pH 2.6 (the expected pH of 0.1 M ascorbic acid solution, Fig. S48†) suggest that all three materials should have a large driving force for proton reduction and reasonable driving force for the overall oxidation of ascorbic acid. The driving force for the initial one-hole oxidation of ascorbic acid is very small for TBAP and TPyP and actually slightly negative for TPhP. Thus, we might expect the oxidation of the scavenger to be rate limiting in these systems, accounting for the increased activity of TBAP-α when changing the pH level from pH 2.6 to pH 7, because the driving force for ascorbic acid oxidation increases. This change also reduces the driving force for proton reduction, however, even at pH 7, the driving force for proton reduction remains large (>1.5 V) and, crucially, larger than that for ascorbic acid oxidation (0.9 V for the 2-hole and 0.5 V for the intermediate one-hole oxidation). This does not explain the decrease in rate of TPyP-α when changing the pH level from pH 2.6 to pH 7, but we note that the partial dissolution of the material and the influence of protonation on the substrate electronics under acidic conditions could play a role.
When 5 vol% triethylamine in water solution was used as the sacrificial system, TBAP-α was found to dissolve and no measurable hydrogen was produced under visible light irradiation. TPyP-α and TPhP-α were both stable in 5% triethylamine but had low HERs of 40 and 6 μmol g−1 h−1, respectively. None of the materials produced hydrogen when tested in a 5 vol% TEA solution adjusted to pH 7.
At first glance, these results could suggest that porosity is the dominant factor for hydrogen evolution, because the porous TBAP-α structure greatly outperforms the non-porous analogues, even though TPyP-α also contains significantly overlapped π-stacking. This may however be an oversimplification. As evident from water sorption isotherms, TBAP-α has no significant water uptake between partial pressures of 0.2–0.6 (Fig. S52†). There is significant water uptake at higher relative pressure, which might suggest water adsorption but could also indicate water molecules condensing on the crystal surfaces or between crystals. We note that this wetting behaviour may be different in the presence of the sacrificial agent, ascorbic acid, and in this respect, water sorption isotherms may not reflect the photocatalysis conditions.
Perhaps more significantly, it is unclear that the platinum cocatalyst, which is possibly the site for proton reduction, actually resides in the pores of the TBAP HOF. This raises additional doubts that porosity alone can account for the superior performance of this material. In the absence of added Pt co-catalyst, TBAP-α had a dramatically reduced rate of 59 μmol g−1 h−1. It is possible that residual palladium from synthesis acts as the active site in this case. Pd levels by ICP-MS were found to be below the 10 ppm detection level of the instrument but we note that very low concentration can be sufficient to give limited photocatalytic activity.63 The addition of 1 wt% Pt was found to give the highest catalytic activity, while increasing the loading to 4 wt% appeared to reduce HER (ESI, Table S6†), perhaps due to reduced light absorption or recombination.64,65 Element analysis (ICP-MS; ESI, Table S6†) and STEM imaging were employed to confirm that in situ photo deposition of platinum had been successful (Fig. 5 and ESI Fig. S54†). Platinum nanoparticles of between 2 and 15 nm formed on both TBAP-α and amorphous TBAP. It was noted that the distribution of Pt in the samples with 1 wt% Pt was generally more even, and that larger, less well-dispersed clusters could be observed at the higher 4 wt% Pt loading (Fig. 5 and ESI, Fig. S54†).
Fig. 5 STEM images of TBAP-α in (A) BF mode and (B) HADF mode, with 1 wt% photodeposited Pt in (C) BF mode and (D) HADF mode and with 4 wt% photodeposited Pt in (E) BF mode and (F) HADF mode. |
Photolysis experiments using 4 wt% Pt were repeated for three different batches of TBAP-α and showed good reproducibility between batches (ESI Fig. S23†). The PXRD patterns and sorption isotherms of the three batches were also very similar (ESI Fig. S19–S22†) with BET surface areas of 2001, 2270 and 2074 m2 g−1. Likewise, static light scattering showed only modest variations in the particle size distribution between batches (average diameters of 9.25–15.52 μm, ESI Fig. S62 and Table S9†).
To investigate further, we conducted photocatalysis experiments using large (13 nm by DLS, ESI Fig. S55†) pre-made platinum nanoparticles as the co-catalyst source, rather than photodeposition of Pt from solution. These pre-made Pt particles are too large to fit within the 1.9 × 2.1 nm TBAP-α nanopores and, hence, if these pores are the main sites for hydrogen evolution, then we might expect to see a very large decrease in rate compared to materials with photodeposited Pt particles, which might be small enough to form within the pore channels. This was not observed. The rate with the pre-formed nanoparticles was 813 μmol g−1 h−1; that is, only around 37% lower than for the photodeposited Pt sample under those catalysis conditions (1293 μmol g−1 h−1). ICP-MS analysis indicated the Pt contents were similar for both the photodeposited and the pre-made Pt materials and STEM imaging showed deposition had occurred in both cases.
Scanning electron microscopy imaging also indicated the inaccessibility of the TBAP-α pores: even the small, photodeposited Pt particles appeared to reside on the crystal surface, rather than within the pore channels (ESI, Fig. S57†). At present, it is not possible to state the precise role that micropores in TBAP-α play in the photocatalytic process, although the location of the Pt cocatalyst instead suggests that the external crystal surface is very important. It is conceivable that hole scavenging can occur within the pore channels, even if the main sites for photocatalytic hydrogen production are on the crystal surface; if so, then charge transport of the holes and electrons in the material would be important, which might explain the benefit of the extended π-stacks in TBAP-α.
Visible light absorption is also an important aspect of photocatalysis, and the three materials do show different absorption onsets in the solid state. These were determined to be 504, 477 and 449 nm for TBAP-α, TPyP-α and TPhP-α, respectively. However, based on experience with other materials, we would not expect this degree of blueshift to give the large, order-of-magnitude, change in photocatalytic activity between TBAP-α and the other two materials.
While the two π-stacked materials do outperform the non-π-stacking TPhP-α, the blue-shifted optical gap (ESI, Fig. S34†), lower wettability (ESI, Fig. S63 and S65†) and smaller driving force (ESI, Fig. S48†) of TPhP-α, relative to TPyP-α, could also contribute to the differences in HER between these two materials. As such, we cannot conclude that TPyP-α outperforms TPhP-α solely due to more optimal π-stacking.
The photocatalytic activity of TBAP-α was also investigated over a longer-term experiment and a slow loss of activity was observed over time (Fig. 6c); the rate reduced to 653 μmol g−1 h−1 after 20 hours, to 369 μmol g−1 h−1 after 60 hours and 156 μmol g−1 h−1 after 110 hours of photolysis. While this points to long-term instability, the sustained photochemical production of hydrogen at these rates over 110 hours is remarkable when one considers that TBAP-α is a molecular HOF where the framework is held together by hydrogen bonds and van der Waals interactions, which are weak in comparison to bonded frameworks such as CMPs and COFs.
The observed loss of activity is accompanied by a gradual decrease in crystallinity, as evident from a drop in X-ray peak intensity (Fig. 7c). No changes in the solution UV-vis spectrum (Fig. 7a) were observed after 110 hours photolysis; also, the photoluminescence spectrum (Fig. 7b) and the solution 1H NMR spectrum (ESI, Fig. S32†) remain unchanged. This suggests that the loss of catalytic activity is connected to changes in the TBAP packing, rather than any chemical degradation. Consequently, amorphous TBAP, isolated from the photocatalysis experiments after filtration, could be used to regenerate TBAP-α, after re-crystallisation from DMF/CHCl3.
It was thought that the loss in crystallinity of TBAP-α we observed during the photocatalysis experiments in aqueous solutions may be reduced by decreasing the polarity of the reaction medium. TBAP-α was therefore also tested using a 9:1 mixture of MeCN:water (Fig. S71†). In this case, the initial HER of 358 μmol g−1 h−1 remained approximately linear over the whole 118 hour experiment, with TBAP-α producing hydrogen at a rate of 347 μmol g−1 h−1 over hours 102–118, very close to the initial rate. This improved performance in the MeCN:water dispersant was also accompanied by a slower drop in the TBAP-α crystallinity over the experiment (Fig. S71†), as compared to the fully aqueous system (Fig. 7). We ascribe the lower rate in the MeCN:water medium to the smaller driving force for scavenger oxidation expected in this system as well as the fact that photodeposition of platinum appeared to occur less efficiently, with only 0.05 wt% Pd measured by ICP-MS (Table S6†).
To test the link between crystallinity and HER further, we deliberately produced amorphous material by the rapid precipitation of TBAP from a basic solution (ESI, general methods and Fig. S33†). We found that this amorphous TBAP had a very poor photocatalytic activity (Table 1 and Fig. 6b); with 1 wt% photodeposited Pt, the rate was 6 μmol g−1 h−1 (after an induction period). This rate is over 200 times lower than the rates observed for crystalline TBAP-α under equivalent conditions. This rationalises the steady loss of activity that is observed for TBAP-α as the crystallinity is reduced during photocatalysis. Amorphous TBAP under equivalent conditions but without added Pt produced no measurable hydrogen over 5 hours of irradiation.
DFT61,62 calculations that consider (i) an isolated TBAP molecule immersed in water, and; (ii) a stacked column of TBAP molecules from the TBAP-α structure suggest that the effect of packing on the TBAP potentials should be small (see ESI, Section 1.5, Table S8, Fig. S49 and S50†), in keeping with the similarity of the absorption spectra for amorphous and crystalline analogues. These DFT calculations also suggest that TBAP-α and its amorphous counterpart should both have the required driving force for proton reduction and ascorbic acid oxidation.
Like TBAP-α, the amorphous TBAP material also showed better activity when tested at pH 7; HER increased to 156 μmol g−1 h−1, (ESI, Fig. S24†), which is consistent with an increased driving-force for ascorbic acid oxidation as discussed above. This rate for amorphous TBAP is still 20-times lower than that of crystalline material under equivalent conditions (3108 μmol g−1 h−1), again suggesting that crystallinity and the packing of TBAP units is the dominant factor in the photocatalytic activity of this systems.
Time-correlated single photon counting (TCSPC) was employed to investigate the lifetime of the excited state in these materials. TBAP-α and amorphous TBAP suspended in water showed almost identical TCSPC spectra, when excited at 405 nm with average lifetimes of 2.12 and 2.13 ns, respectively. When using an aqueous suspension in ascorbic acid (0.1 M) significantly reduced fluorescence lifetimes were observed for both crystalline and amorphous TBAP (ESI, Fig. S47†). Even though TCSPC can only be used to study emissive states, these results suggest that exciton generation is similar for the amorphous and crystalline TBAP and that ascorbic acid is an effective hole scavenger for both materials.
By contrast, TPhP-α in ascorbic acid solution showed no significant quenching of the excited state lifetime (ESI, Fig. S69†). Similarly, neither TPhP-α nor TPyP-α showed a significant reduction in excited state lifetime in the presence of TEA (ESI, Fig. S70†). This indicates that poor interaction with the ascorbic acid or TEA hole scavenger might limit the activity of TPhP-α and TPyP-α for proton reduction in comparison to TBAP-α.
Taking these various experimental observations and DFT calculations together, we suggest that the dramatic difference in catalytic activity for crystalline and amorphous TBAP might be explained, at least in part, by restricted charge/exciton transport in the amorphous material to the active sites on the catalyst surface.47 While the exciton quenching kinetics by the hole scavenger between amorphous TBAP and TBAP-α appear to be similar, it is possible that transfer of a subsequent electron or electron polaron to a Pt active site is aided by interlayer conjugation in the crystalline HOF, which does not exist in the amorphous analogue. Materials with a higher degree of order and closer π–π stacking have been shown previously to have higher charge-carrier mobilities; for example, in conjugated polymers such as poly(thiophene-thieno[3,2-b]thiophene)66 and poly(cyclopentadithiophene-benzothiadiazole).67 Given this apparent benefit of extended π-stacked pyrene units, the low activity of TPyP-α, which has an overlap close to that of TBAP-α (0.86 vs. 0.91) is somewhat surprising. It is possible that in the case of TPyP-α, inefficient hole scavenging (as observed by TCSPC) prevents efficient generation of polarons from excitons reducing the amount of hydrogen that can be produced.54 Thus, any potential increase in charge-carrier mobility due to the packing of TPyP-α in the crystal is less relevant due to the low quantity of polarons generated.
While it is not yet possible to predict photocatalytic activity directly from crystal structure alone, CSP has strong potential to guide the discovery of new molecular building blocks by searching for crystal structure landscapes that contain stable structures with useful features, such as extended columns of π-stacked molecules with good overlap of the molecular faces. This is also useful for other energy materials, such as organic electronics.
Looking forward, CSP and ESF maps could be used to find other molecular systems that combine properties such as microporosity and π–π stacking. By extension, CSP could also be used to search for other packing motifs as our understanding of structure–activity relationships for organic photocatalysts expands. Because the computational search methods require no experimental input, this applies to hypothetical candidate molecules that have not yet been synthesised. Increasingly, these CSP calculations are much faster than experimental techniques such as molecular synthesis and crystallisation screens, allowing us to focus experimental resources on the most promising systems.
Conformers for each molecule were generated using a mixed torsional/low-mode sampling method72,73 implemented in Schrödinger's Maestro package,74 with energies modelling using the OPLS2005 force field.75 All unique conformers were then re-optimised using dispersion-corrected DFT, at the B3LYP-D3/6-311G(d,p) level of theory, leading to 4 conformers for each of TPhP and TPyP, differing in the orientation of phenyl and pyridyl groups with respect to the pyrene core, and 28 conformers TBAP, where orientations of the carboxylic acid groups expand the conformational space.
All conformers were used as starting points for CSP, which was performed with low discrepancy sampling of crystal packing variables, using the Global Lattice Energy Explorer software.76 Crystal structures were generated in the 25 most common space groups with one molecule in the asymmetric unit, then lattice energy minimised with molecular geometries fixed at their gas phase DFT geometries. Intermolecular interactions within the predicted crystal structures were calculated using the FIT atom–atom force field77 combined with atomic multipole electrostatics. Total energies were calculated as a sum of the force field intermolecular energy and the dispersion corrected DFT energy of the molecular conformer. All lattice energy calculations were performed with the DMACRYS crystal structure modelling software,78 and structures up to 100 kJ mol−1 above the global minimum were kept, as highly porous structures can occupy very high energy regions of the lattice energy landscape.36,79 This led to CSP landscapes with relatively large numbers of hypothetical structures; for example, for TBAP, the 28 conformers led to a CSP landscape with over 100000 independent hypothetical crystal structures within 100 kJ mol−1 of the global minimum in the 25 space groups studied (Fig. 2a). Full details of the crystal structure prediction methods are provided in the ESI.†
All CSP data, including predicted crystal structures, energies, accessible surface areas and stacking analysis, can be accessed at https://doi.org/10.5258/SOTON/D1015.
IP = −(E(N) − E(N − 1)) − 4.44 |
EA = −(E(N + 1) − E(N)) − 4.44 |
In the case of TBAP the IP and EA values were also calculated using an alternative strategy, starting from the crystal structure of TBAP-α. Initially, the experimental crystal structure of TBAP-α is energy minimised in a periodic DFT calculation using the B97-3c approach as implemented in Crystal17.89 Subsequently, three cluster models were cut out of the DFT optimised crystal structure, corresponding to one monomer, one monomer (1C) with a molecule above and below it, as well as the phenyl groups of the laterally adjacent molecules (1C+, see Fig. S49†), and an analogous structure with a tetramer in the centre (4C+). The IP and EA values of the three cluster models were calculated in the same way as for the isolated molecules discussed above, other than that in the last two cases we used the ONIOM QM/MM approach90 and described the molecule (fragments) around the monomer and tetramer using the UFF forcefield.91
All DFT optimised cluster models have been uploaded as electronic supplementary information.
Footnotes |
† Electronic supplementary information (ESI) available. CCDC 1919802–1919809. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/d0ta00219d |
‡ These authors contributed equally to this work. |
This journal is © The Royal Society of Chemistry 2020 |