Stefan
Smetaczek
*a,
Andreas
Wachter-Welzl
a,
Reinhard
Wagner
b,
Daniel
Rettenwander
cd,
Georg
Amthauer
b,
Lukas
Andrejs
e,
Stefanie
Taibl
a,
Andreas
Limbeck
a and
Jürgen
Fleig
a
aInstitute of Chemical Technologies and Analytics, Technische Universität Wien, Vienna, Austria. E-mail: stefan.smetaczek@tuwien.ac.at
bDepartment of Chemistry and Physics of Materials, University of Salzburg, Salzburg, Austria
cInstitute for Chemistry and Technology of Materials, Graz University of Technology, Graz, Austria
dChristian Doppler Laboratory for Lithium Batteries, Institute for Chemistry and Technology of Materials, Graz University of Technology, Graz, Austria
eAustrian Institute of Technology, Vienna, Austria
First published on 21st February 2019
The attractiveness of cubic Li7La3Zr2O12 (LLZO) based garnets lies in their high ionic conductivity and the combination of thermal and electrochemical stability. However, relationships between the exact chemical composition and ionic conductivity are still not well understood and samples of very similar composition and processing may show very different properties. In this contribution, we investigated Al stabilized LLZO by employing a combination of local conductivity measurements and 3D local chemical analysis using laser ablation inductively coupled plasma optical emission spectroscopy (LA-ICP-OES). These measurements revealed significant conductivity variations across some samples, with variations up to almost one order of magnitude, as well as an inhomogeneous elemental distribution of Al and Li, largely along the samples. Surprisingly, neither the local Al content nor the local Li content showed a clear correlation with the local Li-ion conductivity. Accordingly, an in-depth understanding of the conduction properties of Al stabilized LLZO has to concentrate on aspects beyond the simple chemical composition. Yet unknown factors (e.g. oxygen vacancies or local cation site occupancies) seem to have a much higher impact on the ionic conductivity than the exact stoichiometry.
A lot of recent research has focused on LLZO and knowledge became available on the ideal synthesis route, effects of water, degradation phenomena, maximizing ionic conductivity, conduction paths, etc.20–40 However, many aspects are still not understood and several challenges remain. The instability of LLZO in an ambient atmosphere (due to moisture and CO2) as well as interfacial issues could hamper the application of LLZO in all-solid-state batteries.41,42 Also resource criticality should be kept in mind, e.g. of La, particularly in the case of large-scale LLZO production.43
Another important aspect that needs clarification is the partly very different ionic conductivity reported for nominally identical samples. In addition to interfacial effects on the local or effective conductivity of samples, also strongly scattering bulk conductivities are reported. For LLZO stabilized by Al, which was the first element found to stabilize the cubic phase at room temperature,44,45 published values for bulk conductivity cover essentially the entire range from 10−4 S cm−1 to 10−3 S cm−1, whereby overall effective conductivities are often lower than the bulk value.46–49 A recent study on 44 nominally identical samples revealed conductivities from about 10−5 S cm−1 to 0.8 × 10−3 S cm−1 for LLZO with 0.2 Al per formula unit (pfu).46 Conductivity variations thus cannot be attributed only to the different preparation routes and degradation effects and must be caused by another (yet unknown) parameter. One likely reason is unintended composition changes during LLZO synthesis,46 which usually involves high temperatures at which lithium oxide loss and associated diffusion processes are to be expected. Investigating the impact of compositional variations and correlating them with the observed conductivity fluctuations is therefore crucial for understanding the conductivity behaviour of LLZO.
Possible effects of local compositional variations beyond the surface region (i.e. in depths of several μm) were already investigated by laser-induced breakdown spectroscopy (LIBS) and even 3D elemental maps of LLZO ceramics became accessible.29,33–35 However, they were not correlated with local conductivities, rather the focus was on the impact of the composition on the overall interfacial resistances of the samples.29,33
In our study, bulk compositional variations within Al stabilized LLZO ceramics are investigated as a possible reason behind the scatter of bulk conductivity values. The relationship between the local bulk ionic conductivity and the corresponding local bulk composition is investigated by combining microelectrode impedance spectroscopy with laser ablation (LA) inductively coupled plasma (ICP) optical emission spectroscopy (OES). Impedance spectroscopic studies on circular microelectrodes with diameters of 100 μm revealed the bulk conductivity of a similarly sized region beneath.50–54 Afterwards, in several depths beneath the electrodes (up to 100 μm deep) the chemical composition was analyzed by LA-ICP-OES. These measurements enabled the analysis of correlations between local Al or Li content and local ionic conductivity. In contrast to many common interpretations of ion conduction in LLZO, neither the Al nor the Li content shows a simple correlation with the ionic conductivity. Thus, further local structural or compositional factors have to play a key role in the ionic conductivity of LLZO.
Local conductivities were measured by means of microelectrodes. Using photolithographic techniques in combination with ion beam etching, circular electrodes with diameters of 100 μm were prepared from the macroscopic Pt/Ti thin films on top of the samples. Microelectrode measurements were performed at ambient temperature (T = 23.5 °C). Tungsten needles were used to contact the microelectrodes under an optical microscope. The position of the needles was adjusted using mechanically controlled micromanipulators. Fig. 1a illustrates the measurement setup and Fig. 1b shows a part of the microelectrode array on top of the sample. Impedance spectroscopic measurements (Novocontrol Alpha Analyzer) were then performed between a microelectrode and a counter electrode on the bottom side (Pt thin film with Pt paste at the sample edges for contact reasons).
![]() | ||
Fig. 1 (a) A scheme of the microelectrode measurement setup. (b) The optical microscope image shows a part of the microelectrode array on top of a sample. |
All LA-ICP-OES experiments were performed using a NWR213 laser ablation system (ESI) equipped with a frequency quintupled 213 nm Nd:YAG laser and a fast-washout ablation cell. The laser ablation device was coupled to an iCAP 6500 RAD ICP-OES instrument (ThermoFisher Scientific) using PTFE tubing with an inner diameter of 0.4 mm. A plasma torch with a corrosion-resistant ceramic injector tube of 1.5 mm inner diameter was utilized for all experiments. Cell washout was performed using He as the carrier gas, which was mixed with Ar make-up gas upon introduction into the plasma. For data acquisition iTEVA software (v.2.8.0.96) provided by the manufacturer of the instrument was used. The two most sensitive and non-interfering emission lines of Al, La, Li and Zr were measured. For each element, only one emission line was used for signal quantification, the other one was used for quality control of the analysis.
Two-dimensional elemental distribution images of whole LLZO pellets were created using line scan ablation patterns with adjoining lines. For all documented experiments a laser beam diameter of 100 μm and a scan speed of 100 μm s−1 were used, which result in a lateral resolution of 100 μm. In order to remove surface contaminations as well as near-surface segregation phases (e.g. Li2CO3 due to air exposure) a pre-ablation step consisting of a similar line scan pattern with a laser beam diameter of 250 μm was used prior to the measurements. This ablation step removes LLZO (or related phases) of about 2 μm thickness. Using the obtained time-resolved signal intensities, image processing was performed using the software ImageLab (v.2.41, Epina GmbH). Detailed information about the parameters used for the LA-ICP-OES experiments is given in the ESI.†
Compositional analysis was performed in sample depths of ca. 5, 50 and 100 μm. This depth (with respect to the original surface) was reached by mechanical polishing of the surface with SiC grinding paper (#2000) for certain times. In additional calibration studies these polishing times were related to the amount of removed material. Hence, a kind of three-dimensional compositional map became accessible. This approach also ensured that the measured “local” electrical conductivities (averaging over depths between 0 and 200 μm) could be truly related to the bulk composition within this specific sample region.
For signal quantification, matrix-matched standards were prepared by isostatically pressing LLZO powders with different Al contents into pellets. To obtain these powders, sintered LLZO pellets with various Al concentrations were crushed, ground, and homogenized using an agate mortar. To determine their chemical composition, a part of each powder was digested using a borax fusion, dissolved in a HF/HCl/HNO3 mixture, and analyzed by conventional liquid ICP-OES.
Measurement of the pressed LLZO standards was performed directly before sample analysis, using the same instrumental parameters and laser settings. To ensure that the ablated material is representative of the whole pellet, several line scan patterns with adjoining lines were used for each standard. Per standard, four ablation patterns each having an area of about 0.4 mm2 were applied on different positions of the pellet. Each pattern was ablated three times.
For the correction of variations in sample ablation and transport, as well as instrumental drifts, an internal standard-independent normalization technique presented by Liu et al. was adapted for LLZO samples.57 Based on the consideration that the sum of all metal oxides in LLZO adds up to 100 wt%, this approach uses the intensities of all measured elements expressed as the corresponding oxides (Al2O3, Li2O, La2O3, ZrO2) for signal normalization.
The accuracy of the LA-ICP-OES analysis was estimated using the prepared LLZO standards for external validation. For the Al/La, Li/La and Zr/La atomic ratio average relative deviations of 9.1%, 1.2% and 0.8% were obtained, respectively. The comparably high deviation of the Al/La ratio is attributable to the relatively low Al concentration.
More detailed information on the quantification of the described LA-ICP-OES measurements will be published in a separate paper.
Thickness (h) cm | Area (A) cm2 | |
---|---|---|
A | 0.41 | 0.34 |
B | 0.46 | 0.35 |
C | 0.10 | 0.45 |
D | 0.21 | 0.39 |
![]() | ||
Fig. 2 SEM image of an Al stabilized LLZO (Al = 0.20) garnet after sintering at 1230 °C for 6 h. Dark spots indicate the existence of some pores. |
![]() | ||
Fig. 3 XRPD patterns of samples A–D as well as the cubic reference pattern (bottom). Very little secondary phases could be found, (*) indicates LiO2 in sample A. |
To quantify the impedance spectrum properly, a resistor in parallel to a constant phase element (R1‖CPE1) is used for the bulk contribution, in series to a constant phase element CPE2, describing the partially blocking electrode. The impedance of a constant phase element includes the two fit parameters Q and n as well as the angular frequency ω according to ZCPE = Q−1(iω)−n. The inductance due to wiring was subtracted from all impedance data. The equivalent circuit is also shown in Fig. 4 and leads to a reliable fit (dashed line) of the impedance spectra. The fit parameters of the spectra in Fig. 4 are given in Table 2. The resulting relative permittivities (εr), calculated from CPE1 in accordance with ref. 58, are in the range of 45 to 78. Such a value confirms the bulk-type character of this part of the spectrum.
R 1, Ω | σ macro, mS cm−1 | Q 1, pF sn−1 | n 1 | C 1, pF | ε r | |
---|---|---|---|---|---|---|
A | 5608 | 0.21 | 31.8 | 0.88 | 3.82 | 52 |
B | 4045 | 0.32 | 8.96 | 0.97 | 5.27 | 78 |
C | 1352 | 0.16 | 336 | 0.83 | 18.0 | 45 |
D | 1112 | 0.48 | 248 | 0.82 | 11.5 | 70 |
From the resistance R1, the effective bulk conductivity σmacro can be calculated using
![]() | (1) |
![]() | (2) |
On each sample an array of circular electrodes with a diameter of 100 μm was applied. By subsequently contacting individual microelectrodes and measuring the impedance spectra between a microelectrode and a macroscopic counter electrode we are thus able to resolve local conductivity variations on the length scale of a few 100 μm. However, one has to keep in mind that these microelectrode measurements still average over a certain (small) sample volume.
Fig. 5 displays typical impedance spectra of such microelectrode measurements for all four samples. The impedance spectra are analyzed using the simple equivalent circuit shown in this figure. The resistive element RSpread describes the ionic charge transport in the probed LLZO volume. The serial constant phase element (CPE2) is attributed to the ionically partially blocking electrode material, as for the macroelectrode measurement in Fig. 4. Parallel to these two elements is a stray capacitance (CPE1) with a value in the range of 200 fF. This is caused by the measurement setup, particularly by the capacitance between the contacting needle and the counter electrode. Actually, such a stray capacitance is also present in macroscopic measurements, but there it can be neglected compared to the other capacitances and thus it is not needed in the circuit of Fig. 4. On the other hand, an additional constant phase element in parallel to RSpread, reflecting the very small LLZO bulk capacitance (in the 100 fF range), leads to an over-parameterization of the microelectrode equivalent circuit and was therefore avoided. More details on the appropriateness of this approach to analyze microelectrode measurements on LLZO are given in ref. 66. The obviously less pronounced separation between the electrode spike and bulk arc in microelectrode measurements (cf.Fig. 4 and 5) is simply caused by the additional stray capacitance (CPE1) and the corresponding peak frequency shift of the bulk arc.
![]() | ||
Fig. 5 Exemplary impedance spectra obtained with microelectrodes of diameter d = 100 μm, for all samples. The corresponding fits (dashed line) are based on the equivalent circuit shown above. |
This demonstrates that parts of the sample exhibit higher conductivities than the effective σmacro (2.1 × 10−4 S cm−1) while others are less conductive. The corresponding electrical conductivity variation might also be the reason behind the reduced n-value (0.88) of the macroscopic experiment. The existence of a gradient is also in agreement with a previous study, where reducing the size of LLZO samples and repeatedly measuring σmacro also revealed conductivity variations.46 Measurements from top to the bottom of the sample (vertical coloured bar in Fig. 6) indicate further conductivity fluctuations but not such a clear trend as found for the horizontal measurements.
Similar local conductivity measurements were done on sample B (Fig. 7). This sample was covered with differently sized microelectrodes (10 to 300 μm), but for the sake of comparability, only electrodes with a diameter of 100 μm were investigated. On average the local conductivities are even slightly higher than σmacro = 3.2 × 10−4 S cm−1, but overall little conductivity variations and good agreement with the macro measurements are found. The rather good agreement of microelectrode and macroelectrode experiments in this sample further supports the interpretation of microelectrode measurements measuring bulk related properties. Hence, we assume that the conductivity gradients observed for sample A truly reflect bulk conductivity inhomogeneities. The rather homogeneous conductivity of sample B might also cause the high n-value of the CPE (0.97, see Table 2).
![]() | ||
Fig. 7 Local conductivity measurements obtained from sample B using 100 μm electrodes. The local conductivity does not show very pronounced variations. |
The macroelectrode measurement of sample C revealed a lower effective conductivity compared to samples A and B (1.6 × 10−4 S cm−1). Microelectrode measurements were performed on both sides of this sample. On side one (Fig. 8 – diamond), measured conductivities are either similar (many spots) or smaller (some spots) than σmacro but did not show any spatial trend. Compared to sample B much more scattering is present. The measurements on side two (Fig. 8 – circles) resulted in very similar conductivities, again some spots were less conductive, the majority, however, was close to the obtained macroelectrode value. Accordingly, we find a very reasonable agreement between the effective macroscopic conductivity (probing the entire depth) and typical local conductivities, probing bulk regions only down to about 200 μm depth. Thus a comparatively low overall ionic conductivity seems to be present across the entire sample and the low σmacro value is not caused by a low conductive layer close to the surface. Such a layer would have strongly affected the local measurements (probing less depth), which is not seen here.
![]() | ||
Fig. 8 Microelectrode measurements on sample C from side one (diamonds) and side two (circles) of the sample. |
Microelectrode measurements on sample D resemble those on sample B with little spatial variation of the local conductivity (not shown). However, the average local conductivity in the investigated sample parts is in the range of 2.0 × 10−4 S cm−1 and thus smaller than the effective macroscopic conductivity, suggesting that at least some conductivity inhomogeneity is present somewhere in the sample. The latter is also in accordance with the lowest n-value of the CPE (0.82).
In order to investigate the correlation of local conductivity and local composition it is essential to ensure that both local measurements refer to the same sample region. Microelectrode measurements using d = 100 μm probe spatial regions in the same range; as an estimate, often a hemisphere with radius 2d is used as the relevant volume.50 The probed sample volume of LA-ICP-OES measurements strongly depends on the analyzed material, utilized laser and chosen laser settings, however, it is typically significantly smaller than that. While analyzed areas (defined by laser beam size) with a diameter of a few hundred μm lie in the accessible range of commonly used instruments, the probed depth hardly surpasses more than a few μm.
As a first step we have to avoid chemical analysis of near surface layers on LLZO; they are well-known to exhibit different chemical compositions and may include different phases.27,29,33–35 Also Li+/H+ exchange caused by water from the air has been reported and should be taken into account.23,25–27,67 In ref. 27 proton exchanged surface regions with thicknesses between 1 and 2 μm are reported. Effects of such surface regions are avoided by two measures: first, before the first chemical analysis by LA-ICP-OES ca. 5 μm LLZO is mechanically removed by polishing (grinding), see the Experimental. Second, any further reaction products resulting from the contact with ambient air after grinding are removed by a pre-ablation step performed directly before the measurement (ca. 2 μm). Hence, already the very first analysis takes place in a depth of several μm (represented by estimated 5 μm in the following) and thus reflects the bulk composition in this depth rather than any surface effect. However, still the composition in 5 μm depth might differ from the average composition in the ca. 200 μm range probed by the microelectrodes. In ref. 34, for example, chemical composition changes in the first few μm of LLZO ceramics were found and characterized by LIBS. Therefore, by further grinding the LLZO samples, regions in 50 μm and 100 μm depth (with respect to the original surface) were investigated and we could thus obtain a kind of (rudimentary) 3D compositional map of some samples. The results of this depth dependent analysis are discussed in the following.
The result of the local chemical analysis on sample A, obtained by performing LA-ICP-OES in a depth of ca. 5 μm, is shown in Fig. 9 for Al (a), Li (e) and Zr (g). The measurement data are not smoothened, each pixel represents the stoichiometry, calculated based on the measured atomic ratios (choosing 3 pfu La as the fixed-point).
Most obvious is the gradient in Al content from left top to right bottom, with values varying by about a factor of two from 0.15 pfu to 0.30 pfu. In accordance with the common model of Al reducing the Li content, this is accompanied by a (scattered) Li gradient with the highest Li-content in the Al poor region. The Zr signal shows very little variation and a homogeneous distribution across the sample, as well as good agreement with the theoretical Zr content of 2.0 pfu. A simple explanation for these Al and Li variations cannot be given yet. During sintering one edge of the sample might have been exposed to a lower local lithium oxide activity in the gas phase. Thus, stronger Li loss might have taken place, which in turn might create a driving force for Al diffusion.
The same experiment was repeated in 50 μm and 100 μm depth and the distribution images of Al are also shown in Fig. 9b and c; Li distributions are given in the ESI.† It is very obvious that the Al gradient is also present in this depth thus indicating the absence of strong compositional variation perpendicular to the surface within about 100 μm. This is confirmed quantitatively by comparing the Al and Li content in the cross-sections (white lines) of the three distribution images, see Fig. 9d and f. In these cross-sections, averaging over six neighboring pixels was performed to better reveal overall trends. Clearly the Al stoichiometry variation is present in all depths.
Also samples B, C and D were chemically analyzed by LA-ICP-OES, for D again in three depths, for C in two depths, and B in a depth of 5 μm only. The Al distribution images of all three samples are shown in Fig. 10a–c; further Al and all Li distribution images are shown in the ESI.† First of all, also the measurements in different depths confirmed that the data found for ca. 5 μm depth already correspond to the bulk composition in much deeper regions. Hence, we can conclude that the chemical composition found by our LA-ICP-OES measurements represents the chemical composition of the local region probed by microelectrodes, irrespective of the depth used for analysis. It is thus meaningful to use these data for looking at correlations between local chemical composition and local conductivity, with “local” meaning average values in regions with a size of ca. 100 to 200 μm in each direction.
Lateral chemical variations of Al are also present in samples B and C (Fig. 10) with B showing an Al enrichment in the centre part of the pellet and C exhibiting a “scattered” minimum in the Al content somewhere in the centre. Sample D was the most homogeneous sample with only moderate scatter of Li and Al contents. In these samples, a correlation between the Al content and Li content is not visible. This, however, might be due to the fact that a very high precision of the Li determination is required to see the small relative changes of the Li content caused by the Al variations. Moreover, lithium oxide loss during sintering with formation of oxygen vacancies68 might affect the cation ratios.
From all measured LA-ICP-OES data of a given sample, i.e. by averaging laterally and for different depths, average chemical compositions can be calculated. They are shown in Table 3.
Li/La | Al/La | Zr/La | Composition | |
---|---|---|---|---|
A | 2.07 | 0.076 | 0.669 | Li6.20Al0.23La3Zr2.01O11.95 |
B | 2.16 | 0.069 | 0.637 | Li6.49Al0.21La3Zr1.91O11.86 |
C | 2.10 | 0.089 | 0.631 | Li6.30Al0.27La3Zr1.89O11.83 |
D | 2.13 | 0.080 | 0.682 | Li6.39Al0.24La3Zr1.92O12.15 |
In sample A the measured Al content was slightly higher and the measured Li content was lower than expected, leading to a total composition of Li6.20Al0.23La3Zr2.01O11.95 when calculating the oxygen content for fixed cation charges. A lower amount of lithium is not too surprising, due to its volatility during sintering. The behaviour of aluminium is less straightforward to interpret. Possibly a gradient in Al not only exists laterally, on the mm scale (cf.Fig. 9), but also across the entire sample with somewhat less Al in the center parts beyond the 100 μm depth probed by LA-ICP-OES. Even though the sample pellets were placed between two additional pellets of pure Li7La3Zr2O12 during sintering to avoid contamination, also incorporation of additional Al3+ from alumina crucible cannot be ruled out. Sample B–D show slightly increased Al content as well, supporting this hypothesis.
While sample C shows a lower Li content than expected, like sample A, the measured Li amounts of sample D and B are close to the nominal value or higher than expected, indicating that different amounts of lithium oxide were lost during sintering.
Fig. 12 displays plots of conductivity versus composition (Li or Al) for the other samples. Samples B and D show only moderate composition variation with Al as well as Li content in a similar range. Also the conductivity variations are comparatively small, however, sample B shows significantly higher local conductivities (about 3.5 × 10−4 S cm−1) compared to sample D (about 2.0 × 10−4 S cm−1). Sample C, finally, includes a substantial variation of the chemical composition while staying at relatively low conductivity. Hence a broad range of Al and Li content leads to very similar conductivity values.
From all these plots it also becomes clear that for the same Al content the conductivity can vary by almost an order of magnitude. The same is true for the Li content. This is exemplified by the conductivity and composition data of four different spots on sample A (Table 4). Two of them exhibit a rather high Li amount (6.35 and 6.37 Li pfu), but conductivities vary from rather low (0.56 × 10−4 S cm−1) to acceptably high (2.42 × 10−4 S cm−1). On the other hand, if the amount of lithium is pretty low (6.12 to 6.13 pfu) the ionic conductivity can still be either acceptably high (2.54 × 10−4 S cm−1) or more than three times lower. Accordingly, there seems to be no optimal Li content (within the range considered here, i.e. between ca. 6.0 and 6.7 pfu) in contrast to suggestions in the literature.
Spots | σ Me, S cm−1 | Composition |
---|---|---|
Li enriched 1 | 2.42 × 10−4 | Li6.37Al0.20La3Zr2.02O12.04 |
Li enriched 2 | 0.56 × 10−4 | Li6.35Al0.18La3Zr2.03O12.03 |
Li poor 1 | 2.54 × 10−4 | Li6.13Al0.28La3Zr2.01O12.01 |
Li poor 2 | 0.79 × 10−4 | Li6.12Al0.26La3Zr2.02O12.00 |
Overall, we thus conclude that there is no clear and simple relationship between ionic conductivity and Al or Li content in Al stabilized LLZO. Hence, other aspects have to be of high relevance for the Li ion conductivity that are not simply reflected by the concentration of the stabilizing cation (Al) or the absolute Li content. Clarifying the true reasons behind varying conductivities requires detailed further measurements and is beyond the scope of this paper. We suggest two possible reasons: first, Li ion conduction in LLZO is a complex phenomenon with different cation sites being involved.37–40,69 The exact site occupancy in a LLZO sample, and thus also the effective mobility of Li, may depend on the preparation and might exhibit gradients within samples. They are not accessible by mapping local compositions only. Accordingly, the detailed local structure might be more important for the ionic conductivity than the exact Al or Li content.
Second, existence of oxygen vacancies and their effect on the lithium ion conduction might be another important aspect to consider. Recent tracer exchange experiments with subsequent secondary ion mass spectrometry (SIMS) analysis revealed the presence of possibly substantial amounts of oxygen vacancies in various LLZO samples.68 Also the results of our study indicate understoichiometric amounts of oxygen: calculating the average oxygen content of our LLZO samples from the measured amounts of cations with fixed charge leads to oxygen contents <12 for most samples, cf.Table 3. These oxygen vacancies may also affect the local occupancy of cation sites as well as activation barriers for Li migration.
Hence, a deeper understanding of the Li ion conductivity in Al stabilized LLZO might require much more information than the Al content; rather, detailed information on oxygen vacancy concentrations and local site occupancies might be key for in-depth knowledge.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c9ta00356h |
This journal is © The Royal Society of Chemistry 2019 |