Home Journals Articles Updates About
CN
image
收藏切换
Switching the O-O Bond Formation Pathways of Ru-pda Water Oxidation Catalyst by Third Coordination Sphere Engineering
收藏切换
PDF
Yingzheng Li1, Shaoqi Zhan2, 3, Lianpeng Tong4, Wenlong Li1, Yilong Zhao1, Ziqi Zhao1, Chang Liu1, Mårten S. G. Ahlquist2, Fusheng Li1, Licheng Sun1, 5, 6
Research. Vol 2021 Article ID 9851231
Less
收藏切换
Research. Vol 2021 Article ID 9851231
• Research Article •
Switching the O-O Bond Formation Pathways of Ru-pda Water Oxidation Catalyst by Third Coordination Sphere Engineering
Full
Yingzheng Li1, Shaoqi Zhan2, 3, Lianpeng Tong4, Wenlong Li1, Yilong Zhao1, Ziqi Zhao1, Chang Liu1, Mårten S. G. Ahlquist2, Fusheng Li1, Licheng Sun1, 5, 6
Affiliations
  • 1State Key Laboratory of Fine Chemicals, Institute of Artificial Photosynthesis, DUT-KTH Joint Education and Research Centre on Molecular Devices, Dalian University of Technology, Dalian 116024, China
  • 2Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry Biotechnology and Health, KTH Royal Institute of Technology, 10691 Stockholm, Sweden
  • 3Department of Chemistry, University of California, Riverside, California 92521, USA
  • 4School of Chemistry and Chemical Engineering/Institute of Clean Energy and Materials, Guangzhou University, No. 230 Wai Huan Xi Road, Higher Education Mega Center, Guangzhou 510006, China
  • 5Department of Chemistry, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, 10044 Stockholm, Sweden
  • 6Center of Artificial Photosynthesis for Solar Fuels, School of Science, Westlake University, 310024 Hangzhou China
About Author:

Authors’ Contributions

F. Li and L. Sun conceived the idea, designed the experiments, and guided the whole project as the corresponding author. Y. Li and S. Zhan contributed equally to this work. Y. Li conceived the idea, conducted the experiments, analyzed the data, and wrote the draft. S. Zhan conducted the DFT calculations and MD simulations, analyzed the data, drew the figures, and composed the manuscript. All authors discussed the results and revised the manuscript.

Published: 2021-04-20 doi: 10.34133/2021/9851231
Outline
收藏切换

Water oxidation is a vital anodic reaction for renewable fuel generation via electrochemical- and photoelectrochemical-driven water splitting or CO2 reduction. Ruthenium complexes, such as Ru-bda family, have been shown as highly efficient water-oxidation catalysts (WOCs), particularly when they undergo a bimolecular O-O bond formation pathway. In this study, a novel Ru(pda)-type (pda2– =1,10-phenanthroline-2,9-dicarboxylate) molecular WOC with 4-vinylpyridine axial ligands was immobilized on the glassy carbon electrode surface by electrochemical polymerization. Electrochemical kinetic studies revealed that this homocoupling polymer catalyzes water oxidation through a bimolecular radical coupling pathway, where interaction between two Ru(pda)–oxyl moieties (I2M) forms the O-O bond. The calculated barrier of the I2M pathway by density-functional theory (DFT) is significantly lower than the barrier of a water nucleophilic attack (WNA) pathway. By using this polymerization strategy, the Ru centers are brought closer in the distance, and the O-O bond formation pathway by the Ru (pda) catalyst is switched from WNA in a homogeneous molecular catalytic system to I2M in the polymerized film, providing some deep insights into the importance of third coordination sphere engineering of the water oxidation catalyst.

