Guozhen Xu,
Jian Zhang*,
Shengping Wang,
Yujun Zhao* and
Xinbin Ma
Key Laboratory for Green Chemical Technology of Ministry of Education, Collaborative Innovation Centre of Chemical Science and Engineering, School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, China. E-mail: zhangjian1985@tju.edu.cn; yujunzhao@tju.edu.cn
First published on 18th May 2016
This work reported a well-fabricated PtSn/SiO2 catalyst prepared by a modified two-step sol–gel (MTSG) method. The adoption of strong electrostatic adsorption (SEA) and sol–gel (SG) method could give both a significantly high Pt dispersion and a large amount of Lewis acid sites. The homogenous distribution of Sn species could not only provide more opportunity for the Pt precursor to disperse on the support by forming Pt–(O–Sn
)y2−y analogous species, but also enhance the synergy between Pt and Sn species. Consequently, an excellent activity was achieved in the hydrogenation of acetic acid (AcOH) with a conversion of 100% and ethanol (EtOH) selectivity of 93%. Investigations on the effect of Sn/Pt molar ratio showed that the dispersion of Pt decreases obviously with the increasing Sn/Pt ratio due to the geometric or electronic effects caused by SnOx species. A balancing effect between Pt active sites and Lewis acid sites was found to be responsible for the superior catalytic performance in AcOH hydrogenation. Moreover, a parallel reaction path model was proposed for the hydrogenation of AcOH over PtSn/SiO2 catalyst, in which ethanol and ethyl acetate (AcOEt) are formed competitively through the adsorbed ethoxy intermediate.
O bond activation. Meanwhile, ethanol is an attractive product that has been widely used as a fuel additive, solvent and antiseptic.5,6
AcOH hydrogenation has been studied since the 1990's. It was found that the direct catalytic hydrogenation of AcOH requires the presence of both metal and metal oxide phases.7 Further researches showed that the adsorption of AcOH on the oxide phases plays a major role in determining the hydro-conversion kinetics and the role of metal is to provide mobile, activated hydrogen atoms. However, the mechanism of AcOH hydrogenation is still on debate. Two types of mechanism involving Mars–Van Krevelen mechanism and Langmuir–Hinshelwood mechanism have been proposed. The difference between the two mechanisms is the definition of active site for AcOH activation. The former attributes the active sites to the surface vacancies, while in the later it is assigned to the Lewis acid sites/metal cations on the catalyst surface.1,2,8,9
Recently, bimetallic supported noble-based catalysts such as Pt–Sn, Pt–Fe, Ru–Fe, Pd–Re etc., are reported to be effective for AcOH hydrogenation.7,9–12 Among these catalysts, PtSn supported bimetallic catalyst acquires a lot of attention since it exhibits good activation for the hydrogenation of AcOH to ethanol. Nevertheless, most of the PtSn catalysts are generally prepared by impregnation (IM) method nowadays, which is quite limited to control the structural parameters. Consequently, the hydrogenation process still suffers from a variety of issues such as lower ethanol selectivity, harsh reaction conditions and high catalyst costs.12,13
PtSn bimetallic catalysts have been widely used in many catalytic processes. Various ways of preparing PtSn bimetallic catalyst have been developed to tailor the electronic and geometric structures for different applications.14–16 Up to now, supported PtSn bimetallic catalysts could be prepared by co-reduction method, co-precipitation method, sol–gel (SG) method, hydrothermal approach, carbonyl route, and a controlled surface reaction method.13,17,18 SG method has been reported to be one of the promising technologies to prepare PtSn catalysts with well-defined physical properties. It increased the homogeneity and purity of the solid, with better microstructural control of the support, well-defined pore size distribution, and higher thermal stability of the catalysts. Two types of sol–gel methods have been reported in preparing PtSn/Al2O3 catalyst for n-heptane reforming.19 One is the one-step sol–gel (OSG) method, which starts from both the support and metal precursors and then go all way to the final catalysts in one batch. The other is two-step sol–gel (TSG) method, where platinum is uploaded by IM method after tin is introduced into alumina support by SG method. Both of them could give catalysts with higher activity and stability.
On another hand, high dispersion of metal nanoparticles is thought to be crucial in the design of supported noble metal catalysts, in order to achieve high performance and minimize the cost. It was reported that Pt-based catalysts prepared by strong electrostatic adsorption (SEA) method exhibited higher dispersion of Pt than those prepared by IM method.20 By controlling the pH value of the impregnation solution, the hydroxyl groups that bonding on oxide surface become protonated or deprotonated, thus become positively charged below a characteristic pH value or negatively charged above the characteristic pH value. Therefore, the metal complex could deposit onto the support surface by SEA method, resulting in a relative higher metal dispersion.21–23
Herein, we use a modified two-step sol–gel (MTSG) method for the synthesis of highly dispersed and homogeneously distributed PtSn catalyst, in which Sn/SiO2 was firstly synthesized by SG method and then Pt was uploaded by SEA method sequentially. The effects of the synthesis methods on the surface structure and catalytic properties of AcOH hydrogenation were studied. It shows that this MTSG strategy could enhance both the catalytic activity and selectivity for the hydrogenation of AcOH into ethanol, due to the enhanced synergy between Pt and Sn species. Moreover, the influence of Sn/Pt molar ratio was investigated and a reaction model for the formation of ethanol and AcOEt was proposed.
:
C2H5OH
:
H2O was 1
:
3.85
:
10.20. TEOS was mixed with ethanol and an alcoholic solution of SnCl2·2H2O. After that, deionized water was added dropwise into the solution with vigorously stirring at room temperature. The pH of the solution was kept to be ca. 4 during the water dropping process. After that, the clear sol was allowed to gel in air for 48 h. Subsequently, the obtained gel was aged for 24 h at 60 °C. The resulting solids were filtered, washed with ethanol until negative chloride test, and then dried in oven at 105 °C. The obtained xerogel was then calcined in air at 500 °C for 4 h. For comparison, Sn precursor was also uploaded by IM method with controlled concentration of SnCl2 solution (1% wt/wt Sn loading), and the acquired substrate was denoted as Sn/SiO2-IM.
Supported Pt catalysts were obtained by depositing the metal precursor onto the support surface by SEA method. Platinum(II) tetra-ammine chloride ([Pt(NH3)4]Cl2, 99%, Sinopharm Chemical Reagent Co., Ltd) was used as the metal precursor. The SiO2 or Sn/SiO2-SG support was firstly suspended in NH4OH solution (aq., pH = 10.6) under stirring for 0.5 h. The solution of [Pt(NH3)4]Cl2 was then added with controlled concentration to obtain 1% wt/wt Pt loading on the support. The suspension was kept under stirring for 10 h and then was separated by filtration under vacuum. The solid was rinsed with water, and then was dried at 105 °C for 5 h. These catalysts were named as PtSnx/SiO2-MTSG and x represents for the theoretical Sn/Pt molar ratio. For comparison, other two types of PtSn/SiO2 catalysts were also prepared. One was obtained by loading Pt onto Sn/SiO2-IM by SEA method, which was denoted as PtSn1.6/SiO2-MIM. The other was prepared by impregnation of [Pt(NH3)4]Cl2 onto Sn/SiO2-SG, and was denoted as PtSn1.6/SiO2-TSG.
The turn over frequency (TOF) of the reaction was measured under the conditions (270 °C, 2.6 MPa, H2/HOAc (mol mol−1) = 20) where AcOH conversion was 60–70%. The TOFPt value was based on exposed Pt atoms, indicating the moles of AcOH converted per hour by per mole of exposed Pt atoms at the catalyst surface.
The in situ DRIFTs experiments were carried out on a Thermo Scientific Nicolet 6700 spectrometer (32 scans, 4 cm−1), equipped with a Harrick Scientific DRIFT cell and a mercury–cadmium–telluride (MCT) detector cooled by liquid N2. All the samples were pretreated at 100 °C under Ar (30 mL min−1) for 1 h, and then FTIR spectra was collected using KBr as background spectrum.
Temperature-programmed reduction (H2-TPR) tests were carried out on a Micromeritics AutoChem 2920 TPD/TPR apparatus using ca. 0.15 g of catalyst. The catalyst was prior treated at 200 °C for 2 h in flowing He. A flow rate of 30 mL min−1 of 10 vol% H2/Ar was used for the reduction with the temperature range from 50 °C to 900 °C at a rate of 10 °C min−1. A thermal conductivity detector (TCD) was employed to determine the amount of hydrogen consumed.
The metal dispersion of the catalysts was determined according to literature12 by pulsed hydrogen titration of oxygen pre-covered sites using a Micromeritics 2920 TPD/TPR apparatus. Prior to the analysis, catalysts (0.15 g) were reduced in flowing 10 vol% H2/Ar at 300 °C for 2 h, followed by purging with He to remove any physical adsorbed hydrogen. After cooling to 50 °C in flowing He, the catalysts were exposed to 10 vol% O2/He for 0.5 h followed by purging with Ar. Titration with pulses of 10 vol% H2/Ar was then employed until no further H2 uptake occurred.
The in situ FTIR of chemisorbed CO on reduced catalysts were collected on a Nicolet 6700 (32 scans, 4 cm−1) spectrometer using a stainless steel cell connected to a gas dosing and evacuation system. The sample powders were prior reduced at 300 °C in H2 for 4 h and then cooled to room temperature. After that, the powder was pressed into a self-supporting disk (ca. 20 mg) under the pressure of 8 MPa and then was introduced into the IR cell. The disk was re-reduced in 10 vol% H2/Ar at 300 °C for 0.5 h and then cooled to 30 °C under a high-purity He flow for collecting the background spectrum. Afterward, the sample was exposed to the flow of CO at 30 °C for 0.5 h. FTIR spectra were collected after purging He for 0.5 h and referenced to the background spectrum of reduced catalysts. For normalization, the intensity of the spectrum was enlarged by the same multiple as the disk weight used to be enlarged to 1 g.
Fourier transform infrared (FTIR) spectra of pyridine adsorption were also collected on a Thermo Scientific Nicolet 6700 (32 scans, 4 cm−1). The pretreatment for the sample was same as that for CO adsorption. After the pretreatment, the disk was re-reduced in 10 vol% H2–Ar at 300 °C for 0.5 h and then was evacuated under reduced pressure at the same temperature. After cooling to 50 °C, the samples were saturated with pyridine vapor for 0.5 h and then evacuated at 150 °C for 0.5 h, respectively. Spectra were recorded at the evacuation temperature in the range of 4000–650 cm−1. The amount of the Lewis acid sites in samples was calculated by peak deconvolution of characteristic band at ca. 1450 cm−1.24
X-ray photoelectron spectrometer (XPS) was carried out on a Thermo ESCALAB 250Xi system, using a monochromatic Al-Kα (1486.6 eV) source operating at 10.8 mA × 15 kV and a pass energy of 50 eV under an operating pressure of 7.1 × 10−5/7.1 × 10−7 Pa. The energy resolution is 1 eV. Catalysts were mounted on a special-designed holders. After reduction in situ in the rapid sampling chamber and pre-treatment chamber, samples were transferred into analysis chamber by vacuumed gradually. Possible deviations caused by electric charge on the samples were corrected through taking the Si 2p line as an internal standard at 103.4 eV. Peak fitting was achieved using XPSPeak 41 software.
)y2−y analogous species could be formed on the catalyst surface.25,26 Consequently, the interaction between Pt and Sn species could be enhanced when Pt is uploaded by SEA method. However, such interaction should be weakened when Pt is uploaded by IM method because of the lack of such interaction between Pt and Sn species.
![]() | ||
| Fig. 1 TEM images of the reduced catalysts. (A) Pt/SiO2-MTSG, (B) PtSn1.6/SiO2-MTSG, (C) PtSn1.6/SiO2-MIM, (D) PtSn1.6/SiO2-TSG. | ||
)y2−y analogous species might be formed by the interaction between [Pt(NH4)3]2+ and SiOH groups.28,29
)y2−y species.31 In the case of Sn/SiO2-SG catalyst, the reduction profile is characterized by a broad peak with a maximal temperature (TM) at 503 °C, and the reduction peak begins at 300 °C and ends at 700 °C, which is due to the reduction of SnO2 species.14
H2-TPR profiles of the PtSn/SiO2 catalysts display features related to the preparation methods. Comparing with the Pt/SiO2-MTSG catalyst, PtSn1.6/SiO2-MTSG catalyst presents a overlapped reduction peak at 282 °C, which could be attributed to the co-reduction of Pt species with weak interaction with the support and Pt–(O–Sn
)y2−y analogous species. A tiny reduction peak at 430 °C could be ascribed to the reduction of SnO2.14 For the PtSn1.6/SiO2-MIM catalyst, two reduction peaks are centered at 272 °C and 360 °C, the former could be assigned to the co-reduction of Pt and Sn species and the later could be the reduction of Pt–(O–Si
)y2−y species. Considering the nonuniform distribution of Sn over PtSn1.6/SiO2-MIM, the presence of Pt–(O–Si
)y2−y species is reasonable. The catalyst (PtSn1.6/SiO2-TSG) in which Pt precursor is uploaded by IM method exhibits three peaks, approximately at 91 °C, 193 °C and 420 °C. The peaks at 91 °C and 193 °C could be associated to the reduction of Pt precursor with no and weak interaction with the support respectively. The peak at 420 °C is attributed to the reduction of SnO2. The reduction temperature of PtSn1.6/SiO2-TSG is lower than that of PtSn1.6/SiO2-MTSG, indicating that the strong interaction between Pt and Sn exists in the MTSG preparation process. Based on the TPR results, it could be proposed that Pt species in PtSn1.6/SiO2-MTSG catalyst have the closest interaction with Sn species.
)y2−y and Pt–(O–Si
)y2−y species. Nevertheless, the dispersion of Pt particles in PtSn1.6/SiO2-TSG is 30.2%, resulting from the aggregation of Pt particles. These results are consistent with the TEM images.
| Catalyst | Loadinga (wt%) | SBET (m2 g−1) | Dpore (nm) | Vpore (cm3 g−1) | Db (nm) | Pt dispersionc (%) | |
|---|---|---|---|---|---|---|---|
| Pt | Sn | ||||||
| a Metal loading determined by ICP-OES.b Particle size determined by TEM images.c Pt dispersion determined by H2–O2 titration method. | |||||||
| Sn/SiO2-SG | — | 0.87 | 975 | 2.91 | 0.40 | — | — |
| Pt/SiO2-MTSG | 0.89 | — | 451 | 2.98 | 0.19 | 1.42 | 70.1 |
| PtSn1.6/SiO2-MTSG | 0.87 | 0.89 | 492 | 3.86 | 0.49 | 1.46 | 55.6 |
| PtSn1.6/SiO2-MIM | 0.90 | 0.81 | 538 | 3.22 | 0.35 | 1.84 | 47.9 |
| PtSn1.6/SiO2-TSG | 0.86 | 0.89 | 582 | 2.81 | 0.18 | 4.16 | 30.2 |
| PtSn1.0/SiO2-MTSG | 0.89 | 0.62 | 395 | 4.01 | 0.38 | 1.43 | 61.9 |
| PtSn3.3/SiO2-MTSG | 0.89 | 1.64 | 449 | 3.30 | 0.30 | 1.70 | 34.4 |
| PtSn8.2/SiO2-MTSG | 0.88 | 4.29 | 483 | 3.31 | 0.33 | 1.98 | 17.6 |
| Catalyst | Conv. (%) | Selectivity (%) | |
|---|---|---|---|
| EtOH | AcOEt | ||
| a Reaction conditions: 270 °C, 2.6 MPa, H2/AcOH (mol ratio) = 20, WHSV(AcOH) = 2.0 h−1. | |||
| Pt/SiO2-MTSG | 11.3 | 76.6 | 21.4 |
| Sn/SiO2-SG | 0.4 | 15.5 | 84.5 |
| PtSn1.6/SiO2-MTSG | 98.3 | 91.8 | 7.6 |
| PtSn1.6/SiO2-MIM | 78.4 | 90.0 | 9.2 |
| PtSn1.6/SiO2-TSG | 52.4 | 73.8 | 25.0 |
Based on these characterization results, it could be concluded that the MTSG method was an effective approach to fabricate the PtSn/SiO2 catalyst with high dispersion of Pt particles. Besides, this method presented a larger amount of Lewis acid sites, which is crucial for AcOH adsorption.8,9,38–40 Moreover, the homogenous distribution of Sn species in the catalyst prepared by MTSG method can give abundant surface hydroxyl, providing more opportunity for Pt precursor to disperse on the support by forming Pt–(O–Sn
)y2−y analogous species. Herein, the synergy between Pt and Sn species was enhanced over PtSn1.6/SiO2-MTSG, so it exhibited excellent catalytic performance. Nevertheless, PtSn catalysts prepared by other methods either have low Pt dispersion or lack enough Lewis acid sites, thus resulting in poor catalytic activity and lower ethanol selectivity.
![]() | ||
| Fig. 6 Effect of Sn/Pt molar ratio on the catalytic performance. Reaction conditions: 270 °C, 2.6 MPa, H2/AcOH (mol mol−1) = 20, WHSV(AcOH) = 1 h−1. | ||
Fig. S3† shows the TEM images and particle size distributions of bimetallic PtSnx/SiO2-MTSG catalysts with various Sn/Pt ratio. It shows that the average particle sizes increase with the increase of Sn/Pt ratio. According to the H2–O2 titration results, the dispersion of Pt particles decreases from 61.9% to 17.6% with the increasing Sn/Pt ratio (listed in Table 1). It may be resulted from the geometric or electronic effect between Pt and Sn species.41 The decrease of Pt dispersion suggests that the hydrogen adsorption ability of PtSn catalysts decreases with the increasing Sn/Pt ratio.42
CO-FTIR profiles of PtSnx/SiO2-MTSG are shown in Fig. 7. The adsorbed CO on the Pt/SiO2-MTSG catalyst clearly presents a sharp peak at 2067 cm−1 with a shoulder peak at 2043 cm−1, which corresponds to CO linearly bonded onto surface Pt atoms (COL, adsorbed atop of Pt atoms). The broad peak in the range of 1800–1900 cm−1 is attributed to the bridge-bound CO (COB, adsorbed on Pt–Pt pair sites).32,43,44 The linear CO peak intensity decreases notably and the bridge-bound CO peak disappears simultaneously with the increasing Sn/Pt ratio. Such decrease could be due to the electronic effect caused by Sn species or the migration of Sn species onto Pt particle surfaces, suggesting the close interaction between Pt and Sn species exists in the prepared PtSn catalysts.32,45 This agrees well with the dispersion results detected by H2–O2 titration method (listed in Table 1).
Pyridine-FTIR spectra of the bimetallic PtSnx/SiO2-MTSG catalysts with different Sn/Pt ratio are shown in Fig. 4B. All the catalysts exhibit obvious Lewis acid sites evidenced by the bands at 1452 cm−1 and 1609 cm−1.34 The catalysts with lower Sn/Pt ratio (below 1.6) show new bands at 1446 cm−1 and 1596 cm−1, which could be attributed to the adsorbed pyridine molecules bonded on the surface hydroxyl groups of silica.35 However, both bands disappear with the increasing Sn/Pt ratio (above 3.3), on account of the coverage of surface hydroxyl groups by Sn species. It is clear that the intensity of the band at 1452 cm−1 increases with the increasing Sn/Pt ratio, indicating the increasing surface acidity of PtSn catalysts. In order to correlate catalytic activity with the Lewis acid sites, the relative amount of Lewis acid sites in the PtSn catalysts are calculated.24 As listed in Table S2,† the relative amount of Lewis acid sites increase from 5.1 μmol g−1 to 56.4 μmol g−1 with the Sn/Pt ratio, resulting from the increment of surface SnOx species.
It was proposed that the reaction rate of acetic acid hydrogenation could be determined by the formation of acyl species or the reaction between H atoms and acyl species, depending on the ratio of the amount of Lewis acid sites versus the amount of exposed Pt active sites (denoted as L/Pt).7–9 However, none of these literatures distinguished the balance problem between Lewis acid sites and the exposed Pt active sites. Thus, in order to elucidate this balancing effect, the TOF based on exposed Pt atoms was correlated with the amount of Lewis acid sites. As shown in Fig. 8, the TOFPt increases with the increasing molar ratio of L/Pt in region I (from 0.18 to 0.54). Molar ratio of L/Pt is very low in this region, thus there is enough exposed Pt to provide dissociative hydrogen. So this increasing trend could originate from the increasing adsorption sites of AcOH. Hence, the determining step is the formation of acetyl. On the other hand, the catalytic activity barely changes with the further increment of L/Pt in region II (from 0.54 to 7.10). Such trend could be deduced that the Lewis acid sites in region II was so abundant that all the exposed Pt on catalyst surface could be involved in the reaction. Owing to the relative insufficient dissociated hydrogen, the hydrogenation is therefore controlled by the reaction between H and acyl species. Thus, conversion of AcOH decreases when the Sn/Pt ratio increases from 1.6 to 8.2. The decrease of ethanol selectivity could be attributed to either the decreasing of Pt dispersion or the increasing surface acidity of catalysts. Explanation for this speculation will be discussed in the next part.
![]() | ||
| Fig. 8 Effect of L/Pt molar ratio on the catalytic activity. Reaction conditions: 270 °C, 2.6 MPa, H2/AcOH (mol mol−1) = 20. | ||
The influence of residence time on the conversion of AcOH as well as the selectivity of ethanol and AcOEt was evaluated by varying the space velocity of AcOH. As shown in Fig. 9, AcOH conversion increases with the increase of residence time and achieve ca. 100% after the residence time exceeds 130 s. When the residence time is not long enough, the selectivity of ethanol remains constant. However, when the residence time is further increased, for example above 130 s, the selectivity of ethanol shows an increasing trend, indicating that partial of AcOEt could be conversed to ethanol through deep-hydrogenation. Similar results were obtained by the study on other kinetic parameters such as temperature, pressure and H2/AcOH molar ratio (listed in Table S3†). According to the variation of the catalytic performance, parallel reactions may exist in AcOH hydrogenation system.
![]() | ||
| Fig. 9 Hydrogenation performance of AcOH over PtSn1.6/SiO2-MTSG catalyst as a function of residence time. Reaction conditions: 270 °C, 2.6 MPa, H2/AcOH (mol mol−1) = 20. | ||
Hence, a schematic model regarding the hydrogenation of AcOH over PtSn1.6/SiO2-MTSG catalyst is proposed (Scheme 1), in which the formation of ethanol and AcOEt are proposed to be parallel reactions. According to literatures, AcOH is firstly dissociated into adsorbed acetyl on Lewis acid sites. Then the acetyl intermediates could combine with H atoms activated by metallic Pt to form adsorbed acetaldehyde. The adsorbed acetaldehyde species then can undergo hydrogenation to form ethoxy groups.9 In the parallel reactions model, ethoxy groups can not only react with activated H atoms but also those adsorbed acetyl, which produced ethanol and AcOEt respectively. Affluent activated H atoms could favor the formation of ethanol, while the limited amount of activated H atoms may lead to the formation of AcOEt. Therefore, the competitive effect of two reactions could lead to the changes of products selectivity. Moreover, the hydrogenation of AcOEt to ethanol could occur by deep-hydrogenation when more H atoms were provided. Herein, one can infer that the synergy effect among Pt and SnOx species is responsible for the hydrogenation of AcOH into ethanol. The parallel reaction mechanism for AcOH hydrogenation proposed by this study may provide valuable help in catalyst design and process development.
![]() | ||
| Scheme 1 Schematic model of AcOH hydrogenation to ethanol and AcOEt over the PtSn1.6/SiO2-MTSG catalyst. | ||
The long-term catalytic behavior results are also shown in Fig. 9. It is demonstrated that the well fabricated PtSn1.6/SiO2-MTSG catalyst remains its excellent activity and ethanol selectivity in the hydrogenation of AcOH into ethanol. It shows a prospective future in the massive synthesis of ethanol by AcOH hydrogenation from an industrial application viewpoint.
Footnote |
| † Electronic supplementary information (ESI) available: Fig. S1 shows the N2 adsorption–desorption isotherms and pore size distribution of reduced catalysts. Fig. S2 displays the HAADF-STEM and EDX elemental-mapping images of the reduced PtSn1.6/SiO2-MTSG catalyst. Fig. S3 shows the TEM images of reduced PtSnx/SiO2-MTSG catalysts. Fig. S4 represents the deconvolution of the fit about FTIR spectra of chemisorbed pyridine on reduced catalysts. Fig. S5 displays the hydrogenation performance of AcOH over PtSn1.6/SiO2-MTSG catalyst as a function of WHSV(AcOH). Table S1 shows the fitted results of H2-TPR experiments of catalysts. Table S2 displays the integral quantity of Lewis acid sites of reduced catalysts. Table S3 illustrates the hydrogenation performance of AcOH over PtSn1.6/SiO2-MTSG catalysts. See DOI: 10.1039/c6ra09199g |
| This journal is © The Royal Society of Chemistry 2016 |