research papers
High-Focusing Bragg-Crystal Polychromator Design for Energy-Dispersive X-ray Absorption Spectroscopy
aEuropean Synchrotron Radiation Facility, BP 220, F-38043 Grenoble CEDEX, France, and bLaboratoire de Magnétisme Louis Néel, CNRS, BP 166, F-38042 Grenoble CEDEX, France
*Correspondence e-mail: sanmigue@dpm.univ-lyon1.fr
The design of a highly focusing profiled Si(111) Bragg crystal polychromator for dispersive
(extended X-ray absorption fine structure) experiments is described. The contour of the crystal has been optimized to give the best focus over the full 4–13 keV energy range. The profile optimization has been improved taking into account all the of the geometry and the results of X-ray tracing simulations. The profile of the crystal has been calculated to take full advantage of the new possibilities given by the undulator source and the optics of beamline ID24 at the European Synchrotron Radiation Facility. Full spot sizes have been measured to be between 20 and 40 µm in the 5–12 keV energy range. These values compare well with X-ray tracing simulations and are the smallest spots ever obtained with energy-dispersive optics, keeping, however, a wide enough energy bandpass for most X-ray absorption spectroscopy experiments.Keywords: EXAFS; focusing; polychromators; X-ray optics.
1. Introduction
X-ray absorption in dispersive mode (DXAS) collects in parallel all the data of ; Flank et al., 1983). This scheme is based on dispersive optics, i.e. a bent perfect crystal which reflects a polychromatic beam out of a white source. After being focused at the sample position, the photons are collected by a position-sensitive detector able to work under high-flux conditions (Tolentino et al., 1988, 1990; Koch et al., 1996). The energy–direction correlation given by the crystal polychromator is transformed into an energy–position correlation at the detector position. Since all the data of a spectrum are acquired simultaneously, the time resolution achievable is in the millisecond range. In addition, since no mechanical movement takes place during the data collection, the stability of this experimental set-up is excellent. These characteristics of the spectrometer are of interest mainly for time-resolved and X-ray experiments. High-pressure experiments using diamond anvil cells also take advantage of the high quality of the focus of the DXAS instrument (Itié et al., 1992). In that case the small size and the of the focus are essential parameters. Consequently, such a focus is required for aberration-free optics. At other synchrotron facilities, up to now, the best achievement of focus sizes on a DXAS spectrometer range from 500 to 100 µm. The reduction of the focus size with the optics described in this paper has allowed the study of smaller samples and consequently opens up new types of experiments. In particular, we can cite the first experiment over megabar pressures (San Miguel, Polian et al., 1998), the first results on on single crystals at high pressure (Itié et al., 1997), the development of nanosecond-resolved XMCD experiments (Bonfim et al., 1997; Pizzini et al., 1998) or an energy-dispersive microXANES experiment (Mosbah et al., 1998). Examples of reference data are given elsewhere (Hagelstein et al., 1997).
(X-ray absorption spectroscopy) spectra over a large spectral range (Matsushita & Phyzackerly, 1981This paper describes the design and the commissioning of a Bragg-contoured polychromator crystal which achieves a focus size between 20 and 40 µm over the 5–12 keV energy range.
In DXAS experiments the crystal polychromator can be set either in transmission geometry (Laue case) or in reflection geometry (Bragg case). In our energetic range, as shown by Hagelstein et al. (1995), the Bragg geometry is much more convenient.
The most commonly used technique to achieve a desired bending of the crystal polychromator in DXAS experiments has been to apply a force to the apex of a triangular-shaped crystal whose base is clamped (Tolentino et al., 1988). The resulting cylindrical shape differs from the ideal elliptical profile (Tolentino et al., 1990) and it gives an aberration proportional to the square of the illuminated length of the crystal (Ice & Sparks, 1984). To achieve the elliptical profile the curvature has to be increased from the base to the apex of the crystal. This is accomplished by tailoring its width, which no longer varies linearly (Tolentino et al., 1990). However, as one end of the crystal must remain fixed, a different contour is needed, ideally for each energy. This may be too expensive and not too practical. An improvement of this method yields additional It consists of using contoured rectangular crystals mounted on four-point benders, which permit almost independent adjustment of two moments (i.e. two curvatures) at two points near each end of the crystal (Allen et al., 1993). This is the principle that we have applied. In addition, we have performed more detailed simulations in order to optimize the contour of the crystal to average the minimization of the spot size over the whole 4–13 keV energy domain.
Therefore we computed the crystal elastic behaviour and performed X-ray tracing simulations in order to find the more effective contour to produce the best focus in the 4–13 keV energy domain. Consistent with the theoretical investigation, the smallest spot ever obtained in an energy-dispersive
spectrometer has been achieved. The optimized bender and the profiling procedure of the polychromator using X-ray tracing is described in the following sections before the presentation of the results of the experimental tests.2. Theory
To decide the contour of a crystal suited to yield ellipse-induced focusing, we need to choose the medium diffracted energy (Eo) and to set the appropriate geometrical factors of the spectrometer.
To find the best-fitted contour for a crystal polychromator working in the 4–13 keV energy domain, we proceed via two steps. First, some geometrical parameters are chosen and the appropriate contoured crystals are calculated for different energies (optimization energies, Eo). Second, the focusing performance of each of these calculated crystals is simulated for other energies (hereafter called calculation energies, Ec) chosen over the 4–13 keV range. This is performed by looking for the best moments that must be applied to the crystal to simulate the new ellipse through a defined figure of merit. The geometrical parameters and the results are compared through this figure of merit. X-ray tracing is used to verify the correctness of the figure of merit and to try to predict the spot size at the focus.
2.1. Geometrical considerations
We define a reference frame (x,y) having its origin in the central point (P) of the bent crystal (Fig. 1). The X-ray source (S) and the sample position (I) define the ellipse that the symmetric crystal polychromator must follow in order to attain focusing aberration-free behaviour. In this reference frame, the shape of the ellipse can be expressed as (see, for example, Howells & Lunt, 1993)
where a2, a3 and a4 are geometrical coefficients depending on the source-to-crystal distance (p), the crystal-to-focus distance (q) and the (θB) at the centre of the crystal. Equation (1) is a good approximation when the radius of curvature of the crystal is larger than its length, as in the present case. As the direction of diffraction of the X-rays depends on the angle between the incident rays and the planes of diffraction, the evaluation of the quality of approximation (1) is given by the differences between the slope derived from (1) and that of the exact formula of the ellipse (the figure of merit, defined later, is based on this calculation). This difference is less than 1%.
Keeping only the first term in (1), the basic equations of cylindrical optics remain, i.e.
where R is the radius of curvature of the crystal at P. Equation (2) defines the energy bandpass, given by
The use of an undulator source, as is the case in ID24, can introduce an additional limitation in the energy bandpass. In most of the cases, this can be avoided by tapering the undulator and/or using the even harmonics of the undulator source (see Hagelstein et al., 1997).
In the absence of optical aberrations, the energy resolution of the polychromator (Tolentino et al., 1988) depends on the spatial resolution of the detector, the size of the X-ray source, the of the incident X-rays into the bent crystal, and the intrinsic Darwin width of the rocking curve of the perfect crystal. The two latter factors deal with the polychromator optics. For a fixed curvature radius of the crystal, the only way to increase the energy resolution will be to increase the of the diffraction planes. Nevertheless, the ultimate limiting factor of the energy resolution in experiments is a given fraction of the core-hole width of the probed atom. A good compromise between energy resolution and energy bandpass, in the energy domain under consideration, is obtained with an Si(111) reflection. As explained by Tolentino et al. (1988), the contribution of the source is minimized if the crystal-to-sample distance, q, is approximately equal to the sample-to-detector distance, d. The resolution decreases slowly if d > q (but the decrease is negligible for small X-ray sources as the undulators at ESRF beamlines) and increases very fast if d < q. Consequently we will set d ≥ q.
Aberrations can contribute to losses in the energy resolution that can be of the order of the core-hole width. Then, any correction of the aberration should improve the energy–position correlation. It is easy to see from the size of the spot which we have measured with the profiled crystal that the remaining aberrations introduce losses of energy resolution smaller than the core-hole width.
2.2. Elastic behaviour and crystal profile
The elastic equation of a beam bent at two points is (see, for example, Landau & Lifschitz, 1986)
where EY is the Young modulus in the compression/dilatation direction, i.e. the x direction in Fig. 2(a). The right-hand side of (4) describes the distribution of internal moments along the beam. M0 and M1 are the internal moments at the points of application of the external forces, which are separated by a distance Lm. The curvature of the crystal is approximated by the second derivative of y with respect to x. I is the section moment around the neutral axis, and is given by I = wt3/12.
Equation (4) is valid for small flexures (the radius of curvature of the beam is large compared with its length) and if the external forces are concentrated at discrete points of the beam, as is the case with a four-point bender.
From (4) one realizes that an x-dependent moment, I(x), enables one to design the real variable curvature along the crystal close to the ideal one. This can be performed by varying either the thickness or, with better control, the width of the crystal. Such a width variation required to obtain a crystal with the elliptic curvature is given by
Any pair of moments M0 and M1 can be used to calculate the crystal profile. However, to move as little as possible away from the unprofiled crystal (rectangular), we choose, following Underwood (1977), the M0 and M1 values that approximate better the rectangular crystal to the ellipse. By introducing these moments in (5) we obtain the profile of the crystal (see Fig. 2b) optimized for the diffracting energy Eo of the central point P.
2.3. Evaluation of the crystal behaviour at new energies
This is the second step of the evaluation before executing the profile cutting. For each crystal contour given by the first iteration we calculate its performances at other energies Ec of the central point. Other energies imply other Bragg angles and consequently other ellipses rotated and translated with respect to the first one.
Under new moments, the crystal with the width given by (5) will take a shape given by (4). It is valuable to estimate how large the deviation is from the ideal conditions in order to obtain an energy range of acceptance. The solution of the differential equation gives rise to a system of equations that may be solved numerically.
To compare the ideal ellipse at energy Ec with the shape given by the curved crystal, a figure of merit, Σ, based on the slope comparison of both shapes, is required (since diffraction behaviour is our concern), where
Equation (6) must be evaluated taking into account that (i) the crystal must be rotated by an angle θB − θB′, (ii) the derivatives must be taken in the appropriate reference system, and (iii) a careful choice of the points of comparison is requested (for approximated ellipses the X-rays coming from the source do not impinge on the crystal at the same point as in the case of an exact ellipse).
We look for the best curvatures at the calculation energy Ec to approach the ideal ellipse at this energy. To accomplish this in the simulation we minimize numerically the figure of merit.
The procedure that we have presented to calculate the profile at Eo features two iterative steps: firstly, moment determination, and secondly, profile optimization with a fixed moment distribution. Finally the figure of merit is derived by an additional adjustment of the moments imposed at the profile calculated in the second step. A second cycle of iteration could be envisaged. For three reasons it seems to us unreasonable to proceed further in that direction: (i) the computational accuracy is limited by the finite size of the sampling of the crystal; (ii) the initial calculation proceeds with approximations of the ideal ellipse and formula; (iii) the figure of merit which has been chosen is strongly correlated to the focus size, but other parameters contribute to it, i.e. the of X-rays in the crystal or the change of position of each individual reflecting point of the surface with respect to the ideal elliptical profile.
2.4. Choice of parameters
Here we consider the geometrical parameters for the ID24 DXAS beamline at the European Synchrotron Radiation Facility (ESRF, Grenoble, France). The source-to-crystal distance is p = 29.7 m. The crystal-to-sample distance, q, can be set between 0.7 and 1.5 m. We have the freedom to change q within this range for different photon energies between 4 and 13 keV. Nevertheless, the choice of q results from a compromise between the creation of a large enough energy bandpass and the need to avoid a too large bending of the crystal.
As yielded by (3), the spectral range limitation is critical at low energies. However, using an undulator as the X-ray source means that there is already a limited bandpass at low energy since the emission derives from the undulator first harmonic. At the ESRF, at low energies, a bandpass ΔE/E = 0.045 matches the undulator bandpass for a typical value of q = 0.8 m. For that q value the energy Eo = 7 keV produces a calculated crystal that optimizes the figure of merit, Σ, for the 4–13 keV range. Fig. 3 shows Σ for crystals contoured for different Eo values with the corresponding experimental q values. Fig. 3 also shows evidence of the presence of a second minima at an energy (E1) higher than the one (Eo) used for the optimization. The gap between Eo and E1 increases rapidly with Eo. The iteration process which is used for the profile evaluation is in fact an approximate route to reach the optimized parameters. Therefore, it may be possible that a profile derived from a second step of the iteration could provide us with a slightly better energy dependence of the figure of merit and possibly with a single minima. To attain that kind of accuracy is probably not possible given all the limits inherent to this calculation.
2.5. X-ray tracing
Before proceeding with the X-ray tracing we used the multilamellar crystal code in the XOP package (Sanchez del Rio et al., 1998) to check that, with the curvatures and the energy range used, no significant asymmetric broadening of the rocking curve occurs with the curved crystal as analysed by Hagelstein et al. (1995). The X-ray penetration at the V, Fe, Cu and Ge K-edges in the relevant experimental conditions are 2.9, 2.3, 3.5 and 4.3 µm, respectively (determined by the extinction length le). Therefore, we can make the simulation by considering the crystal as a mirror. The calculated spot sizes will be then underestimated by a quantity of the order of lecosθBsinθB ≃ 1 µm due to the finite extinction length.
The parameters of the X-ray source of the crystal are obtained by X-ray tracing of the previous elements of the beamline. The effective source observed by the polychromator of the ID24 beamline is located at the focus point of the second mirror of the Kirkpatrick–Baez optics (see Fig. 1 of Hagelstein et al., 1997). This secondary source can be conveniently described by a Gaussian distribution both for the spatial and angular coordinates with parameters σx = 100 µm, σz = 509 µm, σx′ = 410 µrad and σz′ = 16 µrad, where x is the demagnification direction and z is the vertical one. Table 1 summarizes the X-ray tracing results for both the profiled crystal optimized at 7 keV and the rectangular crystal (non-profiled crystal). The energies chosen for the simulations are the same as those used in the test experiments. The spot size is estimated as the spatial gap containing 90% of the X-rays to ensure we take into account properly the tails of the image spot. The correlation between the spot size and the figure of merit is 0.86.
|
3. Experience and discussion
For tests performed at the ESRF energy-dispersive Sc) was set with the help of two motorized and cooled copper slits placed just before the polychromator. A complete description of the mechanical part of the system is given elsewhere (San Miguel, Hagelstein et al., 1998).
beamline (ID24), the crystal is mounted on a specially designed four-point bender that provides independent and opposite couples at two points near each end of the crystal. The crystal is held by cylindrical rollers at only four points, minimizing the strain and consequently the distortions of the The two couples are computer controlled by step-by-step motors. A motorized twist correction has also been implemented combining a differential screw and a piezoelectric transducer. This actuator restores the cylindrical shape of the reflecting surface. This correction has proved to be essential for obtaining a well defined focus. The heat load on the crystal is partially drained by immersing the bottom third of the crystal in a eutectic Ga–In bath which produces a temperature gradient only in the direction perpendicular to the plane of energy dispersion. The eutectic bath is kept at a constant temperature of 292 K. The beam size on the crystal (Two symmetrically cut Si(111) crystals were tested: a profiled crystal (Eo = 7 keV) and one identical but with a rectangular contour.
Three different sets of data were collected to investigate the following:
(a) The focus size at different energies. The focus size was evaluated using a 7.5 µm-wide slit made of 20 µm-thick molybdenum foils moved by 5 µm steps. We tuned the applied moments to optimize the spot size and the uniformity of the distribution of energies within the spot. Fig. 4 compares the spots obtained with the profiled and the rectangular crystals at three different energies, after deconvoluting them by a 7.5 µm-wide slit function. It is worth noting that the focus distance was kept constant at 820 mm for all the studied energies except for 5.465 keV, for which we had to work under vacuum conditions. At this energy the geometrical constrains imposed by the equipment forced us to work at a distance of 1400 mm. Because of these constraints, the bandpass is reduced for this energy.
(b) The field depth of the focus. The focus depth (Fig. 5, Table 2) was measured at the energy of the Fe K-edge (7.1 keV) giving a value between 1 and 2 mm. The two-step slope appearing both in the X-ray tracing and in the experiment can be attributed to an aberration of the type sketched in Fig. 5(b), and can be caused by an error in the slope of the crystal at its middle point.
|
(c) The dependence of the spot size with the illuminated length of the crystal. We optimized the spot size for three different lengths of the footprint of the beam impinging on the crystal (which is the projection of Sc onto the crystal) but we did not observe any significant change in the spot size.
These results, together with the results of the X-ray tracing in the same geometrical conditions, are summarized in Table 2, which also includes the more relevant geometrical parameters.
Fig. 6(a) shows the experimental spot sizes for both crystals compared with the computed figure of merit. Both parameters show a similar energy dependence, proving the validity of our choice of the figure of merit. In Fig. 6(b) the calculated and measured spot sizes for both crystals are shown. First of all we observe the essential improvement yielded by the profiled crystal with respect to the rectangular crystal. This is especially true at high energies. Secondly it also illustrates the proper optimization energy selection that makes it possible to keep a small spot size for the whole energy range.
Despite the reduction of the focus size obtained with the profiled crystal, we observe that the measured spot sizes are systematically larger than the calculated ones (by a factor of ∼1.8). Different origins have to be envisaged:
(i) An underestimation of the source size. The source characteristics provided by the X-ray tracing are generated by idealized mirrors located upstream in the beamline. In comparison, real slope errors increase the size of the secondary source. The slope errors originate very likely from the deformation of the surface of the first mirror owing to the large thermal load deposited by the undulator X-ray source. This point is strengthened by the fact that the correlation between the spot size and the figure of merit (that does not depend on the source size) is 0.87, whereas the correlation between the measured size and the calculated size is only 0.7.
(ii) A residual twist in the crystal.
(iii) A variable thickness t of the crystal. If the thickness of the polished crystal is not uniform, the local curvature generated by the bender can fluctuate following the t−3 dependence of the It is easy to understand that a variation Δt = (t/3)(Δw/w) results in the same ΔI variation as that produced by a variation Δw of the width of the contour. In the present case, the maximum value of Δw/w is equivalent to Δt ≃ 20 µm. The measured thickness error was 10 µm at some points.
(iv) The procedure followed to measure the spot size. With 7.5 µm-wide slits moved in 5 µm steps we are just at the limit for accurate evaluation of image sizes smaller than 30–35 µm (reducing the slit gap leads to too large a sensitivity to artefact effects introduced by coherence). This is clearly appreciated when one takes into account that, with Sc values of 10 and 14 mm, the experimental size is 2.6 times the theoretical one (larger than 1.8, the average factor). Furthermore, in Fig. 6(b) it is clearly seen that the ratio between the measured and calculated spot sizes becomes smaller as the spot size increases.
(v) Finally, moving in the other direction, the sharp energy spectrum from the undulator gives a profile of the spot that also looks sharper and this effect is not taken into account in the calculations.
All these sources of errors can also explain that, for the optimizing energy (7 keV), no clear improvement of the focusing has been observed.
It is appropriate to compare the performances of our optics with that of previous works. In Table 3 we have summarized the results obtained by other authors. Full sizes mean that 90% of intensity is included. Specifically it is interesting to look at the results obtained by Allen et al. (1993), where using a rectangular profiled crystal they are able to focus a 5 mm-wide source (q/p = 0.03) down to 175 µm FWHM spots, which is 25 µm bigger than the theoretical values. The full size is of the order of 600 µm. In our situation, with rectangular crystals and the optimization of relevant parameters we are able to focus a 235 µm FWHM source (q/p = 0.028) to 10–17 µm FWHM spots, with ideal values of 6.5–11 µm. The full sizes are of the order of 20–40 µm. This achievement of one order of magnitude of reduction in size and a spectacular reduction of the tails of the image, which are the most significant of the aberration figures, is due to the quality of the third-generation source, which, consequently, is very demanding in terms of aberration rejection: the optical aberration to be beaten scales with the ideal image size.
4. Conclusions
We have designed a bent polychromator for energy-dispersive optics that takes advantage of the third-generation synchrotron characteristics. We have investigated the effect of the relevant parameters implemented in an extensive simulation, and have tested experimentally the behaviour of the produced contoured crystal mounted in a four-point bender. Considerable improvements with respect to the rectangular crystal were achieved. Spot sizes as small as 20–40 µm at energies ranging from 5 to 12 keV were measured.
Footnotes
‡Present address: Institut de Ciència dels Materials, Universitat de València, Dr. Moliner 50, Ed. Investigació, E-46100 Burjassot, València, Spain. E-mail: julio.pellicer@uv.es.
§Present address: Département de Physique des Matériaux, Batiment 203, Université Lyon I, 43 Bd du 11 novembre 1918, F 69622 Villeurbanne, France.
Acknowledgements
We are grateful to the staff of the ESRF where this work was performed, especially to G. Marot for the bender mechanical design, M. Sanchez Del Rio for his help in the X-ray tracing and M. Hagelstein for his useful comments. One of us (JP) is supported by the Ministerio de Educación y Ciencia of Spain (grant No. PN95 24351836).
References
Allen, P. G., Conradson, S. D. & Penner-Hahn, J. E. (1993). J. Appl. Cryst. 26, 172–179. CrossRef Web of Science IUCr Journals
Bonfim, M., Mackay, K., Pizzini, S., San Miguel, A., Tolentino, H., Giles, C., Neisius, T., Hagelstein, M., Baudelet, F., Malgrange, C. & Fontaine, A. (1997). Highlights of the ESRF 1996/1997, pp. 23–24. ESRF, Grenoble, France.
Flank, A. M., Fontaine, A., Jucha, A., Lemonnier, M., Raoux, D. & Williams, C. (1983). Nucl. Instrum. Methods, 208, 651–654. CrossRef CAS Web of Science
Hagelstein, M., Ferrero, C., Hatje, U., Ressler, T. & Metz, W. (1995). J. Synchrotron Rad. 2, 174–180. CrossRef CAS Web of Science IUCr Journals
Hagelstein, M., San Miguel, A., Fontaine, A. & Goulon, J. (1997). J. Phys. IV, 7(C2), 303–308.
Howells, M. R. & Lunt, D. (1993). Opt. Eng. 32, 1981–1989. CrossRef Web of Science
Ice, G. E. & Sparks, C. (1984). Nucl. Instrum. Methods, 222, 121–127. CrossRef CAS Web of Science
Itié, J. P., Baudelet, F., Dartyge, E., Fontaine, A., Tolentino, H. & San Miguel, A. (1992). High Press. Res. 8, 697–702.
Itié, J. P., Polian, A., Gauthier, M. & San Miguel, A. (1997). Highlights of the ESRF 1996/1997, pp. 57–58. ESRF, Grenoble, France.
Koch, A., Hagelstein, M., San Miguel, A., Fontaine, A. & Ressler, T. (1996). IS&T/SPIE Proceedings of Symposium on Electronic Imaging: Science and Technology, San Jose, CA, USA, Vol. 2416, p. 85.
Landau, L. & Lifschitz, E. (1986). Theory of Elasticity. Oxford: Pergamon.
Lee, P. L., Beno, M. A., Jennings, G., Ramanathan, M., Knapp, G. S., Huang, K., Bai, J. & Montano, P. A. (1994). Rev. Sci. Instrum. 65, 1–6. CrossRef CAS Web of Science
Matsushita, T. & Phyzackerly, R. P. (1981). Jpn. J. Appl. Phys. 20, 2223. CrossRef Web of Science
Mosbah, M., Metrich, N., Dura, J. P., San Miguel, A., Delaney, J. S. & Mccammon, C. (1998). In the press.
Pizzini, S., Bonfim, M., Baudelet, F., Tolentino, H., San Miguel, A., Mackay, K., Malgrange, C., Hagelstein, M. & Fontaine, A. (1998). J. Synchrotron Rad. 5, 1298–1303. Web of Science CrossRef IUCr Journals
Sanchez del Rio, M., Ferrero, C. & Mocella, V. (1998). Proc. SPIE, 3151, 312–323. (http://www.esrf.fr/computing/scientific/xop/.)
San Miguel, A., Hagelstein, M., Borrel, J., Marot, G. & Renier, M. (1998). J. Synchrotron Rad. Submitted. CrossRef IUCr Journals
San Miguel, A., Polian, A., Gauthier, H. & Itié, J. P. (1998). In preparation.
Tolentino, H., Baudelet, F., Dartyge, E., Fontaine, A., Lena, A. & Tourillon, G. (1990). Nucl. Instrum. Methods, A289, 307–316. CrossRef CAS Web of Science
Tolentino, H., Dartyge, E., Fontaine, A. & Tourillon, G. (1988). J. Appl. Cryst. 21, 15–21. CrossRef CAS Web of Science IUCr Journals
Underwood, J. M. (1977). Space Sci. Instrum. 3, 259–270.
© International Union of Crystallography. Prior permission is not required to reproduce short quotations, tables and figures from this article, provided the original authors and source are cited. For more information, click here.