Yingzheng Li, Shaoqi Zhan, Lianpeng Tong, Wenlong Li, Yilong Zhao, Ziqi Zhao, Chang Liu, Mårten S. G. Ahlquist, Fusheng Li, Licheng Sun. Switching the O-O Bond Formation Pathways of Ru-pda Water Oxidation Catalyst by Third Coordination Sphere Engineering[J]. Research, 2021 , 2021 (4) : 9851231 - . DOI: 10.34133/2021/9851231
Full Text
Less
1. Introduction
Less
In a typical solar fuel generation device of either electrochemical or photoelectrochemical driven, it often consists of a fuel-forming cathodic half-reaction, such as hydrogen evolution reaction (HER), carbon dioxide reduction reaction (CO2-RR), or nitrogen reduction reaction (NRR), and acquires protons and electrons. Water oxidation (2H2O → O2 + 4H+ + 4e−), on the other hand, plays the unique role of an anodic half-reaction that provides protons and electrons in the overall reaction of fuel generation [1, 2]. This oxidation process is a four-proton four-electron process, which is not only thermodynamically demanding but also kinetically sluggish. Therefore, developing highly efficient water oxidation catalysts to promote this process is imperative for any practical application of HER, CO2-RR, and NRR [3, 4]. To meet this demand, a comprehensive understanding of the catalytic mechanism of the water oxidation reaction is required [5, 6]. Studying a multielectron and multiproton transfer process is a difficult task where molecular catalysts with tailor-designed organic ligands show their advantages of precise structures and relatively easy identification of intermediates involved in the catalytic processes.
Molecular water oxidation catalysts (MWOCs) containing various transition metal cores have therefore drawn significant attention in recent years, because their structures can be readily tailored for systematic structure-reactivity investigation, and their homogeneous catalytic property allows illustration of reaction mechanisms using classic spectroscopic techniques [7–9]. According to the mechanistic study upon numerous MWOCs, the O-O bond formation is often believed to be the rate-determining step (RDS) with the highest energy barrier through the catalytic water oxidation cycle [8, 10]. Two major approaches have been proposed and widely accepted for the essential O-O bond formation step (Scheme 1): the water nucleophilic attack mechanism (WNA) and the oxygen radical coupling mechanism (I2M). Briefly, proton-coupled redox reaction upon the metal center (Mn) of a MWOC yields a high valence metal-oxo (Mn+2=O) species, which can either be attacked by a water molecule in a nucleophilic manner, forming the O-O bond (WNA), or transform to a corresponding metal-oxyl radical (Mn+1-O•) intermediate, immediately affording a peroxo-bond with another adjacent metal-oxyl radical (I2M) [11].
Among Ru-based MWOCs, the family of mononuclear Ru(bda)L2 (bda2- = 2,2’-bipyridine-6,6’-dicarboxylic acid, L = N–cyclic aromatic ligands, e.g., pyridine) complexes have shown outstanding performances in terms of catalytic turnover frequency (TOF) and overpotential [12–14]. It has been demonstrated that water oxidation by these complexes proceeds through a dinuclear I2M pathway and the O-O bond forms via intermolecular coupling of two formal RuV=O moieties with significant RuIV-oxyl character. The hydrophobic nature of the RuV=O species was believed to be one major force to drive the bimolecular radical coupling of Ru (bda) under aqueous catalytic conditions [15, 16]. The series of Ru(pda)L2 (pda2- =1,10-phenanthroline-2,9-dicarboxylate) complexes have remarkably similar coordination geometry with the Ru(bda)L2 catalysts. However, small modification of the backbone organic ligands from bda to a relatively rigid pda dramatically altered the catalytic pathway of the O-O bond formation. Unlike Ru(bda)L2, Ru(pda)L2 prefers a WNA water oxidation pathway when using ceric ammonium nitrate as a sacrificial oxidant [17].
Density-functional theory (DFT) study displayed that, compared with the more flexible bda backbone ligand, the pda ligand increased the steric repulsion among catalysts and thus hindered two formal (pda)RuV=O species from approaching each other [17]. Meanwhile, the (pda)RuV=O species possesses less RuIV-oxyl radical character than the (bda)RuV=O species, and less likely to proceed through the I2M pathway [18]. For the case study of Ru(bda)L2 as well as its analogs illustrate that the catalytic behavior of Ru MWOCs can be affected by even a small change of the primary coordination environment [17]. While the orientation of hydrophobic/hydrophilic substituent groups of MWOCs is a key factor contributing to the different reaction mechanisms [18], modifying the third coordination sphere of the catalyst is another approach to adjust their local environment and thus facilitate the I2M O-O bond formation pathway. These findings encouraged us to manipulate the catalytic water oxidation pathway of Ru(pda) complexes using combined strategies of direct ligand modification and local catalytic environment design.
Herein, two Ru(pda) MWOCs were synthesized with 4-picoline (1) and 4-vinylpyridyl (2) as axial ligands, respectively (Scheme 2). Kinetic studies revealed that electrochemical-driven water oxidation by 1 went through a WNA O-O bond formation pathway, similarly to what was observed in the previous CeIV-driven water oxidation by 1. Nevertheless, complex 2 after immobilized on the glassy carbon (GC) electrode surface by electropolymerization (this homocoupling polymer is denoted as poly-2) in triggers intermolecular O-O radical coupling between metal-oxyl radicals during the electrocatalytic water oxidation process. Our experimental results clearly revealed that the WNA O-O bond formation pathway of single-site MWOCs can switch to the I2M mechanism by changing the local environment, such as shortening the intermolecular distance of catalyst active center, via a homopolymerization method in this case.
2. Results and Discussion
Less
Electrochemical properties of complex 1 in solution and complex 2 that is immobilized on glass carbon (GC) electrode via electropolymerization were investigated by cyclic voltammograms (CVs) and differential pulse voltammograms (DPVs) (as shown in Figure 1). The voltammograms were recorded in a neutral sodium phosphate buffer (pH=7.0, ionic strength=1 M) containing 10% (v/v) CH3CN and glassy carbon (0.071 cm2) electrodes were employed as the working electrodes for the measurement of complex 1 (1/GC) under homogeneous conditions and polymerized complex 2 (poly-2@GC) under heterogeneous conditions. The CV of poly-2@GC exhibited consecutive redox waves of RuIII/II, RuIV/III, and RuV/IV at of 0.81, 1.03, and 1.42 V vs. normal hydrogen electrode (NHE), respectively. These potential values were very close to the corresponding redox potentials of 1/GC (RuIII/II at 0.81 V, RuIV/III at 1.07 V, and RuV/IV at 1.42 V vs. NHE). Both CVs of 1/GC and poly-2@GC showed an onset potential for catalytic current at approximately 1.25 V vs. NHE, where the electrochemical RuV/IV oxidation started, implying that the [RuV=O]+ species triggered O2 evolution in both cases [17]. Considering the backbone ligand of poly-2 tend to draw Ru(pda) units closer in a face-to-face manner of the active sides than the catalyst 1 in a dilute solution, the O-O bond formation via I2M is more likely to occur in the former situation than the latter one, leading to different catalytic water oxidation pathways by 1/GC and poly-2@GC.
To gain some more deep insights into the different O-O bond formation mechanisms between 1/GC and poly-2@GC, three key kinetic indicators were compared between catalytic water oxidation reactions by these two Ru(pda) catalytic systems [19]: (1) reaction orders of the catalyst, (2) reaction orders of the phosphate proton acceptor, and (3) deuterium kinetic isotope effects (KIEsH/D). In a WNA scenario as shown in Scheme 1, the catalytic reaction rate would show a first-order relationship with respect to the catalyst concentration, and increasing concentration of an effective proton acceptor, which promotes the formation of metal-hydroperoxide species (by atom-proton transfer), would facilitate the O-O bond formation. Meanwhile, the catalytic reaction rate usually displays a primary KIEH/D (>2) if the O-H bond cleavage is involved in the rate-determining step (RDS). In an I2M scenario, by contrast, the catalytic reaction rate would exhibit a second-order relationship depending on the catalyst concentration if the binuclear coupling is the RDS. At the same time, the overall rate of water oxidation should be insensitive to the proton acceptor concentration and deuterium substitution because proton transfer is not involved in the RDS.
Kinetic studies for water oxidation by Ru(pda) were firstly conducted by electrochemical methods in the homogeneous system (1/GC). The scan rate-dependent CVs of complex 1 ranging from 50 to 125 mV/s in aqueous sodium phosphate buffer solution (ionic strength=1 M, Na2SO4) are shown in Figure S2a. The position of RuIII/II redox peak at E1/2=0.55 V is independent of the scan rates; however, its peak current (jp) varies linearly to υ1/2 (Figure S2b). Meanwhile, as shown in Figure S2c and S2d, the peak current shows a linear relationship with the concentration of complex 1. This linear relationship is consistent with the Randles-Sevcik relation (Equation S1), indicating a diffusion-controlled redox process [20, 21].
The concentration-dependent catalytic currents were obtained from LSV curves of the 1/GC at various catalyst concentrations (as shown in Figure 2(a)). The apparent electrode reaction order of the catalyst (ρcatjcat) based on the catalytic current can be calculated according to Equation S2, and all slopes approximately to 1 at different applied potentials indicate a single-site electrocatalytic pathway for the 1/GC system [22, 23]. The reaction rate constant of catalyst (kcat) can be obtained according to Equation S3.
The influence (ρcatkcat) of the concentration of complex 1 on its kcat for the water oxidation reaction can be evaluated (Equation S4). Because the reaction orders of the catalytic current (ρcatjcat) and peak current (ρcatjp) on the concentration of complex 1 are both 1; the calculated ρcatkcat is zero, suggesting the independence of the rate constant (kcat) on the concentration of the catalyst in a dilute solution. This is in agreement with the single-site WNA O-O bond formation mechanism [17]. In the case of a heterogeneous system catalyzed by the polymerized Ru(pda) (poly-2@GC) as described above, the high local catalyst concentration in the polymer (poly-2) backbone may shorten the intermolecular distance of Ru units and afford a Ru face-to-face environment, which is beneficial to the dinuclear I2M pathway. The Faraday efficiency of the complex 2 modified electrode (poly-2@GC) was measured (Equation S5) before the kinetic study [19]. As shown in Figure S3, the high Faraday efficiency of poly-2@GC indicated that the cumulative charges passing through the electrochemical system were almost quantitatively depleted in the water oxidation reaction, and therefore, the corresponding catalytic currents of the electrodes were directly applied for the kinetics analysis. Meanwhile, because this electrode reaction for the poly-2@GC system is not a catalyst diffusion-controlled process, the catalytic current density of poly-2@GC can represent the rate constant as well as the turnover frequencies of the catalyst on the surface. By adding the polarization times of electrochemical polymerization, poly-2@GC with various catalyst coverages (Γ) were obtained (Figure S4 and Equation S6). As shown in Figure 2(b), the catalytic current density (jcat) raised with an increasing Γ at the certain applied potential. According to Equation S2, S7, and S8, the reaction orders on the catalyst concentration (ρcatjcat) were calculated to be 1.8~1.91 regarding applied potentials. Furthermore, the ρcatkcat is equal to the ρcatjcat as discussed above. The approximate value 2 of ρcatkcat indicates that the water oxidation reaction catalyzed by poly-2@GC is a second-order reaction on Ru(pda) unit. This observation is consistent with the I2M O-O bond formation pathway. Contrary to the homogeneous catalytic conditions, the kcat of Ru(pda) on the electrode surface shows a strong correlation with the catalyst coverage. The difference of reaction kinetics between the homogenous and heterogeneous conditions confirms the promotion of I2M pathway by reducing intermolecular distance of Ru(pda) moieties on the electrode surface.
As described previously, atom-proton transfer (APT) can dramatically promote the rate of a WNA process [24]. As shown in Figure 2(c), the catalytic currents of 1/GC increase significantly with the increment of phosphate concentration; correspondingly, the electrode reaction order on phosphate (ρpijcat) can be defined and calculated to be 0.3 (Equation S9). The influence of phosphate concentration on the catalytic reaction rate constant ρpikcat was evaluated in a similar manner (Equation S10). As shown in Figure S5a, the peak current (jp) does not change in the presence of various concentrations of phosphate, corresponding to a zero-order kinetic process, and ρpikcat was then calculated as twice of ρpijcat (approximately 0.6). Both ρpijcat and ρpikcat demonstrated the occurrence of phosphate-assisted atom proton transfer process during the electrocatalytic water oxidation reaction by complex 1 under homogeneous conditions. These findings are in line with the proposed WNA O-O bond formation mechanism. The [Pi]-dependent LSV curves of poly-2@GC were recorded by setting the catalyst coverage Γ and other parameters to constants and changing only the concentration of phosphate buffer (Na2SO4 salt is added to maintain the ionic strength to 1.0 M). As shown in Figure 2(d), the reaction rate of poly-2@GC did not rise by increasing the phosphate concentration, corresponding to a zero-order reaction on [Pi]. This suggests no involvement of the atom-proton transfer process in the RDS of water oxidation by poly-2@GC and a I2M pathway for Ru(pda) in the homopolymer form.
The H2O/D2O kinetic isotope effects (KIEsH/D) were further investigated to verify the proposed reaction pathways. As shown in Figure 3(a), KIEsH/D based on the current of the electrode reaction (KIEsH/Djcat) was found to be ~2.3 according to EquationsS11~S14 for 1/GC. The KIEsH/D based on the catalyst rate constant (KIEsH/Dkcat) can also be calculated (Equation S15) to be larger than 4 (Figure S6). The primary kinetic isotope effect on deuterium revealed that the O-H bond cleavage is involved in the O-O bond formation step for complex 1 in the 1/GC system, which is consistent with the WNA O–O bond formation mechanism. Under the heterogeneous catalytic conditions, LSV curves of poly-2@GC in H2O and D2O electrolytes were compared (Figure 3(b)), and no apparent H2O/D2O kinetic isotope effects (KIEsH/D=1.05~1.16) were observed in the potential range from 1.4 V to 1.6 V vs. NHE. These KIEsH/D values suggest that the O-H bond cleavage is not involved in the RDS of water oxidation catalysis by poly-2 and is consistent with the features of a radical coupling (I2M) mechanism.
Based on the experimental data, Figures 3(c) and 3(d) depicted the schemes for the O-O bond formation pathway of water oxidation by complex 1 in a homogeneous system and complex 2 on the electrode surface. As shown in Figure 3(c), the key points for the catalytic pathway of 1/GC are summarized as follows: (i) it is a single site electrochemical reaction in solution, (ii) water oxidation is assisted by the APT process using phosphate as the proton acceptor, and (iii) O-H bond cleavage is involved in the rate-determining O-O bond formation step. All these findings point to water oxidation mechanism via a single Ru site, and the WNA is the dominant O-O bond formation pathway for the 1/GC system. On the other hand, the O-O bond formation pathway of Ru(pda) in poly-2@GC system is proposed as Figure 3(d), which is a binuclear reaction with neither the assistance of phosphate nor the O-H bond cleavage in the rate-determining O-O bond formation step. The O-O radical coupling interaction of two metal-oxyl radicals is the dominant reaction pathway for Ru(pda) in poly-2@GC. The compact intermolecular distance of Ru centers on the electrode surface makes the approaching of two formal (pda)RuV=O species overcome the steric hindrance of the pda ligand and benefits the I2M mechanism reaction rather than the WNA pathway of pristine Ru(pda)-type MWOC.
DFT calculations and molecular dynamics (MD) simulations were performed to further illustrate the reaction pathway of electrochemical driven water oxidation by Ru(pda) in poly-2@GC and properties of the catalyst. Our previous study has shown that the face-to-back configuration is also more favored than the face-to-face configuration of two (pda)RuV=O, due to the hydrophobic pda backbone and the hydrophobic oxo, which results that Ru(pda) complex disfavored the binuclear I2M pathway for homogeneous water oxidation in aqueous solution [18]. The kinetic study on poly-2@GC showed a binuclear I2M reaction pathway. Hence, a supramolecular dimer with face-to-face configuration linked by two propyl groups was revisited as a reasonable initial structural model for the calculations. The calculated redox potentials of RuII−OH2/RuIII−OH2+,RuIII−OH2/RuIV−OH+,RuIV−OH/RuV=O+ are 0.63, 0.73, and 1.22 V at pH 7.0, respectively, which are close to the experimental values (Figure S7). For the key intermediate, [RuV=O]+, MD simulations were performed to understand solvation properties of the oxo in the supramolecular dimer, by using the previously parameterized Ru(pda) model with new partial charge parameters [18]. The H-bond analysis and the radial distribution function (RDF) analysis (Figure S8) of oxygen atoms with water molecules showed that the oxo of the supramolecular dimer is also hydrophobic, the same as that in the molecular Ru(pda) catalysts. The high spin density of the oxo (0.69) from DFT calculations is similar to that of the oxo in molecular Ru(pda) (0.71). With the O-O bond formation step of two (pda)RuV=O species, as shown in the relative free energy profile (Figure 4), the reaction is exergonic with a free energy of -21.16 kcal mol-1 from the prereactive to the product, accompanied by an activation free energy of 0.53 kcal mol-1 from the prereactive to the transition state. The low activation free energy of the (pda)RuV=O coupling process, comparing to the activation free energy of 20.6 kcal mol-1 by WNA pathway [18], supports the experimental results that the Ru(pda) in poly-2@GC catalyzed water oxidation through the I2M mechanism. Once the catalysts overcoming the arrangement step and yielding the prereactive dimer state, the O-O bond is favored to form via coupling of two linked Ru(pda)V=O species. The theoretical evidence implied the possibility of switching the O-O bond formation pathway of the single-site Ru(pda) catalyst from the WNA mechanism to the I2M by manipulating the third coordination sphere, that is the hydrophobic/hydrophilic directionality of the catalysts. For instance, forcing the Ru(pda) catalyst to be a prereactive face-to-face dimer configuration via structural modification is a good strategy for I2M O-O bond formation pathway.
3. Discussion
Less
Recently, several MWOCs based on first-row transition-metal such as Mn, Fe, Co, Ni, and Cu have been developed [25–35]. Although the catalytic mechanism and structure-activity correlation of these WOCs are still not known as much as that of Ru-based WOCs, they are very attractive for opening up the way to develop WOCs based on earth-abundant elements. The O-O bond formation is often the RDS for most WOCs; meanwhile, for most single-site molecular WOCs, the O2 is produced through WNA. How to improve the activity of single-site WNA WOCs attracts attention [11]. Among these, designing binuclear complexes from single-site catalysts has attracted much interest, owning to the potential possibility to obtain complexes that undergo the fast I2M pathway [7]. However, there is no clear experimental evidence that the reported binuclear complexes react via I2M mechanistic pathway. So far, most structure-activity studies of MWOCs have focused on the first and second coordination sphere surrounding the metal center of the catalysts (the chemical structure of the catalyst and proton transfer relay) [36–38]. A significant influence of intermolecular environments (the third coordination sphere) over the O-O bond formation pathway has rarely been reported. This is the first successful study on the switching of a WNA O-O bond formation pathway of a catalyst to a binuclear I2M pathway via manipulating the third coordination sphere of a catalyst. Our approaches might inspire related research on nonnoble metal WOCs-based water splitting devices in the future.
In summary, we studied the O-O bond formation mechanisms of electrochemical water oxidation by Ru(pda) type catalysts. Ru(pda) bearing 4-picoline axial ligands (1/GC) electrochemically catalyzed water oxidation through the water nucleophilic attack (WNA) pathway in a homogeneous system. When the intermolecular distance of the Ru(pda) type catalyst was forced to diminish and configurated in a face-to-face manner, by immobilizing complex 2, a Ru(pda) with 4-vinylpyridyl axial ligands, on glassy carbon electrode surface via electropolymerization (poly-2@GC), the O-O bond formation mechanism switched from WNA to intermolecular radical coupling (I2M). Kinetic studies, including the phosphate concentration and deuterium kinetic isotope effects upon reaction orders, revealed that Ru(pda) in the homocoupling polymer (poly-2) triggers two metal-oxyl radicals coupling interaction during water-oxidation catalysis. This I2M reaction pathway was also supported by DFT calculations. Our results suggest that the WNA O-O bond formation pathway of single-site MWOCs can switch to a binuclear I2M mechanism by manipulating the third coordination spheres of catalysts through structural modification. Our results provide a novel perspective for the design of more advanced water-oxidation catalysts by regulating the catalyst environment without changing the structures of catalysts.
4. Materials and Methods
Less
4.1. Materials and Instrumentation
All commercial chemical reagents were used as received, and the water applied in this work was deionized using the Milli-Q technique. cis-Ru(DMSO)4Cl2 was prepared based on published works [39]. NMR spectra were measured by Bruker Advance 500 spectrometer. Mass spectra were collected on a Finnigan LCQ Advantage MAX mass spectrometer.
4.2. Synthesis
Complex 1 was prepared via a simple one-pot reaction. In a microwave vial (25 mL), 268 mg of pda =1,10-phenanthroline-2,9-dicarboxylate (pda), 1 mL of triethylamine, 484 mg of cis-Ru(DMSO)4Cl2, and 1 mL of 4-methyl pyridine were added together with MeOH (15 mL). The reaction mixture was heated at 100°C for 30 min in a microwave reactor (Biotage Initiator+). The solvent was then removed in the vacuum, and the residues were washed with diethyl ether. The crude product was purified by column chromatography over alumina using DCM : MeOH (100 : 3) as eluents, and complex 1 was obtained as a dark brown solid (567 mg, 38% yield) (38% yield). The 1H-NMR, MS, and elemental analyst were consistent with the work published previously [17]. Complex 2 was prepared via a simple one-pot reaction. In a microwave vial (25 mL), 268 mg of pda =1,10-phenanthroline-2,9-dicarboxylate (pda), 1 mL of triethylamine, 484 mg of cis-Ru(DMSO)4Cl2, and 1 mL of 4-vinyl pyridine were added together with MeOH (15 mL). The reaction mixture was heated at 100°C for 30 min in a microwave reactor (Biotage Initiator+). Methanol was then removed, and the residues were purified by chromatography using an Al2O3 column with DCM : MeOH (10 : 1) as an eluent, yielding a dark red powder as the desired product (468 mg, 30% yield) (30% yield). 1H-NMR (500 MHz, d4-CD3OD, Figure S1): 8.98 (d, j=10, 1H), 8.91 (d, j=10, 1H); 8.62 (s, 2H), 8.47 (d, j=10, 2H), 8.27 (d, j=10, 4H), 7.05 (dd, j=10, 15, 2H), 6.43 (d, j=20, 2H), 5.92 (d, j=7.5, 2H); MS (ESI): calcd for 601.0633 (M+Na+), found m/z+=601.0633. Elemental analysis calculated to 52.57% C, 8.98% N, and 4.222% O.
4.3. Preparation of Poly-2@GC Electrodes
Before electropolymerization, all solutions were deaerated with Ar for 10 min. The preparation of poly-2@GC was based on the previous method [40]. Glassy carbon electrodes (diameter=3 mm) were polarized from 0 to -2.4 V (vs. Ag/AgNO3) in successive cyclic scans in an acetonitrile solution containing 0.1 M TBAPF6 and 0.5 mM complex 2. The resulting poly-2@GC electrodes were rinsed with ethanol and dried under N2 flow.
4.4. Electrochemistry
Electrochemical measurements were carried out with a CH Instruments 660E electrochemical workstation at room temperature. A three-electrode configuration, where a glassy carbon was used as the working electrodes, a Ag/AgCl (3 M KCl) was used as the reference electrode, and a platinum wire was used as the counterelectrode, in a single compartment cell. All the measured potentials were converted to NHE according to a previous report [41]. All GC electrodes were polished with 1  μm alumina powder before all experiments.
4.5. Computational Details
All DFT calculations for the estimation of Gibbs free energies were carried out with the Jaguar 8.3 program package by Schrödinger LLC [42]. Molecular geometries were optimized using Becke’s three-parameter hybrid functional and the LYP correlation functional (B3LYP) with D3 correction of Grimme et al. with the LACVP∗∗ basis set [43–46]. To identify the transition states for O-O bond formation, we searched the potential energy surface by scanning the terminal O-O bond distance [RuV=O•••O=RuV] of the antiferromagnetic open-shell singlet. The thermochemical corrections for estimations of the Gibbs free energy barrier from the prereactive dimers were calculated at the B3LYP-D3/LACVP∗∗ level for both the prereactive dimers and transition state structures. Single-point energy corrections were performed with the B3LYP-D3 functional using the LACV3P∗∗++ basis set augmented with two f-functions on the metal. Based on the gas-phase optimized geometries, the implicit solvation energies were estimated by single-point calculations using the Poisson-Boltzmann reactive field implemented in Jaguar (PBF) in water. The Gibbs free energy was defined by the following equation G=E B3LYP−D3/LACV3P∗∗++2f on Ru+Gsolv+ZPE+H298−TS298+1.9 kcal mol−1 (the value 1.9 kcal mol-1 is a concentration correction to the free energy of solvation, which by default is calculated at 1 M (g) to 1 M (aq) in Jaguar).
MD simulations were performed with the Gromacs 5.1.4 MD software package [47]. A 20 ns MD run was performed in a 44×41×38 Å3 periodic box filled with TIP3P water molecules and two chloride ions to neutralize the charge [48]. In MD simulations, the resulting systems were subject to 100 000 steps of steepest descent minimization. The periodic boundary condition was applied in the simulation. The cutoff radius for the Lennard-Jones and electrostatic interactions were set to be 10 Å. For an accurate evaluation of the long-range Coulomb interactions, Particle Mesh Ewald (PME) summation method is used for electrostatic interactions beyond the cutoff [49]. The system was heated to 300 K in 100 ps by using a v-rescale thermostat for the canonical ensemble (NVT) simulations [50]. During this process, the Linear Constraint Solver (LINCS) algorithm was used to constrain all the bond lengths [51]. The isothermal isobaric ensemble (NPT) was used in the subsequent simulations, with the pressure set to 1 bar in 100 ps, controlled using a v-rescale thermostat and Parrinello-Rahman barostat [52]. Thereafter, the systems were simulated for 20 ns. Three repeated MD simulations with different initial velocities were also performed. A time step of 2.0 fs was used throughout the simulations.
Funding
Less
  • K & A Wallenberg Foundation(KAW 2016.0072)
  • Vetenskapsrådet(2017-00935)
  • National Natural Science Foundation of China(21120102036)
  • Fundamental Research Funds for the Central Universities(DUT19LK16)
