Linwei
He‡
a,
Shengtang
Liu‡
a,
Long
Chen‡
a,
Xing
Dai
a,
Jie
Li
a,
Mingxing
Zhang
ab,
Fuyin
Ma
a,
Chao
Zhang
c,
Zaixing
Yang
*a,
Ruhong
Zhou
a,
Zhifang
Chai
a and
Shuao
Wang
*a
aState Key Laboratory of Radiation Medicine and Protection, School of Radiation Medicine and Protection, Collaborative Innovation Center of Radiological Medicine of Jiangsu Higher Education Institutions, Soochow University, Suzhou 215123, China. E-mail: shuawang@suda.edu.cn; zxyang@suda.edu.cn
bShanghai Institute of Applied Physics, Chinese Academy of Sciences, No. 2019 Jialuo Rd., Jiading Dist., Shanghai, 201800, China
cSchool of Materials Science and Engineering, Anhui University of Science and Technology, Huainan 232001, China
First published on 19th February 2019
99Tc is one of the most problematic fission products in the nuclear fuel cycle owing to its large inventory in used nuclear fuel, long half-life, potential radiation hazard, high environmental mobility of its major species 99TcO4−, and its redox-active nature. Ideally, 99TcO4− should be removed at the first stage, when the used fuel rods are dissolved in highly concentrated nitric acid solution, which can substantially reduce its interference with the solvent extraction process through catalytic redox reactions with the key actinides and diminish the chance of discharge into the environment as the volatile species during the waste vitrification process. However, this task cannot be achieved by any of the reported anion-scavenging materials including traditional polymeric anion-exchange resins, inorganic cationic framework materials, and recently developed cationic metal–organic framework materials, because they either are not stable under the extreme conditions of the combined high acidity and strong radiation field or do not possess the required uptake selectivity towards 99TcO4− in the presence of a huge excess of competing anions such as NO3− and SO42−. Herein, we present the first study of 99TcO4− removal under extreme conditions by a two-dimensional conjugated cationic covalent organic framework material, SCU-COF-1. This material exhibits ultrahigh acid stability, great resistance towards both large-dose β and γ irradiation and unprecedented 99TcO4− uptake capabilities including extremely fast sorption kinetics (sorption equilibrium can be reached within 1 min), ultrahigh uptake capacity (702.4 mg g−1 for the surrogate ReO4− at a slightly elevated temperature), and good anion-exchange selectivity towards 99TcO4−. These excellent features endow SCU-COF-1 with the practical capabilities of separating 99TcO4− from both simulant highly acidic fuel reprocessing solutions (3 M nitric acid) and low-activity waste streams at the US legacy nuclear site. The anion-exchange mechanism and the 99TcO4− uptake selectivity are further demonstrated and clearly visualized by the molecular dynamics simulation investigations.
99Tc is predominately present in the nuclear fuel cycle and in aqueous environments as the pertechnetate anion (99TcO4−). Unlike its congener MnO4−, 99TcO4− is not a strong oxidant, giving rise to its high stability under a variety of environmental conditions. The high symmetry (Td) and low charge density of 99TcO4− lead to its non-complexing nature and extremely high solubility (11.3 M at 20 °C).9 Therefore, 99TcO4− can facilely migrate in the natural water system and its transportation velocity is not affected by the majority of natural materials. Currently, vitrification followed by geological disposal has been established to be the only acceptable strategy for treating the high levels of nuclear waste in several countries.10 However, this strategy is not applicable to 99Tc because it would easily escape from the vitrification glass as the volatile 99Tc2O7.11 This issue along with improper waste management further results in the identified contamination of 99Tc in the groundwater at a variety of legacy nuclear sites, with 99Tc being the major radioactivity contributor in the vadose zone.12 Another significant issue for 99Tc is its catalytic redox capability with uranium, neptunium, and plutonium during used fuel reprocessing (e.g. the plutonium–uranium redox extraction (PUREX) process), hindering the precise valence state control of these key actinides, even when 99Tc is present in very low concentrations.13 Therefore, it would be highly desirable for 99TcO4− to be first separated when used nuclear fuel rods are dissolved in concentrated nitric acid prior to the PUREX process. This would be beneficial not only for the effective extraction of actinides but also in the elimination of 99Tc discharge into the environment during the following waste disposal processes.
Up to now, there have only been a handful of materials, including inorganic cationic framework materials,14 polymeric networks and anion-exchange resins,15 as well as cationic metal–organic frameworks (MOFs),16 reported to exhibit TcO4− separation capability from nuclear waste solutions. Among these, crystalline inorganic cationic framework materials, such as layered double hydroxides (LDHs),14h,i Y2(OH)5Cl,14f and NDTB-114a suffer from the clear demerits of poor selectivity, low capacity and slow uptake kinetics. Especially, poor selectivity makes these materials inapplicable because, either in nuclear waste solutions or natural water systems, coexisting anions such as NO3−, SO42−, CO32−, and PO43− that can induce anion-exchange competition with TcO4− are all present in 100- to 10000-fold of excess. Commercially available polymeric anion-exchange resins show a notable advantage in 99TcO4− uptake selectivity, but their poor radiation resistance impedes their practical applications in the radiological field especially in used fuel repossessing, where strong radiation fields (high fluxes of α, γ, β, and neutron irradiations) are present at extremely high dose rates. It was reported that the 99TcO4− uptake capacity would gradually decrease as the receiving radiation doses increase.15f,17 This critical issue is further amplified given that the 99TcO4− uptake kinetics by these materials are rather slow, leading to a significantly higher dosage received during the anion-exchange process. The recently developed systems based on cationic MOFs do offer solutions to the aforementioned drawbacks and find extensive application potentials in the remediation of 99TcO4− from contaminated water systems under relatively neutral pHs.16 However, none of these cationic MOFs can survive the highly acidic conditions (e.g. 3 M nitric acid) that are required for used fuel reprocessing.
Covalent organic frameworks (COFs) are a class of emerging crystalline polymers with periodic structures and well-defined nanopores.18 Importantly, all these structural features are predictable and highly tunable and can be tailored for specific applications, similar to the MOF system. Different from the common COFs with neutral skeletons, ionic COFs (iCOFs) are a subgroup of COFs with charged linkages in the framework or attachments in the side chain, which would promote particular properties such as selective ion-exchange, compared with neutral COFs.19 We initially proposed that the system of iCOFs may overcome the challenge for TcO4− separation under the combined extreme conditions of high acidity, strong radiation field, and high ionic strength for several reasons. First, compared with the majority of MOF materials, COFs exhibit clear elevations in the stabilities under different types of extreme conditions such as high acidity, alkalinity, and salinity, owing to the increased robustness of the covalent bond-based linkages for building the structure.18d,e,g Second, many COFs are highly conjugated and can substantially stabilize the radical-based intermediates generated by radiation, a merit not possessed by traditional polymeric anion-exchange resins. This may lead to excellent radiation resistance for these COF materials, one of the most important prerequisites for nuclear applications. Third, the porous nature and ordered channels may aid in the rapid transportation of 99TcO4− in the COF structure, which is advantageous for accelerating the anion-exchange process, in comparison to amorphous polymeric materials. Finally, the local hydrophobic environment in the COFs is expected to enhance the uptake selectivity for relatively hydrophobic 99TcO4−, compared to coexisting anions with higher charge densities, a design rule that has already been applied in selective polymeric anion-exchange resins.15a Recently, two 3D cationic COF materials 3D-ICOF-1 and 3D-ICOF-2 were tested for the removal of MnO4− under neutral conditions by Qiu et al.,19f which initially suggests the promise of 99TcO4− remediation using these types of materials, although the capability may be overestimated owing to the partial chemical reduction of MnO4− during the sorption experiment. Despite that, to the best of our knowledge, there has been no study on the separation and remediation of 99TcO4− by COFs through the direct radiological test.
iCOFs are in general very challenging to build because the charge repulsion between the adjacent building blocks impedes their ordered connection, often leading to poor crystallinity or an amorphous structure.19g Herein, we provide a new strategy to constitute cationic COFs with good crystallinity based on positively charged linkers and neutral nodes via irreversible tautomerism. The synthesized two-dimensional viologen-based iCOF (denominated as SCU-COF-1) displays an unusual ABC-staggered stacking model and unprecedented 99TcO4− uptake capability including large uptake capability, fast sorption kinetics, high anion-exchange selectivity towards 99TcO4−, great acid stability, and decent radiation resistance. These excellent features endow SCU-COF-1 with the promising ability to remove 99TcO4− from 3 M nitric acid solution, a key step for potentially establishing a new process aiming at solving the “technetium issue” in the nuclear fuel cycle.
The powder X-ray diffraction (PXRD) pattern of SCU-COF-1 (Fig. 1e) identifies an ordered structure, where a strong low angle peak at 2θ = 4.02° is observed. Two other prominent peaks at 8.23° and 25.59° were also identified. The peak at 4.02° corresponds to the (100) plane of the structure. The broad feature at 8.23° is recognized as a superposition of peaks at 7.85° and 8.60°, which are attributed to the (110) and (200) planes, respectively. Additionally, the presence of another broad peak at 25.59°, which should be assigned to the (001) plane, suggests highly ordered π–π stacking existing between the adjacent layers and the good planarity of the 2D layers.21
To obtain the simulated crystal structure of SCU-COF-1, theoretical simulations were conducted using Materials Studio. Three stacking models were generated and optimized: the eclipsed stacking (AA-stacking) model, the AB-staggered stacking model, and the ABC-staggered stacking model (see ESI† for details). The simulated PXRD patterns corresponding to the three different models are shown in Fig. 1e for comparison. The experimental pattern matches well with the simulated pattern derived from the unique ABC-stacking model. The features of stacking are shown in the top and side views of SCU-COF-1 (Fig. 1d). In order to distinguish the three distinct layers more clearly, we marked out layer A, layer B, and layer C with different colours in a periodic unit cell, respectively (Fig. 1b). After DFT optimization and Pawley refinement with good agreement factors (Rp = 1.55%, Rwp = 1.98%), we verified the structure in the R3 space group with the unit cell parameters of a = b = 43.47 Å, c = 10.77 Å, α = β = 90°, γ = 120°. In contrast, the simulated patterns of the AA-stacking model and the AB-staggered stacking model in the P6 space group clearly deviate from the experimental pattern.
Overall, the spatial structure of SCU-COF-1 can be described as parts of two reacted ligands being connected sequentially into a flat cationic 2D layer with a large hexagonal channel (∼4.34 nm). These cationic layers are further stacked via π–π interaction in the ABC-staggered manner, which in fact minimizes the charge repulsion between the adjacent layers. A series of shrunken triangular channels (∼1.44 nm) along the c-axis can be still observed after the stacking. We further simulated the electrostatic potential (ESP) distribution of the cationic layer (Fig. 1c). The results show that the regions close to the N atoms in the pyridine rings provide the most positive ESP (blue areas in the hexatomic ring), which are appropriate for accommodating the charge-balancing Cl− anion, whose existence can be identified through EDS analysis (Fig. S3†). In addition, the positive charges are delocalized on the pyridine ring, which is expected to play a leading role in the anion-exchange process.
The porosity of SCU-COF-1 is probed by the nitrogen gas adsorption isotherm measured at 77 K (Fig. 2a, red dots and lines). The calculated Brunauer–Emmett–Teller (BET) surface area of 22.98 m2 g−1 indicates moderate porosity, which partially originates from the staggered stacking of layers that may hinder the efficient transportation of nitrogen in the structure. A similar phenomenon was observed in other 2D COF materials with the staggered stacking.19h,j,21 On the other hand, the presence of charge-balancing Cl− in the pore may further block the entrance of guest gas molecules. However, this relatively small surface area is not an issue for rapid anion-exchange as discussed below. In addition, the nonlocal density functional theory (NLDFT) analysis gives an average pore width of 1.37 nm, which is in good agreement with the simulated crystal structure of SCU-COF-1 (Fig. 2b).
To elevate the thermal stability of SCU-COF-1, thermogravimetric analysis (TGA) was measured from 25 °C to 900 °C in N2 atmosphere. The weight loss before 100 °C could be assigned to the loss of lattice water, and SCU-COF-1 was thermally stable up to 225 °C (Fig. S4†). The chemical stabilities of SCU-COF-1 were investigated by dispersing fresh samples in various solvents or acids including H2O, DMF, THF, aqueous HCl (1 M), aqueous HNO3 (1 M and 3 M), and aqueous NaOH (1 M) solutions for 48 h. In various solvents (H2O, DMF, THF) or 1 M HCl, the PXRD patterns were almost identical to that of the original sample, indicating that SCU-COF-1 can retain its crystalline structure in these solvents. Even in 1 M HNO3 and 3 M HNO3, the crystallinity of SCU-COF-1 slightly decreased but maintained the ordered structure (Fig. 2d). The unchanged characteristic vibrational peaks in the IR spectra (Fig. 2c) and a BET surface area of 18.25 m2 g−1 (Fig. 2a, blue dots and lines) for samples treated with 3 M HNO3 further underscored the acid stabilities of SCU-COF-1. Note that the slightly reduced surface area may partially originate from the larger ionic radii of NO3− exchanged into the structure compared with the original Cl−, which may further block the transportation pathway of gas molecules. SCU-COF-1 is unstable in highly basic solutions (1 M NaOH). This can be explained by the presence of the cationic pyridine hexatomic rings, which are easily attacked by the hydroxyl as the nucleophile in solution.22 More impressively, the highly conjugated structure of SCU-COF-1 protects the functional cationic pyridinium rings from high doses (200 kGy, 400 kGy, and 600 kGy) of β-irradiation provided by an electron accelerator (1.5 MeV, dose rate of 2250 kGy h−1) or γ-irradiation from a 60Co source (92.42 PBq, dose rate of 3.125 kGy h−1). The PXRD patterns of the irradiated samples confirm that SCU-COF-1 fully retains its crystallinity after irradiation (Fig. 2e and S5†). In addition, the nearly unchanged IR spectra (Fig. 2c) and BET surface areas after being exposed to β- (20.01 m2 g−1) or γ- (19.00 m2 g−1) irradiations (Fig. 2a, green and purple dots) of SCU-COF-1 synergistically demonstrate its excellent radiation resistance. Moreover, the radiation resistance can be further confirmed by the anion-exchange capacity experiment using the irradiated sample (see discussion below). In the Fourier-transform infrared (FT-IR) spectrum, the strong peaks at 1580 and 1263 cm−1 originate from the stretching of CC and C–N, respectively, confirming the existence of the cis-keto form (Fig. 2c). Furthermore, the 13C nuclear magnetic resonance (NMR) spectrum exhibits four characteristic peaks. The signal at ∼185 ppm, which is assigned as the C atoms of ketone, demonstrates the formation of the cis-keto structure of SCU-COF-1 (Fig. 2f).20b
To investigate the sorption kinetics of SCU-COF-1 towards 99TcO4− or ReO4− as the non-radioactive surrogate of 99TcO4− owing to their almost identical charge densities and therefore similar anion-exchange properties,4 20 mg of solid materials was added into 20 mL of aqueous solution containing 28 ppm ReO4− or 14 ppm 99TcO4− at 300 K. As shown in Fig. 3a, the concentrations of 99TcO4− as a function of contact time were measured by its characteristic peak at 290 nm in the UV-vis spectra. In addition, the kinetics curve can be also confirmed by measuring the time-dependent radioactivity of the solution through a liquid scintillation counting (LSC) technique (Fig. 3b). The results show that SCU-COF-1 can almost quantitatively remove 99TcO4− from solution within an extremely short of time of ∼30 s, which is, in fact, the shortest time required for experimental operation. Therefore, it is difficult to establish a sorption kinetics model, as the first data point has already reached the sorption equilibrium. The results of ReO4− removal kinetics experiment are identical with those of 99TcO4− (Fig. 3c). This ultrafast anion-exchange kinetics sets SCU-COF-1 as the record-holding material along with the recently reported SCU-CPN-1 and is significantly faster than those of other types of anion-exchange materials (Table S1†). For inorganic cationic framework materials and commercial polymeric anion-exchange resins, the anion-exchange equilibrium would require at least hours (up to days) to establish. For example, the first purely inorganic 3D cationic framework NDTB-1 can remove approximately 72% of ReO4−/99TcO4− in more than 36 h.14b The anion-exchange resins Purolite A532E and A530E resins15a require at least 2.5 h to reach equilibrium. Only 32% and 38% of ReO4−, respectively, can be removed from the solution in a short time (5 min).15c Recently developed cationic MOFs show progress in this regards, but still require at least 10 min to quantitatively remove 99TcO4− from solution.16b We propose that this excellent feature is a result of ordered channels in the structure of SCU-COF-1 that allows for the efficient transportation of guest anions and the hydrophobic pore that leads to the fast exchange-out of relatively hydrophilic Cl− back into the aqueous solution (more profound discussions are provided in the section of the sorption mechanism). The advantage in anion-exchange kinetics is the benefit of the rapid remediation of anionic radionuclide contamination during an emergency event as well as reducing the received radiation dosage when dealing with highly radioactive systems such as used fuel reprocessing to minimize the radiation-induced degradation of materials. Moreover, the distribution coefficient (Kd) of SCU-COF-1 calculated from the ReO4− sorption kinetics experiment is 3.89 × 105 mL g−1, which is higher than most reported materials (Table S2†), exhibiting a promising removal depth for the remediation task. FT-IR (Fig. 2c) and EDS mapping measurements (Fig. S3†) also confirm the successful substitution of Cl− by ReO4− in SCU-COF-1. The sorption isotherm experiment was carried out at 300 K to gain insight into the anion-exchange capacity of SCU-COF-1 for trapping ReO4−/99TcO4−. From the curve, the saturated sorption capacity calculated based on the Langmuir model is 367.6 mg ReO4− per g sorbents (Fig. 3d and Table S3†). This value is notably lower than the theoretic sorption capacity (872.7 mg g−1) calculated by assuming Cl− anions are completely substituted, suggesting the presence of a kinetic barrier for the complete anion-exchange, which can be confirmed by the theory study discussed below. We therefore checked the sorption capacity at an elevated temperature of 373 K. Upon heating and refluxing for 24 h, the maximum ReO4− uptake capacity reached 702.4 mg g−1, as the actual experimental value, being quite close to the theoretic value. This sets SCU-COF-1 as a significant anion-exchange material with one of the highest reported exchange capacities. In comparison with inorganic materials such as Yb3O(OH)6Cl (48.6 mg g−1),14d,16c NDTB-1 (49.4 mg g−1)14a,16c and LDH (130.2 mg g−1),14h,15a polymeric anion-exchange resins such as Purolite A532E (446 mg g−1),16a,c and cationic MOFs including UiO-66-NH3+ (159 mg g−1),16e SCU-100 (541 mg g−1),16c and SCU-101 (247 mg g−1),16b SCU-COF-1 exhibits a notably larger sorption capacity towards ReO4−.
Since anions such as NO3− and SO42− often coexist in huge excess with 99TcO4− both in highly radioactive waste solutions (e.g. nitric acid solution) and under environmentally relevant conditions (such as in groundwater), the anion-exchange selectivity of sorbents is critically required for practical applications. Because those anions with higher charge density often exhibit stronger electrostatic interactions with the cationic host, they are often captured with priority by typical anion-exchange materials, which is known as Hofmeister bias selectivity, unless other types of driving forces are present.23 SO42−, CO32− and PO43− (one equivalent) were therefore chosen to evaluate the selectivity for ReO4− removal (Fig. 3e). In competing with different anions, SCU-COF-1 exhibits excellent uptake selectivity, and the removal percentages remain in the range of 85% to 99%. Even in the presence of a high concentration of NO3− (NO3−:ReO4− = 100:1, molar ratio), the removal percentage of ReO4− still reached 60% (Fig. 3f). The high selectivity for ReO4− is likely owing to the hydrophobic nature of the COF skeleton, which provides a remarkable affinity for relatively hydrophobic anions with lower charge densities.15c Some commercial resins, for example SuperLig-639,24 with a high uptake capacity towards 99TcO4−/ReO4−, have obvious drawbacks including low selectivity in the presence of a large excess of NO3− anions.3 More impressively, experiments for removing radioactive 99TcO4− were conducted in a simulated Hanford low-activity waste (LAW) melter recycle stream (Table S4†), which contains far higher amounts of NO3−, NO2− and Cl− (at least 300-fold in excess for each anion) than 99TcO4−. The results show that SCU-COF-1 can remove 20.9% and 62.8% of 99TcO4− at a solid/liquid ratio of 1 and 5, respectively (Table S5†).
To probe the uptake capability of SCU-COF-1 in acidic and alkaline conditions, solutions containing 400 ppm ReO4− in various pH values were used in the experiments. The uptake capacities of SCU-COF-1 were almost unchanged through a wide pH range from 4 to 10 (Fig. 4a). More significantly, adsorption experiments in extreme conditions were conducted to estimate its potential application in highly acidic conditions that is relevant to used fuel reprocessing. We further prepared a 1 M HNO3 solution containing ReO4− (c0 = 28 ppm or c0 = 338 ppm) and a 3 M HNO3 solution containing ReO4− (c0 = 253 ppm, NO3−:ReO4− = 2922, molar ratio). The latter system simulates the real situation of used fuel reprocessing. Impressively SCU-COF-1 exhibits superior extraction performance, as approximately 40% of ReO4− are removed at a solid–liquid ratio of 1:20 in the solution containing 338 ppm ReO4− (Fig. 4b). Under the conditions of 3 M HNO3, SCU-COF-1 can still remove 52.3% of ReO4− at a solid–liquid ratio of 1:40 in the solution (Fig. 4c). These significant results indicate the SCU-COF-1 are indeed feasible for 99Tc separation at the first stage of used fuel reprocessing prior to the PUREX process.
After treatment in different doses of β- or γ-irradiation, the uptake for ReO4− by SCU-COF-1 is uninfluenced. The removal capacity at room temperature for ReO4− was completely unchanged after irradiation (Fig. 4d and e). In comparison, two commercial resins Purolite A532E and A530E have reductions of 9% and 18% of sorption capacities after exposure to 600 kGy β-irradiation (or γ-irradiation), respectively.15c Additionally, SCU-COF-1 is also a reversible anion-exchange material for multiple sorption/desorption trials. More than 95% of ReO4− could be desorbed back to solution after treatment with 3 M NaCl solution at 300 K (Table S6†). Even after 4 cycles of sorption/desorption, the sorption performance of SCU-COF-1 was not affected (Fig. 4f). More impressively, the COF can keep almost the same regeneration property after removing ReO4− in 3 M HNO3 (Fig. S6†), which can significantly reduce the cost for practical applications.
We further explored the adsorption dynamics of 99TcO4− anions into the structure of SCU-COF-1 and the underlying uptake mechanism using all-atom molecular dynamics (MD) simulations (Fig. 5a). Impressively, after placing the solid structure of SCU-COF-1 into the aqueous 99TcO4− solution, a large number of 99TcO4− anions initially diffused in aqueous solution were rapidly sorbed into the interior of SCU-COF-1. A video revealing the anion-exchange process is provided as a ESI file.† During this process, many Cl− anions originally dwelling in the pores of SCU-COF-1 were desorbed from the solid structure (Fig. 5b, also see the Movie S1 in ESI† for more intuitive illustration). To quantitatively describe this dynamics process, the adsorption rate of 99TcO4− anions from the external solution environment (Fig. 5c red curve) and the residence rate of original Cl− anions in SCU-COF-1 (Fig. 5c green curve) along with the simulation time were calculated. As shown in Fig. 5c, a very clear crossover occurs at ∼26 ns. After this time point, the amount of 99TcO4− anions entering SCU-COF-1 exceeds the amount of Cl− anions that still remain in SCU-COF-1. After ∼40 ns, a dynamic equilibrium is reached, i.e., the 99TcO4− adsorption rate is nearly stabilized at a constant value of (60.6 ± 2.3)%. Meanwhile, the reserving rate of Cl− anions is nearly unchanged at (34.1 ± 2.0)%. In addition, another control simulation in pure water demonstrates that in the absence of 99TcO4−, nearly (87.2 ± 1.4)% of Cl− anions could be retained in SCU-COF-1 throughout the full simulation period (Fig. 5c blue curve). This suggests that the sorption of 99TcO4− anions initiates and further complicates the following anion-exchange process.
To uncover the driving forces for the anion-exchange process, the time evolution of the nonbonded interaction energies of COF-99TcO4−, COF-Cl−, and 99TcO4−-Cl− are calculated and analyzed, and they were further deconvoluted into the contributions from the electrostatic and van der Waals (vdW) interactions as well (Fig. 5d). During the main 99TcO4− uptake process (<40 ns; after 40 ns, simulations are converged), the nonbonded interaction energy between 99TcO4− anions and SCU-COF-1 sharply decreases by ∼13000 kJ mol−1 (decreases from 0 to −13000 kJ mol−1; Fig. 5d solid red square) (i.e., more favorable). The contribution from the electrostatic part is ∼10000 kJ mol−1 (Fig. 5d hollow red inverted triangle), which is about 3.3-fold of that from the vdW part (∼3000 kJ mol−1) (Fig. 5d hollow red triangle). During the same period, many Cl− anions were desorbed from SCU-COF-1. The nonbonded interaction energy between Cl− anions and SCU-COF-1 increases by ∼8000 kJ mol−1 (increases from −13000 to −5000 kJ mol−1; Fig. 5d solid green square) (less favorable), of which the electrostatic part contributes ∼7400 kJ mol−1 (Fig. 5d hollow green inverted triangle) and the vdW part contributes ∼600 kJ mol−1 (Fig. 5d hollow green triangle). Meanwhile, the nonbonded interaction energy between the Cl− anions and 99TcO4− anions increased by ∼850 kJ mol−1 (increases from 0 to 850 kJ mol−1; Fig. 5d solid blue square) (less favorable), with the electrostatic part increasing by 900 kJ mol−1 (repulsive) (Fig. 5d hollow blue inverted triangle), whereas the vdW part decreases by ∼50 kJ mol−1 (attractive) (Fig. 5d hollow blue triangle). These results suggest that the strong direct nonbonded interaction (particularly the electrostatic interactions) between 99TcO4− anions and SCU-COF-1 intensively drives this vigorous 99TcO4− uptake process. During this process, the 99TcO4− anions gained considerable energy, which can overcompensate for the energetic penalty induced by the desorption of Cl− anions. Moreover, the direct electrostatic repulsive interaction between 99TcO4− anions and Cl− anions further promotes the desorption of Cl− anions. To gain a more microscopic insight into the substituting procedure of Cl− anions by 99TcO4− anions, we further compared the binding free energy of a single 99TcO4− anion to the binding site that is most energetic favorable for the Cl− anion. This binding site is determined by the analysis of the Cl− anions residing probability distribution in the control simulation (for a better sampling). It was demonstrated that Cl− ions prefer to be located near the pyridine ring on the inner wall of SCU-COF-1 (Fig. 6a), where the most positive charges are distributed (Fig. 1f). The potential of mean force (PMF) data demonstrates that the binding free energies are −9.9 kJ mol−1 and −13.7 kJ mol−1 for the Cl− anion and 99TcO4− anion, respectively, indicating that the 99TcO4− anion can successfully outcompete the Cl− anion in binding behavior (Fig. 6b). In other words, when one Cl− anion (at it most energetic favorable binding site) is exchanged with one 99TcO4− anion, it can gain 3.8 kJ mol−1 of energy total, which can be attributed to the preferential binding of 99TcO4− anions to SCU-COF-1 over Cl− anions. Therefore, an entire picture for the anion-exchange process can be revealed: on one hand, the much stronger binding free energy between the 99TcO4− anion and SCU-COF-1 leads to successful substitution of Cl− by taking over the anion binding site in SCU-COF-1 as time progresses; on the other hand, the occupation of the 99TcO4− anion significantly weakens the direct interactions between the SCU-COF-1 and Cl− anion and gradually excludes Cl− anions out of the SCU-COF-1 due to the direct electrostatic repulsive interactions between the 99TcO4− anion and Cl− anion.
All MD simulations were performed with the GROMACS software package30 (version 4.6.7). VMD31 software (version 1.9.3) was applied for trajectory visualization and analysis. The SPC/E water model32 was used for water. In the production run, a time step of 2.0 fs was used, and the data were collected every 10 ps. For each system, three independent trajectories of 50 ns were generated, which resulted in a total aggregate simulation time of 300 ns. All production runs were under isothermal-isochoric (NVT) ensemble; the temperature was maintained at 300 K using a υ-rescale thermostat.33 The particle-mesh Ewald (PME) method34 was used for long-range electrostatic interactions. A cutoff distance of 1.2 nm was used for the calculation of van der Waals (vdW) interactions. All solute bonds between heavy atoms and hydrogens were constrained with the LINCS algorithm.35 Periodic boundary conditions were applied in all directions, and the COF were frozen throughout the simulation process. We used umbrella-sampling simulations to investigate the potential of mean force (PMF) of a single 99TcO4− and Cl−. The reaction path for the PMF calculation was defined by pulling center of mass of the anion (99TcO4− or Cl−) from the most favourable binding site of Cl− along the upper surface of COF membrane to the bulk water (from ∼0.0 nm to ∼1.2 nm of the COF upper surface) with a harmonic force of 2000 kJ mol−1 nm−2. The reaction coordinate was divided into 12 windows, with each window simulated for 3 ns and thus accumulating a total simulation of 36 ns. The weighted histogram analysis method36 (WHAM) was applied to calculate the free energy. All PMF profiles were normalized to zero at a distance of 1.2 nm.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c9sc00172g |
‡ These authors contribute equally. |
This journal is © The Royal Society of Chemistry 2019 |