Lihui
Ou
*a and
Junxiang
Chen
b
aHunan Province Cooperative Innovation Center for the Construction & Development of Dongting Lake Ecologic Economic Zone, College of Chemistry and Materials Engineering, Hunan University of Arts and Science, Changde, 415000, China. E-mail: oulihui666@126.com
bCAS Key Laboratory of Design and Assembly of Functional Nanostructures, Fujian Provincial Key Laboratory of Nanomaterials, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou, 350002, China
First published on 25th November 2019
A recently proposed CO coverage-dependent electrochemical model combined with the calculation of electronic structure is applied for the first time to study the effect of the overpotential on Cu-catalyzed CO electroreduction mechanisms by changing the coverage of surface adsorbed CO. The results show that the presently defined CH2O and CHOH pathways may be able to occur parallelly under different overpotentials. However, high overpotentials will facilitate CO electroreduction, thus explaining why a high overpotential is required during CO2 electroreduction in experiments on Cu electrodes. The potential-limiting step may be further CO electroreduction into CHO, which is considered as the origin of the experimentally observed high overpotentials. The analyses of electronic structure show that an adsorbed COδ− species is formed on the Cu electrodes, validating the previous experimental speculations on electron transfer between CO and Cu electrodes. More and more electrons are transferred into the π antibonding orbitals of the adsorbed CO with increasing surface CO coverage, leading to increasing overpotential and weaker and weaker CO bonding with the Cu surface. Thus, the significantly lower barrier of further CO electroreduction at higher overpotential can be correlated with lower CO adsorption energy. Interestingly, it is found that there is greater localization of electrons around the C than the O atom in the adsorbed CO molecule, explaining why the hydrated proton prefers to reach the C atom to form intermediate CHO rather than COH.
In situ infrared spectroscopy has offered some insights into the nature of adsorbed CO2, CO and formate. Moreover, the observed CH2 species is able to lead to the formation of hydrocarbons CH4 and C2H4 as shown by in situ X-ray photoelectron spectroscopy and Auger electron spectroscopy on Cu electrodes.5,23,24 However, experimental understanding of CO2 electroreduction mechanisms has been challenged and hindered by difficulties in performing spectroscopic characterizations of the various possible intermediates that are crucial to the production of the observed electroreduction products at electrochemical interfaces. Most recently, the feeding technique conducted by Schouten et al. has been used to analyze experimentally CO2 electroreduction mechanisms on Cu by feeding postulated intermediates into reaction systems and the subsequent product distribution is taken as evidence that the intermediates may indeed play a crucial role in the electrochemical reduction mechanisms.7 However, this approach has significant limitations since only stable intermediates can be studied. Thus, only rough and incomplete CO2 electroreduction pathways can be speculated based on early experiments. In recent decades, theoretical calculations have become a powerful tool for investigating various key electrocatalytic reactions and numerous theoretical studies have contributed to our understanding of CO2 electroreduction mechanisms on Cu electrodes by using various theoretical models,8,15,20–22,25–29 which can provide mechanistic insights that are not accessible from experiments alone and identify favored intermediates. Taking CO2 electroreduction into the CH4 product on Cu(111) as an example with the application of various theoretical models, Nørskov et al. reported the first full studies of CO2 electroreduction into CH4 on the Cu surface by using the computational hydrogen electrode (CHE) model with a hexagonal H2O overlayer.8 They concluded that CO, CHO, CH2O and CH3O are key intermediates during CH4 formation along which CO electroreduction into CHO is the potential-limiting step. The CHE model could calculate conveniently the reaction free energy changes of the proton-coupled electron transfer (PCET) steps at applied potentials. However, the activation barriers of all PCET steps are neglected by the CHE model and thus accurate kinetics are not available. The kinetic barriers for each PCET step are subsequently considered with incorporation of applied potentials at the electrode/aqueous interfaces. In contrast to the study of the CHE model, Asthagiri et al. reported the complete CO2 electroreduction pathways by calculating the applied potential-dependent barriers on Cu(111) in the aqueous phase based on Butler–Volmer theory,15,20,27 in which an adsorbed H atom is used to represent a proton–electron couple under the equilibrium potential and the COH species is predicted to be the kinetically dominant key intermediate, sequentially resulting in CHx (x = 0–4) species. The difference in electroreduction pathways between both of these models may show the importance of considering kinetic barriers. However, the generality of the assumption that the adsorbed H direct transfer represents the kinetics of a PCET step may be still in question. Moreover, the inclusion of only one or two explicit H2O molecules or a fixed H2O bilayer would be insufficient for the modeling of the aqueous-phase environment. By changing the number of H atoms added to the bilayer H2O structure to tune the electrode potentials, a DFT-based extrapolation scheme was applied to studies of CO2 electroreduction into CH4 on Cu(111) by extrapolating the calculated barriers in model systems of increasing cells.28,29 Skúlason et al. speculated that the COH, CHOH, CH2OH, CH2 and CH3 species are the reduction intermediates after CO formation. The major issue for this scheme is that it is extremely costly computationally due to extrapolation of an increasingly larger size. However, inconsistent conclusions were obtained by Nørskov et al. by using a similar electrochemical model,21 in which CHO formation is predicted to be more favorable than COH on Cu(111) thermodynamically and kinetically, subsequently leading to the formation of CHOH and CHx (x = 2–4), in contradiction to the report from Asthagiri and Skúlason. By using the implicit solvation model, CO2 electroreduction pathways were studied based on a DFT methodology that predicted onset potentials at constant electrochemical potentials by Goddard et al.22 They proposed that forming COH is more favorable than forming CHO although the difference in the barrier is subtle on Cu(111). The subsequent COH electroreduction results in the formation of CHOH and CHx (x = 1–4), which is different from the pathways reported by Asthagiri et al. The drawback of this model may be that the implicit solvation model cannot describe hydrogen bonds and interactions between adsorbates and H2O molecules. Despite numerous theoretical efforts, the mechanism for CH4 formation on Cu(111) is still under debate, and the difference in the CO2 electroreduction pathways may lie in the different electrochemical models. Furthermore, most previous theoretical investigations were done at constant charge, resulting in changes of the work function/electrode potential during a PCET step due to the limited cell size, whereas the experimental PCET steps take place at a constant electrode potential. Thus, the modeling of electrode/aqueous interfaces remains a subject of ongoing discussion to obtain sufficiently accurate computational models for describing electrochemical systems and fill in the lacking intermediate in the uncertain reaction mechanisms.30–32
The previous experimental reports from the groups of Frese, Ito and Hori felt that adsorbed CO was a crucial intermediate during CO2 electroreduction into CH4, with further CO electroreduction as the slowest step in the overall reduction pathways, which can dominate most of the Cu electrode surface.11,33–36 Thus, it can be speculated that the CO2 electroreduction kinetics is dependent on the coverage of CO. Based on this, an improved CO coverage-dependent electrochemical model can be used to predict CO coverage-dependent equilibrium potentials and elucidate CH4 production mechanisms during CO2 electroreduction on Cu(111) by calculating potential-dependent free energy profiles at the predicted equilibrium potentials. The solvent effect is essential in modeling electrode/aqueous interfaces and the previous study from Nørskov et al. validated that even a single H2O bilayer is sufficient to describe the electrical double layer.37 An explicit solvation model containing two relaxed H2O bilayers is included in the present model which allows us to better simulate the interactions between H2O molecules and adsorbates. Furthermore, the interactions between solvent and adsorbents have been elaborated in our most recent work.38,39 It is aim of this paper to offer in depth insight focused on specifically the influence of the overpotential on the CO2 reduction mechanisms on Cu(111) by applying the CO coverage-dependent electrochemical model, which may be able to result in better understanding of the possible electroreduction pathways.
CO2 + * + 2(H+ + e−) → CO* + H2O | (1) |
ΔG(θ) = ΔE(θ) + 2eU − TΔS(θ) + ΔZPE + 2kBTln(θ/1 − θ) | (2) |
(3) |
Considering that the entropies of the adsorbed molecules are small when compared with the entropies of gases and that the zero-point energy of the adsorbed CO molecules on Cu(111) is little based on our present calculations, the contribution from the changes of entropy and ZPE together to ΔG(θ) has been estimated to be −0.42 eV for standard temperature (298 K) according to the available data from Nørskov et al.8 Thus, eqn (3) can be converted to the following form of eqn (4).
(4) |
The values of , EH2O, ECO2 and EH2 are directly available via DFT calculations. By changing the coverage of adsorbed CO on Cu(111), a series of values of can be obtained. The values of ΔE(θ) at various coverages can be obtained according to eqn (3) by differentiating the plots of against θ.
Subsequently, further CH2O and CHOH electroreduction are considered. Further hydrogenation of intermediate CH2O may be able to form CH2OH and CH3O species on Cu(111). Barriers of 0.32 and 0.23 eV are required for CH2OH and CH3O formation at the overpotential of 0.13 V, respectively, being slightly higher than those of the corresponding processes (0.19 and 0.16 eV) at the relatively higher overpotential of 0.37 V, as shown in Fig. S3 (ESI†). Thus, the results show that further CH2O electroreduction is more favorable at the higher overpotentials. However, both of these intermediates can be easily formed in the range of low overpotentials due to the surmountable barriers at room temperature. Our presently calculated barriers for further CH2O electroreduction into CH2OH and CH3O are significantly lower than those reported by Asthagiri et al. at low overpotentials,20 in which barriers of ca. 0.90 and 0.60 eV need to be overcome at an applied potential of ca. −0.10 V (vs. RHE). The difference in the used solvation models may lead to the inconsistencies, where only a single H2O molecule or fixed H2O bilayer structure was employed in previous theoretical studies. Due to the favorable CHOH formation, the CH2OH species can be also formed by CHOH hydrogenation with almost identical barriers of ca. 0.50 eV at the low overpotentials of 0.13 and 0.37 V (see Fig. S4, ESI†). Simultaneously, we also note that the CH species can be formed by CHOH hydrogenative dissociation along with H2O formation with extremely low barriers of ca. 0.20 eV in the range of low overpotentials, which are significantly lower than those of CHOH hydrogenation into CH2OH, thereby suggesting that CH formation is more favorable during further CHOH electroreduction at low overpotentials. Considering the almost identical and extremely low barriers between further CH2O electroreduction into CH2OH and CH3O species and further CHOH electroreduction into CH species, we can speculate that parallel pathways may occur through the CH2O and CHOH intermediates during CO2 electroreduction in the range of low overpotentials.
Further electroreduction of intermediate CH2OH into CH3OH through direct hydrogenation and into the CH2 species through hydrogenative dissociation along with H2O formation are considered on Cu(111). The calculated barriers for further CH2OH electroreduction into CH3OH and CH2 are ca. 0.36 and 0.07 eV at the overpotential of 0.13 V, which is in good agreement with the previously reported values of ca. 0.29 and 0.05 eV from Goddard et al. for the formation of product CH3OH and intermediate CH2,46 further validating the rationality of our presently employed electrochemical model. However, they are 0.34 and 0.31 eV when the overpotential is increased to 0.37 V, respectively, as shown in Fig. S5 and S6 (ESI†). The results suggest that CH2 formation is more and more difficult, whereas CH3OH formation becomes more and more favorable with increasing overpotential. However, it is observed that the CH3OH and CH2 species may be able to be formed simultaneously due to the surmountable barriers at room temperature in the range of low overpotentials. Although the CH3O species can be easily formed, a significantly higher barrier of ca. 0.70 is required for its further hydrogenation into CH3OH compared to CH2OH electroreduction at the overpotentials of 0.13 and 0.37 V (see Fig. S5 and S6, ESI†). Thus, we can speculate that the CH3O species is only a spectator during CO2 electroreduction at low overpotentials. Additionally, the CH2 intermediate can be also produced with a low barrier of ca. 0.35 eV at low overpotentials by further hydrogenation of the CH species formed by CHOH hydrogenative dissociation. Subsequently, the final CH4 product can be formed by CH2 serial hydrogenation, in which almost identical barriers of ca. 0.32 and 1.30 eV need to be overcome for CH2 hydrogenation into CH3 and then CH4 at the overpotentials of 0.13 and 0.37 V, respectively, as shown in Fig. S7 (ESI†). Almost identical barriers for CH serial hydrogenation into CH4 are also obtained at low overpotentials based on the electrochemical model proposed by Asthagiri et al.20 Thus, the present results suggest that the barriers for CH serial hydrogenation into CH4 will be almost not affected in the range of low overpotentials. The corresponding energies of various possible elementary steps during CO2 electroreduction into CH4 and CH3OH on Cu(111) under conditions of low overpotentials are listed in Tables S1 and S2 (ESI†). Based on the energetics, the presently defined CH2O and CHOH pathways on Cu(111) may be able to be proposed at low overpotentials, as can be seen in Fig. 2 and 3.
Fig. 2 The overall energy diagram of CO electroreduction into CH4 and CH3OH on Cu(111) at the low overpotential of 0.13 V: (a) CH2O pathway; (b) CHOH pathway. |
Once CH2O is formed, the further hydrogenation of CH2O into CH3O occurs more easily with a barrier of ca. 0.29 eV compared to that into CH2OH with a barrier of ca. 0.54 eV on Cu(111) under conditions of high overpotentials (see Fig. S10, ESI†). In fact, based on the electrochemical model previously proposed by Asthagiri et al.,20 activation barriers of ca. 0.25 and 0.60 eV for CH3O and CH2OH formation on Cu(111) at an applied potential of ca. −0.80 V (vs. RHE) can be obtained, respectively. Thus, the results calculated by using our present and previous electrochemical models show that CH3O formation is more favorable during further CH2O electroreduction on Cu(111) at high overpotentials. However, CH2OH may be able to be easily formed by further CHOH hydrogenation with a surmountable barrier of ca. 0.24 eV at room temperature at high overpotentials, as shown in Fig. S11 (ESI†). Furthermore, we find that the barrier for further CHOH electroreduction into the CH2OH species is considerably decreased on Cu(111) compared with that at low overpotentials, again indicating that the effect of the overpotential on the barriers is significant. A barrier of ca. 0.29 eV is required for CHOH hydrogenative dissociation into the CH species along with H2O formation at high overpotentials, which is almost identical to CHOH hydrogenation into CH2OH (see Fig. S11, ESI†), suggesting that intermediates CH2OH and CH may be able to be formed simultaneously through further CHOH electroreduction. Due to more favorable formation of CH3O and CH2OH at high overpotentials, CH3OH may be able to be produced by further CH3O and CH2OH hydrogenation with surmountable barriers of ca. 0.01 and 0.40 eV at room temperature, respectively. Moreover, hydrogenative dissociation of CH2OH into CH2 and H2O needs to overcome a higher barrier of ca. 0.56 eV than into CH3OH, as shown in Fig. S12 (ESI†). Compared with those under conditions of low overpotentials, the barriers of CH3OH production are considerably reduced, suggesting that CH3OH formation may be kinetically favorable at high overpotentials. Intermediate CH2 may be also able to be formed by CH direct hydrogenation. The calculated barrier is 0.36 eV under conditions of high overpotentials, which is lower than that of CH2 formation through CH2OH hydrogenative dissociation, showing that the CH2 species is mainly formed by CH direct hydrogenation. Finally, the product CH4 can be formed by CH2 serial hydrogenation, in which barriers of ca. 0.44 and 0.01 eV need to be overcome for CH2 hydrogenation into CH3 and then its further hydrogenation into the CH4 product at high overpotentials (see Fig. S13, ESI†), respectively. We note that the hydrogenation barriers for CH into CH2 and then CH3 will be almost not affected with increasing overpotential. However, compared with that under conditions of low overpotentials, the barrier for CH3 hydrogenation into CH4 is considerably reduced to barrierless at high overpotentials. In fact, based on the electrochemical model from Asthagiri et al.,20 the hydrogenation barrier of CH3 to CH4 on Cu(111) also goes from ca. 0.90 eV at low overpotentials to barrierless at high overpotentials, confirming the rationality of the presently proposed electrode/aqueous interface model. The energies of the various possible elementary steps during CO2 electroreduction into CH4 and CH3OH on Cu(111) under conditions of high overpotentials are listed in Table S3 (ESI†). The corresponding CH2O and CHOH pathways on Cu(111) are given in Fig. 4.
By comprehensive consideration of the overall energy diagrams during CO2 electroreduction into the CH4 product at the overpotentials of ca. 0.13 and 0.37 V (see Fig. 4 and 5), we speculate that the presently defined CH2O pathway can occur through the CHO, CH2O, CH2OH and CHx (x = 2–4) species in the range of low overpotentials with barriers of ca. 1.10 and 0.88 eV, respectively. Considering that there are two possibilities for further CHOH electroreduction into CH4, the overall energy diagrams are scrutinized in the CHOH pathways. It is found that the CHOH pathway can parallelly occur through CHOH electroreduction into the CH and CH2OH species due to almost identical barriers of ca. 1.10 eV at the overpotential of ca. 0.13 V, whereas it may be only able to take place through CHOH electroreduction into the CH species due to the lower required barrier than that of CHOH electroreduction into CH2OH (1.00 eV vs. 1.22 eV) at the relatively higher overpotential of 0.37 V. In the meantime, we note that the presently defined CH2O and CHOH pathways may be able to parallelly occur due to the subtle difference in the barriers in the range of low overpotentials. However, higher overpotentials will facilitate CO2 electroreduction based on the above difference in the barriers. The optimal CO electroreduction pathways on Cu(111) can be proposed under conditions of low overpotential by the above analysis, as can be seen in Fig. 5.
To further determine the effect of the overpotential on CO2 electroreduction into the CH4 product, the overall energy diagrams at the higher overpotential of ca. 0.80 V are scrutinized (see Fig. 5). The CH2O pathway can occur through the CHO, CH2O, CH2OH and CHx (x = 2–4) species at high overpotentials with barriers of ca. 0.51 eV. In the CHOH pathways, two possibilities via further electroreduction of the CHOH species into CH along with H2O formation and the CH2OH species may be able to parallelly occur due to the identical barriers of ca. 0.51 eV. Thus, it is also speculated that the CH2O and CHOH pathways can parallelly occur and the optimal CO electroreduction pathways on Cu(111) can be proposed at high overpotentials, as shown in Fig. 5. The present results are found to be in excellent agreement with the available experimental data. For example, Dewulf et al. presented that adsorbed CHO and CH2 species are crucial intermediates for CH4 formation at Cu electrodes using in situ X-ray photoelectron spectroscopy and Auger electron spectroscopy during electroreduction of CO2,5 and the most recent experiment from Koper et al. studied CO2 electroreduction into CH4 on two basal planes, Cu(111) and Cu(100), using online electrochemical mass spectrometry, in which the electrochemical reduction of the various species can be measured online when the electrode potentials are changed by the tip-based sampling technique and the formation and consumption of intermediates during CO2 reduction can be followed, and showed that CHO is possibly the key intermediate toward the formation of the CH2 and CH3 species and CH4 product through the cleavage of the C–O bond,7,14 further validating the rationality of the presently employed methodology. However, the presence of high overpotentials notably reduces the barriers of the overall energy pathway. Simultaneously, it is observed that CH4 formation via further CH3 electroreduction requires a high barrier of ca. 1.40 eV at low overpotentials, whereas it is an almost barrierless process (ca. 0.01 eV) under conditions of high overpotentials. Thus, the present results can conclude that high overpotentials will be in favor of CO2 electroreduction into CH4, again explaining why a high overpotential is required during CO2 electroreduction into hydrocarbons in experiments on Cu electrodes.3–6
Additionally, we observed that an extremely low barrier of ca. 0.20 eV is required for the inverse process of further CO electroreduction into CHO, which is lower than that of further CHO electroreduction under conditions of low overpotentials. Thus, it can be speculated that intermediate CHO may be not stable in the range of low overpotentials, and can easily transform back to CO by C–H bond cleavage, leading to not easy formation of key intermediate CHO and preventing CO2 electroreduction on Cu electrodes. However, the barrier for the inverse process of further CO electroreduction into CHO is increased to ca. 0.39 eV at the high overpotential of ca. 0.80 V, which is notably higher than that of further CHO electroreduction (ca. 0.10 eV). Simultaneously, the barrier for further CO electroreduction into CHO is also significantly lowered at high overpotentials. Thus, we can speculate that the potential-limiting step may be further CO electroreduction into CHO during CO2 electroreduction into the CH4 product, which is considered as the origin of the observed experimentally high overpotentials. In fact, the previous experiments from Hori et al. and Frese et al. also showed that adsorbed CO electroreduction is the key rate-determining step for the overall reduction reaction of CO2 into hydrocarbons.11,35 Thus, the present finding on the potential-limiting step is also in excellent agreement with the available experimental data.
To further determine the origin of the overpotential, partial density of states (PDOS) analyses of an adsorbed CO molecule that can be further electrochemically reduced into hydrocarbons are performed on Cu(111) with different surface CO coverages in the present study since it can describe the number of electronic states near the Fermi energy level that are available to be occupied, as shown in Fig. 6. A high intensity of the PDOS near the Fermi energy level means high localization of electrons. For comparison, the PDOS of free CO is also included. As can be seen in Fig. 6(a) and (c), almost no obvious peaks exist in the PDOS of s orbitals of C and O atoms of adsorbed CO compared with that of the free CO molecule near the Fermi level, suggesting that maybe no significant electron transfer occurs into the s state of adsorbed CO. However, obvious peaks are observed in the p orbitals of C and O atoms near the Fermi level (ca. 0–3 eV), as shown in Fig. 6(b) and (d), indicating that electrons from the surface and electrochemical interfaces are transferred into π antibonding orbitals of the adsorbed CO molecule. Thus, an adsorbed COδ− species is formed on Cu(111), validating the previous experimental speculations on electron transfer between CO and Cu electrodes.5,11 To ascertain the effect of the CO coverage on electron transfer, the PDOS of p orbitals near the Fermi level in Fig. 6 is amplified, as can be seen in Fig. 7. The calculated value of the Fermi level is ca. 3.45, 3.49, 3.72 and 3.79 eV in the present research range of CO coverage, respectively, which gradually increases with increasing CO coverage. Furthermore, it is clearly observed that the intensity of the peaks becomes larger and larger with increasing CO coverage, indicating more and more electron transfer into π antibonding orbitals of the adsorbed CO molecule and higher and higher localization of electrons around C and O atoms. The higher peak intensities of π antibonding orbitals near the Fermi level (ca. 0–3 eV) indicate weaker CO bonding with the Cu surface, as can be confirmed by the calculated less and less negative CO adsorption energies of −2.86, −1.15, −0.76 and −0.59 eV with increasing surface CO coverage at the presently simulated electrochemical interfaces, suggesting that higher overpotentials will reduce the CO adsorption strength. Thus, the abovementioned significantly lower activation barrier of further CO electroreduction into CHO on Cu(111) at higher overpotential can be correlated with lower CO adsorption energy. Simultaneously, higher peak intensities in the PDOS of p orbitals of C atoms are observed compared to those of O atoms (see Fig. 7), suggesting greater localization of electrons around C atoms, which will be further confirmed by the following Löwdin population analysis.
The quantative analysis of electronic structures will facilitate the further understanding of the origin of the overpotentials and the effect of surface CO coverage on electron transfer during CO2 electroreduction. The Löwdin charge of an adsorbed CO molecule that can be further electrochemically reduced into hydrocarbons on Cu(111) with different surface CO coverages can be obtained from Löwdin population analysis based on the projected electron densities of states. As shown in Table 1, the net electron gains of the p orbitals of the C and O atoms in the adsorbed CO molecule that can be further electrochemically reduced into hydrocarbons gradually increased with increasing surface CO coverage, which are filled into π antibonding orbitals based on the above PDOS analysis. Furthermore, the total net electron gain of adsorbed CO also gradually increased, indicating that the net electron gain of the s orbital is negligible and further conforming the results obtained by the above PDOS analysis. Interestingly, it is noted that the C atom obtains more electrons than the O atom in the adsorbed CO molecule, which quantatively confirms the abovementioned greater localization of electrons around C atoms and explains why the hydrated proton tends to reach the C atom to form intermediate CHO rather than COH as indicated in the above MEP analysis. Simultaneously, the greater and greater electron transfer into the π antibonding orbitals of the adsorbed CO molecule with increasing CO coverage may lead to increasing overpotential and easier and easier occurrence of further CO electroreduction.
Difference of electron/atom (Δq)a | |||||
---|---|---|---|---|---|
Orbital | 1CO | 2CO | 3CO | 4CO | |
a Positive values of Δq imply an electron gain by the component. | |||||
C | s | −0.2467 | −0.2363 | −0.2356 | −0.2308 |
p | +0.7740 | +0.7809 | +0.7833 | +0.8039 | |
Total | +0.5237 | +0.5446 | +0.5477 | +0.5730 | |
O | s | +0.0208 | +0.0192 | +0.0169 | +0.0179 |
p | +0.0452 | +0.0619 | +0.0778 | +0.0920 | |
Total | +0.0660 | +0.0811 | +0.0947 | +0.1099 |
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c9cp05043d |
This journal is © the Owner Societies 2020 |