References
Less
1
Vilan A., and Cahen D., “Chemical modification of semiconductor surfaces for molecular electronics,” Chemical Reviews, vol. 117, no. 5, pp. 4624–4666, 2017
2
Jiang C., Moniz S. J. A., Wang A., Zhang T., and Tang J., “Photoelectrochemical devices for solar water splitting – materials and challenges,” Chemical Society Reviews, vol. 46, no. 15, pp. 4645–4660, 2017
3
Li W., Li F., Yang H., Wu X., Zhang P., Shan Y., and Sun L., “A bio-inspired coordination polymer as outstanding water oxidation catalyst via second coordination sphere engineering,” Nature Communications, vol. 10, no. 1, p. 5074, 2019
4
Song J., Wei C., Huang Z.-F., Liu C., Zeng L., Wang X., and Xu Z. J., “A review on fundamentals for designing oxygen evolution electrocatalysts,” Chemical Society Reviews, vol. 49, no. 7, pp. 2196–2214, 2020
5
Haschke S., Mader M., Schlicht S., Roberts A. M., Angeles-Boza A. M., Barth J. A. C., and Bachmann J., “Direct oxygen isotope effect identifies the rate-determining step of electrocatalytic OER at an oxidic surface,” Nature Communications, vol. 9, no. 1, p. 4565, 2018
6
Seh Z. W., Kibsgaard J., Dickens C. F., Chorkendorff I., Nørskov J. K., and Jaramillo T. F., “Combining theory and experiment in electrocatalysis: insights into materials design,” Science, vol. 355, no. 6321, article eaad4998, 2017
7
Kärkäs M. D., Verho O., Johnston E. V., and Åkermark B., “Artificial photosynthesis: molecular systems for catalytic water oxidation,” Chemical Reviews, vol. 114, no. 24, pp. 11863–12001, 2014
8
Shaffer D. W., Xie Y., and Concepcion J. J., “O-O bond formation in ruthenium-catalyzed water oxidation: single-site nucleophilic attack vs. O-O radical coupling,” Chemical Society Reviews, vol. 46, no. 20, pp. 6170–6193, 2017
9
Meyer T. J., Sheridan M. V., and Sherman B. D., “Mechanisms of molecular water oxidation in solution and on oxide surfaces,” Chemical Society Reviews, vol. 46, no. 20, pp. 6148–6169, 2017
10
Ghosh S., and Baik M.-H., “The mechanism of O-O bond formation in Tanaka’s water oxidation catalyst,” Angewandte Chemie International Edition, vol. 51, no. 5, pp. 1221–1224, 2012
11
Matheu R., Neudeck S., Meyer F., Sala X., and Llobet A., “Foot of the wave analysis for mechanistic elucidation and benchmarking applications in molecular water oxidation catalysis,” Chem Sus Chem, vol. 9, no. 23, pp. 3361–3369, 2016
12
Duan L., Bozoglian F., Mandal S., Stewart B., Privalov T., Llobet A., and Sun L., “A molecular ruthenium catalyst with water-oxidation activity comparable to that of photosystem II,” Nature Chemistry, vol. 4, no. 5, pp. 418–423, 2012
13
Daniel Q., Wang L., Duan L., Li F., and Sun L., “Tailored design of ruthenium molecular catalysts with 2,2′-bypyridine-6,6′-dicarboxylate and pyrazole based ligands for water oxidation,” Dalton Transactions, vol. 45, no. 37, pp. 14689–14696, 2016
14
Wang L., Duan L., Wang Y., Ahlquist M. S. G., and Sun L., “Highly efficient and robust molecular water oxidation catalysts based on ruthenium complexes,” Chemical Communications, vol. 50, no. 85, pp. 12947–12950, 2014
15
Zhan S., Zou R., and Ahlquist M. S. G., “Dynamics with explicit solvation reveals formation of the prereactive dimer as sole determining factor for the efficiency of Ru(bda)L2 catalysts,” ACS Catalysis, vol. 8, no. 9, pp. 8642–8648, 2018
16
Zhan S., Mårtensson D., Purg M., Kamerlin S. C. L., and Ahlquist M. S. G., “Capturing the role of explicit solvent in the dimerization of RuV(bda) water oxidation catalysts,” Angewandte Chemie International Edition, vol. 56, no. 24, pp. 6962–6965, 2017
17
Tong L., Duan L., Xu Y., Privalov T., and Sun L., “Structural modifications of mononuclear ruthenium complexes: a combined experimental and theoretical study on the kinetics of ruthenium-catalyzed water oxidation,” Angewandte Chemie International Edition, vol. 50, no. 2, pp. 445–449, 2011
18
Zhan S., Zhang B., Sun L., and Ahlquist M. S. G., “Hydrophobic/hydrophilic directionality affects the mechanism of Ru-catalyzed water oxidation reaction,” ACS Catalysis, vol. 10, no. 22, pp. 13364–13370, 2020
19
Zhuo Q., Zhan S., Duan L., Liu C., Wu X., Ahlquist M. S., Li F., and Sun L., “Tuning the O–O bond formation pathways of molecular water oxidation catalysts on electrode surfaces via second coordination sphere engineering,” Chinese Journal of Catalysis, vol. 42, no. 3, pp. 460–469, 2021
20
Bard A. J., and Faulkner L. R. Electrochemical methods: fundamentals and applications, Wiley, New York, vol. 2, 1980
21
Bamford C. H., Tipper C. F. H., and Compton R. G. Electrode Kinetics: Principles and Methodology: Principles and Methodology, Elsevier, 1986
22
Chen Z., Concepcion J. J., Jurss J. W., and Meyer T. J., “Single-site, catalytic water oxidation on oxide surfaces,” Journal of the American Chemical Society, vol. 131, no. 43, pp. 15580–15581, 2009
23
Coggins M. K., Zhang M. T., Chen Z., Song N., and Meyer T. J., “Single-site copper(II) water oxidation electrocatalysis: rate enhancements with HPO42− as a proton acceptor at pH 8,” Angewandte Chemie, International Edition in English, vol. 53, no. 45, pp. 12226–12230, 2014
24
Chen Z., Vannucci A. K., Concepcion J. J., Jurss J. W., and Meyer T. J., “Proton-coupled electron transfer at modified electrodes by multiple pathways,” Proceedings of the National Academy of Sciences, vol. 108, no. 52, pp. E1461–E1469, 2011
25
Du P., and Eisenberg R., “Catalysts made of earth-abundant elements (Co, Ni, Fe) for water splitting: recent progress and future challenges,” Energy & Environmental Science, vol. 5, no. 3, pp. 6012–6021, 2012
26
Poulsen A. K., Rompel A., and McKenzie C. J., “Water oxidation catalyzed by a dinuclear Mn complex: a functional model for the oxygen-evolving center of photosystem II,” Angewandte Chemie International Edition, vol. 44, no. 42, pp. 6916–6920, 2005
27
Dismukes G. C., Brimblecombe R., Felton G. A. N., Pryadun R. S., Sheats J. E., Spiccia L., and Swiegers G. F., “Development of bioinspired Mn4O4−cubane water oxidation catalysts: lessons from photosynthesis,” Accounts of Chemical Research, vol. 42, no. 12, pp. 1935–1943, 2009
28
Brimblecombe R., Swiegers G. F., Dismukes G. C., and Spiccia L., “Sustained water oxidation photocatalysis by a bioinspired manganese cluster,” Angewandte Chemie International Edition, vol. 47, no. 38, pp. 7335–7338, 2008
29
Ma L., Wang Q., Man W.-L., Kwong H.-K., Ko C.-C., and Lau T.-C., “Cerium(IV)-driven water oxidation catalyzed by a manganese(V)–nitrido complex,” Angewandte Chemie International Edition, vol. 54, no. 17, pp. 5246–5249, 2015
30
Lee W.-T., Muñoz S. B., Dickie D. A., and Smith J. M., “Ligand modification transforms a catalase mimic into a water oxidation catalyst,” Angewandte Chemie International Edition, vol. 53, no. 37, pp. 9856–9859, 2014
31
Wasylenko D. J., Ganesamoorthy C., Borau-Garcia J., and Berlinguette C. P., “Electrochemical evidence for catalyticwater oxidation mediated by a high-valent cobalt complex,” Chemical Communications, vol. 47, no. 14, pp. 4249–4251, 2011
32
Dogutan D. K., McGuire R., and Nocera D. G., “Electocatalytic water oxidation by cobalt(III) hangman β-octafluoro corroles,” Journal of the American Chemical Society, vol. 133, no. 24, pp. 9178–9180, 2011
33
Fillol J. L., Codola Z., Garcia-Bosch I., Gomez L., Pla J. J., and Costas M., “Efficient water oxidation catalysts based on readily available iron coordination complexes,” Nature chemistry, vol. 3, no. 10, pp. 807–813, 2011
34
Ellis W. C., McDaniel N. D., Bernhard S., and Collins T. J., “Fast water oxidation using iron,” Journal of the American Chemical Society, vol. 132, no. 32, pp. 10990–10991, 2010
35
Zhang M.-T., Chen Z., Kang P., and Meyer T. J., “Electrocatalytic water oxidation with a copper(II) polypeptide complex,” Journal of the American Chemical Society, vol. 135, no. 6, pp. 2048–2051, 2013
36
Zhang B., and Sun L., “Ru-bda: unique molecular water-oxidation catalysts with distortion induced open site and negatively charged ligands,” Journal of the American Chemical Society, vol. 141, no. 14, pp. 5565–5580, 2019
37
Blakemore J. D., Crabtree R. H., and Brudvig G. W., “Molecular catalysts for water oxidation,” Chemical Reviews, vol. 115, no. 23, pp. 12974–13005, 2015
38
Zeng Q., Lewis F. W., Harwood L. M., and Hartl F., “Role of ligands in catalytic water oxidation by mononuclear ruthenium complexes,” Coordination Chemistry Reviews, vol. 304-305, pp. 88–101, 2015
39
Alessio E., Mestroni G., Nardin G., Attia W. M., Calligaris M., Sava G., and Zorzet S., “Cis- and trans-dihalotetrakis (dimethyl sulfoxide) ruthenium (II) complexes (RuX2(DMSO)4; X = Cl, Br): synthesis, structure, and antitumor activity,” Inorganic Chemistry, vol. 27, no. 23, pp. 4099–4106, 1988
40
Ashford D. L., Sherman B. D., Binstead R. A., Templeton J. L., and Meyer T. J., “Electro-assembly of a chromophore-catalyst bilayer for water oxidation and photocatalytic water splitting,” Angewandte Chemie International Edition, vol. 54, no. 16, pp. 4778–4781, 2015
41
Daniel Q., Duan L., Timmer B. J. J., Chen H., Luo X., Ambre R., Wang Y., Zhang B., Zhang P., Wang L., Li F., Sun J., Ahlquist M., and Sun L., “Water oxidation initiated by in situ dimerization of the molecular Ru(pdc) catalyst,” ACS Catalysis, vol. 8, no. 5, pp. 4375–4382, 2018
42
Bochevarov A. D., Harder E., Hughes T. F., Greenwood J. R., Braden D. A., Philipp D. M., Rinaldo D., Halls M. D., Zhang J., and Friesner R. A., “Jaguar: a high-performance quantum chemistry software program with strengths in life and materials sciences,” International Journal of Quantum Chemistry, vol. 113, no. 18, pp. 2110–2142, 2013
43
Christiansen P. A., Lee Y. S., and Pitzer K. S., “Improved ab initio effective core potentials for molecular calculations,” The Journal of Chemical Physics, vol. 71, no. 11, pp. 4445–4450, 1979
44
Grimme S., Antony J., Ehrlich S., and Krieg H., “A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu,” The Journal of Chemical Physics, vol. 132, no. 15, article 154104, 2010
45
Goerigk L., and Grimme S., “A thorough benchmark of density functional methods for general main group thermochemistry, kinetics, and noncovalent interactions,” Physical Chemistry Chemical Physics, vol. 13, no. 14, pp. 6670–6688, 2011
46
Zhao Y., and Truhlar D. G., “The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals,” Theoretical Chemistry Accounts, vol. 120, no. 1, pp. 215–241, 2008
47
Abraham M. J., Murtola T., Schulz R., Páll S., Smith J. C., Hess B., and Lindahl E., “GROMACS: high performance molecular simulations through multi-level parallelism from laptops to supercomputers,” Software X, vol. 1-2, pp. 19–25, 2015
48
Jorgensen W. L., Chandrasekhar J., Madura J. D., Impey R. W., and Klein M. L., “Comparison of simple potential functions for simulating liquid water,” The Journal of Chemical Physics, vol. 79, no. 2, pp. 926–935, 1983
49
Essmann U., Perera L., Berkowitz M. L., Darden T., Lee H., and Pedersen L. G., “A smooth particle mesh Ewald method,” The Journal of Chemical Physics, vol. 103, no. 19, pp. 8577–8593, 1995
50
Bussi G., Donadio D., and Parrinello M., “Canonical sampling through velocity rescaling,” The Journal of Chemical Physics, vol. 126, no. 1, article 014101, 2007
51
Hess B., Bekker H., Berendsen H. J. C., and Fraaije J. G. E. M., “LINCS: a linear constraint solver for molecular simulations,” Journal of Computational Chemistry, vol. 18, no. 12, pp. 1463–1472, 1997
52
Parrinello M., and Rahman A., “Polymorphic transitions in single crystals: a new molecular dynamics method,” Journal of Applied Physics, vol. 52, no. 12, pp. 7182–7190, 1981
Appendix
Less
Year 2021 volume 2021 Issue 4
PDF
72
40
Cite this Article
BibTeX
Article Info
doi: 10.34133/2021/9851231
  • Receive Date：2021-01-18
  • Online Date：2026-10-27
  • Published：2021-04-20
