Open Access Article
Catherine Mollart
*,
Ellena Sherrett,
Patrick Heasman,
Michael J. G. Peach
,
Adam Rowling,
Lewis J. Beck,
Ellie Varley,
David Seed and
Abbie Trewin
*
Department of Chemistry, Lancaster University, Lancaster, LA1 4YB, UK. E-mail: a.trewin@lancaster.ac.uk; c.l.mollart@bham.ac.uk
First published on 24th November 2025
An alternative route to synthesise TCNQ-CTF by using trifluoromethanesulfonic (TFMS) acid catalysis is presented, in comparison to a previously reported ZnCl2-catalysed synthesis. The new synthetic route yields a polymer with additional structural diversity compared to the previously reported material. The composition of the framework is rationalised by ‘artificial’ acid-catalysed synthesis of TCNQ-CTF, together with a novel approach to structural feature identification, with a range of alternative structural features appearing that were not present in the previously reported polymer formed by ZnCl2 catalysis. These results will inform the design of new CTF materials with additional functionality and broader applications.
One class of such functional polymer materials are covalent triazine frameworks (CTFs), with applications in areas such as carbon dioxide capture,11,12 drug delivery,13 heterogeneous catalysis,14 and photocatalysis.11,15 These are microporous polymer materials with a covalent, π-conjugated skeleton, first reported in 2008 by Kuhn and co-workers.16 The original CTF material, CTF-1, was formed using a high-temperature ionothermal synthesis of 1,4-dicyanobenzene, catalysed using one equivalent of molten ZnCl2.16 The resulting material was highly porous (with a Brunauer–Emmett–Teller [BET] surface area of 791 m2 g−1), chemically and thermally stable, and crystalline, with an analogous stacked topology to covalent organic framework (COF)-1.16
Whilst molten ZnCl2 is a suitable catalyst for CTF formation, it proves challenging to separate the catalyst from the product.11 This led Ren and co-workers to investigate more facile synthetic routes to CTF synthesis; operating at room temperature or under microwave-assisted conditions in the presence of trifluoromethanesulfonic (TFMS) acid catalyst, with a much-reduced reaction time and no reduction in chemical stability.11
However, unlike the original ZnCl2-catalysed CTF materials, which were synthesised under thermodynamic control, the acid catalysed CTFs are kinetic products, exhibiting little to no order within the structures. Additionally, the TFMS-catalysed CTFs synthesised at room temperature and under microwave-assisted conditions using the same monomer as CTF-1 (named P1 and P1M, respectively) display reduced porosity in comparison to the crystalline CTF-1 material, demonstrating the synthetic route employed in CTF preparation is a key influence on the structure and properties of the polymer.17 A comparison of the synthetic routes utilised to make CTF-1, P1M and P1, and the properties of the resulting materials, are given in Scheme S1 and Table S1.
CTF materials synthesised from tetracyanoquinodimethane (TCNQ) by Li and co-workers (Scheme S2), named TCNQ-CTFs, showed intriguing properties as electrodes for supercapacitor materials due to their high surface areas and nitrogen content (with observed values of greater than 3600 m2 g−1 and 8%, respectively).18 Yet despite being synthesised using a high temperature molten ZnCl2-catalysed reaction, which has previously been the default synthetic route to crystalline CTFs, the resulting TCNQ-CTFs were amorphous as judged by powder X-ray diffraction (PXRD, Fig. S2).18
We have ourselves previously rationalised the FT-IR spectrum of the amorphous P1 material by generating snapshots of its structure using our in-house code Ambuild, which is designed to model amorphous microporous polymers under kinetic control.17 By generating P1 structures using Ambuild, a range of structural diversity beyond the expected triazine ring thermodynamic product was observed, arising from the kinetic control of the reaction. IR spectra were calculated for fragments observed in the generated structure; it was found the peaks, in combination with adsorbed CO2, H2O and H3O+ guests, matched the additional peaks in the P1 FT-IR spectrum compared to the spectrum of the crystalline CTF-1 product.
In this work, we have synthesised TCNQ-CTF using the TFMS-catalysed reaction (Fig. 1). We rationalise the structure of TCNQ-CTFs using 100 models of TCNQ-CTF clusters, each grown independently to contain 10 TCNQ building blocks. By generating this series of small clusters, we can assess the range of structural diversity present and the prevalence of each observed structural feature across our dataset. This approach is a variation on a protocol we recently established in ref. 19 for the study of amorphous pyrene-based conjugated microporous polymers.19 Further we establish an automated protocol for the substructure counting process by searching for the SMILES (simplified molecular line entry system) strings of each of the identified structural features within the clusters using RDKit in Python. This allows our models to be rapidly analysed in a more high-throughput manner using a method readily adaptable to the study of new materials. Using fragments representing the identified structural features, we calculated the infrared (IR), nuclear magnetic resonance (NMR) and UV–Vis spectra of each structural feature observed within the models, to produce an overall spectrum representative of the experimental material. The increased functionality of our materials in comparison with the original synthesis route results in a structure that more closely resembles covalent adaptable networks (CANs)20–22 and thus increases their range of functionalities. These include but are not limited to shape tuneability, recyclability and self-healing properties due to their smart response to thermal or photochemical stimuli.23 This is particularly important in the current global strive towards reduction, reuse and recycling of waste materials, in line with the United Nations Sustainable Development Goals.24
Previously, we have successfully used Ambuild to generate P1 network models.17 Here we can apply the same strategy with details of the mechanism given in Fig. S1. Full details of the Ambuild network generation procedure can be found in section S2, however a brief description of the process is as follows. Firstly, one TCNQ building block is seeded into a cell and a growBlocks step used to bond a second block to the first, thus mimicking the next step in the mechanism. A free end group is randomly identified by the Ambuild code, and a further building block is attached, following bonding rules that a carbon end group can only bond to a nitrogen end group and vice versa. The orientation of the attached building block is determined by the cap atom vectors and any close contacts with other atoms in the cell. The hydrogen cap atoms are removed once a successful bond is formed, leading to intermediate 2. This process is repeated to produce intermediate 3. A zipBlocks test is undertaken to form a triazine ring from intermediate 3, where each end group is tested for the possibility of forming a bond between it and neighbouring end groups by measuring the distance and angle between them against pre-determined criteria. If the zipBlocks test is successful, the two cap atoms are removed and a bond is formed between the end groups. This is repeated until 10 building blocks have been grown onto the network, and any further bonding has been tested, shown in Fig. 2.
We have previously determined the parameters used within Ambuild to describe the bonds, angles and dihedrals using a simplified test mechanism, calibrating against structures obtained from the polymer consistent forcefield (PCFF). Full details of this parameterisation can be found in ref. 17. All structural feature analysis was performed using Gaussian 09,25 and full details of the model chemistries, calculation methods and spectra generation can be found in section S6.
Fig. 3 outlines the structural features identified via our computational synthesis approach and their prevalence across the 100 clusters produced. A full list of the structural features identified and their prevalence in the cluster models can be found in section S5. Amorphous networks are very challenging to analyse due to the limited characterisation techniques available. Some chemical functionality can be identified through IR spectroscopy or solid-state NMR, but it is often difficult to distinguish between similar groupings of chemical functionalities or structural features. Therefore, it is possible to miss functionality that may not be immediately obvious or is unexpected from the complete reaction taking place. This is a particular challenge in amorphous networks, where intermediates, otherwise not expected to be observed, may remain as they are fixed into their position within the polymer structure and so cannot react further, particularly in later stages of the network formation. Although triazine rings were one of the most common features present, with 19 counted across the clusters, the distribution was dominated by the closely related triazocine feature alongside macrocycles triangle C and 12-ring (42, 39 and 25 observed, respectively). Previously identified in our P1 study,17 4-ring and 8-ring features were also observed in our TCNQ-CTF clusters. This range of structural features contributes to the amorphicity of the material. We expect that this broader structural diversity compared to P1 arises because of the difference in structure and functionality of the monomer. The 1,4-dicyanobenzene monomer of P1 possesses two cyano groups para to one another, whereas the TCNQ monomer has four cyano groups, with two at each end separated by a single carbon atom. This reduction in distance, and variation in angle between cyano groups in the TCNQ monomer, leads to a wider range of bonding opportunities between the building blocks, introducing much greater diversity in macrocycle formation.
This contrasts with the structural composition of P1, where it was found that triazine rings were most prevalent as one would expect for this reaction pathway.17 Additionally, several ‘pre-rings’ were identified in the clusters, including pre-4-rings, pre-triazines and pre-triazocines. The pre-4-ring was the most prevalent feature overall, with 314 counted across the models. Notably, no pre-8-rings were observed, a feature identified in our P1 study.17
The experimental IR spectrum of TCNQ-CTF synthesised using a TFMS catalysed synthesis, shown in Fig. 4(a) in black, exhibits peaks at ∼1200 cm−1, ∼1500 cm−1, ∼2100 cm−1, ∼2200 cm−1, ∼2350 cm−1 and a broad peak with features at ∼3200 cm−1, ∼3350 cm−1, ∼3450 cm−1 and ∼3600 cm−1. The peaks at ∼2200 cm−1 and ∼2350 cm−1 can be attributed to unreacted cyano groups. We accounted for the presence of atmospheric CO2 by calculating its IR spectrum (Fig. S46) and obtained a peak overlapping a cyano group peak at ∼2350 cm−1. Although the peaks at ∼1500 cm−1 and ∼1200 cm−1 have previously been assigned to triazine rings, and indeed our calculated spectra in this study support this, all ring and pre-ring calculated spectra possess peaks in these regions originating from the C
C and C
N bonds produced upon macrocycle formation. This demonstrates the difficulty in rationalising a diverse, amorphous structure using only IR spectroscopy—the complexity of the spectrum cannot be rationalised by triazine rings alone. Upon analysis of the calculated spectra (Fig. S15–45(a)), the peak at ∼2100 cm−1 can be attributed to charged pre-ring intermediates present in the network, specifically stretches between the carbocation and neighbouring carbon and nitrogen atoms.17 The peaks at ∼3200 cm−1 and ∼3350 cm−1 arise from C–H stretches in the phenyl rings. It is well-known that porous materials readily adsorb water and carbon dioxide,17 and although TCNQ-CTF was found to be largely non-porous, these species are still able to adsorb onto the material's surface. We calculated the IR spectrum of a water molecule (Fig. S47(a)), giving a peak at ∼3650 cm−1 which likely explains the experimental feature at ∼3600 cm−1.
![]() | ||
| Fig. 4 The (a) IR, (b) NMR, and (c) UV–Vis spectra of our Ambuild-derived computed models (blue), weighted according to the statistics in Fig. 3, and reduced models (orange), compared to the experimental spectra of the synthesised TFMS acid catalysed TCNQ-CTF (black). See Fig. S52 for an expanded version of the IR spectrum. | ||
Due to the presence of acid from the TFMS synthesis, it is reasonable to expect hydronium ions to form from water, and its calculated IR spectrum (Fig. S48(a)) exhibits a strong peak at ∼3450 cm−1 and further peaks at ∼3400 cm−1, ∼1600 cm−1 and ∼900 cm−1. Each of these peaks is observed in the experimental spectrum. The combined weighted IR spectrum (Fig. 4(a), blue) derived by statistically combining the individual spectra of 25 identified fragments (Fig. S16–40) is a good overall match to experiment. The additional peaks arising from carbon dioxide, water and the hydronium ion can be found in Fig. S46–48.
The weighted calculated NMR spectrum (Fig. 4(b), blue) is a good match to experiment (the full experimental spectrum is shown in Fig. S8) between 70–200 ppm, suggesting the structural features identified in our clusters dominate within the synthesised material.
As none of our calculated models show peaks at chemical shifts larger than 200 ppm, we choose to plot our calculated NMR spectra over the range of ∼0–200 ppm. The most intense peak of the calculated spectrum at 120 ppm can be attributed to unreacted cyano groups present in the material, which will naturally be overestimated in the fragment approach considered here. The peaks found between 90–125 ppm arise from the carbon atom adjacent to the cyano groups, becoming more deshielded when one of the cyano groups becomes an imine when participating in bonding within a macrocyclic ring.
Peaks from 150–165 ppm belong to phenyl group carbons, becoming more deshielded the closer to the cyano groups they are. In the calculated spectra, the only peak unique to one specific feature is that at around 200 ppm from the carbon atoms participating in the 4-ring. The experimentally observed feature at 230 ppm is suggested to be a spinning side band of the intense peak at 130 ppm, with an assigned matching band observed as part of the dominant peak at ∼30 ppm.
The experimental UV–Vis spectrum (Fig. S6) has a broad high intensity peak from 1.5–3.5 eV and further lower-intensity peaks at 4.1 eV and 4.8 eV. This suggests the material has an optical gap of around 1.5 eV. The weighted calculated UV–Vis spectrum (Fig. 4(c), blue) shows a broad set of peaks between ∼1.5 and 3.4 eV in-line with the experimental UV–Vis spectrum. The peaks at lower values between 1.5 and 2.2 eV appear strongly in the 4-,8- and 10-ring fragments, whereas peaks between 2.2 and 2.8 eV are common throughout all fragments but are particularly dominant in triazocine fragments (Fig. S27 and S28) and triangles A (Fig. S29 and S34), and triangle E (Fig. S23). The peak at ∼3.1 eV is particularly dominant in rectangles A and B (Fig. S22 and S23), and triangles C and D (Fig. S31 and S32). The small peak at ∼4.1 eV is not prominent in any fragment but does feature in rectangle A (Fig. S22), triangle A (Fig. S29), and macrocycle B (Fig. S39).
We note that five of our identified fragments yield very low (≪1 eV) excitation energies, with calculations suggesting that these fragments may have triplet (i.e., diradical) character. Given the small timescales of our Ambuild artificial synthesis calculations compared to experiment, this “instability” may be representative of an incomplete artificial reaction. In the fully reacted system, we would therefore expect to see fewer pre-rings than we currently observe in the clusters. To avoid artificially biasing the resultant weighted spectra, we have only included contributions from the 25 identified “stable” fragments—the spectra in Fig. 4 do not include any contributions from these five “unstable” fragments.
We emphasise however that here, we are considering isolated fragments, which do not account for any of the longer-range structure or interactions of the real-world polymer, which may potentially stabilise such motifs.
We also note that we are considering around 40 fragments identified from a subset of 100 small clusters, so we may not have captured all possible structural environments. In addition to this, several of the studied fragments are similar in composition and functionality. It is incredibly challenging to assign specific activities to single sites with absolute certainty, as similar structures will give rise to similar characterisation spectra. We also acknowledge that as our material is amorphous, there may be different local structural environments and functionalities. This may increase the complexity in establishing structure–property relationships for target applications with our current knowledgebase.
Assessing the experimental and weighted calculated spectra as a whole, we can see that the set of peaks in the NMR at ∼30 ppm and peaks between 2.6 and 3.0 eV and above ∼3.8 eV in the UV–Vis spectra are not fully rationalised. In particular, whilst part of the peaks at ∼30 ppm in the NMR can be attributed to the spinning side band for the intense peak at 130 ppm, the peak is too large to be solely due to that side band. Our models could not identify the remaining contributions to this peak. This suggests that either additional features are produced in the experimental synthesis that we are yet to identify, or that adsorbed species, such as residual solvent, are present.11 To explore the possibility of as-yet unidentified structural features being the origin of these peaks, we consider the potential for alternative reactions. We note the peak is within the region we would expect aliphatic carbons to appear. It is known that TFMS is able to catalyse the hydrogenation of alkanes,26 therefore we consider this possibility.
Two models were considered based on a fragment of the idealised CTF polymer shown in Fig. 1(a). In the first model, one of the two alkene bonds is reduced, and in the second model both alkene bonds are reduced. The NMR spectrum was then calculated for both fragments. Full details are given in section S7. We see that both models have a dominant peak at ∼30 ppm (Fig. S50(b) and S51(b)), which when combined and compared to the experimental spectrum (Fig. 4(b), orange) are a good match to the previously unidentified set of peaks. Therefore, we assign these as due to the presence of aliphatic groups caused by the reduction (full or partial) of the 1,4-cyclohexadiene rings. The UV–Vis spectrum was also calculated for both fragments (Fig. S50(c) and S51(c)) which when combined is compared to the experimental spectrum in Fig. 4(c) – orange line. A good match is again observed, with a peak at ∼2.9 eV (originating from reduced model 1) and a peak starting at ∼3.9 eV (originating from reduced model 2). For completeness, the IR spectrum was also calculated for both fragments, shown in Fig. S50(a) and S51(a), which were again combined and compared to the experimental spectrum in Fig. 4(a) – orange line. We can see that the combined spectrum is again consistent with the experimental spectrum.
Considering the experimental spectra and the calculated spectra together, they can only be rationalised through our approach of mimicking the synthetic method to identify the full range of structural features that can form.
Our approach highlights the breadth of structural features that form when a TFMS-acid catalysed approach is used to generate TCNQ-CTF. A number of these structural features identified are difficult to synthesise in isolation but have a wide range of potential uses, and so are exciting target chemistries for incorporation into polymer materials. For example, macrocyclic C
N groups can be classed as secondary ketimines, which are popular in coordination chemistry as salen ligands and Schiff bases.27–31 Additionally, the larger macrocyclic rings have the potential for applications available to polyimines, such as their use as CANs to achieve thermoset-like properties, recycling and self-healing.20–22 Thermoset-like properties, which are obtained from the highly crosslinked polymer backbone, help to increase rigidity and mechanical strength within the materials.22 Due to the strong covalent bonds present, framework rigidity is maintained even at high temperatures, giving various engineering applications. Self-healing and recycling properties are also advantageous when considering the global requirement to reduce the quantity of waste generated, which is beneficial from both an environmental and economic viewpoint. Self-healing polymers, which utilise the framework structure to ‘heal’ the effects of degradation over time, have been used in biomedicine, electronics and coatings, extending the lifetime of the material.32 The free primary ketimines produced when the cyano carbon bonds to another nitrogen atom provides the possibility for post-functionalisation, further broadening the range of potential applications of this diverse material.
Alongside the experimental synthesis, ‘artificial synthesis’, growing a series of 100 cluster models of TCNQ-CTF, is used to probe potential structural features generated in this kinetically controlled reaction. To help interpret the resultant clusters, we have developed an efficient categorisation protocol, allowing the prevalence of different structural features to be automatically identified. Using these statistical data and a series of structural fragments derived from the clusters, we have efficiently rationalised the experimentally observed IR, NMR, and UV–Vis spectra of the resultant material. Importantly, the spectra can only be rationalised assuming a diverse set of structural features form experimentally, highlighting the importance of allowing for deviation from the thermodynamic product when modelling these exciting materials.
The structural diversity in the synthesised TCNQ-CTF thus provides additional chemical functionality, leading to possible applications across a broad range of areas, including within coordination chemistry, and as covalent adaptable networks.
The code for Ambuild can be found at https://github.com/linucks/ambuild.
| This journal is © The Royal Society of Chemistry 2025 |