Chemical Science View Article Online
Open Access Article. Published on 07 January 2016. Downloaded on 11/09/2017 09:12:49. This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.
EDGE ARTICLE
Cite this: Chem. Sci., 2016, 7, 1773
View Journal | View Issue
Single-chain polybutadiene organometallic nanoparticles: an experimental and theoretical study† Inbal Berkovich,a Sudheendran Mavila,a Olga Iliashevsky,b Sebastian Kozuchac and N. Gabriel Lemcoff*a High molecular weight polybutadienes and rhodium complexes were used to produce single chain organometallic nanoparticles. Irradiation of high cis-polybutadiene in the presence of a photosensitizer
Received 25th November 2015 Accepted 7th January 2016
isomerised the double bonds to produce differing cis/trans ratios within the polymer. Notably, a higher cis percentage of carbon–carbon double bonds within the polymer structure led to faster binding of
DOI: 10.1039/c5sc04535e
metal ions, as well as their faster removal by competing phosphine ligands. The experimental results
www.rsc.org/chemicalscience
were supported and rationalized by DFT computations.
Introduction The creation of new compounds and the study of their properties is one of the most rewarding endeavours for experimental scientists. A new area of macromolecular research deals with the study and preparation of organic nanoparticles (ONPs).1–4 An emerging technique for the facile generation of tailor-made well-dened ONPs is the intramolecular cross-linking of single polymer chains; these ONPs are called single chain nanoparticles (SCNPs).5–11 Moreover, the insertion of metals to these polymeric structures poses promising prospects for various applications, such as catalysis, drug delivery and sensing.12–18 In part, this potential arises because metallopolymers may combine the physical and chemical properties of both polymer and the embedded metal. Thus, the study of the factors affecting the formation, folding and architecture of metal containing SCNPs is an appealing goal.19–25 The concept of single chain collapse driven by metal coordination has been introduced recently by several research groups. Pomposo and coworkers reported on the preparation of organometallic SCNPs by the intramolecular cross linking of copolymers bearing bketoester chelating groups with Cu(II).19 The use of these nanoparticles in oxidative coupling reactions revealed an enhanced catalytic specicity compared to the classical catalysts such as CuCl2, Cu(OAc)2 and Cu(acac)2. In another recent a
Department of Chemistry, Ben-Gurion University of the Negev, Beer-Sheva 84105, Israel. E-mail: lemcoff@bgu.ac.il
b
Department of Chemical Engineering, Ben-Gurion University of the Negev, Beer-Sheva 84105, Israel
c
Lise Meitner – Minerva Center for Computational Quantum Chemistry, Israel
† Electronic supplementary information (ESI) available: Detailed experimental procedures and characterization of ONPs including: 1H and 13C NMR, SEC, DLS, DSC, TGA, UV-vis, GC-MS spectra and computational calculations. See DOI: 10.1039/c5sc04535e
This journal is © The Royal Society of Chemistry 2016
work, Pomposo et al. have extended this methodology for the preparation of water-soluble globular SCNPs with polymerase mimicking activity.20 Paik et al. detailed a cyclotetramerization cross-linking reaction within polymers bearing phthalonitrile pendant groups with CuCl resulting in SCNPs containing copper phthalocyanines.21 The Barner-Kowollik group recently used telechelic phosphine ligands to fold polymers with Pd(II).22 Moreover, Pd-SCNPs formed by a similar strategy have shown catalytic activity for a benchmark Sonogashira coupling reaction.23 The carbon–carbon double bond p-electron coordination to metals is one of the most ubiquitous bonding types in organometallic chemistry.26–28 A recent example by Manners et al. shows how Pt ions may be used to reversibly cross-link a corona of polyisoprene to control the overall architecture of the macrostructure.29 Our group has recently developed an intramolecular chain collapse approach for the synthesis of well-dened organometallic ONPs by a direct exchange of labile ligands of chlorobis(ethylene)rhodium(I) dimer by 1,5-hexadiene units present in poly-1,5-cis-cyclooctadiene (PCOD) produced by ring-opening metathesis polymerization (ROMP).24 The method was also expanded for the preparation of SCNPs with Ir(I) and Ni(0), showing promising applications in catalysis.25 The polymers used were about 50 kDa in size and had 20% cis content; typical values for ROMP reactions using Grubbs initiators.30,31 While the production of PCOD by this method is fairly straightforward, we sought for a way to achieve higher molecular weight polymers and to control the cis/trans ratio, as it may widen the scope of applications for this family of SCNPs and in turn afford new insights on the intricacies of metal–p bonding in these materials. 1,4-Polybutadiene (PBD), which has the same repeat unit as PCOD, is the second largest synthetic rubber product by volume and the one with the best performance to price ratio (1.1$ kg1, 18.2 performance/price ratio).32 In addition, it is available in high molecular weights
Chem. Sci., 2016, 7, 1773–1778 | 1773
View Article Online
Chemical Science
Edge Article
Open Access Article. Published on 07 January 2016. Downloaded on 11/09/2017 09:12:49. This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.
and various cis/trans contents (depending on the catalyst used). Herein, we report on the preparation and characterization of rhodium(I) organometallic nanoparticles prepared by single chain collapse of high molecular weight polybutadienes of differing cis/trans ratios and the experimental and computational study of the effect of the double bond stereochemistry on the polymers' metal-binding properties.
Results and discussions Synthesis of Rh(I)-ONPs from PBD with varying cis content 1,4-Polybutadienes with varying cis/trans content were obtained by irradiating toluene solutions of commercial PBD (3 105 g mol1, PDI ¼ 2.5, 95% cis, 5% vinyl) with 365 nm UV light in the presence of diphenyl disulde as a photosensitizer for different periods of time (Fig. S1†).33 The polymers obtained had a nal cis content of 95%, 72%, 51% and 20% (determined by 1H-NMR) without a signicant change in their molecular weights. Having obtained the series of PBDs with varying cis/trans ratios, their single chain collapse was attempted via the addition of chlorobis(ethylene)rhodium(I) dimer (1), as rhodium donor, under dilute conditions. To our satisfaction, the PBDs quickly reacted at room temperature to afford polymer bound complexes 2 (Scheme 1), as determined by UV-vis, NMR, DLS and GPC analyses. The size of the organometallic SCNP was controlled by systematically varying the percentage of rhodium(I) introduced onto the polymer chain. As in the PCOD examples, a linear dependence between the amount of rhodium(I) (cross-linker) added and the reduction in the hydrodynamic radius of the Rh(I)-ONP was observed. The degree of chain collapse was
Overlay of SEC plots obtained for ONPs prepared with varying concentrations of Rh(I).
Fig. 1
Table 1
Rh(I)-ONPs with varying rhodium(I) contenta
Entry
Rh(I) [mol%]b
Mwc [105 g mol1]
PDI
1 2 3 4
0 2 5 8
3.14 3.21 2.92 3.11
2.47 1.97 1.91 1.62
c
Intrinsic viscosityc (mL g1)
Rhd (nm)
168.4 150.6 110.8 85.5
27.0 23.5 16.6 12.8
Conditions: solvent ¼ THF, [PBD] ¼ 0.5 mg ml1 (1.5 mM), T ¼ 35 C. Condition: Rh(I) ¼ [RhCl(C2H4)2]2, mol% relative to 1,5-hexadiene units. c Condition: determined by triple-detector SEC in THF. d Condition: determined by DLS in THF. a b
monitored by both SEC and dynamic light scattering (DLS) in THF (Fig. 1 and Table 1, Fig. S9–S12†). Kinetic study of rhodium(I) complexation
Scheme 1
Schematic illustration of the preparation of Rh(I)-ONPs.
1774 | Chem. Sci., 2016, 7, 1773–1778
The observation that 1,4-polybutadienes with different cis content readily react with 1 to produce well-dened ONPs gave us the opportunity to study the effect of the double bond stereochemistry on the metal binding propensity. Thus, the rate of rhodium complexation was monitored with stopped ow kinetics by observing the change in the characteristic UV-vis absorption spectrum. ROMP-derived PCOD (50 kDa, PDI 1.43, 20% cis, 80% trans) was used as a control under the same conditions. In each experiment, 10 mol% of 1 was mixed with a dilute solution of polymer and a plot of the absorbance at 360 nm versus time was recorded. Interestingly, the rate constants increased almost linearly as the cis content of the polymer was raised (Fig. 2). In addition, both 20% cis PBD and 20% cis PCOD showed similar rate constants, suggesting that the rate of the complexation is independent of the polymer size or origin. The inuence of the double bond stereochemistry on binding was also studied with two 4,8-dodecadienes, one all cis and the other 75% trans, as additional model compounds. The reaction of the high cis
This journal is © The Royal Society of Chemistry 2016
View Article Online
Open Access Article. Published on 07 January 2016. Downloaded on 11/09/2017 09:12:49. This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.
Edge Article
Fig. 2 Binding rate constants versus cis content of PBD and PCOD for the formation of Rh(I)-ONPs.
content 4,8-dodecadiene with 1 was about 7 times faster than the reaction of 75% trans 4,8-dodecadiene in accordance with the polymer experiments (Fig. S19†). These results suggest that the intrinsic stereochemistry of the carbon–carbon double bond inuences the metal binding, as a similar trend was observed for both the small molecule model compounds and the polymers. Computational study on the effect of double bond stereochemistry on rhodium binding Quantum mechanical computations (at the M06/def2tzvpp(THF)//b97d3/def2-svp level,34–38 see the ESI†) shed some light on the reason for the faster association of the cis conformers of the 1,4-polybutadienes. cis- and trans-2-butene were used as models for the rst step of the complexation process with complex 1. The rst reaction with 2-butene, as shown in Fig. 3, involves four elementary steps. The rst step is the formation of an
Fig. 3 Mechanism for the first ligand substitution of 2-butene with ethene. The numbers are Gibbs relative energies in kJ mol1 for the cis-2-butene, in a THF continuum solvent.
This journal is © The Royal Society of Chemistry 2016
Chemical Science
adduct stabilized only by van der Waals interactions (1b), which – although exothermic – must pay a high entropy penalty, thus resulting endergonic. The second step is the insertion of the butene to form an 18e complex (1c through transition state TSbc). This is followed by the dissociation of the ethene (TScd and 1d), which results in the determining transition state for the complete ligand substitution. Finally, the ethene moiety is completely released to form 1e. The complete reaction should result endergonic (at standard concentrations), but it is driven by the release of ethene as a gas. Several conformers were considered (including the inversion of the bent structure between the two coordination planes),39 but only the most relevant conformers will be discussed here (for higher energy structures see the ESI†). Fig. 3 shows the relative Gibbs energies in THF for the pathway including cis-2-butene. All the other pathways have very similar energy proles, and therefore only the relation between the DG of the determining transition states40–42 (TScd) of the various mechanisms is relevant to understand the experimentally observed kinetics for the cis and trans species. The most probable mechanism for the trans-2butene substitution presents a relative DG of 76.1 kJ mol1 for TScd, only 3.1 kJ mol1 higher than the equivalent state in the cis reaction. At room temperature this corresponds to 3.5 times faster reaction for the ethene substitution with a cis conformer compared to the trans one, agreeing with the experimental results (considering the computational and experimental differences in the reactant models). The reason for the higher energy of the trans pathway is shown in Fig. 4: the repulsion between the methyl group of the trans conformer and the ethene neighbour ligand produces a distortion of the verticality of the butene, which translates in a slightly higher energy all along the reaction (see the ESI†). While in the cis conformer the Cl–Rh–C–C dihedral angle is virtually perpendicular (89 ), the two possible trans conformers have angles of 77 and 109 (Fig. 4). Note that the chloride ligand in the opposite side of the ethene also exerts some steric hindrance, pushing out the methyl groups of the butene. However, in the cis case the chloride ligand pushes both methyl groups simultaneously, without affecting the dihedral angle of the butene and maintaining the back-bonding character of the metal–olen bond. Although the second ligand substitution may be of kinetic importance in the case of 2-butene, for chelating dienes there will be a small entropy penalty on the second substitution, as
Fig. 4 (A) cis conformer of the monosubstituted dirhodium complex 1e, where the butene double bond is perpendicular to the complex plane, optimizing the back-bonding interaction. (B) trans conformers; the curves show steric impediments which diminish the p backbonding of the metal to the ligand.
Chem. Sci., 2016, 7, 1773–1778 | 1775
View Article Online
Chemical Science
Edge Article
the second double bond is already in the vicinity of the metal. In the overall reaction there will be a net entropy gain when a diene species displaces two ethene molecules, and as a result the only critical step of the complete mechanism will be at the rst substitution.
Open Access Article. Published on 07 January 2016. Downloaded on 11/09/2017 09:12:49. This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.
Kinetic and computational study of rhodium(I) release Having shown that this type of organometallic ONPs are reversible in nature,24 the release of Rh from PBDs was also studied by addition of a competing ligand, tricyclohexylphosphine (PCy3, Scheme 2). Thus, excess PCy3 was added to 10 mol% Rh(I)-ONPs and the absorbance at 360 nm was monitored using stopped ow kinetics. Addition of the phosphine led to a reduction in the absorption at 360 nm, indicating the disruption of the organometallic bond. A plot of the rate constants versus the cis content of the polymer showed that the reaction of 95% cis Rh(I)-ONPs with PCy3 was the fastest among the series of PBDs used; while the slowest rate constants were obtained with 20% cis-PCOD and PBD (Fig. 5). Also rhodium complexes of 4,8-dodecadienes were mixed with PCy3 and the absorption at 360 nm was monitored (Fig. S20†). The control experiment with cis,cis-4,8-dodecadiene–rhodium complex was also faster than that with the trans isomer. Thus, experiments showed that not only the association, but also the dissociation of the alkenes – driven by substitutions with PCy3 – are faster in the cis system, a fact that may seem to be at odds with the greater stability of the h2 bond in intermediate complex 1e with cis-2-butene compared to trans-2-butene. To gain a computational insight into this occurrence, 2,6-octadiene with a double cis or double trans conguration was used as a model. Computations showed that the cis-2,6-octadiene forms a complex 36 kJ mol1 higher in energy than the double trans case, thus being more reactive to ligand substitution.
Rate constants versus cis content of PBD and PCOD for the disruption of Rh(I)-ONPs.
Fig. 5
cis (A) and trans (B) 2,6-octadienes as chelating double h2 ligands. In the cis species the double bonds cannot achieve the ideal vertical position that maximizes the back-bonding interaction, and therefore results in a higher energy, more reactive complex. Fig. 6
Again, this is mainly due to strain arising when the olen double bonds cannot achieve a vertical position (that is, perpendicular to the complex plane where the back-bonding interaction is maximized). As can be seen in Fig. 6, the double trans diene can be positioned with both p bonds perpendicular to the Cl–Rh–Cl plane, a preferred geometry, but the double cis conformer (Fig. 6A) is forced to rotate one of the double bonds. In summary, the formation of complexes based on cis-dienes is faster than with the trans ligands because the rst substitution (which determines the rate of association) has less steric impediments in the cis case. The dissociation of cis-dienes is also faster than trans, since the chelating diene in the cis conformer cannot achieve a stable geometry with both double bonds perpendicular to the complex plane, partially breaking the back bonding stabilization and making this a more reactive complex prone to a faster ligand substitution, in accordance to the experimental observations.
The effect of the intermolecular cross-linking of PBD on the thermal properties
Scheme 2
Reaction of 2 with PCy3.
1776 | Chem. Sci., 2016, 7, 1773–1778
Having shown that dilute solutions of polymers containing the 1,5-hexadiene unit have the ability to form organometallic complexes, we sought to investigate whether the bulk polymer could also act as a metal sponge and change its properties when ions are absorbed. It has been shown that the coordination of transition metals to olen containing polymers, such as Pd(II)
This journal is © The Royal Society of Chemistry 2016
View Article Online
Edge Article Table 2
Tg of 95% cis-PBD with varying amount of Rh(I)
Entry
Rh(I) [mol%]
Tg ( C)
1 2 3
0 5 10
101.9 88.5 71.2
a
Open Access Article. Published on 07 January 2016. Downloaded on 11/09/2017 09:12:49. This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.
Chemical Science
a
mol% relative to 1,5-hexadiene units.
with atactic 1,2-polybutadiene, has a signicant effect on the glass transition temperature (Tg) of the composite.43 It has also been reported that Rh containing PBD copolymers (SBS) could be prepared by heating a THF or benzene solution in the presence of RhCl3.44 However, the inuence of different amounts of metal coordination on the Tg of 1,4-polybutadiene has not been studied. Thus, samples of PBD were soaked overnight in THF solutions with varying concentrations of rhodium. The effect of this procedure was studied by analyzing the change in the glass transition temperature of the polymer sample and the shi in the UV-vis absorption spectrum (Fig. S21†). Indeed, as shown in Table 2, the metal led to a signicant change in the polymer's physical properties as expected from the intermolecular crosslinking caused by the metal bridges.
Conclusions Polybutadiene, an abundant starting material, was used for the rst time as the precursor single polymer chain for the formation of organometallic SCNPs. Because these polymers are much larger than those previously used, they may bind many more metal ions per nanoparticle; an important requisite for future applications (e.g. templates for inorganic nanoparticles). Interestingly, the introduction of rhodium ions to neat PBD produced increased glass transition temperatures as a consequence of the restricted local chain mobility caused by the metal bridges. By investigating polymers with differing cis/trans ratios, we have experimentally and computationally shown that there is a stark difference in the rate of metal sequestration and release when polymers with different microstructures are used to form organometallic nanoparticles. Indeed, polymers with higher ratio of cis double bonds bind rhodium ions faster and, when competitive ligands are added, also release them faster. To put this in perspective, under our experimental conditions the rhodium ions could be bound in less than ten seconds by a 300 kDa all cis-PBD and they could be completely removed in about seven minutes by the addition of excess phosphine, while with a high trans-PBD it would take thrice as long. The availability of cis-selective olen metathesis catalysts45,46 may allow the formation of novel 1,5-cyclooctadiene polymers or copolymers with high cis/trans ratios which may be useful for quick binding and release of selected metal ions for future applications of this method. DFT based computational studies supported the experimental results and revealed the mechanistic and steric reasons for the behaviour observed. Taking into account the great importance of PBD cross-linking in the rubber industry, the study of so, reversible, cross-linking
This journal is © The Royal Society of Chemistry 2016
methodologies by metals and the way the polymer structure affects the metal intake is of signicant value.
Acknowledgements The United States – Israel Binational Science Foundation (Grant No. 2014116) is gratefully acknowledged for partial funding of this research.
Notes and references 1 J. P. Rao and K. E. Geckeler, Prog. Polym. Sci., 2011, 36, 887– 913. 2 K. Li and B. Liu, Chem. Soc. Rev., 2014, 43, 6570–6597. 3 M. Elsabahy and K. L. Wooley, J. Polym. Sci., Part A: Polym. Chem., 2012, 50, 1869–1880. 4 M. Elsabahy, G. S. Heo, S.-M. Lim, G. Sun and K. L. Wooley, Chem. Rev., 2015, 115, 10967–11011. 5 S. Mavila, O. Eivgi, I. Berkovich and N. G. Lemcoff, Chem. Rev., 2015, DOI: 10.1021/acs.chemrev.5b00290. 6 C. K. Lyon, A. Prasher, A. M. Hanlon, B. T. Tuten, C. A. Tooley, P. G. Frank and E. B. Berda, Polym. Chem., 2015, 6, 181–197. 7 O. Altintas and C. Barner-Kowollik, Macromol. Rapid Commun., 2012, 33, 958–971. 8 M. Gonzalez-Burgos, A. Latorre-Sanchez and J. A. Pomposo, Chem. Soc. Rev., 2015, 44, 6122–6142. 9 M. Aiertza, I. Odriozola, G. Caba˜ nero, H.-J. Grande and I. Loinaz, Cell. Mol. Life Sci., 2012, 69, 337–346. 10 T. Mes, R. van der Weegen, A. R. A. Palmans and E. W. Meijer, Angew. Chem., Int. Ed., 2011, 50, 5085–5089. 11 P. J. M. Stals, M. A. J. Gillissen, T. F. E. Paffen, T. F. A. de Greef, P. Lindner, E. W. Meijer, A. R. A. Palmans and I. K. Voets, Macromolecules, 2014, 47, 2947–2954. 12 K. A. Williams, A. J. Boydston and C. W. Bielawski, Chem. Soc. Rev., 2007, 36, 729–744. 13 M. Ramanathan and S. B. Darling, Polym. Int., 2013, 62, 1123–1134. 14 G. R. Whittell, M. D. Hager, U. S. Schubert and I. Manners, Nat. Mater., 2011, 10, 176–188. 15 J. Xiang, C.-L. Ho and W.-Y. Wong, Polym. Chem., 2015, 6, 6905–6930. 16 X. Wang and R. McHale, Macromol. Rapid Commun., 2010, 31, 331–350. 17 Y. Liu, T. Pauloehrl, S. I. Presolski, L. Albertazzi, A. R. A. Palmans and E. W. Meijer, J. Am. Chem. Soc., 2015, 137, 13096–13105. 18 M. Artar, E. R. J. Souren, T. Terashima, E. W. Meijer and A. R. A. Palmans, ACS Macro Lett., 2015, 4, 1099–1103. 19 A. Sanchez-Sanchez, A. Arbe, J. Colmenero and J. A. Pomposo, ACS Macro Lett., 2014, 3, 439–443. 20 A. Sanchez-Sanchez, A. Arbe, J. Kohlbrecher, J. Colmenero and J. A. Pomposo, Macromol. Rapid Commun., 2015, 36, 1592–1597. 21 J. Jeong, Y.-J. Lee, B. Kim, B. Kim, K.-S. Jung and H.-J. Paik, Polym. Chem., 2015, 6, 3392–3397.
Chem. Sci., 2016, 7, 1773–1778 | 1777
View Article Online
Open Access Article. Published on 07 January 2016. Downloaded on 11/09/2017 09:12:49. This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.
Chemical Science
22 J. Willenbacher, O. Altintas, P. W. Roesky and C. BarnerKowollik, Macromol. Rapid Commun., 2014, 35, 45–51. 23 J. Willenbacher, O. Altintas, V. Trouillet, N. Knofel, M. J. Monteiro, P. W. Roesky and C. Barner-Kowollik, Polym. Chem., 2015, 6, 4358–4365. 24 S. Mavila, C. E. Diesendruck, S. Linde, L. Amir, R. Shikler and N. G. Lemcoff, Angew. Chem., Int. Ed., 2013, 52, 5767– 5770. 25 S. Mavila, I. Rozenberg and N. G. Lemcoff, Chem. Sci., 2014, 5, 4196–4203. 26 C. Huber, F. Bangerter, W. R. Caseri and C. Weder, J. Am. Chem. Soc., 2001, 123, 3857–3863. 27 J. B. Johnson and T. Rovis, Angew. Chem., Int. Ed., 2008, 47, 840–871. 28 M. A. Bennett, Chem. Rev., 1962, 62, 611–652. 29 P. A. Rupar, G. Cambridge, M. A. Winnik and I. Manners, J. Am. Chem. Soc., 2011, 133, 16947–16957. 30 C. W. Bielawski and R. H. Grubbs, Angew. Chem., Int. Ed., 2000, 39, 2903–2906. 31 B. Allaert, N. Ledoux, N. Dieltiens, H. V. Mierde, C. V. Stevens, P. van Der Voort and F. Verpoort, Catal. Commun., 2008, 9, 1054–1059. 32 H. Brandt, Ullmann's encyclopedia of industrial chemistry, Wiley-VCH, 2012. 33 M. A. Golub, J. Polym. Sci., 1957, 25, 373–377.
1778 | Chem. Sci., 2016, 7, 1773–1778
Edge Article
34 Y. Zhao and D. Truhlar, Theor. Chem. Acc., 2008, 120, 215– 241. 35 S. Grimme, J. Comput. Chem., 2006, 27, 1787–1799. 36 S. Grimme, S. Ehrlich and L. Goerigk, J. Comput. Chem., 2011, 32, 1456–1465. 37 A. V. Marenich, C. J. Cramer and D. G. Truhlar, J. Phys. Chem. B, 2009, 113, 6378–6396. 38 F. Weigend and R. Ahlrichs, Phys. Chem. Chem. Phys., 2005, 7, 3297–3305. 39 G. Aull´ on, G. Ujaque, A. Lled´ os, S. Alvarez and P. Alemany, Inorg. Chem., 1998, 37, 804–813. 40 J. R. Murdoch, J. Chem. Educ., 1981, 58, 32–36. 41 S. Kozuch and S. Shaik, Acc. Chem. Res., 2011, 44, 101–110. 42 S. Kozuch, in Understanding Organometallic Reaction Mechanisms and Catalysis, Wiley-VCH Verlag GmbH & Co. KGaA, 2014, pp. 217–248, DOI: 10.1002/9783527678211.ch9. 43 L. A. Belore and M. P. McCurdie, J. Polym. Sci., Part B: Polym. Phys., 1995, 33, 105–124. 44 L. M. Bronstein and P. M. Valetsky, J. Inorg. Organomet. Polym., 1994, 4, 415–424. 45 B. K. Keitz, A. Fedorov and R. H. Grubbs, J. Am. Chem. Soc., 2012, 134, 2040–2043. 46 M. B. Herbert and R. H. Grubbs, Angew. Chem., Int. Ed., 2015, 54, 5018–5024.
This journal is © The Royal Society of Chemistry 2016
Open Access Article. Published on 07 January 2016. Downloaded on 11/09/2017 09:12:49. This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.
EDGE ARTICLE
Cite this: Chem. Sci., 2016, 7, 1773
View Journal | View Issue
Single-chain polybutadiene organometallic nanoparticles: an experimental and theoretical study† Inbal Berkovich,a Sudheendran Mavila,a Olga Iliashevsky,b Sebastian Kozuchac and N. Gabriel Lemcoff*a High molecular weight polybutadienes and rhodium complexes were used to produce single chain organometallic nanoparticles. Irradiation of high cis-polybutadiene in the presence of a photosensitizer
Received 25th November 2015 Accepted 7th January 2016
isomerised the double bonds to produce differing cis/trans ratios within the polymer. Notably, a higher cis percentage of carbon–carbon double bonds within the polymer structure led to faster binding of
DOI: 10.1039/c5sc04535e
metal ions, as well as their faster removal by competing phosphine ligands. The experimental results
www.rsc.org/chemicalscience
were supported and rationalized by DFT computations.
Introduction The creation of new compounds and the study of their properties is one of the most rewarding endeavours for experimental scientists. A new area of macromolecular research deals with the study and preparation of organic nanoparticles (ONPs).1–4 An emerging technique for the facile generation of tailor-made well-dened ONPs is the intramolecular cross-linking of single polymer chains; these ONPs are called single chain nanoparticles (SCNPs).5–11 Moreover, the insertion of metals to these polymeric structures poses promising prospects for various applications, such as catalysis, drug delivery and sensing.12–18 In part, this potential arises because metallopolymers may combine the physical and chemical properties of both polymer and the embedded metal. Thus, the study of the factors affecting the formation, folding and architecture of metal containing SCNPs is an appealing goal.19–25 The concept of single chain collapse driven by metal coordination has been introduced recently by several research groups. Pomposo and coworkers reported on the preparation of organometallic SCNPs by the intramolecular cross linking of copolymers bearing bketoester chelating groups with Cu(II).19 The use of these nanoparticles in oxidative coupling reactions revealed an enhanced catalytic specicity compared to the classical catalysts such as CuCl2, Cu(OAc)2 and Cu(acac)2. In another recent a
Department of Chemistry, Ben-Gurion University of the Negev, Beer-Sheva 84105, Israel. E-mail: lemcoff@bgu.ac.il
b
Department of Chemical Engineering, Ben-Gurion University of the Negev, Beer-Sheva 84105, Israel
c
Lise Meitner – Minerva Center for Computational Quantum Chemistry, Israel
† Electronic supplementary information (ESI) available: Detailed experimental procedures and characterization of ONPs including: 1H and 13C NMR, SEC, DLS, DSC, TGA, UV-vis, GC-MS spectra and computational calculations. See DOI: 10.1039/c5sc04535e
This journal is © The Royal Society of Chemistry 2016
work, Pomposo et al. have extended this methodology for the preparation of water-soluble globular SCNPs with polymerase mimicking activity.20 Paik et al. detailed a cyclotetramerization cross-linking reaction within polymers bearing phthalonitrile pendant groups with CuCl resulting in SCNPs containing copper phthalocyanines.21 The Barner-Kowollik group recently used telechelic phosphine ligands to fold polymers with Pd(II).22 Moreover, Pd-SCNPs formed by a similar strategy have shown catalytic activity for a benchmark Sonogashira coupling reaction.23 The carbon–carbon double bond p-electron coordination to metals is one of the most ubiquitous bonding types in organometallic chemistry.26–28 A recent example by Manners et al. shows how Pt ions may be used to reversibly cross-link a corona of polyisoprene to control the overall architecture of the macrostructure.29 Our group has recently developed an intramolecular chain collapse approach for the synthesis of well-dened organometallic ONPs by a direct exchange of labile ligands of chlorobis(ethylene)rhodium(I) dimer by 1,5-hexadiene units present in poly-1,5-cis-cyclooctadiene (PCOD) produced by ring-opening metathesis polymerization (ROMP).24 The method was also expanded for the preparation of SCNPs with Ir(I) and Ni(0), showing promising applications in catalysis.25 The polymers used were about 50 kDa in size and had 20% cis content; typical values for ROMP reactions using Grubbs initiators.30,31 While the production of PCOD by this method is fairly straightforward, we sought for a way to achieve higher molecular weight polymers and to control the cis/trans ratio, as it may widen the scope of applications for this family of SCNPs and in turn afford new insights on the intricacies of metal–p bonding in these materials. 1,4-Polybutadiene (PBD), which has the same repeat unit as PCOD, is the second largest synthetic rubber product by volume and the one with the best performance to price ratio (1.1$ kg1, 18.2 performance/price ratio).32 In addition, it is available in high molecular weights
Chem. Sci., 2016, 7, 1773–1778 | 1773
View Article Online
Chemical Science
Edge Article
Open Access Article. Published on 07 January 2016. Downloaded on 11/09/2017 09:12:49. This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.
and various cis/trans contents (depending on the catalyst used). Herein, we report on the preparation and characterization of rhodium(I) organometallic nanoparticles prepared by single chain collapse of high molecular weight polybutadienes of differing cis/trans ratios and the experimental and computational study of the effect of the double bond stereochemistry on the polymers' metal-binding properties.
Results and discussions Synthesis of Rh(I)-ONPs from PBD with varying cis content 1,4-Polybutadienes with varying cis/trans content were obtained by irradiating toluene solutions of commercial PBD (3 105 g mol1, PDI ¼ 2.5, 95% cis, 5% vinyl) with 365 nm UV light in the presence of diphenyl disulde as a photosensitizer for different periods of time (Fig. S1†).33 The polymers obtained had a nal cis content of 95%, 72%, 51% and 20% (determined by 1H-NMR) without a signicant change in their molecular weights. Having obtained the series of PBDs with varying cis/trans ratios, their single chain collapse was attempted via the addition of chlorobis(ethylene)rhodium(I) dimer (1), as rhodium donor, under dilute conditions. To our satisfaction, the PBDs quickly reacted at room temperature to afford polymer bound complexes 2 (Scheme 1), as determined by UV-vis, NMR, DLS and GPC analyses. The size of the organometallic SCNP was controlled by systematically varying the percentage of rhodium(I) introduced onto the polymer chain. As in the PCOD examples, a linear dependence between the amount of rhodium(I) (cross-linker) added and the reduction in the hydrodynamic radius of the Rh(I)-ONP was observed. The degree of chain collapse was
Overlay of SEC plots obtained for ONPs prepared with varying concentrations of Rh(I).
Fig. 1
Table 1
Rh(I)-ONPs with varying rhodium(I) contenta
Entry
Rh(I) [mol%]b
Mwc [105 g mol1]
PDI
1 2 3 4
0 2 5 8
3.14 3.21 2.92 3.11
2.47 1.97 1.91 1.62
c
Intrinsic viscosityc (mL g1)
Rhd (nm)
168.4 150.6 110.8 85.5
27.0 23.5 16.6 12.8
Conditions: solvent ¼ THF, [PBD] ¼ 0.5 mg ml1 (1.5 mM), T ¼ 35 C. Condition: Rh(I) ¼ [RhCl(C2H4)2]2, mol% relative to 1,5-hexadiene units. c Condition: determined by triple-detector SEC in THF. d Condition: determined by DLS in THF. a b
monitored by both SEC and dynamic light scattering (DLS) in THF (Fig. 1 and Table 1, Fig. S9–S12†). Kinetic study of rhodium(I) complexation
Scheme 1
Schematic illustration of the preparation of Rh(I)-ONPs.
1774 | Chem. Sci., 2016, 7, 1773–1778
The observation that 1,4-polybutadienes with different cis content readily react with 1 to produce well-dened ONPs gave us the opportunity to study the effect of the double bond stereochemistry on the metal binding propensity. Thus, the rate of rhodium complexation was monitored with stopped ow kinetics by observing the change in the characteristic UV-vis absorption spectrum. ROMP-derived PCOD (50 kDa, PDI 1.43, 20% cis, 80% trans) was used as a control under the same conditions. In each experiment, 10 mol% of 1 was mixed with a dilute solution of polymer and a plot of the absorbance at 360 nm versus time was recorded. Interestingly, the rate constants increased almost linearly as the cis content of the polymer was raised (Fig. 2). In addition, both 20% cis PBD and 20% cis PCOD showed similar rate constants, suggesting that the rate of the complexation is independent of the polymer size or origin. The inuence of the double bond stereochemistry on binding was also studied with two 4,8-dodecadienes, one all cis and the other 75% trans, as additional model compounds. The reaction of the high cis
This journal is © The Royal Society of Chemistry 2016
View Article Online
Open Access Article. Published on 07 January 2016. Downloaded on 11/09/2017 09:12:49. This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.
Edge Article
Fig. 2 Binding rate constants versus cis content of PBD and PCOD for the formation of Rh(I)-ONPs.
content 4,8-dodecadiene with 1 was about 7 times faster than the reaction of 75% trans 4,8-dodecadiene in accordance with the polymer experiments (Fig. S19†). These results suggest that the intrinsic stereochemistry of the carbon–carbon double bond inuences the metal binding, as a similar trend was observed for both the small molecule model compounds and the polymers. Computational study on the effect of double bond stereochemistry on rhodium binding Quantum mechanical computations (at the M06/def2tzvpp(THF)//b97d3/def2-svp level,34–38 see the ESI†) shed some light on the reason for the faster association of the cis conformers of the 1,4-polybutadienes. cis- and trans-2-butene were used as models for the rst step of the complexation process with complex 1. The rst reaction with 2-butene, as shown in Fig. 3, involves four elementary steps. The rst step is the formation of an
Fig. 3 Mechanism for the first ligand substitution of 2-butene with ethene. The numbers are Gibbs relative energies in kJ mol1 for the cis-2-butene, in a THF continuum solvent.
This journal is © The Royal Society of Chemistry 2016
Chemical Science
adduct stabilized only by van der Waals interactions (1b), which – although exothermic – must pay a high entropy penalty, thus resulting endergonic. The second step is the insertion of the butene to form an 18e complex (1c through transition state TSbc). This is followed by the dissociation of the ethene (TScd and 1d), which results in the determining transition state for the complete ligand substitution. Finally, the ethene moiety is completely released to form 1e. The complete reaction should result endergonic (at standard concentrations), but it is driven by the release of ethene as a gas. Several conformers were considered (including the inversion of the bent structure between the two coordination planes),39 but only the most relevant conformers will be discussed here (for higher energy structures see the ESI†). Fig. 3 shows the relative Gibbs energies in THF for the pathway including cis-2-butene. All the other pathways have very similar energy proles, and therefore only the relation between the DG of the determining transition states40–42 (TScd) of the various mechanisms is relevant to understand the experimentally observed kinetics for the cis and trans species. The most probable mechanism for the trans-2butene substitution presents a relative DG of 76.1 kJ mol1 for TScd, only 3.1 kJ mol1 higher than the equivalent state in the cis reaction. At room temperature this corresponds to 3.5 times faster reaction for the ethene substitution with a cis conformer compared to the trans one, agreeing with the experimental results (considering the computational and experimental differences in the reactant models). The reason for the higher energy of the trans pathway is shown in Fig. 4: the repulsion between the methyl group of the trans conformer and the ethene neighbour ligand produces a distortion of the verticality of the butene, which translates in a slightly higher energy all along the reaction (see the ESI†). While in the cis conformer the Cl–Rh–C–C dihedral angle is virtually perpendicular (89 ), the two possible trans conformers have angles of 77 and 109 (Fig. 4). Note that the chloride ligand in the opposite side of the ethene also exerts some steric hindrance, pushing out the methyl groups of the butene. However, in the cis case the chloride ligand pushes both methyl groups simultaneously, without affecting the dihedral angle of the butene and maintaining the back-bonding character of the metal–olen bond. Although the second ligand substitution may be of kinetic importance in the case of 2-butene, for chelating dienes there will be a small entropy penalty on the second substitution, as
Fig. 4 (A) cis conformer of the monosubstituted dirhodium complex 1e, where the butene double bond is perpendicular to the complex plane, optimizing the back-bonding interaction. (B) trans conformers; the curves show steric impediments which diminish the p backbonding of the metal to the ligand.
Chem. Sci., 2016, 7, 1773–1778 | 1775
View Article Online
Chemical Science
Edge Article
the second double bond is already in the vicinity of the metal. In the overall reaction there will be a net entropy gain when a diene species displaces two ethene molecules, and as a result the only critical step of the complete mechanism will be at the rst substitution.
Open Access Article. Published on 07 January 2016. Downloaded on 11/09/2017 09:12:49. This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.
Kinetic and computational study of rhodium(I) release Having shown that this type of organometallic ONPs are reversible in nature,24 the release of Rh from PBDs was also studied by addition of a competing ligand, tricyclohexylphosphine (PCy3, Scheme 2). Thus, excess PCy3 was added to 10 mol% Rh(I)-ONPs and the absorbance at 360 nm was monitored using stopped ow kinetics. Addition of the phosphine led to a reduction in the absorption at 360 nm, indicating the disruption of the organometallic bond. A plot of the rate constants versus the cis content of the polymer showed that the reaction of 95% cis Rh(I)-ONPs with PCy3 was the fastest among the series of PBDs used; while the slowest rate constants were obtained with 20% cis-PCOD and PBD (Fig. 5). Also rhodium complexes of 4,8-dodecadienes were mixed with PCy3 and the absorption at 360 nm was monitored (Fig. S20†). The control experiment with cis,cis-4,8-dodecadiene–rhodium complex was also faster than that with the trans isomer. Thus, experiments showed that not only the association, but also the dissociation of the alkenes – driven by substitutions with PCy3 – are faster in the cis system, a fact that may seem to be at odds with the greater stability of the h2 bond in intermediate complex 1e with cis-2-butene compared to trans-2-butene. To gain a computational insight into this occurrence, 2,6-octadiene with a double cis or double trans conguration was used as a model. Computations showed that the cis-2,6-octadiene forms a complex 36 kJ mol1 higher in energy than the double trans case, thus being more reactive to ligand substitution.
Rate constants versus cis content of PBD and PCOD for the disruption of Rh(I)-ONPs.
Fig. 5
cis (A) and trans (B) 2,6-octadienes as chelating double h2 ligands. In the cis species the double bonds cannot achieve the ideal vertical position that maximizes the back-bonding interaction, and therefore results in a higher energy, more reactive complex. Fig. 6
Again, this is mainly due to strain arising when the olen double bonds cannot achieve a vertical position (that is, perpendicular to the complex plane where the back-bonding interaction is maximized). As can be seen in Fig. 6, the double trans diene can be positioned with both p bonds perpendicular to the Cl–Rh–Cl plane, a preferred geometry, but the double cis conformer (Fig. 6A) is forced to rotate one of the double bonds. In summary, the formation of complexes based on cis-dienes is faster than with the trans ligands because the rst substitution (which determines the rate of association) has less steric impediments in the cis case. The dissociation of cis-dienes is also faster than trans, since the chelating diene in the cis conformer cannot achieve a stable geometry with both double bonds perpendicular to the complex plane, partially breaking the back bonding stabilization and making this a more reactive complex prone to a faster ligand substitution, in accordance to the experimental observations.
The effect of the intermolecular cross-linking of PBD on the thermal properties
Scheme 2
Reaction of 2 with PCy3.
1776 | Chem. Sci., 2016, 7, 1773–1778
Having shown that dilute solutions of polymers containing the 1,5-hexadiene unit have the ability to form organometallic complexes, we sought to investigate whether the bulk polymer could also act as a metal sponge and change its properties when ions are absorbed. It has been shown that the coordination of transition metals to olen containing polymers, such as Pd(II)
This journal is © The Royal Society of Chemistry 2016
View Article Online
Edge Article Table 2
Tg of 95% cis-PBD with varying amount of Rh(I)
Entry
Rh(I) [mol%]
Tg ( C)
1 2 3
0 5 10
101.9 88.5 71.2
a
Open Access Article. Published on 07 January 2016. Downloaded on 11/09/2017 09:12:49. This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.
Chemical Science
a
mol% relative to 1,5-hexadiene units.
with atactic 1,2-polybutadiene, has a signicant effect on the glass transition temperature (Tg) of the composite.43 It has also been reported that Rh containing PBD copolymers (SBS) could be prepared by heating a THF or benzene solution in the presence of RhCl3.44 However, the inuence of different amounts of metal coordination on the Tg of 1,4-polybutadiene has not been studied. Thus, samples of PBD were soaked overnight in THF solutions with varying concentrations of rhodium. The effect of this procedure was studied by analyzing the change in the glass transition temperature of the polymer sample and the shi in the UV-vis absorption spectrum (Fig. S21†). Indeed, as shown in Table 2, the metal led to a signicant change in the polymer's physical properties as expected from the intermolecular crosslinking caused by the metal bridges.
Conclusions Polybutadiene, an abundant starting material, was used for the rst time as the precursor single polymer chain for the formation of organometallic SCNPs. Because these polymers are much larger than those previously used, they may bind many more metal ions per nanoparticle; an important requisite for future applications (e.g. templates for inorganic nanoparticles). Interestingly, the introduction of rhodium ions to neat PBD produced increased glass transition temperatures as a consequence of the restricted local chain mobility caused by the metal bridges. By investigating polymers with differing cis/trans ratios, we have experimentally and computationally shown that there is a stark difference in the rate of metal sequestration and release when polymers with different microstructures are used to form organometallic nanoparticles. Indeed, polymers with higher ratio of cis double bonds bind rhodium ions faster and, when competitive ligands are added, also release them faster. To put this in perspective, under our experimental conditions the rhodium ions could be bound in less than ten seconds by a 300 kDa all cis-PBD and they could be completely removed in about seven minutes by the addition of excess phosphine, while with a high trans-PBD it would take thrice as long. The availability of cis-selective olen metathesis catalysts45,46 may allow the formation of novel 1,5-cyclooctadiene polymers or copolymers with high cis/trans ratios which may be useful for quick binding and release of selected metal ions for future applications of this method. DFT based computational studies supported the experimental results and revealed the mechanistic and steric reasons for the behaviour observed. Taking into account the great importance of PBD cross-linking in the rubber industry, the study of so, reversible, cross-linking
This journal is © The Royal Society of Chemistry 2016
methodologies by metals and the way the polymer structure affects the metal intake is of signicant value.
Acknowledgements The United States – Israel Binational Science Foundation (Grant No. 2014116) is gratefully acknowledged for partial funding of this research.
Notes and references 1 J. P. Rao and K. E. Geckeler, Prog. Polym. Sci., 2011, 36, 887– 913. 2 K. Li and B. Liu, Chem. Soc. Rev., 2014, 43, 6570–6597. 3 M. Elsabahy and K. L. Wooley, J. Polym. Sci., Part A: Polym. Chem., 2012, 50, 1869–1880. 4 M. Elsabahy, G. S. Heo, S.-M. Lim, G. Sun and K. L. Wooley, Chem. Rev., 2015, 115, 10967–11011. 5 S. Mavila, O. Eivgi, I. Berkovich and N. G. Lemcoff, Chem. Rev., 2015, DOI: 10.1021/acs.chemrev.5b00290. 6 C. K. Lyon, A. Prasher, A. M. Hanlon, B. T. Tuten, C. A. Tooley, P. G. Frank and E. B. Berda, Polym. Chem., 2015, 6, 181–197. 7 O. Altintas and C. Barner-Kowollik, Macromol. Rapid Commun., 2012, 33, 958–971. 8 M. Gonzalez-Burgos, A. Latorre-Sanchez and J. A. Pomposo, Chem. Soc. Rev., 2015, 44, 6122–6142. 9 M. Aiertza, I. Odriozola, G. Caba˜ nero, H.-J. Grande and I. Loinaz, Cell. Mol. Life Sci., 2012, 69, 337–346. 10 T. Mes, R. van der Weegen, A. R. A. Palmans and E. W. Meijer, Angew. Chem., Int. Ed., 2011, 50, 5085–5089. 11 P. J. M. Stals, M. A. J. Gillissen, T. F. E. Paffen, T. F. A. de Greef, P. Lindner, E. W. Meijer, A. R. A. Palmans and I. K. Voets, Macromolecules, 2014, 47, 2947–2954. 12 K. A. Williams, A. J. Boydston and C. W. Bielawski, Chem. Soc. Rev., 2007, 36, 729–744. 13 M. Ramanathan and S. B. Darling, Polym. Int., 2013, 62, 1123–1134. 14 G. R. Whittell, M. D. Hager, U. S. Schubert and I. Manners, Nat. Mater., 2011, 10, 176–188. 15 J. Xiang, C.-L. Ho and W.-Y. Wong, Polym. Chem., 2015, 6, 6905–6930. 16 X. Wang and R. McHale, Macromol. Rapid Commun., 2010, 31, 331–350. 17 Y. Liu, T. Pauloehrl, S. I. Presolski, L. Albertazzi, A. R. A. Palmans and E. W. Meijer, J. Am. Chem. Soc., 2015, 137, 13096–13105. 18 M. Artar, E. R. J. Souren, T. Terashima, E. W. Meijer and A. R. A. Palmans, ACS Macro Lett., 2015, 4, 1099–1103. 19 A. Sanchez-Sanchez, A. Arbe, J. Colmenero and J. A. Pomposo, ACS Macro Lett., 2014, 3, 439–443. 20 A. Sanchez-Sanchez, A. Arbe, J. Kohlbrecher, J. Colmenero and J. A. Pomposo, Macromol. Rapid Commun., 2015, 36, 1592–1597. 21 J. Jeong, Y.-J. Lee, B. Kim, B. Kim, K.-S. Jung and H.-J. Paik, Polym. Chem., 2015, 6, 3392–3397.
Chem. Sci., 2016, 7, 1773–1778 | 1777
View Article Online
Open Access Article. Published on 07 January 2016. Downloaded on 11/09/2017 09:12:49. This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.
Chemical Science
22 J. Willenbacher, O. Altintas, P. W. Roesky and C. BarnerKowollik, Macromol. Rapid Commun., 2014, 35, 45–51. 23 J. Willenbacher, O. Altintas, V. Trouillet, N. Knofel, M. J. Monteiro, P. W. Roesky and C. Barner-Kowollik, Polym. Chem., 2015, 6, 4358–4365. 24 S. Mavila, C. E. Diesendruck, S. Linde, L. Amir, R. Shikler and N. G. Lemcoff, Angew. Chem., Int. Ed., 2013, 52, 5767– 5770. 25 S. Mavila, I. Rozenberg and N. G. Lemcoff, Chem. Sci., 2014, 5, 4196–4203. 26 C. Huber, F. Bangerter, W. R. Caseri and C. Weder, J. Am. Chem. Soc., 2001, 123, 3857–3863. 27 J. B. Johnson and T. Rovis, Angew. Chem., Int. Ed., 2008, 47, 840–871. 28 M. A. Bennett, Chem. Rev., 1962, 62, 611–652. 29 P. A. Rupar, G. Cambridge, M. A. Winnik and I. Manners, J. Am. Chem. Soc., 2011, 133, 16947–16957. 30 C. W. Bielawski and R. H. Grubbs, Angew. Chem., Int. Ed., 2000, 39, 2903–2906. 31 B. Allaert, N. Ledoux, N. Dieltiens, H. V. Mierde, C. V. Stevens, P. van Der Voort and F. Verpoort, Catal. Commun., 2008, 9, 1054–1059. 32 H. Brandt, Ullmann's encyclopedia of industrial chemistry, Wiley-VCH, 2012. 33 M. A. Golub, J. Polym. Sci., 1957, 25, 373–377.
1778 | Chem. Sci., 2016, 7, 1773–1778
Edge Article
34 Y. Zhao and D. Truhlar, Theor. Chem. Acc., 2008, 120, 215– 241. 35 S. Grimme, J. Comput. Chem., 2006, 27, 1787–1799. 36 S. Grimme, S. Ehrlich and L. Goerigk, J. Comput. Chem., 2011, 32, 1456–1465. 37 A. V. Marenich, C. J. Cramer and D. G. Truhlar, J. Phys. Chem. B, 2009, 113, 6378–6396. 38 F. Weigend and R. Ahlrichs, Phys. Chem. Chem. Phys., 2005, 7, 3297–3305. 39 G. Aull´ on, G. Ujaque, A. Lled´ os, S. Alvarez and P. Alemany, Inorg. Chem., 1998, 37, 804–813. 40 J. R. Murdoch, J. Chem. Educ., 1981, 58, 32–36. 41 S. Kozuch and S. Shaik, Acc. Chem. Res., 2011, 44, 101–110. 42 S. Kozuch, in Understanding Organometallic Reaction Mechanisms and Catalysis, Wiley-VCH Verlag GmbH & Co. KGaA, 2014, pp. 217–248, DOI: 10.1002/9783527678211.ch9. 43 L. A. Belore and M. P. McCurdie, J. Polym. Sci., Part B: Polym. Phys., 1995, 33, 105–124. 44 L. M. Bronstein and P. M. Valetsky, J. Inorg. Organomet. Polym., 1994, 4, 415–424. 45 B. K. Keitz, A. Fedorov and R. H. Grubbs, J. Am. Chem. Soc., 2012, 134, 2040–2043. 46 M. B. Herbert and R. H. Grubbs, Angew. Chem., Int. Ed., 2015, 54, 5018–5024.
This journal is © The Royal Society of Chemistry 2016