Article Data
Affiliations
History
  • Received：2021-01-18
  • Accepted：2021-03-19
Funding
K & A Wallenberg Foundation(KAW 2016.0072)
Vetenskapsrådet(2017-00935)
National Natural Science Foundation of China(21120102036)
Fundamental Research Funds for the Central Universities(DUT19LK16)
Affiliations
    1State Key Laboratory of Fine Chemicals, Institute of Artificial Photosynthesis, DUT-KTH Joint Education and Research Centre on Molecular Devices, Dalian University of Technology, Dalian 116024, China
    2Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry Biotechnology and Health, KTH Royal Institute of Technology, 10691 Stockholm, Sweden
    3Department of Chemistry, University of California, Riverside, California 92521, USA
    4School of Chemistry and Chemical Engineering/Institute of Clean Energy and Materials, Guangzhou University, No. 230 Wai Huan Xi Road, Higher Education Mega Center, Guangzhou 510006, China
    5Department of Chemistry, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, 10044 Stockholm, Sweden
    6Center of Artificial Photosynthesis for Solar Fuels, School of Science, Westlake University, 310024 Hangzhou China
References
Share
https://castjournals.cast.org.cn/joweb/research/EN/10.34133/2021/9851231
Share to
QR

Scan QR to access full text

Cite this article
BibTeX
Citations
表12种不同金属材料的力学参数
科
Family		 属数
Number of
genus		 种数
Number of
species		 占总种数比例
Percentage of
total species (%)		 属
Genus		 种数
Number of
species		 占总种数比例
Percentage of total
species (%)
鹅膏菌科Amanitaceae 2 11 5.26 鹅膏菌属 Amanita 10 4.78
小菇科 Mycenaceae 2 12 5.74 丝盖伞属 Inocybe 5 2.39
多孔菌科 Polyporaceae 8 14 6.70 蜡蘑属 Laccaria 5 2.39
红菇科 Russulaceae 3 23 11.00 小皮伞属 Marasmius 6 2.87
小菇属 Mycena 11 5.26
光柄菇属 Pluteus 5 2.39
红菇属 Russula 17 8.13
栓菌属 Trametes 5 2.39
关闭全屏
  • BibTeX
  • EndNote
  • RefWorks
  • TxT
Links
Articles
Latest Articles
Most Read
Collections
Updates
Events
News
Multimedia
About
About Us
Contact
No. 86 Xueyuan South Road, Haidian District, Beijing
100081
010-62662699(Operation and Technical Support)
qkjq@cast.org.cn

Copyright © 2025 China Association for Science and Technology. All rights reserved. For all open access content, the relevant licensing terms apply.
Sponsored by the Office of the Leading Group for Cybersecurity and Informatization of CAST, and supported by Science and Technology Review Publishing House