Photocatalytic proton reduction by a computationally identi ﬁ ed, molecular hydrogen-bonded framework †

We show that a hydrogen-bonded framework, TBAP - a , with extended p -stacked pyrene columns has a sacri ﬁ cial photocatalytic hydrogen production rate of up to 3108 m mol g (cid:1) 1 h (cid:1) 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-ﬂ exible pyrene derivatives. Speci ﬁ cally, we show that energy – structure – function (ESF) maps can be used to identify molecules such as TBAP that are likely to form extended p -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.


Introduction
The de novo design of solid-state energy materials is challenging because function is dened by features that span multiple length scales.2][3][4][5][6][7][8] Organic materials are promising candidates for photocatalytic hydrogen production, but predicting the best combination of properties is difficult because the underlying structure-activity rules are poorly understood.So far, most studies involving heterogeneous organic photocatalysts have been conducted on carbon nitride materials or amorphous conjugated polymers, 9 where insolubility and lack of long-range order make thorough structural characterization difficult.Consequently, it is hard to deconvolute the structure-activity relationships in organic photocatalysts where the extended packing is poorly dened.
Molecular organic crystals have highly ordered structures that can be prepared in a modular way using solutionprocessable 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 oen 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 sacricial hole scavenger. 12ecently, 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 y 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,8tetrapyridin-4-yl pyrene (TPyP), and 1,3,6,8-tetra(4 0 -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 a-form, pyrene has a herringbone-type crystal packing arrangement, and this a-polymorph transforms to the structurally-related b-polymorph below 163 C. 15 By crystallising pyrene from CH 2 Cl 2 at pressures > 0.3 GPa, it is possible to obtain the denser g-polymorph, which features aligned pstacked columns of pyrene molecules. 16However, this structure was reported to dissolve in the crystallisation solvent during decompression, making it impossible for us to access this pstacked polymorph for photochemical experiments.The pyrene derivative, TPhP, crystallises from CH 2 Cl 2 /n-hexane, 17 to produce a solid that also lacks p-p 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.9][20] Similarly, labile C-H/N bonding interactions in molecular crystals, comprising pyridyl functionalised molecules, have been used to stabilise low density crystal packings. 21Hence, TPyP and TBAP both feature functional groups that have the potential to direct lowdensity, porous crystal packings with extended p-stacks, which is attractive for organic photocatalysts.
3][24] A crystal structure of TBAP was subsequently published while this study was in progress. 25olecular crystals are not subject to the same intuitive design rules associated with materials such as MOFs [26][27][28] and COFs, [29][30][31] because the crystallisation of organic molecules is governed by an interplay of many weak intra-and intermolecular interactions.It therefore remains a signicant challenge to design, from rst principles, new organic molecular crystals that feature porosity and/or extended p-conjugation.Porosity is also much less common in organic crystals because of their strong preference for close packing. 32At the atomistic level, discovering a new material involves substantiating a stable minimum on the free energy surface dened by chemical composition and the relative positions of atoms within an extended solid. 33Traditionally, 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 postanalysis of new materials to genuinely predictive methods that can be applied in advance of, or in tandem with, experiment.Such computationally-led approaches [34][35][36] promise to accelerate experimental searches and to unearth materials that might otherwise have remained undiscovered. 377][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,51These 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 congurational space of available crystal packings; this is a high dimensional problem including molecular positions, orientation, and unit cell dimensions, as well as exible 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][51][52] where a lack of conformational freedom leads to a reduced search space.TPhP, TPyP and TBAP are conformationally exible and in this study, we have implemented a CSP method that accounts for this exibility during structure searching, thus expanding the scope of CSP in the de novo design of functional materials.Specically, we performed CSP calculations with TPhP, TPyP and TBAP, to determine whether these molecules were likely to form porous or p-stacked structures that could increase their photocatalytic activity for hydrogen production from water.

Crystal structure landscapes
The crystal structure landscapes of TPhP and TPyP (Fig. 2b and  c) are similar and show the usual strong correlation between energetic stability and crystal density.We calculated the accessible surface area for all predicted crystal structures, and this showed a lack of porous structures in stable regions of the energy landscape for these two molecules.Thus, we can decide a priori that neither molecule is a promising candidate for the formation of porous solids.
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 energydensity 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 innite 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 p-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 MOFs 22-24 but without the inclusion of metals.During the course of this study, although aer these calculations, structure 1 was in fact reported independently by another research group, 25 conrming our prediction.

