Alberto
de Juan
a,
Alejandro
López-Moreno
a,
Joaquín
Calbo
b,
Enrique
Ortí
*b and
Emilio M.
Pérez
*a
aIMDEA Nanociencia, C/Faraday 9, Ciudad Universitaria de Cantoblanco, 28049, Madrid, Spain. E-mail: emilio.perez@imdea.org
bInstituto de Ciencia Molecular, Universidad de Valencia, 46980 Paterna, Spain. E-mail: enrique.orti@uv.es
First published on 7th September 2015
Single-walled carbon nanotubes (SWNTs) are one of the most promising nanomaterials and their supramolecular chemistry has attracted a lot of attention. However, despite well over a decade of research, there is no standard method for the quantification of their noncovalent chemistry in solution/suspension. Here, we describe a simple procedure for the determination of association constants (Ka) between soluble molecules and insoluble and heterogeneous carbon nanotube samples. To test the scope of the method, we report binding constants between five different hosts and two types of SWNTs in four solvents. We have determined numeric values of Ka in the range of 1–104 M−1. Solvent effects as well as structural changes in both the host and guest result in noticeable changes of Ka. The results obtained experimentally were validated through state-of-the-art DFT calculations. The generalization of quantitative and comparable association constants data should significantly help advance the supramolecular chemistry of carbon nanotubes.
Here, we describe a simple method for the determination of association constants between insoluble and heterogeneous nanotube samples and soluble molecules. To prove its validity, we have determined the association constants of five molecules towards two types of SWNTs in four different solvents.
The experimental procedure is described graphically in Scheme 1 and can be summarized as follows: SWNTs (1 mg mL−1, unless stated otherwise) are suspended in a solution of known concentration of the host molecule in a given solvent, and the mixture stirred for 2 hours to allow it to reach equilibrium. After this time, the suspension is filtered through a 0.2 μm-pore polytetrafluoroethylene membrane, retaining the host–SWNT complex. The solid is analysed through TGA (N2, 50 °C min−1) to quantify the amount of host in the complex, from which the concentration of free species is calculated by subtraction. Specifically, we measure the weight loss up to 600 °C, where all of the associated host has been desorbed and the nanotubes are still intact. From the degree of functionalization and the mass of the sample analysed, we calculate the total mass of host in the complex, from which its initial concentration in the equilibrium is immediate. Alternatively, the concentration of free species can be directly measured in the filtrate.‡ The same procedure is repeated for several initial concentrations of the host molecule, ranging from 0 to near saturation in the solvent under study. A blank experiment to determine the adsorbed/encapsulated solvent was run in all cases, and the data subtracted. All the experiments were performed at room temperature.
Scheme 1 Procedure for the measurement of [H]bound and [H]free. A known concentration of host molecule H and SWNTs are allowed to reach equilibrium, and then complexed and free species are physically separated through filtration. The concentration of [H]bound is measured by TGA (typical results for a titration experiment are shown, see the ESI† for full data set). The concentration of [H]free can then be calculated by subtraction or directly measured in the filtrate. |
The binding isotherms are obtained by plotting the degree of functionalization against the concentration of free host, and were analyzed using a standard 1:1 isotherm:23
Pyrene is by far the most widely used supramolecular partner for SWNTs, so we based our investigations on pyrene and its derivatives. Fig. 1 shows the chemical structure of the hosts for SWNTs used in the present work. First, we titrated 1 against plasma-purified SWNTs (pp-SWNTs, 98% purity, 0.8–1.6 nm in diameter) in tetrahydrofuran (THF), dimethylformamide (DMF), tetrachloroethane (TCE), and methanol (MeOH) at room temperature. Fig. 2 shows results of these titrations, where each data point is the average of three separate experiments.
As a first test of the experimental validity of our approach, we decided to get data for titrations with significant variations in the concentration of SWNTs. In particular, we used 0.1, 1 and 10 mg mL−1 of nanotubes in THF (Fig. 2a), which afforded Ka = 16.4 ± 0.8 M−1,§ 24 ± 6 M−1, and 21 ± 4 M−1, respectively. We were pleased to find that all values for Ka are identical within experimental error. The main variability comes from the degree of functionalization at saturation, which is significantly larger for the more dilute sample. This reflects a more efficient disaggregation of the nanotubes, which in turn results in an increase in the availability of binding sites for 1. Therefore, the method works correctly for samples with significantly different degrees of aggregation of the SWNTs.
With regards to the effect of the solvent, the association constants increase with decreasing ability to solvate SWNTs, showing that solvophobic interactions play a relevant role in the binding event. In DMF and TCE, solvents commonly used to disperse SWNTs, the binding constants are very small: Ka = 9 ± 3 and 4.5 ± 0.9 M−1, respectively (Fig. 2b and c). In THF there is an increase of one order of magnitude, to Ka = 24 ± 6 M−1 (Fig. 2a), which is further amplified in MeOH, a notoriously bad solvent for SWNTs, to reach millimolar affinity with Ka = (2.6 ± 0.2) × 103 M−1 (Fig. 2d).
In order to investigate whether the method is sufficiently sensitive to detect small changes in the structure of the nanotubes, we carried out titrations of 1vs. (6,5)-enriched SWNTs (93% purity, 0.7–0.9 nm in diameter) in THF, DMF, TCE and MeOH at room temperature. The association constant towards (6,5)-SWNTs are: Ka = 41 ± 8 M−1 in THF, 1.6 ± 0.4 M−1 in DMF, 1.6 ± 0.1 M−1 in TCE, and (1.0 ± 0.1) × 103 M−1 in MeOH (Fig. 3). Therefore, with the only exception of THF, in which some unexpected solvent effect takes place, the association constants are smaller than those towards pp-SWNTs. Considering the planar geometry of pyrene, it is expected to establish stronger van der Waals interactions with nanotubes of larger diameter, a tendency that is corroborated by DFT calculations (see below). These results confirm that the method is sensitive enough to such subtle differences in the structure of the nanotube as a decrease in the average diameter of the sample.
The method is also sensitive towards the structure of the host. To get experimental evidence, we designed a collection of hosts composed by 1,6-diaminopyrene (2), the benzoic and isophthalic esters of pyrene-1-methanol (3 and 4, respectively), and bis-pyrene U-shape molecule 5, which we have used in the synthesis of mechanically interlocked derivatives of SWNTs.25–28 Hosts 3–5 were titrated vs. pp-SWNTs in THF at room temperature, while we used DMF for the titration of 2 for solubility reasons.
Both electron-rich conjugated compounds and amines are known to interact strongly with SWNTs, so we expected 2 to show a significantly larger association constant compared to pyrene. This is indeed the case, as we calculated Ka = (2.2 ± 0.5) × 102 M−1 for the 2·pp-SWNTs associate (Fig. 4a), which is more than two orders of magnitude larger than the Ka of 1 in the same solvent. Addition of an extra aromatic ring in 3 also results in a significant increase in binding constant with respect to pyrene, reaching Ka = (9 ± 3) × 10 M−1 in THF (Fig. 4b). Bivalent tweezers-like hosts are a particularly popular design for the supramolecular association of SWNTs and fullerenes, as they typically show very good affinity at a relatively low synthetic cost.29 Indeed, 4 shows Ka = (6.5 ± 0.6) × 103 M−1 towards pp-SWNTs in THF (Fig. 4c). Finally, we decided to get an insight into the association of U-shaped molecule 5, which associates pp-SWNTs with Ka = (7 ± 2) × 103 M−1 (Fig. 4d), slightly larger than that of 4.
Note that 4 and 5 feature two pyrene binding motifs each, and might show multivalency and/or cooperativity phenomena.30 Since our method is based on measuring the concentration of the complex by desorbing it completely, it would not be valid to determine stepwise association constants. A possible approach to investigate such issues would be to utilize the Hill equation.31–33 Considering the clearly hyperbolic shape of the binding isotherms, we have determined average binding constants only.
To validate our experimental results, theoretical calculations were performed for the list of host–guest nanotube assemblies under the density functional theory (DFT) framework. The atom pair-wise Grimme's dispersion correction in its latest version (D3)34 was coupled to the hybrid density functional of Perdew–Burke–Hernzerhof (PBE0)35 through the Becke-Johnson damping function36 and including the three-body dispersion correction (EABC).37 The double-zeta Pople's 6-31G** basis set38 was employed throughout and the basis set superposition error (BSSE) was corrected according to the counterpoise (CP) scheme of Boys and Bernardi.39 The intensity of the interaction between host and guest was calculated by means of two different quantities. The interaction energy (Eint) is defined as the energy difference between the host–guest complex (HG) and the individual moieties separately (H and G), with all of them at the geometry of the complex:
Eint = EHGHG − EHGH − EHGG |
Otherwise, the binding (or association) energy (Ebind) was calculated taking into account the relaxation of the separate monomers and, therefore, considering the deformation energy required to transform the host/guest moieties from their minimum-energy geometries to the geometry acquired in the assembly:
Ebind = Edef + Eint |
Edef = (EHGH − EHH) + (EHGG − EGG) |
As a general model for the pp-SWNTs we have utilized a fragment of a zig-zag (10,0)-SWNT. The effect of the length of the nanotube into the intermolecular interaction was assessed by increasing the SWNT size in a 1·SWNT complex, showing that the association energy is nearly converged with sizes slightly larger than the host length (see Table S2 and Fig. S21†).
Fig. 5 displays the minimum-energy geometries for the 1–5 hosts assembled with the pp-SWNT model of C160H20 computed at the PBE0-D3/6-31G** level of theory in gas phase. Among the different closely energetic conformations of 1 over pp-SWNT, the diagonal arrangement is found to be the most stable, with close π–π contacts in the range of 3.2–3.5 Å. The interaction energy of 1·pp-SWNT is computed at −15.24 kcal mol−1, which is slightly reduced to −14.84 kcal mol−1 for the binding energy as a consequence of the deformation energy penalty (0.59 kcal mol−1). Moving from the pyrene system to 1,6-diaminopyrene (2), additional n–π interactions arise from close nitrogen⋯nanotube contacts (approximately at 4.0 Å). The Eint of 2·pp-SWNT is calculated 1.3 kcal mol−1 larger than for 1·pp-SWNT, but this difference is not maintained in the binding energy (Table 1). The deformation energy, calculated to be 2.83 kcal mol−1 for 2·pp-SWNT, explains this trend. The inclusion of an extra aromatic ring in 3 results in a significant increase of the interaction energy up to −23.68 kcal mol−1, with close π–π benzene⋯SWNTs (3.5 Å) and CO⋯SWNTs (3.2 Å) contacts. Bivalent tweezers-like hosts further improve the supramolecular affinity vs. pp-SWNT with Eint as large as −38.78 and −63.23 kcal mol−1 in 4 and 5, respectively. The binding energy in host 4 (−36.42 kcal mol−1) is indeed approximately the sum of Eint for its constituting moieties 1 and 3 (−14.84 + (−21.52) = −36.36 kcal mol−1), which supports the theoretical approach undertaken. Whereas the Edef of 4 is computed similar to 2 and 3, it amounts 20.46 kcal mol−1 for 5 due to the accommodation of the alkoxy chains around the nanotube (Fig. 4). This disposition confers 5·pp-SWNT an increased Ebind of −42.78 kcal mol−1 due to close CH⋯π contacts calculated in the range of 2.7–3.2 Å, which contribution to the total binding energy amounts 6 kcal mol−1.
Fig. 5 Minimum-energy geometry for the supramolecular assemblies formed by hosts 1–5vs. the pp-SWNTs model calculated at the PBE0-D3/6-31G** level of theory. |
System | E int | E def | E bind | CAa (Å2) |
---|---|---|---|---|
a The intermolecular contact area (CA) was calculated using the UCSF Chimera 1.7 software according to the formula: (area of the host + area of the guest − area of the complex)/2, where the areas used refer to solvent-excluded molecular surfaces, composed of probe contact, toroidal, and reentrant surface. | ||||
1·(6,5)-SWNT | −13.85 | 0.81 | −13.04 | 42.20 |
1·pp-SWNT | −15.24 | 0.59 | −14.84 | 42.70 |
2·pp-SWNT | −16.53 | 2.83 | −13.70 | 47.25 |
3·pp-SWNT | −23.68 | 2.16 | −21.52 | 75.30 |
4·pp-SWNT | −38.78 | 2.36 | −36.42 | 126.85 |
5·pp-SWNT | −63.23 | 20.46 | −42.78 | 188.55 |
Finally, the influence of the structure of the nanotube in the stability of the host–guest assembly was assessed by comparing the associates of pp-SWNT and (6,5)-SWNT with pyrene 1. The Eint of 1·(6,5)-SWNT was computed at −13.85 kcal mol−1, which is 1.4 kcal mol−1 smaller than the Eint of 1·pp-SWNT. The minimum energy structures calculated for the supramolecular complexes between pyrene and the two types of nanotubes (Fig. S22†) reveal subtle differences in terms of intermolecular contacts. The diameter of (6,5)-SWNT is computed at 7.5 Å, slightly smaller than for the pp-SWNT model (7.9 Å), which provokes a less efficient supramolecular assembly with pyrene. The deformation energy of 1·(6,5)-SWNT is computed somewhat larger than 1·pp-SWNT (Table 1), suggesting that the pyrene core is required to have a larger deformation to accommodate over the more-curved nanotube surface of (6,5)-SWNT. Moreover, the intermolecular contact area for 1·(6,5)-SWNT is calculated to be 0.5 Å2 smaller than in 1·pp-SWNT.
Most remarkably, the calculated Ebind energies and the experimentally determined Ka values show excellent quantitative agreement, despite the fact that desolvation and solvation energies are not included in our calculations. A plot of the lnKavs. −Ebind for molecules 1, 3, 4, and 5, towards pp-SWNTs in THF at room temperature, the largest set for which we have extracted comparable Ka data, is shown in Fig. 6. Fixing the intercept to 0, the data fit well (r2 = 0.984) to a straight line of slope 0.22 ± 0.01. Therefore, our analysis shows that the ΔGbind determined experimentally is proportional to the calculated Ebind.
To illustrate the scope and limitations of this methodology, we have tested five different hosts and two types of SWNTs in four different solvents for a total of 17 binding constant determinations. The data fit well to the binding isotherm in all cases, with a minimum r2 of 0.937 (Table S1†). The method is sensitive to solvent effects, as well as to small structural changes in both the SWNT and the host. The numeric values of Ka span over approximately four orders of magnitude, showing that the method is valid both for very small and large binding constants. Our data were validated by DFT calculations, which correctly reproduce the trends observed experimentally.
Although the main objective of the present work was to develop a standard method for the determination of binding constants towards carbon nanotubes, several interesting observations can be made with the present data set of Ka. Perhaps the most relevant conclusion is that our results back the utilization of a single unit of pyrene as a noncovalent anchor to SWNTs in polar protic solvents,40 but caution against assuming that it will “adsorb irreversibly” to the nanotubes in any organic solvent,41 as the association constants can be as low as 1 M−1. This is particularly relevant in cases where the pyrene-SWNT supramolecular construct will be subjected to further modifications after association. In this respect, using two pyrene units connected to form a tweezers-like receptor seems a valid alternative.
Taking into account the simplicity of the methodology described, we sincerely hope that the determination of association constants will become routine for anyone interested in the supramolecular chemistry of carbon nanotubes. The generalization of such quantitative data will undoubtedly produce a significant leap in our understanding of their noncovalent chemistry.
The techniques and methods described here should also be applicable to other insoluble nanomaterials, such as few-layer graphene. We are currently working towards the extension of this method to such nanomaterials.
Footnotes |
† Electronic supplementary information (ESI) available: Synthetic details and characterization, thorough description of the method to determine the concentration of the relevant species, full set of TGA data, isotherms not shown in the main text, table with all the association constants, computational details and theoretical calibration of the nanotube length. See DOI: 10.1039/c5sc02916c |
‡ Somewhat counterintuitively, we find the direct measurement of [H]free in the filtrate more problematic experimentally. We believe the main reasons behind this are variations in the volume of solvent during the filtration process and spectral overlap of the hosts with the carbonaceous impurities present in the filtrate. |
§ Considering that each data point consists of three separate experiments, we have reported errors directly as obtained from the fitting software. Based on our previous experience determining association constants, an experimental error within 20% can be expected. |
This journal is © The Royal Society of Chemistry 2015 |