research papers
Local structural studies of the cubic Cd1–xCaxO system through Cd K-edge extended X-ray absorption spectroscopic studies
aSurface Physics Division, Saha Institute of Nuclear Physics, Kalkata 700064, India, bCondensed Matter Physics Division, Materials Science Group, Indira Gandhi Centre for Atomic Research, Kalpakkam 603102, India, and cPhoton Factory, Institute of Materials Structure Science, High Energy Accelerator Research Organization, Oho 1-1, Tsukuba 305-0801, Japan
*Correspondence e-mail: varadu@igcar.gov.in
Cd K-edge extended X-ray absorption fine-structure spectroscopic studies were carried out on Cd1–xCaxO (0 ≤ x ≤0.9) solid solutions and the first and second nearest neighbour (NN) distances and their mean square relative displacement σ2 were estimated. The first NN distance, dCd–O(x), was found to be smaller than its expected value, a(x)/2, obtained from the X-ray diffraction measurements. It increases monotonically and non-linearly with a negative curvature, comparable with that of the a(x) value variation. The variation σ2 of the 1NN with x is consistent with a disordered model. The 2NN distances dCd–Cd(x) and dCd–Ca(x) are found to follow the average values obtained by X-ray diffraction with dCd–Ca(x) > dCd–Cd(x). From detailed analysis it is argued that the exhibits a of the 1NN distances, dCd–O(x) and dCa–O(x), and that the system belongs to a persistent type.
Keywords: EXAFS; CdO; CaO; ternary solid solution; oxide semiconductors.
1. Introduction
The effect of local atomic rearrangements on the properties of semiconductor alloys has been emphasized recently on the basis of both experimental and theoretical results (Kuzmin et al., 1993; Tsai et al., 1990; Aldrich et al., 1994; Koteski et al., 2004). Such studies received further impetus with the availability of synchrotron radiation facilities to carry out extended X-ray absorption fine-structure (EXAFS) (Stern, 1974; Lytle et al., 1975; Stern et al., 1975; Bunker, 2010) and atomic pair distribution function (PDF) (Egami & Billinge, 2003) studies. These studies confirm that the local atomic arrangements show substantial deviation from the average structure, though the translational symmetry prevails over much larger length scales (Kuzmin et al., 1995; Lebedev et al., 2001). The effects of such deviation include the reduced point-group symmetry, lattice distortions and the modification in the bonding. These in turn affect the physical, optical and electronic properties, which are of fundamental importance for the semiconducting materials. Hence, detailed local structural studies have become of importance to understand not only the average structural changes on alloying but also optical and electronic properties.
Recently we reported a wide and tunable band gap II–VI semiconducting oxide Cd1–xCaxO system. Although the end members CdO and CaO crystallize in the NaCl structure with a lattice parameter mismatch of 2.5%, the lattice parameter (a) variation of their deviates from a linear variation (Vegard's law). The variation exhibits a negative curvature characterized by a bowing parameter b = 0.03088 (Srihari et al., 2011). Employing a Williamson–Hall type plot (Williamson & Hall, 1953) obtained from the profile analysis of powder X-ray diffraction (XRD), the strain as a function of Ca substitution was calculated. It is found that the strain is negligibly small, a generic feature of the oxides. This implies the absence of strain as a field over the coherently diffracting volume. This clearly indicates that substitution results in local structure deviation from the average cubic symmetry, warranting a detailed structural study at the atomic level. spectroscopy is a widely used tool for such studies and is preferred over PDF studies since it selectively probes the local structure around an absorbing species. In this paper we report Cd K-edge studies on the Cd1–xCaxO system. Our studies indicate that the first nearest neighbour (1NN) distance dCd–O(x) is smaller than a(x)/2 and exhibits a negative deviation with a bowing parameter same as that of a(x)/2. On the other hand, the second nearest neighbour (2NN) distance dCd–Cd/Ca(x) closely follows the variation of a(x)/21/2. It is shown that the 1NN distances dCd–O(x) and dCa–O(x) are different and the optical properties of this system would exhibit characteristics of a persistence mode system.
2. Experimental details
Solid solutions of the Cd1–xCaxO system with x = 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8 and 0.9 were prepared by solid state reaction described in an earlier work (Srihari et al., 2011). Cd K-edge measurements on the Cd1–xCaxO system were carried out on beamline NW10A of the Photon Factory Advanced Ring (PF-AR), Tsukuba, Japan (Nomura, 2001; Nomura et al., 2007). The energy was varied from 500 eV below the Cd K-absorption edge, E0 (= 26711 eV), to 1100 eV above it using a Si (311) double-crystal monochromator (d = 1.63747 Å). Measurements were carried out on finely ground powders of the Cd1–xCaxO system. The optimum sample quantity for each composition was estimated such that the effective thickness of the sample corresponds to approximately one absorption length for that composition. The required quantity of a diluting medium, boron nitride powder, was thoroughly mixed with the sample and pressed under moderate pressure into pellets of diameter 10 mm for measurement. Measurements under ambient conditions were carried out in transmission mode using ionization chambers for I0 and It. Energy calibration was carried out using a standard Cd foil.
3. Data processing and analysis
The absorption spectra were processed using the Athena software package (Ravel & Newville, 2005) for edge alignment and pre-edge treatment (−500 to −150 eV before E0) and background parameters were carefully chosen so that the data suffer the least from artefacts of such a processing. The post-edge treatment corresponding to removal of the K-edge absorption of free Cd2+ ions in the range 100–1100 eV above E0 was carried out to obtain the signal, χ(k). The χ(k) signal was then Fourier transformed (FT) in the range 2 Å−1 ≤ k ≤ 14.5 Å−1 to obtain the real, Re[χ(R)], imaginary, Im[χ(R)] and absolute magnitude, |χ(R)|, of χ(R). The obtained χ(R) was weighted by k3 and fitted to a structural model using Artemis software (Newville, 2001) to obtain local structural information; interatomic distances or the effective path length and mean square relative displacement, σ2, with reference to the absorbing Cd2+ ion. In the analysis, single-scattering (SS) and multiple-scattering (MS) paths with up to four legs within a cluster size Rmax = 6 Å were included. In the present study, the fitting was carried out in the range 2 Å ≤ Reff ≤ 6 Å with the Hanning window and dR = 1 Å; Reff is the effective path length. Although our primary interest is limited to the first and second NN, paths with higher Reff need to be included since they form `leakage' paths into the low-R range.
The following strategy was adopted to incorporate atom substitution, which is not a standard feature of the Artemis package. It is assumed that the substitution is random, a fair assumption for the NaCl-type structure. Additional Ca2+ backscattering paths were included such that for every path (SS or/and MS paths) involving the Cd2+ ion a corresponding identical path involving the Ca2+ ion is also present. The degeneracy N (coordination number in the case of SS paths or multiplicity in the case of MS paths) of such paths was proportioned between Cd2+ and Ca2+ paths conforming to the nominal composition of the sample: NCd = Nx; NCa = N(1 − x). This methodology provides a scope to extract the parameters of both the 2NN paths, Cd–Cd and Cd–Ca. In the case of random alloys it is known that the distribution of NN atoms around the associated average coordination shell is non-Gaussian (Crozier et al., 1988), especially so in the case of the Cd1–xCaxO system, wherein the difference between Cd and Ca is large. Hence, the cumulant expansion method, which does not assume any specific form for the distribution, was adopted to fit the χ(R) signal (Bunker, 2010; Crozier et al., 1988). Accordingly, the first (C1), second (C2), third (C3) and fourth (C4) cumulants correspond to the NN distance or effective (scattering) path length Reff, the mean square relative displacement (MSRD) σ2 of the NN atoms with reference to the absorbing atom, skewness in the distribution and its weightage, respectively. While fitting, correlated parameters were consistently floated and it was seen that the number of floating parameters is always less than the number of independent parameters Nind, as given by the Nyquist criteria (Stern, 1993).
Although we have considered 2–6 Å for obtaining structural parameters, the analysis was broken to overlapping windows in R-space. For a given fitting, only a few paths were included such that the total number of floated parameters obeys the Nyquist criteria. In all the fitting procedures the pairs of correlated parameters, S0 and σ2 and E0 and ΔR, floated simultaneously. Subsequently, this window range (in R) was sifted so as to include at least one path of the previous fitting session and its parameters were fixed. The entire process was iterated so that good fitting and meaningful errors in the floated parameter were obtained.
4. Results and discussion
Fig. 1 shows the Cd K-edge absorption spectra of the Cd1–xCaxO The spectra did not show perceptible changes in the XANES range (26611–26775 eV), except for a shift in the position of the first peak to higher energies with x. This shift implies changes in bond lengths and in the nature of bonding, corroborating our earlier inference deduced from the electron density distribution from the powder XRD data (Srihari et al., 2011). On the other hand, substantial and systematic changes are observed in the the k3-weighted signal χ(k) (Fig. 2), especially in the range 4 Å−1 ≤ k ≤ 12.5 Å−1, indicating changes in the environment of Cd with Ca substitution. In passing, it is remarked that the absence of noise in k3χ(k), even up to 16 Å−1, indicates that the samples are free of pin holes. The absolute magnitude of χ(R) of the in the range 1–6 Å is shown in Fig. 3. It exhibits several peaks and substructures: a peak at ∼2.2 Å corresponding to the 1NN distance, viz. dCd–O ≃ a/2, estimated from powder XRD studies (Srihari et al., 2011). The shift in the position of this peak to higher R values with x indicates an increase of dCd–O(x) with x, conforming to the lattice dilation with Ca substitution. The second peak is centered at about ∼3.37 Å and exhibits a shoulder at R ≃ 3.05 Å; the shoulder to the second peak in CdO may be due to non-linearity in backscattered photoelectron intensity from the Cd atom. Although the position of the second peak corresponds approximately to the 2NN distance, viz. dCd–Cd/Ca ≃ a/21/2, as will be seen below, the paths with still higher effective scattering length also contribute to this peak. Figs. 4 and 5 show a structural model fit to the real part, imaginary part and absolute magnitude of χ(R), respectively, for x = 0.3. Contributions from a few representative scattering paths up to 1–3NN are shown. It is seen from these figures that the 1NN signal is not contaminated much by higher-order paths. The contribution of the 2NN, Cd–Cd and Cd–Ca paths are not in phase [Figs. 4(a) and 4(b)] and their contributions are located at about 3.37 Å and 3.05 Å, respectively, in |χ(R)|. The contribution of the Cd–Ca path to the second peak also falls near to the shoulder of the Cd–Cd contribution. From Fig. 5 it is observed that the Cd–Ca contribution to the second peak is at a lower R value compared with that of the Cd–Cd contribution; this may be due to the difference in the phase shift between the scattered electron wave from the Cd and Ca atoms. Also, it can be seen from Fig. 5 that contributions from the 3NN and higher-order paths strongly overlap with those of the 2NN contributions; e.g. contributions from the MS path, Cd(absorbing)–Cd(at face-centered positions)–O1(at edge-centered positions), and SS path, Cd–O2 (body-centered position). This necessitates one to include the paths with Reff = 6 Å during fitting. From the analysis of |χ(R)|, the NN distances dCd–O(x), dCd–Cd(x) and dCd–Cd(x) were estimated and are tabulated in Table 1 along with σ2, S02, E0 and the R-factor, a measure of the goodness of fit.
|
In Fig. 6 the variation of dCd–O(x) (red circles) is compared with the corresponding distance a(x)/2 (black squares). It is seen that dCd–O(x) monotonically increases and exhibits a negative deviation from a linear behavior. The variation of dCd–O(x) was fitted to a second-order polynomial equation (Srihari et al., 2011) (red line, `Quadratic Fit') and the bowing parameter was estimated to be dEXAFS = 0.026 Å. This value compares well with the bowing parameter for the lattice parameter variation, dXRD = 0.03088 Å (Srihari et al., 2011). In contrast, the 2NN distances, dCd–Cd(x) and dCd–Ca(x), increase almost linearly with x (Fig. 7) and are comparable with a(x)/21/2. It is also seen that dCd–Cd(x) is smaller than dCd–Ca(x) for all values of x, consistent with the larger ionic radius of Ca2+ compared with that of Cd2+. Variation of , a measure of both thermal and structural disorder, is shown in Fig. 8. The 1NN of undoped CdO itself has a finite value. For x = 0.5, the value of is found to be smaller even than that for the end members. Barring this, the overall variation of exhibits a more or less symmetric variation with a positive curvature having a maximum value of ∼0.0043 Å2 for x ≃ 0.5 (Fig. 8). On the other hand, the disorder associated with the 2NN is comparable, ≃ ≃ 0.01 Å. This is much higher compared with that of 1NN and did not exhibit any systematic variation with Ca substitution. The third cumulant of the 1NN distribution for the CdO (x = 0) is ∼2 × 10−4 Å, implying a small skewness. It increases with armchair structure to ∼3 × 10−4 Å for x = 0.9. While the variation in the skewness in the distribution of dCd–Ca(x) was found to be negligible with no systematic variation, the variation in the skewness in the distribution of dCd–Ca(x) exhibits a bell shape: it is either zero or negligibly small for x = 0, 0.1, 0.8 and 0.9 and has a maximum value of ∼6 × 10−5 Å3 for x = 0.5. Thus, the skewness in the 2NN distribution is estimated to be one order less compared with that of 1NN.
In the past, local structural studies through PDF/EXAFS have been carried out on semiconducting alloys, which can be broadly divided into two groups. For the alloys belonging to the first group, the variation of the lattice parameter follows Vegard's law: Ga1–xInxAs (Egami & Billinge, 2003), Hg1–xCdxTe (Pong et al., 1989), Rb1–xKxBr (Boyce & Mikkelsen, 1985), Th1–xUxO2 and Th1–xPuxO2 (Hubert et al., 2006). For the alloys belonging to the second group, the variation of the lattice parameter deviates from Vegard's law: Mg1–xNixO (Kuzmin et al., 1995), Mg1–xCoxO (Kuzmin et al., 1993) and Mg1–xFexO (Waychunas et al., 1994). In the first group of alloys the 1NN cation–anion distances (e.g. dGa–As and dIn–As for Ga1–xInxAs and dTh–O and dPu–O for Th1–xPuxO2) are either almost the same or are comparable with those of the respective end-member cation–anion distances or they vary linearly with a substantially small slope comparable with the lattice parameter variation. Among the members belonging to the second group, Kuzmin & Mironova (1998) have carried out detailed local structural investigations on Mg1–xNixO (space group ) employing Mg- and Ni-edge Although the lattice parameter of this system exhibits a negative deviation from Vegard's law, with bXRD ≃ 0.01 Å (Kuzmin & Mironova, 1998), both the 1NN distances dMg–O(x) and dNi–O(x) are reported to vary linearly. These authors have used the decrease of 2NN distances dMg–Mg(x) and dNi–Ni(x) to explain the lattice parameter variation. This is to be contrasted with the Cd1–xCaxO system under investigation. For the Cd1–xCaxO system the 1NN distance dCd–O(x) exhibits a negative deviation with a bowing parameter comparable with that of the lattice parameter variation. Considering that dCd–O(x) is always smaller than a(x)/2, Cd1–xCaxO mimics the behavior of the former group of alloys with the exception that both the lattice parameter and the dCd–O(x) variation deviates from Vegard's law. Compared with semiconducting alloys with wurzite or chalcopyrate structure, the difference between a(x)/2 and dCd–O(x) for the Cd1–xCaxO system is smaller, as has been observed for purely ionic systems like K1–xRbxBr (Boyce & Mikkelsen, 1985). Since the wurzite structure allows larger changes in the bond angles rather than in the bond lengths, to minimize the strain, the end-members' 1NN distances are more or less preserved, leading to a larger difference between a(x)/2 and 1NN distances. On the other hand, the NaCl-type structure with the to which the Cd1–xCaxO system belongs, does not allow for large bond angle variations and results in larger changes in the bond lengths. This leads to a smaller difference between a(x)/2 and dCd–O(x) compared with that of the systems belonging to the wurzite structure. Additionally, for systems where the lattice parameter variation obeys Vegard's law, the Phillips ionicity (Phillips, 1970, 1973) of the end members are closely matched (e.g. 0.31 for GaAs and 0.357 for InAs) and changes in the nature of the chemical bonding are not expected. The Phillips ionicity of CdO and CaO are widely different: 0.785 for CdO and 0.913 for CaO. Alloying of such end members brings out change in the nature of the chemical bonding; from more covalent-like for smaller values of Ca substitution to ionic bonding for higher Ca substitution (Srihari et al., 2011). This, coupled with the rather rigid cubic symmetry, is possibly responsible for the observed bowing in the variation of dCd–O(x) and hence in a(x).
To quantitatively explain the variation of dCd–O(x), it is necessary to consider the neighborhood of Ca. Owing to the proximity of the Cd LI-edge (4018 eV) to the Ca K-edge (4038.5 eV), Ca L-edge could not be carried out and dCa–O(x) could not be estimated. However, the local structure averaged over a much larger length scale should conform to the as obtained from the diffraction studies. In other words, the sum of the compositional weighted 1NN distances of Cd and Ca should give the lattice parameter as (Azoulay et al., 1982)
In the above equation, the fitted values of a(x) and dCd–O(x) were used to estimate dCd–O(x), and its variation is plotted as the solid blue line (`expected dCd–O') in Fig. 6. Thus the studies indicate a for 1NN distances.
As seen from Fig. 7, the variation of the 2NN distance dCd–Cd/Ca(x) closely follows that of the average 2NN distance a(x)/2, in contrast to the 1NN distance dCd–O(x) variation. Such a behavior has been reported for both purely covalent and ionic systems like the Ga1–xInxAs (Mikkelsen & Boyce, 1983) and K1–xRbxBr (Boyce & Mikkelsen, 1985) systems, respectively. This is explained by Boyce & Mikkelsen (1985) by invoking a radial force model, wherein bond-stretching radial forces alone are considered while weak bond-bending forces are omitted (Shih et al., 1985). Alternatively, one should also take into account (x) to explain this (Frenkel et al., 1993). Upon lattice dilation, dCd–Cd/Ca(x), as expected, increases by 0.06 Å, from 3.3292 Å (x = 0) to 3.3892 Å (x = 0.9). However, no appreciable changes in (x) are observed indicating that the `packing' around 2NN Cd/Ca is preserved over the entire composition range, as has been reported for the K1–xRbxBr system (Frenkel et al., 1993). Hence, an increase in dCd–Cd/Ca(x) is interpreted as an increase in the size of the ionic species with x. With this in the background, variations in (x) clearly indicate that `packing' around 1NN O is rather loose and it undergoes substantial local structural deviation to accommodate the substitution of the larger Ca2+ for the smaller Cd2+.
5. Conclusions
Cd K-edge measurements at room temperature were carried out on Cd1–xCaxO (0 ≤ x ≤ 0.9). Employing the cumulant expansion method, the first [dCd–O(x)] and second nearest neighbour distances [dCd–Cd(x) and dCd–Ca(x)] and mean square relative displacement, σ2, were estimated. It was found that the dCd–O(x) variation exhibits a negative deviation from linearity with a curvature quite close to that of the lattice parameter variation a(x) and is smaller than a(x)/2. By analyzing the variation of dCd–O(x) and a(x), it is inferred that the 1NN distance for Ca, dCa–O(x), is larger than dCd–O(x), implying a for the nearest neighbour distances. From this, it is expected that the optical properties of the Cd1–xCaxO system should conform to a persistent type system. From the linear increase in 2NN distances dCd–Cd(x) and dCd–Ca(x), with the associated σ2 being almost constant, it is reasoned that the ionic sizes of the species are concentration-dependent.
Acknowledgements
The authors would like to acknowledge the NW10A beam staff and PF staff for their support in carrying out the
measurement. This work was undertaken by financial support from the Department of Science and Technology, India, under the project `Utilization of International Synchrotron Radiation and Neutron Scattering Facilities'. V. Srihari would like to acknowledge financial support from the Department of Atomic Energy.References
Aldrich, D. B., Nemanich, R. J. & Sayers, D. E. (1994). Phys. Rev. B, 50, 15026–15033. CrossRef CAS Web of Science Google Scholar
Azoulay, J., Stern, E. A., Shaltiel, D. & Grayevski, A. (1982). Phys. Rev. 25, 5627–5630. CrossRef CAS Google Scholar
Boyce, J. B. & Mikkelsen, J. C. Jr (1985). Phys. Rev. B, 31, 6903–6905. CrossRef CAS Web of Science Google Scholar
Bunker, G. (2010). Introduction to XAFS: A Practical Guide to X-ray Absorption Fine Structure Spectroscopy. Cambridge University Press. Google Scholar
Crozier, E. D., Rehr, J. J. & Ingalls, R. (1988). X-ray Absorption, edited by D. C. Koningsberger and R. Prins, pp. 373–442. New York: John Wiley and Sons. Google Scholar
Egami, T. & Billinge, S. J. L. (2003). Underneath the Bragg Peaks: Structural Analysis of Complex Materials. Oxford: Pergamon. Google Scholar
Frenkel, A., Stren, E. A., Voronel, A., Qian, M. & Newville, M. (1993). Phys. Rev. Lett. 72, 3485–3488. CrossRef Web of Science Google Scholar
Hubert, S., Purans, J., Heisbourg, G., Moisy, P. & Dacheux, N. (2006). Inorg. Chem. 45, 3887–3894. Web of Science CrossRef PubMed CAS Google Scholar
Koteski, V., Haas, H., Holub-Krappe, E., Ivanovic, N. & Mahnke, H.-E. (2004). J. Alloys Compd. 371, 138–141. Web of Science CrossRef CAS Google Scholar
Kuzmin, A. & Mironova, N. (1998). J. Phys. Condens. Matter, 10, 7937–7944. Web of Science CrossRef CAS Google Scholar
Kuzmin, A., Mironova, N., Purans, J. & Rodionov, A. (1995). J. Phys. Condens. Matter, 7, 9357–9368. CrossRef CAS Web of Science Google Scholar
Kuzmin, A., Mironova, N., Purans, J. & Sazonov, A. (1993). Phys. Status Solidi A, 135, 133–141. CrossRef CAS Web of Science Google Scholar
Lebedev, A., Sluchinskaya, I. & Munro, I. (2001). J. Synchrotron Rad. 8, 800–802. Web of Science CrossRef CAS IUCr Journals Google Scholar
Lytle, F. W., Sayers, D. E. & Stern, E. A. (1975). Phys. Rev. B, 11, 4825–4835. CrossRef CAS Web of Science Google Scholar
Mikkelsen Jr, J. C. & Boyce, J. B. (1983). Phys. Rev. B, 28, 7130–7140. Google Scholar
Newville, M. (2001). J. Synchrotron Rad. 8, 96–100. Web of Science CrossRef CAS IUCr Journals Google Scholar
Nomura, M. (2001). Handbook of EXAFS Experimental Stations. KEK Internal 2001–5, KEK, Japan. Google Scholar
Nomura, M., Yuichiro, K., Masato, S., Atsushi, K., Yasuhiro, I. & Kiyotaka, A. (2007). AIP Conf. Proc. 882, 896–898. CrossRef CAS Google Scholar
Phillips, J. C. (1970). Rev. Mod. Phys. 42, 317–356. CrossRef CAS Google Scholar
Phillips, J. C. (1973). Bonds and Bands in Semiconductors. New York: Academic Press. Google Scholar
Pong, W. F., Mayanovic, R. A. & Bunker, B. A. (1989). Physica B, 158, 617–618. CrossRef CAS Web of Science Google Scholar
Ravel, B. & Newville, M. (2005). J. Synchrotron Rad. 12, 537–541. Web of Science CrossRef CAS IUCr Journals Google Scholar
Shih, C. K., Spicer, W. E., Harrison, W. A. & Sher, A. (1985). Phys. Rev. B, 31, 1139–1140. CrossRef CAS Web of Science Google Scholar
Srihari, V., Sridharan, V., Chandra, S., Sastry, V. S., Sahu, H. K. & Sundar, C. S. (2011). J. Appl. Phys. 109, 013510. Web of Science CrossRef Google Scholar
Stern, E. A. (1974). Phys. Rev. B, 10, 3027–3037. CrossRef CAS Web of Science Google Scholar
Stern, E. A. (1993). Phys. Rev. B, 48, 9825–9827. CrossRef CAS Web of Science Google Scholar
Stern, E. A., Sayers, D. E. & Lytle, F. W. (1975). Phys. Rev. B, 11, 4836–4846. CrossRef CAS Web of Science Google Scholar
Tsai, M.-H., Dow, J. D. & Newman, K. E. (1990). Phys. Rev. B, 41, 7744–7748. CrossRef CAS Web of Science Google Scholar
Waychunas, G. A., Dollase, W. A. & Ross, C. R. II (1994). Am. Miner. 79, 274. Google Scholar
Williamson, G. K. & Hall, W. H. (1953). Acta Metall. 1, 22–31. CrossRef CAS Web of Science Google Scholar
This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.