Using ESF maps to search for candidate photocatalysts
Although porosity has been linked to increased photocatalytic hydrogen evolution activity for conjugated organics polymers and polymeric carbon nitride, 47,53 many of the most active organic photocatalysts are in fact non-porous. 5,54This indicates that photocatalytic hydrogen evolution activity is not determined by a single factor: it is dependent on many variables.There is, however, compelling evidence that certain recurrent structural motifs are more desirable than others.In particular, studies involving COFs 13,55 indicate that extended p-p stacking is linked to increased photocatalytic hydrogen evolution activity.ESF maps offer a method to search for such structural motifs systematically, particularly as predictions are now tractable for more complex molecules and can oen be performed on a timescale that is much faster than the associated experiments (i.e., synthesis and crystallisation of a new candidate molecule coupled with photocatalytic measurements).Fig. 4 shows ESF maps that summarize the propensity of TPhP, TPyP, and TBAP to crystallise to form extended p-p stacks with signicant overlap between the pyrene cores.In this analysis, we dene an extended p-p 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 innite 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 p-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 p-p stacking; that is, these regions of enhanced stability are explicitly linked to p-p stacking.We would therefore predict a priori that TBAP is a more promising candidate for p-stacked HOFs than either of its two structural analogues, TPhP or TPyP.
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 p-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 t (Fig. 4d and ESI †).][58] From a predictive point of view, the distribution of p-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 pp overlap (overlap z 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 signicantly 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 p-stacks with strong p-p 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 p-stacking with a reasonable degree of overlap (0.863 and 0.857, respectively), an experimentally known structure of TPhP 17 that sits well above the global energy minimum displays no p-stacking, suggesting that these more 'optimal' structures may be kinetically difficult to access, at least for TPhP.

Crystallisation experiments
Crystallisation of TPhP.We sublimed the P2 1 2 1 2 1 a-polymorph of TPhP under vacuum at 425 C, which is the same polymorph as that reported from recrystallisation using CH 2 Cl 2 / n-hexane. 17The crystal packing in TPhP-a closely matches a predicted structure from the computed crystal structure landscape, 2 (Fig. 2b).The geometric agreement between these structures is excellent (ESI, Fig. S2 †), but this structure is predicted to lie above the global minimum energy structure.The energetic ranking of this structure might reect limitations of the force eld and rigid-molecule approach adopted in the CSP, but it might also indicate that more thermodynamically stable crystal structures of TPhP exist.
Crystallisation of TPyP.Due to its poor solubility, we were unable to crystallise TPyP from solvent.We therefore sublimed TPyP under vacuum at 450 C to afford a solvent-free structure, TPyP-a, which crystallised in the monoclinic space group P2 1 /n (ESI, Fig. S3 †).TPyP-a is reproduced accurately by the predicted global minimum energy crystal structure, 3 (Fig. 2c and ESI, Fig. S4 †) and features close-packed, interdigitated, columnar stacks of TPyP molecules, which are packed 3.95 A apart along the crystallographic a axis.
Crystallisation of TBAP.TBAP was crystallised by vapour diffusion of CHCl 3 into a saturated N,N-dimethylformamide (DMF) solution.The crystallisation solvent was then exchanged with acetone, and the acetone was subsequently removed from the crystal pores at 120 C under dynamic vacuum (ESI, Fig. S5-S15 and Table S5 †).The desolvated crystal structure, TBAP-a, which has C2/m symmetry is in excellent agreement with the predicted low-density TBAP structure, 1 (Fig. 2a and 3c).The crystal structure of TBAP-a, as determined at 25 C, had a packing distance between pyrene cores of 3.93 A, which is slightly longer than 3.85 A reported for the g-polymorph of pyrene recorded at 0.3 GPa and 25 C. 16 We carried out sorption measurements to estimate the surface area of TBAP-a, which was found to match our ESF map predictions (Brunauer-Emmett-Teller surface area, SA BET ¼ 2001-2270 m 2 g À1 , ESI, Fig. S19-S22; † ESF map prediction ¼ 2273 m 2 g À1 ).This is also in agreement with the report by Yin et al. 25 These three crystalline solids, TPhP-a, TPyP-a, TBAP-a, display different p-p stacking modes of pyrene cores, and the TBAP-a 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.

Photocatalysis experiments
TPhP-a, TPyP-a, TBAP-a have light absorption onsets well into the visible region (Fig. S34 †).We therefore tested activity for sacricial photocatalytic hydrogen production under visible light (l > 420 nm, 300 W Xe light source) with 1 wt% platinum co-catalyst, using ascorbic acid as the sacricial hole-scavenger.TBAP-a was by far the most active material under these conditions with an initial hydrogen evolution rate (HER) of 1293 mmol g À1 h À1 (Table 1).This is the rst example of a crystalline porous HOF that shows appreciable photocatalytic hydrogen evolution from water under sacricial conditions, challenging the paradigm that covalent networks such as COFs or extended polymer chains are required.These photocatalytic rates are signicantly higher than for amorphous pyrene-based conjugated microporous polymers (CMPs) (HER l>420 nm ¼ 174 mmol g À1 h À1 using diethylamine as the hole-scavenger). 12The external quantum efficiency (EQE) of TBAP-a was estimated to be 4.1% at 420 nm using monochromatic light from a LED source, which is higher than many conjugated polymer catalysts such as a benzodithiophene-bipyridine CMP, PCP4e (EQE 350 nm ¼ 0.34% using triethylamine/water mixtures) 8 and a tricyano-benzene-centered phenylenevinylene-co-terphenylene polymer network, OB-POP-3, (EQE 420 nm ¼ 2.0% using water triethanolamine mixtures), 59 but lower than dibenzo[b,d]thiophene-CMP S-CMP3 (EQE 420 nm ¼ 13.2%, using water/methanol/triethylamine mixtures), 47 and certain linear conjugated polymers, such as poly(dibenzo[b,d] thiophene) P10 (EQE 420 nm ¼ 11.6%, using water/methanol/ triethylamine mixtures). 54At 470 nm, the EQE for TBAP-a is reduced to 1.2% and no signicant activity is observed at 595 nm.As such, the photocatalytic efficiency follows the absorption prole of TBAP-a (ESI, Fig. S28 †).
No hydrogen production is observed in the dark or in the absence of the TBAP-a.When D 2 O was used as the proton source, D 2 production was mostly observed (ESI, Fig. S29 †), with a small amount of H 2 , probably due to H-D exchange with protons in the non-deuterated ascorbic acid. 60Essentially 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-a was much less active and produced hydrogen at a rate of just 18 mmol g À1 h À1 , although the basic pyridyl groups in TPyP meant that a signicant proportion of the TPyP-a catalyst (>50 wt%) dissolved in the acidic medium during this measurement.Consequently, a direct comparison between the HER activity of TPyP-a and TBAP-a in ascorbic acid is not possible.TPhP-a is stable in ascorbic acid but had an even lower HER of 2 mmol g À1 h À1 .
To allow for a direct comparison between TBAP-a and TPyPa, 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-a produced no measurable amount of hydrogen over a ve-hour period.By contrast, TBAPa was even more active with a HER of 3108 mmol 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 onehole 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-a 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 TPyPa when changing the pH level from pH 2.6 to pH 7, but we note that the partial dissolution of the material and the inuence of protonation on the substrate electronics under acidic conditions could play a role.
When 5 vol% triethylamine in water solution was used as the sacricial system, TBAP-a was found to dissolve and no measurable hydrogen was produced under visible light irradiation.TPyP-a and TPhP-a were both stable in 5% triethylamine but had low HERs of 40 and 6 mmol g À1 h À1 , respectively.None of the materials produced hydrogen when tested in a 5 vol% TEA solution adjusted to pH 7.
At rst glance, these results could suggest that porosity is the dominant factor for hydrogen evolution, because the porous TBAP-a structure greatly outperforms the non-porous analogues, even though TPyP-a also contains signicantly overlapped p-stacking.This may however be an over-simplication.As evident from water sorption isotherms, TBAPa has no signicant water uptake between partial pressures of 0.2-0.6 (Fig. S52 †).There is signicant 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 sacricial agent, ascorbic acid, and in this respect, water sorption isotherms may not reect the photocatalysis conditions.
Perhaps more signicantly, 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-a had a dramatically reduced rate of 59 mmol 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. 63The 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,65Element analysis (ICP-MS; ESI, Table S6 †) and STEM imaging were employed to conrm 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-a 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 †).Photolysis experiments using 4 wt% Pt were repeated for three different batches of TBAP-a 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 m 2 g À1 .Likewise, static light scattering showed only modest variations in the particle size distribution between batches (average diameters of 9.25-15.52mm, 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 t within the 1.9 Â 2.1 nm TBAP-a 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 mmol g À1 h À1 ; that is, only around 37% lower than for the photodeposited Pt sample under those catalysis conditions (1293 mmol 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-a 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-a 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 benet of the extended p-stacks in TBAP-a.
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-a, TPyP-a and TPhP-a, respectively.However, based on experience with other materials, we would not expect this degree of blueshi to give the large, order-ofmagnitude, change in photocatalytic activity between TBAPa and the other two materials.
While the two p-stacked materials do outperform the non-pstacking TPhP-a, the blue-shied optical gap (ESI, Fig. S34 †), lower wettability (ESI, Fig. S63 and S65 †) and smaller driving force (ESI, Fig. S48 †) of TPhP-a, relative to TPyP-a, could also contribute to the differences in HER between these two materials.As such, we cannot conclude that TPyP-a outperforms TPhP-a solely due to more optimal p-stacking.
The photocatalytic activity of TBAP-a 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 mmol g À1 h À1 aer 20 hours, to 369 mmol g À1 h À1 aer 60 hours and 156 mmol g À1 h À1 aer 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-a 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 aer 110 hours photolysis; also, the photoluminescence spectrum (Fig. 7b) and the solution 1 H 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 aer ltration, could be used to regenerate TBAP-a, aer re-crystallisation from DMF/CHCl 3 .
It was thought that the loss in crystallinity of TBAP-a we observed during the photocatalysis experiments in aqueous solutions may be reduced by decreasing the polarity of the reaction medium.TBAP-a was therefore also tested using a 9 : 1 mixture of MeCN : water (Fig. S71 †).In this case, the initial HER of 358 mmol g À1 h À1 remained approximately linear over the whole 118 hour experiment, with TBAPa producing hydrogen at a rate of 347 mmol 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-a 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 mmol g À1 h À1 (aer an induction period).This rate is over 200 times lower than the rates observed for crystalline TBAP-a under equivalent conditions.This rationalises the steady loss of activity that is observed for TBAP-a 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.
DFT 61,62 calculations that consider (i) an isolated TBAP molecule immersed in water, and; (ii) a stacked column of TBAP molecules from the TBAP-a 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-a and its amorphous counterpart should both have the required driving force for proton reduction and ascorbic acid oxidation.
Like TBAP-a, the amorphous TBAP material also showed better activity when tested at pH 7; HER increased to 156 mmol g À1 h À1 , (ESI, Fig. S24 †), which is consistent with an increased drivingforce 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 mmol 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-a 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) signicantly reduced uorescence 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-a in ascorbic acid solution showed no signicant quenching of the excited state lifetime (ESI, Fig. S69 †).Similarly, neither TPhP-a nor TPyP-a showed a signicant 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-a and TPyP-a for proton reduction in comparison to TBAP-a.
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. 47While the exciton quenching kinetics by the hole scavenger between amorphous TBAP and TBAPa 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 p-p stacking have been shown previously to have higher charge-carrier mobilities; for example, in conjugated polymers such as poly(thiophenethieno[3,2-b]thiophene) 66 and poly(cyclopentadithiophenebenzothiadiazole). 67Given this apparent benet of extended p-stacked pyrene units, the low activity of TPyP-a, which has an overlap close to that of TBAP-a (0.86 vs. 0.91) is somewhat surprising.It is possible that in the case of TPyP-a, inefficient hole scavenging (as observed by TCSPC) prevents efficient generation of polarons from excitons reducing the amount of hydrogen that can be produced. 54Thus, any potential increase in charge-carrier mobility due to the packing of TPyP-a in the crystal is less relevant due to the low quantity of polarons generated.

Conclusions
In summary, we report here the rst example of a molecular crystalline HOF with high photocatalytic activity.The porous, hydrogen-bonded structure of TBAP-a was anticipated using CSP methods and ESF maps.TBAP-a has high photocatalytic activity for proton reduction under sacricial conditions, while the amorphous TBAP material has an activity that is up to 200 times lower over the rst 5 hours of photolysis.To our knowledge, this is the rst example where two different solid phases of the same organic materials have been studied for photocatalysis, i.e.where the solids are chemically identical and differ only in their solid-state packings, and it allows us to deconvolute the effects of molecular structure and crystal structure.The huge difference in catalytic activity for crystalline and amorphous forms of the same molecule shows unambiguously that solid-state crystal packing can have a large effect on photocatalytic activity.Detailed comparisons of the activities and crystal structures of TBAP, TPyP and TPhP suggest that pstacked columns with strong p-p overlap are a desirable packing motif.The specic role of porosity in the photocatalytic activity of TBAP-a is less clear, but porosity is, in general, a desirable feature for heterogeneous catalysis and it also opens up other opportunities, such as doping or dye sensitization. 68hile 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 pstacked 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 nd other molecular systems that combine properties such as microporosity and p-p 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.

Experimental
Crystal structure prediction TPhP, TPyP and TBAP all have multiple low-energy conformers, which required us to expand our CSP and ESF map methodology to consider the ensemble of possible conformers to ensure that the crystal packing arrangements are comprehensively sampled and that experimentally realizable structures are not missed. 69,70For conformationally exible molecules, it is not obvious which form will crystallise experimentally; indeed, observed crystal structures are sometimes based on molecular conformers that have high energies in the gas phase. 71This means that a wide range of conformers must be considered during CSP.
Conformers for each molecule were generated using a mixed torsional/low-mode sampling method 72,73 implemented in Schrödinger's Maestro package, 74 with energies modelling using the OPLS2005 force eld. 75All 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 soware. 76Crystal structures were generated in the 25 most common space groups with one molecule in the asymmetric unit, then lattice energy minimised with molecular geometries xed at their gas phase DFT geometries.Intermolecular interactions within the predicted crystal structures were calculated using the FIT atom-atom force eld 77 combined with atomic multipole electrostatics.Total energies were calculated as a sum of the force eld intermolecular energy and the dispersion corrected DFT energy of the molecular conformer.All lattice energy calculations were performed with the DMACRYS crystal structure modelling soware, 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,79This 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 100 000 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.

DFT potential calculations
The vertical ionization potential (IP) and electron affinity (EA) of TBAP, TPhP and TPyP were calculated using a DDFT approach.First the ground state geometry of each as an isolated molecule was optimised using the B97-3c approach by Grimme and coworkers. 62Next the energy of each of the molecules in its neutral (E(N)), cation (E(N À 1)) and anionic (E(N + 1)) state were obtained from single-point calculations using the B3LYP functional [80][81][82][83] and the 6-31G** basis-set. 84,85Finally, IP and EA were calculated from: where all energies are in eV and the subtraction by 4.44 converts the calculated IP and EA from the vacuum to the standard hydrogen electrode scale.The B3LYP single-point calculations were performed using Gaussian16 (ref.86) and employed the PCM solvation model 87 to describe the aqueous environment of the molecules near the molecular solid-solution interface.The B97-3c calculations were performed using Turbomole 7.3 (ref.88) and employed no solvation model as to as closely as possible match the computational set-up of the periodic DFT calculations on TBAP-a.
In the case of TBAP the IP and EA values were also calculated using an alternative strategy, starting from the crystal structure of TBAP-a.Initially, the experimental crystal structure of TBAPa is energy minimised in a periodic DFT calculation using the B97-3c approach as implemented in Crystal17. 89Subsequently, 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 approach 90 and described the molecule (fragments) around the monomer and tetramer using the UFF forceeld. 91ll DFT optimised cluster models have been uploaded as electronic supplementary information.

Synthesis
Synthesis of 1,3,6,8-tetra(4 0 -carboxyphenyl)pyrene (TBAP).TBAP was synthesised according to literature routes. 22 Synthesis of 1,3,6,8-tetrapyridin-4-yl pyrene (TPyP).1,3,6,8-Tetrabromopyrene (1.04 g, 2 mmol), 4-pyridinylboronic acid (983 mg, 8 mmol), N,N-dimethylformamide (200 mL) and K 2 CO 3 (50 mL) were added to a ask and degassed by N 2 bubbling for 30 minutes.[Pd(PPh 3 ) 4 ] (40 mg, 0.035 mmol) was added and the solution was degassed for a further 10 minutes before heating to 145 C for 48 hours.Aer cooling the mixture was poured into water (1 L) and stirred for 30 minutes.The precipitate was collected by ltration and washed with, water (100 mL), methanol (mL) and dichloromethane (100 mL) before drying under vacuum.The product was obtained as a green solid (986 mg, 1.92 mmol, 96%). 1  Crystallisation, solvent exchange and activation of TBAP-a.500 mg of as-synthesised TBAP was covered by DMF (40 mL)   in a large vial.The mixture was sonicated for 10 min and le overnight so all remaining undissolved material settled at the bottom of the vial.20 mL of the supernatant solution was the ltered into two 40 mL vial using a 0.45 mm PTFE syringe lter to remove any particulates.For batches 1 and 3, the 40 mL vials were capped with a septum that had been pierced using a needle and these vials were placed in a sealed chamber containing chloroform.For batch 2, the 40 mL sample vials were le uncapped and placed in a sealed chamber containing chloroform.Vapour diffusion of chloroform into the TBAP solution was carried out for until the vials containing the TBAP material were nearly full of solvent.Most of the solvent was then removed via syringe, until the level of solvent was just above the TBAP crystals.Acetone (10 mL) was then injected into the solution and subsequently removed via syringe.This process was repeated twice more, and enough acetone was then added to ll the vial completely.This solvent was removed and then replenished every 12 hours for 5 days, which yielded the solvent exchanged the TBAP$x (acetone) solvate (Fig. S9 and S10 †), the yellow crystalline material was ltered off and allowed to dry under ambient conditions.The solid was then evacuated at 120 C for 14 hours to give the activated material (Fig. S12 and S18 †).Complete removal of the solvent was conrmed by the absence of resonances relating to DMF, CHCl 3 , or acetone in the 1 H NMR spectra of the activated sample (Fig. S11 †).
Generation of amorphous TBAP phase.As synthesised TBAP was fully dissolved 2 M KOH (aq.).The aq. solution was then placed in an ice bath for 10 minutes, and excess concentrated HCl (aq.) was added in one portion to neutralise the solution and rapidly precipitate the TBAP material.The TBAP material was collected by ltration and washed with copious amounts of water.The ltrate was then suspended in water (20 mL), sonicated for 1 hour, and re-ltered to ensure all remaining salts were washed out.This process was repeated three times in total.The amorphous TBAP sample was dried under vacuum at 120 C and a PXRD pattern of the amorphous material was recorded, see Fig. S33.† Sublimation of TPhP and TPyP.Crystals of TPhP were obtained by sublimation at 425 C and pressure of 5 Â 10 À4 hPa.Crystals of TPyP were obtained by sublimation at 450 C and pressure of 5 Â 10 À4 hPa.

Fig. 2
Fig. 2 Energy-density distributions of predicted crystal structures (a) TBAP, (b) TPhP and (c) TPyP.Each point corresponds to a distinct crystal structure, coloured by calculated accessible surface area (m 2 g À1 ), and calculated using a 1.2 A probe radius.The observed crystal structures are labelled.Fig. 3 Predicted porous structure of TBAP.(a) Viewed down the channels with voids with the surface indicating the solvent accessible surface (probe radius 1.2 A) and (b) viewed perpendicular to the channels, showing the pyrene p-stacking.(c) An overlay of the predicted crystal structure 1 (red) and the experimental crystal structure TBAP-a (blue).

Fig. 3
Fig. 2 Energy-density distributions of predicted crystal structures (a) TBAP, (b) TPhP and (c) TPyP.Each point corresponds to a distinct crystal structure, coloured by calculated accessible surface area (m 2 g À1 ), and calculated using a 1.2 A probe radius.The observed crystal structures are labelled.Fig. 3 Predicted porous structure of TBAP.(a) Viewed down the channels with voids with the surface indicating the solvent accessible surface (probe radius 1.2 A) and (b) viewed perpendicular to the channels, showing the pyrene p-stacking.(c) An overlay of the predicted crystal structure 1 (red) and the experimental crystal structure TBAP-a (blue).

Fig. 4
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.

Fig. 5
Fig. 5 STEM images of TBAP-a 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.

Fig. 7 (
Fig.7(A) UV-vis and (B) PL (l exc ¼ 360 nm) in DMSO solution and (C) PXRD of TBAP-a loaded with 1 wt% Pt, after photocatalysis experiments performed in ascorbic acid (0.1 M), illuminated with a 300 W Xe light source fitted with a l > 420 nm cut-off filter.After photocatalysis experiments, the samples were collected by filtration and air-dried before analysis.

Table 1
Photocatalytic activity of the materials a Catalyst (25 mg) loaded with 1 wt% Pt, from in situ photodeposition of H 2 PtCl 6 , suspended in water and scavenger (25 mL), irradiated with a 300 W Xe light source tted with a l > 420 nm lter.The HER was determined over ve hours.b HER calculated over 20 hours.c Material fully or partially dissolved under these conditions.