Ramin
Rabbani‡
a,
Sima
Saeedi‡
a,
Md
Nazimuddin
a,
Héctor
Barbero
ab,
Nathalie
Kyritsakas
c,
Travis A.
White
a and
Eric
Masson
*a
aDepartment of Chemistry and Biochemistry, Ohio University, Athens, Ohio 45701, USA. E-mail: masson@ohio.edu
bGIR MIOMeT, IU CINQUIMA/Química Inorgánica, Facultad de Ciencias, Universidad de Valladolid, Valladolid, E47011, Spain
cMolecular Tectonics Laboratory, University of Strasbourg, UMR UDS-CNRS 7140, Institut le Bel, F-67000 Strasbourg, France
First published on 4th November 2021
A cucurbit[8]uril (CB[8])-secured platinum terpyridyl chloride dimer was used as a photosensitizer and hydrogen-evolving catalyst for the photoreduction of water. Volumes of produced hydrogen were up to 25 and 6 times larger than those obtained with the corresponding free and cucurbit[7]uril-bound platinum monomer, respectively, at equal Pt concentration. The thermodynamics of the proton-coupled electron transfer from the Pt(II)–Pt(II) dimer to the corresponding Pt(II)–Pt(III)–H hydride key intermediate, as quantified by density functional theory, suggest that CB[8] secures the Pt(II)–Pt(II) dimer in a particularly reactive conformation that promotes hydrogen formation.
As Pt complexes are prone to dimerization caused by favorable dz2–dz2 orbital overlap of the Pt centers,15–17 Sakai questioned whether sensitization is triggered by complex 1a as a monomer or dimer. A dimerization constant of 4 (±2) × 103 M−1 and the quadratic dependence of the initial rate of H2 production as a function of the concentration of complex 1a (below 0.4 mM) indicate that either (a) photosensitization is mostly achieved by the 1a2 dimer – via a proton-coupled electron transfer (PCET) to the corresponding Pt(II)–Pt(III)–H hydride (see Fig. 1, pathway A), or less likely (b) the overall rate of H2 formation is governed by the coupling of two Pt hydrides, a necessarily bimolecular process (see Fig. 1, pathway B).3
Fig. 1 Two possible pathways for the Pt-catalyzed photoreduction of water.13 Pt(tpy) chloride complexes 1a3 and 1b and amidate-bridged dimer 214 as photosensitizers and/or H2-evolving catalysts. |
Sakai also showed that amidate-bridged Pt dimers (see complex 2, for example, in Fig. 1) are efficient H2-evolving catalysts14,18,19 in a three-component system consisting of the tris(bipyridine)ruthenium(II) cation as photosensitizer and methylviologen as electron relay. Furthermore, Yamashita and coworkers showed that grafting complex 1a onto solid supports (layered niobate20 K4Nb6O17 or mesoporous silica21,22) improves its catalytic activity by packing the complexes into arrangements with short Pt–Pt distances. This led us to consider the cucurbit[8]uril (CB[8])-secured Pt tpy dimers prepared recently in our group as a possible new family of single component catalysts for water reduction.23–26 In the past few years, we have shown that Pt tpy complexes substituted at the 4′-position with CB[8]-binding aryl units form 2:1 head-to-head “stacked” assemblies, with both Pt tpy units sitting on top of each other at the same CB[8] portal, i.e. with two positively charged guests interacting with the same CB[8] portal while leaving the opposite one void of any host/guest interaction (see Fig. 2). The driving forces of this new recognition motif are dispersive interactions between the planar tpy ligands, and the very Pt–Pt interactions through dz2–dz2 orbital overlap that are suspected to enhance hydrogen formation in water reduction processes. We remind the reader that the overlap of the pair of filled dz2 orbitals is only favourable because of symmetry-allowed mixing with the Pt 6pz orbitals. This returns four orbitals: a strongly bonding filled dσ, a weakly anti-bonding filled dσ* (HOMO), a weakly bonding pσ LUMO and a strongly anti-bonding pσ*. The stabilization of the dσ* HOMO confers a net energy gain to the interaction.27,28
Fig. 2 1H NMR spectra of complexes (a) CB[7]·1b, (b) 1b and (c) CB[8]·1b2 in pure D2O; chemical shifts in ppm; see Fig. 1 for numbering. (d) Single-crystal X-ray diffraction structure of complex CB[8]·1b2 showing the two conformations within the unit cell. |
The X-ray crystal structure of assembly CB[8]·1b2 showed two different conformations within the unit cell (see Fig. 2d). This arrangement allows an optimal interaction between the carbonyl rim of one CB[8] unit with the positive Pt dimer and the outer wall of the other CB[8] macrocycle. Both CB[8]·1b2 conformations present slip-stacked Pt tpy units, with geometries primarily governed by the slip angle between the Pt tpy planes29–32 (i.e. the angle between the Pt–Pt axis and the normal line to one tpy NNN plane), and the Caryl–Cpyridyl torsional angles ϕ1 and ϕ2 between the CB[8]-binding tolyl group and the tpy scaffolds. In one assembly, both set of angles rotate in the same direction (θ = +17°, ϕ1 = +21° and ϕ2 = +22°), and in the other, the directions are opposite (θ = +11°, ϕ1 = −24° and ϕ2 = −32°, see Fig. 2). This results in two different Cl–Pt–Pt–Cl dihedral angles (+3° and −21°, respectively) and two Pt–Pt distances (3.63 Å and 3.44 Å, respectively). Complex 1b and its CB[n]-bound assemblies display typical absorption features (see Fig. S4a†). As shown by Che,33,34 McMillin35–39 and Gray,40,41 absorption bands between 280 and 350 nm, typically with vibronic structure and extinction coefficients of approximately 104 M−1 cm−1, are associated with intra-ligand 1[π(tpy) → π*(tpy)] transitions (the π*(tpy) orbital being better described as a 6pz(Pt)/π*(tpy) hybrid). Absorption bands above 350 nm are broad with lower extinction coefficients; they correspond to 1[5d(Pt) → π*(tpy)] metal-to-ligand charge transfers (1MLCT).
Photolysis experiments were carried out in MES buffer (0.10 M, pH 5) in the presence of EDTA as sacrificial electron donor (30 mM). Solutions (6.0 mL in 23 mL vials equipped with PTFE septa) were irradiated with 10 royal blue LEDs (λ = 447 ± 10 nm; 250 mW/LED). Aliquots were removed from the headspace and analyzed by gas chromatography as a function of time (see ESI† section for details). Irradiation was intentionally carried out using the least energetic wavelengths required to trigger the photoreduction. The fact that solar irradiance at sea level peaks at approximately 460 nm (ref. 42) also guided our decision.
Fig. 3 shows the poor performance of free complex 1b at promoting water reduction (up to 0.10 mL H2/10 mL solution at 0.50 mM Pt after 4 h, and a TON decreasing steadily from 2.5 to 0.1 as Pt concentration increases from 0.20 to 2.0 mM). CB[7] encapsulation of the tolyl tpy substituent does not significantly enhance performance, with up to 0.20 mL H2/10 mL solution at 1.0 mM Pt. Increasing the catalyst load above 0.50 mM and 1.0 mM is detrimental to H2 production, in the case of complexes 1b and CB[7]·1b, respectively (with CB[7]·1b, The TON decreases again from 2.6 to 0.4 as Pt concentration increases from 0.20 to 2.0 mM). To the contrary, H2 production is increased by 25- and 6-fold when the CB[8]-secured Pt dimer is used as catalyst, compared to complexes 1b and CB[7]·1b (1.1 mL/10 mL solution at 2.0 mM Pt, compared to 0.04 and 0.17 mL, respectively). A pronounced increase in hydrogen production is observed as the catalytic load increases (see Fig. 3e and f in particular). The TON remains rather steady with Pt concentrations up to 1.0 mM (4.1, 3.7 and 3.5 at 0.20, 0.50 and 1.0 mM, respectively); a mild TON decrease to 2.4 is measured for the 2.0 mM Pt series after 4 h of irradiation. To rule out catalysis by colloidal Pt particles (if present at all), photolysis was carried out again with 1.0 mM Pt as free, CB[7] and CB[8]-bound complexes, in the presence of a few mercury drops. The latter would extract Pt(0) colloids and hamper hydrogen production if those were present and had catalytic activity.43 No significant changes in H2 formation were observed, hence ruling out this possibility. In an earlier study, Yamashita and coworkers indeed showed the absence of Pt(0) formation using X-ray absorption fine-structure spectroscopy (XFAS),44 after irradiation of complex 1a (4.0 mM) and EDTA (30 mM) in a sodium acetate buffer at pH 5.
Photolysis of a 0.50 mM CB[8]·1b2 solution was also monitored by electrospray ionization mass spectrometry, with 1.5 mM EDTA as sacrificial electron donor in water. Before irradiation, assembly CB[8]·1b2, its binary complex CB[8]·1b fragment, its free guest fragment, and an EDTA bridged Pt(tpy) dimer were detected (see Fig. S3†). As observed by Sakai,3 Pt/EDTA adducts underwent extensive decomposition upon irradiation; a noteworthy degradation product was the Pt(tpy) hydride cation (see Fig. S3 and Table S1†).
A linear correlation between the initial rate of H2 formation (up to 40 min) and the Pt concentration might have been detected when complex CB[8]·1b2 is used as the photocatalyst (see Fig. 3f). This result, in light of the possible quadratic relationship observed by Sakai with free complex 1a,3 could suggest that the enhanced catalytic activity of the CB[8] assembly is caused by the macrocycle securing two Pt tpy complexes together, regardless of Pt concentration. However, free Pt complexes are prone to dimerization anyway, even at low concentration, but hydrogen production in the absence of CB[8] remains low (see Fig. 3). As mentioned above, Sakai and coworkers extracted a dimerization constant of 4 (±2) × 103 M−1 for complex 1a. UV-Vis absorption spectra of Pt complex 1b returned concentration-dependent extinction coefficients (see Fig. S4b†). The absorbances between 360 and 480 nm at various concentrations were fit with a simple dimerization model and returned a dimerization constant of 5.7 (±0.2) × 103 M−1 in the MES/EDTA buffer, in excellent agreement with Sakai's measurements. The low solubility of free Pt complex 1b and free CB[8] in aqueous medium precluded the extraction of the thermodynamic parameters of the dimerization and CB[8] binding processes by isothermal titration calorimetry (ITC). However, the more soluble Pt thiolate 1c obtained after exchange of the chloride ligand with L-cysteine afforded a dimerization constant KPt–Pt of 4.0 (±0.4) × 104 by ITC in water, in very good agreement with the 7-fold lower constant measured for complex 1b in the high ionic strength buffer environment. A ternary binding constant β for complex 1c towards CB[8] of 1.9 (±0.6) × 1013 M−2 (see equilibrium 1) was extracted. In other terms, dimer 1c2 (and most likely dimer 1b2) can be considered as a standalone guest forming a very tight 1:1 complex with CB[8], with a binding affinity K′ of 5 (±2) × 108 M−1, obtained from eqn (2), in excellent agreement with one of our recent studies26 on a parent system.
1c + 1c + CB[8] ⇄ CB[8]·1c2 | (1) |
(2) |
The dimerization constants of complexes 1b and 1c return free energy terms of −5.1 (±0.1) and −6.3 (±0.1) kcal mol−1, respectively, in excellent agreement with the typical strength of Pt–Pt interactions.45 They also indicate that approximately 52, 66, 74 and 81% of Pt complex 1b is already present as a dimer at concentrations of 0.2, 0.5, 1.0 and 2.0 mM, respectively, even in the absence of CB[8], in the MES/EDTA buffer. While a simple dimerization process is sufficient to fit enthalpograms and changes in UV-Vis absorption, further aggregation into trimers or larger oligomers cannot be ruled out even if, in our case, considering those processes leads to over-parametrization.
To identify the underlying cause of the catalytic enhancement in the presence of CB[8], we used density functional theory (DFT) to calculate the free energy of Pt(III)–H hydride formation from the corresponding Pt(II) complex (see Fig. 1), using Pt complex 1b as (1) a free monomer, (2) a CB[7]-bound monomer, (3) a dimer in the absence of CB[8], and (4) a CB[8]-secured dimer. The sheer complexity of the PCET process prevents us from calculating free energies of activation: it would necessarily involve proton transfer to the Pt center, most likely from protonated MES at pH 5, and the oxidation of a carboxylate group of EDTA to the corresponding neutral radical, both mechanisms operating with various degrees of concertedness. As these processes are unlikely to be affected by the presence of CB[n], we postulate that more thermodynamically favorable PCET processes would result in lower activation energies and faster hydrogen production.
Sakai notes that Pt–Pt distances obtained by DFT optimization (B3LYP/LANL2DZ)46–49 are often significantly longer than those measured by X-ray crystallography (0.3 Å or longer).11 Here, we can report that optimization in the gas phase with Grimme's newly developed B97-3c functional50 and def2-mTZVP basis sets51,52 returns accurate Pt–Pt distances and geometries highly similar to those measured by X-ray crystallography. The two CB[8]·1b2 arrangements (see Fig. 2d) were reoptimized at this level of theory, and afforded Pt–Pt distances of 3.71 Å (vs. 3.63 Å in the X-ray crystal structure) and 3.35 Å (vs. 3.44 Å), respectively. Cl–Pt–Pt–Cl dihedral angles are +3° and −21°, respectively, in perfect agreement with the X-ray crystal structure. Grimme also showed that for a set of reactions involving transition metals, the B97-3c functional not only returns realistic geometries, but also more accurate free energies compared to the B3LYP functional, even when the latter is supplemented with a dispersive term.50 The enthalpic and entropic contributions at 25 °C, as well as the free energies of solvation, were calculated using the tight-binding semi-empirical method GFN2-xTB currently developed by Grimme and coworkers53–55 (see ESI† section for details; GFN2-xTB reaches the accuracy of the B97-3c functional when calculating the rotational and vibrational contributions to reaction free energies, including for transition metal complexes, at a fraction of computation time).53,54 The HOMO in the arrangement with the short Pt–Pt distance is the expected dz2–dz2 antibonding orbital (see Fig. S6†). To the contrary, the HOMO in the other arrangement is a dxz–dxz antibonding orbital (see Fig. S7a†); dz2–dz2 antibonding is relegated to the HOMO-5 energy level (see Fig. S7b†). Remarkably, both assemblies have similar stabilities, with just a slight preference for the short Pt–Pt distance (0.4 kcal mol−1), thereby indicating an even competition between the Coulombic repulsion of the positively charged Pt centres and the dz2–dz2 orbital overlap. This is also consistent with the presence of both arrangements in the unit cell of the X-ray crystal structure.
Free Pt complex 1b, its CB[7]-bound assembly and dimer 1b2 were optimized at the same level of theory. The Cl–Pt–Pt–Cl dihedral angle in the most stable conformation of dimer 1b2 is 124°, and the Pt–Pt distance 4.99 Å (see Section 11 in the ESI† section). Dimerization is favoured by dispersive interactions between the planar ligands, while Pt–Pt orbital overlap is insignificant.
We then optimized the corresponding open-shell Pt(III)–H hydrides at the same level of theory. As noted by Sakai, the Pt–Pt bond is reinforced compared to Pt(II)–Pt(II) interactions, and can be seen as a three-center two-electon Pt(II)–Pt(III)–H system11,56–58 with a Pt–Pt bond order of 0.5 (further oxidation to a d7–d7 Pt(III)–Pt(III) system would afford a σ bond between both metals, as the dz2–dz2 anti-bonding orbital would be empty). A unique CB[8]·1b2–H assembly is obtained with a pronounced shortening of the Pt–Pt bond (from 3.35 or 3.71 Å to 2.99 Å, see Fig. 4a). The 0.36 Å shortening is very consistent with that observed after oxidation of tetranuclear Pt(II) complexes into mixed-valence Pt(II)2Pt(III)2 complexes.59 The Pt–Pt axis also crosses the tpy planes almost perpendicularly (the slip angle between both tpy planes is 2° vs. 11° and 17° in the pair of CB[8]·1b2 assemblies, see Fig. 2d). The SOMO of assembly CB[8]·1b2 H is indeed a dz2(Pt)-dz2(Pt)-s(H) hybrid (see Fig. 4b). A similar 2.96 Å Pt–Pt distance in dimer 1b2–H (with a Cl–Pt–Pt–Cl dihedral angle of 136°) was also obtained (see ESI† section).
The following observations can be made from the electronic contributions (ΔE) and free energies (ΔG) of the PCET processes (see Table 1): (1) the latter is more favourable by 4.8 kcal mol−1 when applied to a Pt dimer compared to a Pt monomer, even in the absence of CB[8], in strong support of the mechanism proposed by Sakai and coworkers (ΔG = −18.0 vs. −13.2 kcal mol−1). (2) CB[n] encapsulation does not provide any electrostatic contribution to the PCET process, as it is separated from the Pt centers by the tpy ligands (ΔG = −13.6 vs. −13.2 kcal mol−1 in the presence and absence of CB[7], respectively). (3) Aggregation of free Pt complex 1b must severely hamper the catalytic process, as calculations indicate a more favorable PCET process for dimer 1b2 compared to complex CB[7]·1b (ΔG = −18.0 vs. −13.6 kcal mol−1), whereas hydrogen production is even poorer. (4) Most importantly, in addition to efficiently deaggregating Pt complex 1b (and likely preventing the formation of higher order oligomers), CB[8] secures the Pt dimer into a conformation that facilitates the PCET process. If dimerization of Pt complex 1b were quantitative in the presence and absence of CB[8], a 2.2 kcal mol−1 gain for the PCET process would be obtained upon addition of the macrocycle (ΔG = −20.2 vs. −18.0 kcal mol−1 for dimer 1b2, in the presence and absence of CB[8], respectively, see Table 1). However, dimerization is only partial in the absence of CB[8]; to allow a proper comparison, the calculated free energy ΔG of the PCET process must be corrected by −RTlnx, where x is the fraction of Pt complex 1b present as a dimer (74% at 1.0 mM, for example, i.e. a correction of −0.17 kcal mol−1). The energy gain upon addition of CB[8] therefore ranges from 1.8 to 2.1 kcal mol−1 depending on Pt concentration. Also, removing CB[8] from complexes CB[8]·1b2 and CB[8]·1b2–H and carrying out single point energy calculations returns an electronic contribution ΔE of −31.4 kcal mol−1 for the PCET process, a 4.6 kcal mol−1 improvement compared to dimer 1b2 in its most stable conformation in the absence of CB[8] (ΔE = −26.8 kcal mol−1). In addition to being computed using adequate functionals and semi-empirical methods, we note that the quality of these calculated energy gains must be further enhanced by error compensation, as they represent relative stabilities of very similar assemblies. One can thus conclude that CB[8] encapsulation destabilizes dimer 1b2 by at least 1.8 kcal mol−1, enhances its reactivity, and promotes hydrogen production.
Footnotes |
† Electronic supplementary information (ESI) available: Preparation and characterization of Pt/CB[n] assemblies; photocatalytic H2 production details; chemical actinometry; mass spectrometry analysis during photolysis; UV-Vis absorption analysis of Pt/CB[n] assemblies; ITC experiments; single-crystal X-ray crystallographic studies and computational details. CCDC 2094663. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/d1sc03743a |
‡ These authors contributed equally. |
This journal is © The Royal Society of Chemistry 2021 |