Tao Zheng‡
a,
Mo Zhu‡a,
Muhammad Waqasa,
Ahmad Umaira,
Muhammad Zaheera,
Jinxian Yanga,
Xiaozheng Duan*b and
Lianwei Li*a
aDepartment of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China. E-mail: llw@mail.ustc.edu.cn
bState Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, Jilin 130022, China. E-mail: xzduan@ciac.ac.cn
First published on 19th November 2018
In this work, we analyze the catalytic mechanism of P4VP–RuII(bda) polyelectrolyte–metal complex (PMC) as a water oxidation catalyst and elucidate how the unique slow diffusion and multi-charge properties of the polyelectrolyte ligand dominate the catalytic process. Four poly(4-vinyl pyridine)–Ru(bda) (P4VP–Ru) PMCs with different chain lengths and controlled Ru loading amounts were prepared and used as catalysts for catalytic water oxidation. These catalysts present excellent catalytic performance with turnover numbers (TON) from ∼1200 to ∼1700 because of the good hydration properties. Surprisingly, the combined catalysis kinetics and kinetic isotope effect (KIE) studies for P4VP–Ru PMCs confirm the single-site water nucleophilic attack (WNA) mechanism in catalysis, rather than the interaction between two metal oxide units (I2M). A combination of dynamic light scattering characterization, zeta-potential measurement and molecular dynamics simulation reveals that the slow diffusion and multi-charge properties of the polyelectrolyte ligand are responsible for the observed mechanism difference between the P4VP–Ru PMC system and small-molecule multi-nuclear system, though the two systems actually own a similar structural feature (flexible linkages between Ru centers). Our experimental and simulation results highlight the fact that though the existence of flexible linkages between Ru centers could provide large conformation entropy for the occurrence of Ru-dimerization in small-molecule and neutral polymer systems, the entropy elasticity could not overcome the electrostatic interaction energy in the PMC system. Clearly, this work unambiguously clarified why both intra-chain and inter-chain Ru-dimerization (I2M) are prohibited for the PMC system from a perspective of macromolecular chemistry and physics.
Due to the strong electron-donating ability of the pyridine ring, P4VP shows great potential in coordination with various transition metal ions to form polymer–metal complexes. To date, a number of examples of P4VP–metal complexes have been reported to show high activities for catalyzing organic reactions.10–14 However, in most of the previous reports, the used P4VP–metal complexes were normally cross-linked and heterogeneous in nature, i.e., they were neither soluble nor ionized in water. Generally, the heterogeneous catalytic system could provide a convenience for easier separation and recycling of catalysts from the reaction system, but simultaneously sacrifices the opportunity for the deep understanding of the structure–property correlation at a molecular level. For heterogeneous catalysis, the thermal motion of polymer segments is greatly prohibited due to the “frozen” state of chain segments, which ultimately results in the prohibition of possible cooperative catalytic effect originating from neighboring centers. In addition, if the catalysis product is poorly soluble in the reaction medium, the heterogeneous catalyst might be encapsulated in the product and eventually lose its catalytic activity. Such a situation has indeed been previously reported in the study of hydrogenation of 4-nitrophenol catalyzed by a heterogeneous catalyst.15 It is worth pointing out that the P4VPs were typically used as polymer ligands in previous related studies,10–14 whereas how the chain length dependent properties affect the catalysis process has never been explored experimentally.
So far, the interaction between the polymer ligand and metal center remains unclear in most of previous studies, and how the polymer ligand effect influences the catalytic performance as well as the catalytic mechanism has been scarcely investigated. By a careful comparison of small-molecule ligands and polymer ligands, it is obvious that the greatest difference lies in the chain length-dependent properties, such as the unique slow-diffusion and multi-charge properties of polymer ligand.16 However, how these polymer effects influence the catalytic performance and mechanism for a given polymer–metal catalytic system has never been seriously considered.
In this work, we illustrate how to utilize the unique slow diffusion and multi-charge properties of polyelectrolyte ligand to regulate the catalytic performance and mechanism. Specifically, we report the preparation and catalysis study of four P4VP–RuII(bda) polyelectrolyte–metal complexes (PMCs) with different chain lengths as polymer catalysts for homogeneous water oxidation (Scheme 1), where bda is 2,2′-bipyridine-6,6′-dicarboxylate. The P4VP–Ru catalyst was designed based on the following rationales: (1) [Ru(bda)(L)2] (L: donor ligand) is one of the most widely studied water oxidation catalysts (WOCs) showing high activities;17–19 (2) the chain length of P4VP can be tailored via controlled polymerization to regulate the polymer coil size as well as the diffusion coefficient of P4VP–Ru complex in aqueous solution; (3) P4VP–Ru with different loading amounts of Ru can be easily prepared by the efficient ligand exchange reaction between pyridine unit of P4VP and precursor [Ru(bda)(pic)(DMSO)] (pic = 4-picoline); (4) the acidic aqueous solution (pH = 1.0) is typically used for Ru(bda) based catalytic studies, which endows protonated P4VP ligand with multi-charge property.
Scheme 1 Schematic synthesis of P4VP–Ru polymer–metal complex by a combination of RAFT polymerization and pyridine/DMSO exchange reaction. |
The above advantages facilitate the present study on the influence of slow diffusion and multi-charge effects of PMC catalysts on the catalytic mechanism and performance in water oxidation reaction. Our previous studies of polymer-based catalysts for heterogeneous water reduction and oxidation reactions lay the foundation of present work.20–22 This work has revealed, for the first time, that the slow diffusion and multi-charge effects of macromolecular ligand play important roles in determining the catalytic mechanism and performance for the PMC-catalyzed water oxidation reaction.
Fig. 1a shows the SEC characterization result of P4VPs in DMF. The result demonstrates that all the RAFT polymerizations were conducted in a controlled manner, reflected in the small polydispersity indexes (Mw/Mn) ranging from ∼1.05 to ∼1.16. The number average molar mass (Mn) determined by SEC is ∼5.20 × 103 g mol−1, ∼9.10 × 103 g mol−1, ∼1.35 × 104 g mol−1 and ∼3.68 × 104 g mol−1 for P4VP21, P4VP42, P4VP74 and P4VP198, respectively, which are significantly larger than those determined by 1H-NMR (Fig. S2†). The inconsistency is supposed to be attributed to the difference in hydrodynamic sizes of P4VP and polystyrene standard for a given molar mass. The chain length effect can be also observed from the UV-vis spectra shown in Fig. 2b, where the peak associated to tri-thiocarbonate (at ∼310 nm) gradually increases as Mn decreases, and the absorption peak located at ∼260 nm is assigned to absorption of pyridine rings.
Fig. 1 (a) SEC curves of polymer ligands P4VPs with different chain lengths in DMF at T = 35 °C. (b) UV-vis spectra of polymer ligands P4VPs with different chain lengths in MeOH. |
Highly efficient ligand exchange reaction was utilized for the preparation of P4VP–Ru polymer–metal complexes. [Ru(bda)(pic)(DMSO)] (pic = 4-picoline) precursor was prepared according to the reported procedures (Fig. S3†).23,24 The ligand exchange reaction was typically carried out at 45 °C for 12 h. The reaction solution was passed through a filter (450 nm) to remove the trace amount of insoluble substance before precipitation purification. Considering the limited solubility of P4VP–Ru complexes in MeOH and aqueous solution, the mass fraction (wt%) of Ru in P4VP–Ru complexes were controlled to be less than ∼5%. In the initial stage of exchange reaction, [Ru(bda)(pic)(DMSO)] precursor was insoluble in MeOH, while the reaction mixture became more and more clear as the reaction proceeded, indicating a much better solubility of P4VP–Ru complex because of the stronger solvation property of P4VP ligand. The successful preparation of P4VP–Ru complexes was first confirmed by FTIR measurement. As shown in Fig. 2b, P4VPs–Ru present the feature vibration bands of carboxylic ester (∼1650 cm−1, blue band), signifying the successful installation of bda groups onto P4VP polymer backbone; moreover, the skeleton stretching vibration bands of pyridine rings (∼1600 cm−1, red band) is clear for P4VPs–Ru. The ratio (I1/I2) of peak intensities of carboxylic ester (I1, blue band) and pyridine (I2, red band) in FTIR spectra can be used to semi-quantitatively estimate the exchange process (Fig. S4 and Table S1†). Namely, the I1/I2 ratios are typically 0.50–0.70 for P4VPs, but 0.90–0.95 for P4VPs–Ru. Fig. 2a further shows the comparative plots of UV-vis spectra of P4VPs and P4VPs–Ru at a fixed polymer concentration in MeOH. As shown, two newly emerging absorption bands located at 350–400 nm and 450–550 nm are obvious for P4VPs–Ru, which can be assigned to the absorption feature of RuII. The intensities of the two bands seem to be similar for P4VPs–Ru with different chain lengths, indicating similar Ru loading amounts. Meanwhile, P4VP–Ru shows a good stability of RuII oxidation state in methanal (Fig. S5a†). It is worth noting that the initial attempt by using HNO3 (aq, pH = 1.0) as a test solvent failed, as shown in Fig. S5b,† probably due to the oxidation of RuII into RuIII in acidic environment. After the addition of reducing agent ascorbic acid, the oxidized species RuIII can be reversibly reduced to be RuII (Fig. S5b†).
Further, a combination of 1H-NMR and inductively coupled plasma mass spectrometry (ICP-MS) was used to quantitatively analyze the loading amount of Ru in P4VPs–Ru. Compared with the bare P4VP, P4VP–Ru is supposed to own a higher mass fraction of aromatic protons due to the introduction of bda unit, which can be quantitatively reflected in the larger ratios (Aaro/Aali) of integral areas of aromatic protons (Aaro) and aliphatic protons (Aali) in 1H-NMR spectra (Fig. S6†). Overall, the Aaro/Aali ratio increases significantly after ligand exchange reaction (Fig. 3). The experimentally determined and theoretically calculated Aaro/Aali ratios for P4VP and P4VPs–Ru are summarized in Table S2.† Finally, ICP-MS measurement quantitatively determined the mass fractions (wt%) of Ru element of P4VPs–Ru (Fig. 3b and Table S3†). Overall, wt% varies from 2.8–3.5%, slightly smaller than the theoretical values (3.5–3.9%) calculated based on the feeding ratios, probably due to the incompletion of the ligand-exchange reaction. The redox potentials of P4VP42–Ru and P4VP198–Ru were determined by cyclic voltammetry (CV) in aqueous solution (pH = 1.0). The larger ΔE values is probably attributed to the slower electron transfer kinetics for polymer catalyst system. As shown in Fig. 4, when the scan was carried out from 0.2 to 1.2 V, the RuIII/II couple at ∼0.70 V can be clearly observed as well as the electro-catalytic waves associated with the RuIV/III (1.0 V). We noted that a small shift of ±0.02 V is observed for P4VP42–Ru and P4VP198–Ru, which is within the reasonable range. Moreover, the slight difference might also be related with the difference of degree of inter-chain coupling for different P4VP polymers. For comparison, the CV curve was also measured for reference complex [Ru(bda)(pic)2]. Clearly, the result shows that the active Ru sites were successfully attached onto P4VP backbones.
By using these well characterized P4VP–Ru complexes as catalysts, we further explored how the polyelectrolyte ligand effect influences the catalytic behavior. From green chemistry perspective, the aqueous solution without organic co-solvent was used as the reaction medium. The catalytic water oxidation was performed by using cerium(IV) ammonium nitrate (CAN) as sacrificial oxidant. The reaction kinetics was monitored with an oxygen-sensor (Unisense), and the evolved amount of O2 in gas space was analyzed by GC (Fig. S7†). In a typical experiment, 0.5 mL of P4VP–Ru solution (pH = 1.0) was degassed and then injected into O2-free HNO3 solution (9.5 mL, pH = 1.0, [CeIV] = 0.105 M) in a sealed vial with rubber septum under constant stirring rate. Each P4VP–Ru solution was prepared by diluting a mother solution with polymer concentration of 1.0 mg mL−1, and then pre-filtered through a 0.45 μm PTFE filter before injection.
Fig. 5a–d show the reaction time and catalyst concentration dependent oxygen evolution curves by using P4VP–Ru polyelectrolyte–metal complexes (PMCs) as catalysts. As expected, the oxygen evolution rate increases with the molar concentration of Ru (CRu), where CRu is the molar concentration of Ru in the final reaction solution with a total volume of 10 mL. Experimentally, the turnover frequency (TOF) was determined to be ∼6.85 s−1 (P4VP21–Ru, CRu = 14.00 μM), ∼5.79 s−1 (P4VP42–Ru, CRu = 13.86 μM), ∼11.94 s−1 (P4VP74–Ru, CRu = 14.00 μM) and ∼3.01 s−1 (P4VP198–Ru, CRu = 17.33 μM), respectively. According to Henry's law, the produced molecular oxygen in solution is supposed to successively undergo the nucleation, fusion and diffusion processes. Finally, the formed nano bubbles diffuse into gas phase to reach the gas–liquid phase equilibrium on thermodynamics. Thus, the interesting turning points from increase to decrease in Fig. 5 actually reflect the transition points where the O2 diffusion rates start to exceed the O2 evolution rates for P4VPs–Ru. In addition, the existence of turning points also indicates that the widely adopted pressure-sensor method by monitoring the internal pressure in gas phase might be not the best method for the real-time monitoring of the catalysis kinetics, due to the hysteresis of the oxygen diffusion process.
Fig. 5 (a)–(d) Reaction time (t) dependence of the oxygen concentration ([O2]) in situ measured in solution by using P4VPs–Ru as catalysts, where [CeIV] = 100 mM and pH = 1.0. |
To obtain reliable kinetics data, the initial oxygen evolution rates (vi) were calculated on the basis of data collected within the first few seconds (Fig. S8†), and the CRu values were calculated according to the ICP-MS results. As shown in Fig. 6, the quasi-linear relationships between vi and CRu are established in log-scale, i.e., vi ∼ CRuα. The fitted value for scaling exponent α is ∼1.00, ∼1.13, ∼1.00 and ∼1.12 for P4VP21, P4VP42, P4VP74 and P4VP198, respectively. Surprisingly, the constant α shows no dependence on the chain length, which is different from our initial anticipation. However, these values are much smaller than those for homogeneous catalyst [Ru(bda)(pic)2] (α ∼ 2.0) reported in literature.18,25 Theoretically, the exponent α is related with the catalytic mechanism (vide infra).
Fig. 6 Catalyst concentration (CRu) dependence of the initial oxygen evolution rate (vi) for P4VPs–Ru and small-molecule [Ru(bda)(pic)2], where [CeIV] = 100 mM and pH = 1.0. |
In addition, the turnover numbers (TONs) for P4VP–Ru PMCs were measured by GC. As shown in Fig. 7, TON is ∼1450 (P4VP21–Ru, CRu = 14.00 μM), ∼1620 (P4VP42–Ru, CRu = 13.86 μM), ∼1710 (P4VP74–Ru, CRu = 14.00 μM) and ∼1170 (P4VP198–Ru, CRu = 17.33 μM), respectively, which are significantly larger than ∼670 determined for [Ru(bda)(pic)2] in our previous work.25 The reason should be mainly attributed to the much stronger hydration ability of hydrophilic P4VP polymer ligand, which facilitates the diffusion and collision for water molecules around the catalytic centers of P4VPs–Ru. Numerically, P4VP198–Ru shows the smallest TON, which is due to its much poorer solubility in HNO3 aqueous solution. Overall, the above combined TOF and TON measurements clearly demonstrate that the polymer catalysts not only show high catalysis rates but also present excellent performances.
Generally, there are two widely accepted mechanisms for water oxidation at the metal centers, including: (1) single-site water nucleophilic attack (WNA) and (2) interaction between two metal oxide units (I2M).17,26 It is well-known that the Ru(bda) based homogeneous catalysts generally oxidize water through bimolecular radical coupling process, i.e., I2M mechanism, which typically corresponds to the second-order kinetics, i.e., vi ∼ CRuα (α = 2.0).18,25,27 In contrast, only a few studies reported the first-order kinetics for Ru(bda) based catalysis, which can be roughly divided into the following three specific cases: (1) the ligand rigidity induced mechanism change;18,19,27 (2) the restriction of motion of Ru catalytic centers located at interface of heterogeneous system;28–30 (3) the intramolecular dimerization of two neighboring RuVO species to form –O–O– bond.31
For example, Sun et al. reported that the replacement of bda with 1,10-phenanthroline-2,9-dicarboxylate (pda) could lead to the change of catalytic mechanism from I2M to WNA because of the flexibility difference of ligand structure.27 Würthner et al. reported the preparation and catalysis behavior of Ru(bda) based supra-molecular macro-cycles as water oxidation catalysts, similarly, the supra-molecular catalysts were found to oxidize water via WNA mechanism because of the rigid ligand structure of supra-molecular assemblies.18,31 Accordingly, the first-order reaction kinetics was observed in these studies. On the other hand, Ru(bda) based heterogeneous catalysts were also reported to catalyze water oxidation via WNA mechanism, which is generally explained by the restriction effect on motion of Ru catalytic centers located at interface of heterogeneous system. Namely, via kinetic isotope effect (KIE) study, Sun et al. found that the Ru(bda) supported electrodes could achieve electro-catalytic water oxidation via WNA mechanism.29 Our recent work on heterogeneous Ru(bda) based polymer catalyst also revealed the key role of motion of catalytic centers in determining the catalytic mechanism in heterogeneous system.25 In reality, the first-order kinetic could not only exist in mononuclear catalytic process, but also in bi-nuclear catalytic process in some cases. Namely, the recent work by Sun et al. has demonstrated that the Ru(bda) based tri-nuclear complex with flexible carbon–carbon single bond linkage, as homogeneous catalyst, could also present an apparent first-order kinetics, which was attributed to the occurrence of intramolecular radical coupling of two neighboring Ru centers, i.e. I2M pathway.31
The above discussion clearly clarifies that the determined reaction order is not a reliable criteria for the judgment of catalytic mechanism. Overall, for a given homogeneous catalysis system, if the first-order kinetics and the WNA mechanism are experimentally confirmed, the Ru(bda) based catalyst has to be composed of rigid ligands without torsional flexibility; in contrast, if the first-order kinetics and the I2M mechanism are experimentally confirmed, the Ru(bda) based catalyst has to consist of multiple Ru centers with flexible linkages between neighboring centers, favoring intramolecular radical coupling process.
The question left here is to confirm the catalytic mechanism for our observed first-order kinetics in P4VP–Ru PMC system, i.e., whether catalysis is related with the WNA or I2M process. Intuitively, we supposed that our P4VP–Ru system should obey the I2M mechanism. This is because: (1) P4VP–Ru PMC catalyst is a typical type of multi-Ru catalyst; (2) the P4VP backbone is made of rotatable/flexible carbon–carbon single bonds. These two features make P4VP–Ru PMC system very similar with the tri-nuclear catalyst system reported by Sun et al., where the backbone of tri-nuclear catalyst was also made of flexible carbon–carbon single bonds, and the catalytic mechanism was shown to be I2M mechanism.19
Thus, the kinetic isotope effect (KIE) study was further performed to provide clues to the catalytic mechanism. D2O was used as reference solvent to explore whether O–H bond cleavage was involved in the rate-determining step. As shown in Fig. 8, for P4VP21–Ru, P4VP42–Ru and P4VP74–Ru, the oxygen evolution rates in D2O are significantly slower than those measured in H2O. Quantitatively, the KIEH/D values are ∼4.84, ∼4.99 and ∼3.43, where the KIEH/D value is the ratio of the oxygen evolution rates in H2O and D2O. These values are much larger than the reported values for [Ru(bda)(pic)2] (KIEH/D ∼ 0.9).17–19 However, for high molar mass P4VP198–Ru, its solubility in D2O was much poorer than that in H2O due to the weaker hydration ability, and an apparent KIEH/D ∼ 23 was observed in the D2O solution. Theoretically, the KIEH/D value is >2.0 if the proton-coupled process is involved (WNA); while it is normally between 0.7–1.5 if the catalytic mechanism is through I2M.32 Clearly, the present result clarifies that both the inter-chain and intra-chain dimerization processes between RuVO species are prohibited in P4VP–Ru system, though the flexible linkages are supposed to provide enough conformation entropy (entropy elasticity) for the occurrence of dimerization process. The mononuclear WNA is responsible for the observed first-order reaction kinetics.
The combined kinetics and KIE studies clearly indicate that, even if the neighboring Ru centers are chemically linked by flexible carbon–carbon single bonds, it is still possible to achieve water oxidation via the mononuclear WNA mechanism for a given multi-Ru catalyst system. From a perspective of polymer physics, we attributed the observed phenomenon to the unique slow diffusion and multi-charge properties of P4VP–Ru PMCs under catalysis condition, which are generally ignored in the study of small-molecule catalytic systems.
It is well-known that, for a given type of macromolecule in solution, the average diffusion coefficient (〈D〉) generally decreases as the molecular size increases. For small molecules, no significant difference of diffusion coefficients is anticipated; however, it is not the case for polymer ligands due to their significant size effect. Thus, the diffusion behavior of P4VPs and P4VPs–Ru were comparatively investigated by dynamic light scattering (DLS) in MeOH. The obtained diffusion coefficient distributions (f(D)) as well as hydrodynamic radius distributions (f(Rh)) are shown in Fig. S9 and S10.† Notably, three features could be extracted from the results in Fig. S9:† (i) the 〈D〉 values for P4VP–Ru complexes are nearly one order of magnitude smaller than those for P4VP precursors; (ii) P4VP–Ru complexes show much broader diffusion coefficient distributions compared with their P4VP precursors; (iii) 〈D〉 decreases as the chain length increases for both P4VP–Ru complexes and P4VP precursors. Accordingly, an opposite trend was observed for f(Rh) (Fig. S10†).
Quantitatively, Fig. 9a shows the number average molar mass Mn dependence of the average diffusion coefficient 〈D〉 for these samples. The fitted results show that 〈D〉 is approximately scaled to Mn as 〈D〉 ∼ Mnα, where the exponent is close to 0.50 ± 0.05 for P4VPs and P4VPs–Ru, indicating the coil conformations for P4VPs and P4VPs–Ru chains in MeOH. More specifically, 〈D〉 for P4VPs–Ru are ∼0.49 (P4VP21–Ru), ∼0.37 (P4VP42–Ru), ∼0.31 (P4VP74–Ru) and ∼0.36 (P4VP198–Ru) times smaller than their corresponding precursors. It should be noted that only ∼5% of the total pyridine units contain Ru centers, which is not supposed to significantly increase the molar masses of P4VP–Ru PMCs. Moreover, MeOH is a very good solvent for both P4VP backbone and Ru(bda) moiety, which precludes the possibility of inter-chain aggregation induced by solvophobic/solvophilic interaction. Reasonably, the significant decrease in 〈D〉 and increase in polydispersity could be attributed to the chain extension during the DMSO/pyridine ligand exchange reaction. In recent work by Meyer et al. on the interfacial deposition of Ru(bda) complexes,33 they found that the pyridine/pyridine exchange could occur even at room temperature, though with a slow exchange rate. Thus, we expect that some of Ru(bda) centers might be simultaneously coordinated with pyridine units from two different P4VP precursor chains during exchange reaction (Fig. 9b). Overall, the determined 〈D〉 values for P4VP–Ru PMCs are 2–3 orders of magnitude smaller than those for small molecules (typically 4.0 × 10−10 to 9.0 × 10−10 m2 s−1). It is important to note that the slow diffusion of polymer catalyst can greatly enhance the suppression for the inter-chain collision and Ru dimerization processes, which satisfactorily explains why the WNA mechanism was observed in this study. Note that the viscoelastic property of polymer ligand may also play a critical role in preventing the Ru dimerization process.
Reasonably, the slow diffusion feature of P4VPs–Ru can explain why the inter-chain collision is prohibited well, but it still does not help explain why the intra-chain dimerization is also prohibited. Note that the inter-Ru distance for Ru centers on the same P4VP backbone is very short, because the intra-chain Ru centers are separated averagely by only ∼20 monomer units. For classic polyelectrolytes, the multi-charge electrostatic interaction is another key factor to be considered. It is well known that the pKb of pyridine is ∼8.8, which means that most of free pyridine units on polymer ligand should be in the protonated form during catalysis (pH = 1.0). The charged nature of P4VP–Ru PMCs was confirmed by zeta-potential measurement in aqueous solution at pH = 2.0. Fig. S11† shows that the measured zeta-potential is ∼31.4, ∼30.4, ∼31.5 and ∼29.2 for P4VP21–Ru, P4VP42–Ru, P4VP74–Ru and P4VP198–Ru, respectively; in contrast, the zeta-potential is ∼0.0 ± 1.0 for [Ru(bda)(pic)2], indicating a neutral state of small molecule catalyst. It should be noted that, for a suspending nanoparticle in solution, typically, it can be stabilized by the electrostatic repulsion interactions when the interfacial zeta-potential exceeds 15–20. Physically, it is not difficult to understand that the intra-chain electrostatic repulsion interaction can greatly suppress intra-chain Ru-dimerization process due to the enthalpy penalty, which satisfactorily explains why the I2M mechanism is not observed for P4VP–Ru PMC system. The result also reveals that though the existence of flexible linkages between Ru centers could provide enough conformation entropy for the occurrence of Ru-dimerization in small-molecule multi-nuclear system, the entropy elasticity could not overcome the electrostatic interaction energy in PMC system.
To further deepen our understanding of how the intra-chain electrostatic repulsion interaction affects the catalytic process, molecular dynamics simulations were carried out to quantitatively calculate the intra-chain collision probability (Pc) and the electrostatic potential (Ue) for two neighboring Ru centers (molecular dynamics simulation is only for Pc, Ue is calculated by Debye–Hückel Theory). Fig. 10a schematically illustrates the segment conformations for single polymer chains in the charged and neutral states. The simulation result clearly shows that the electrostatic repulsive interaction causes the conformational extension of the single chain and strongly suppresses the collision of intra-chain Ru centers; while a dynamic equilibrium between the states of collision and separation can be established for a neutral polymer chain. Fig. 10b quantitatively shows how the collision probability Pc changes with the dielectric constant (εr) of the system. The result illustrates that when εr = 80 (without screening the electrostatic interaction), Pc = 0, which means that the collision for the two Ru units is completely inhibited in the charged state. As εr increases, the intra-chain electrostatic repulsion interaction decreases, and we observe that Pc exhibits an apparent non-linear increase. Eventually, Pc approaches to the limit of charge-neutral polymer. To highlight the energy barrier caused by the electrostatic repulsion, we further calculated the electrostatic potential (Ue) between two arbitrary segments i and j on polyelectrolyte backbone as a function of the separation distance rij at the ionic concentration of 0.1 M.
Clearly, Ue is approximately 1.5–2.0kBT at rij ∼ 0.45 nm, indicating that the thermal motion of PMC chains (∼kBT) could not overcome the electrostatic repulsive interaction to achieve a collision event, quantitatively explaining why the I2M mechanism is forbidden for PMC system, from a perspective of thermodynamics. Scheme 2 illustrates how the slow diffusion and multi-charge effects (electrostatic interaction) of polyelectrolyte ligand play a synergistically role in determining the catalytic mechanism and make the PMC system different from the small molecule tri-nuclear system.
Scheme 2 Schematic illustration of how the slow diffusion and multi-charge effects (electrostatic interaction) of polyelectrolyte ligand play a synergistically role in determining the catalytic mechanism and make the P4VP–Ru system different from the small molecule tri-nuclear system.19 |
(1) |
(2) |
We set the size of simulation box to be L = 16σ. The polymer is composed of Np = 62 segments, with the 21th and the 42nd ones characterized as segments containing Ru centers, respectively. The system is composed of a total of 3276 particles, and therefore, the concentration of the charged segments in the real units is consistent with our experiments. We vary the dielectric constant of the system from εr = 80 (which corresponds to water) to εr = 200000. The velocity-Verlet algorithm is used to integrate the equation of motion, and the simulations are performed in the NVT ensemble with periodic conditions in three dimensions. We set the temperature to be 1.0, which corresponds to 300 K in the real units, and the time step in solving the equation of motion is t = 0.005. All the simulations are performed using Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) developed by Sandia National Laboratories.36 The final data are taken from 103 to 104 statistically independent samples, each after an equilibration time of the order of 106 MD steps.
Footnotes |
† Electronic supplementary information (ESI) available: Details of the 1H-NMR spectra of all samples, the original schematic diagrams and data of some calculation results, the UV-vis spectra of one P4VP–Ru, the experimental setups of GC and oxygen sensor measurement and the f(Rh) of P4VP–Ru complexes and their corresponding P4VP precursors. See DOI: 10.1039/c8ra08012g |
‡ Tao Zheng and Mo Zhu contributed equally to this work. |
This journal is © The Royal Society of Chemistry 2018 |