Research Article  Open Access
F. N. Dultsev, E. A. Kolosovsky, "Application of Ellipsometry to Control the Plasmachemical Synthesis of Thin TiONx Layers", Advances in Condensed Matter Physics, vol. 2015, Article ID 709308, 8 pages, 2015. https://doi.org/10.1155/2015/709308
Application of Ellipsometry to Control the Plasmachemical Synthesis of Thin TiONx Layers
Abstract
Ellipsometry is often used to determine the characteristics of films. Ellipsometric studies may turn out to be ineffective because several solutions correspond to the same polarization angles. It is demonstrated that the ambiguity is not due to the physical limitations of the method but it has a purely mathematical character. So, additional information about the film is necessary to determine the absolute values of refractive index, attenuation, and thickness.
1. Introduction
Preparation of thin layers with the given stoichiometry and minimal size of granules is important for the development of nanostructures. For example, ultrathin layers of titanium nitride (TiN) are widely used in microelectronics as adhesion layers or diffusion barriers [1, 2], as metalgate materials for metal oxide semiconductor devices. They possess good mechanical properties and thermal stability [3, 4]. In the present work, we used ellipsometry to study the mechanism of nitridation of the titanium matrix deposited on silicon surface [5]. Thus we obtained the layers of titanium nitride or oxynitride at the substrate temperature of 4–6°C. The advantage of this method in comparison with that described in [6–9] is, first of all, low temperature at which nitridation is carried out; second, this method allows us to obtain finestructured layers of required thickness.
On the other hand, reliable control devices are needed to prepare such layers. It is necessary to apply highprecision measuring means intended for efficient diagnostics and tests of promising materials and structures under development. Ellipsometers belong to the measuring means of this kind [10–13]. A broad range of the applications of ellipsometry (from nanoelectronics to biotechnologies), as well as the availability of uptodate ellipsometric instrumentation, make optical ellipsometry a very attractive and available method. However, it is the high sensitivity of this technique to small changes in the parameters of the medium that can cause in some cases incorrect interpretation of measurement results. This is especially so if the formed layers are not transparent. In our work we use ellipsometry to observe the formation of thin films having different compositions and thickness from several monolayers to nm. For the case when a film grows without cluster formation, in order to understand and reject incorrect interpretations of ellipsometric measurements, we will consider the reasons of misinterpretation. We will demonstrate below that misinterpretation of thin film parameters and error recovery are not connected with the physical nature of measuring technique but they have purely mathematical character.
It is known that the solution of inverse ellipsometric problem is ambiguous in the case of discrete light angles. For the ellipsometric examination of the growing TiONx films, we asked ourselves the question whether the ambiguity of the determination of film parameters is conserved if the number of measurements is increased many times embracing the whole possible range of light angles , or an increase in the number of measurements will cause the statistical adjustment of the same parameters. The second question is as follows: is this situation conserved for spectral ellipsometry? And the last question is as follows: is it possible to escape from ambiguity or decrease its extent in the ellipsometric monitoring of growing TiONx?
2. Experimental Section
The metal Ti layer was deposited onto the 4inch epiready Si (001) wafers in high vacuum (10^{−5}–10^{−7} Pa), using the ebeam sputtering technique. The Ti layers were deposited with the molecular beam epitaxy unit. The rate of metal layer deposition was 10 nm/min. Layer thickness was 5 to 50 nm.
After deposition, we used plasmachemical nitridation of Ti nanolayers to obtain TiN nanolayers. Plasmachemical treatment was performed in nitrogen gas plasma at a pressure of 0.8 Torr, with RF power 150 Watt and radio frequency 13.56 MHz (“MATRIX” unit) at a temperature of 4–6°C.
The metaltodielectric transition, that is, the formation of the dielectric film, was studied with the ellipsometer. Ellipsometric measurements were carried out according to multiangle technique with a LEF4M ellipsometer at the wavelength of 632.8 nm. Measurements were carried out at the light angle within the range 45–80 angular degrees. The ellipsometric setup was assembled according to the PCSA scheme (polarizercompensatorfilmanalyzer); beam blanking method (null ellipsometry) was used. Ellipsometer allows carrying out measurements over 4 angle zones. In addition to the ellipsometric measurements, UV absorption spectra of the films were recorded with the Shimadzu UV2401PC spectrometer within the wavelength region 190–500 nm. Spectroscan ellipsometer manufactured at the ISP SB RAS was used for spectral ellipsometry.
The highresolution transmission electron microscopy (HRTEM) images of TiNSiO_{2}Si interfaces were taken with JEM 4000 EX. Measurements were carried out at 300 kV. The surface roughness of films was determined with AFM Solver P47H atomic force microscope (AFM) in the resonant mode operating at 350 kHz and 38 nm amplitude. The size of the tip was 10 nm.
Equations used in ellipsometric calculations are simple and readily testable in the case of the direct problem (by means of the forward substitution). To solve these problems, one can use any mathematical program package: MATLAB, Mathcad, Excel, and so forth. To solve the inverse problem, we used our own original algorithm within the framework of the same packages; the system of equations was solved over the complex numbers field, similarly to the method described elsewhere [14].
3. Results and Discussion
The kinetics of the transformation of metal to dielectric during plasma treatment is shown in Figure 1. These results were obtained from ellipsometric data. Measurements were carried out with the help of ellipsometer at different light angles, as well as at different wavelengths (Spectroscan). Absorption was measured by means of UV spectroscopy. The figure shows how the film thickness changes actually versus nitridation time; during nitridation, we observe total increase in film thickness.
The initial film thickness was 90 Å. Thickness was calculated according to the procedure proposed by us.
We determined the thickness of Ti film before nitridation with the help of AFM. Using AFM, we also observed the change of film relief as a result of nitridation. One can see that nitridation is accompanied by surface smoothing. We demonstrated previously [15] with the help of highresolution transmission electron microscopy that nitridation proceeds layer by layer, so we used the multilayer model. The ellipsometric data show that layer thickness increases during nitridation almost by a factor of 2 (Figure 1). It was shown previously that the resulting layers are uniform, and a model for ellipsometric calculations was proposed [15]. Now we will consider the solution of the inverse problem and demonstrate why a set of solutions arises.
3.1. General Solution of the Ellipsometric Equations
The fundamental equation of ellipsometry connecting ellipsometric parameters delta () and psi () with the complex absolute reflectance coefficients and for  and polarized light (see Figure 3) is written as [16–19] In the case of a multilayered structure, the absolute reflectance coefficients , coincide with coefficient calculated in the recurring manner using (2) and (3) for the corresponding polarizations [14–18]: where are the coefficients of reflection from the system of the first layers; are the coefficients of reflection from the interfaces between the layers; are phase incursions in layers.
Parameters and are determined on the basis of Fresnel’s formulas: Here for polarized waver, and for polarized wave, while angles (the angles at which the wave spreads over the layers) are determined from Snellius law: , is the wavenumber of light in vacuum. For absorbing media, a refractive index of each medium is also represented as a complex factor: . The direct problem of ellipsometry involves determination of parameters delta () and psi (); it is easily solved for the geometrically flat interface. In a multilayered structure with a set of layers this determination, according to (1), holds true unambiguously. It should be noted that (1)–(3) imply isotropy of layers.
When solving the inverse ellipsometric problem, the optical parameters of sample surface are calculated on the basis of measured and values relying on the chosen model. The coincidence of the calculated and with the experimental values is the evidence in favour of the correctness of the chosen optical model. As a rule, it is necessary to solve the inverse ellipsometric problem, that is, to determine the optical characteristics of the reflecting system on the basis of the measured set of and values under different conditions: different light angles at which the light falls on the surface from different media and different light frequencies (socalled spectral ellipsometry) [14, 17, 19, 20].
3.2. Multiangle Measurements at One Wavelength
At first, we will consider a onelayer structure shown in Figure 2: a model metal film with refractive index and thickness Å on silicon substrate (refractive index for silicon at the light wavelength of Å is ). Assigning the angles at which the light falls from the air () within the range from 45 to 80 with a step of 1 degree, we will obtain the solution of the direct ellipsometric problem in the form of a plot in Figure 3, where the numbers above points correspond to the light angle . The side dash lines indicate the area of plot shift with a 10% change of film thickness . One can see that the change of film thickness by only ±1 nm leads to a visible substantial shift at the plot, which confirms the high sensitivity of the ellipsometric technique to film parameters during film formation.
Now we will pass to the inverse problem of restoration of film parameters. Experimental () and () values measured for one specific light angle allow us to calculate one of the film parameters—thickness , or its complex refractive index following a general rule: one equation—one unknown quantity. Usually not two but several measurements for different light angles are used, so that matching of the results with the help of the least squares method could help us to eliminate the errors in angle measurement.
It is known that the inverse problem when solving (1) in a onelayer structure can give several solutions [17, 19] (this fact is not connected with the harmonic exponent in (2) because ). It was demonstrated that the total number of possible solutions can be determined beforehand applying the theory of the functions of complex variables [14]. To confirm this statement, the model data are listed in Table 1 depicting several solutions: different values of refractive index for the films of different thickness values. All of them satisfy (1) and go well into the same central curve of Figure 3 within the whole range of angles and not only at some values of this light angle.

Now it is necessary to answer the question whether the ambiguity of parameter determination in two and more layered films remains if we increase the number of measurements embracing the whole possible range of light angles , or an increase in the number of measurements would cause statistical refinement of the same parameters.
We specially carried out a thorough investigation of this question concerning the solution of the inverse problem. For two and threelayered structures, at first we chose the same total film thickness equal to Å and the same attenuation as those for the onelayer structure (shown in bold characters in Tables 1, 2, and 3). Thus, the starting dependencies of parameters () and () were reduced to one plot coinciding with curve 2 in Figure 3. Then thickness was varied within rather broad range, and a special algorithm was used to reveal the new meanings of refractive index to give coincidence with the starting one by solving the direct problem within the whole range of angles . The mathematical precision of this coincidence exceeded the possible errors of experimental angle measurements by several orders of magnitude. As a result, all the one, two, and threelayered structures marked in Tables 1–3, in spite of essential differences in their parameters, belong to the same ellipsometric curve in Figure 3. We see that an increase in the number of measurements made with a singlewavelength ellipsometric data (in multiangle ellipsometry) conserves the ambiguity of film parameter determination. The experimental values for ellipsometric angles delta and psi are shown in Table 4. These values were used to perform verification.



3.3. Spectral Ellipsometry
As far as spectral ellipsometry is concerned, with the measurements performed at different light wavelengths, a similar situation is still actual. Similarly to the previous case, Figure 4 depicts the solution of the direct problem of spectral ellipsometry for three layers 100, 200, and 300 nm thick, deposited on silicon (line). Then the thickness of two films was varied by ~50 nm, and the new film parameters were determined by solving the inverse problem. One can see in this figure that the new solution (points) graphically coincides with the primary solution. Figures on the plot mark the positions of measurements depending on the light wavelength , which was varied from 250 nm to 850 nm with a step of 12 nm, which corresponds to several hundred measured data values. The light angle is equal to 70 angular degrees. Film parameters for both solutions are indicated in Table 5. One can see that the refractive index for layer changed insignificantly and refractive index and attenuation in layer decreased much more substantially, while the refractive index of the very first layer increased from to 1.8983, and losses for light extinction disappeared completely for this layer.

As the objects for investigation, we prepared three samples of TiONx films with thickness nm on silicon according to the technology described in [15]. The data for titanium layer deposited on silicon substrate for different nitridation times are shown in Figure 5. In view of the fact that the TiOx films were deposited onto silicon, one should take into account the changes of its optical properties versus the wavelength. The data shown in Figure 6 demonstrate that the refractive index changes by a factor of 7, while absorption changes by a factor of several hundreds. This change is not monotonous but, unlike for quartz [21], it has a complicated character due to differences in the mechanism of light absorption in Si (electronhole and latticerelated) [22]. Because of this, the restoration of the parameters of TiONx films was carried out as follows. On each plot of measured data in Figure 5, we could choose any three successive points at a small distance from each other. Because measurements are performed with a small step over wavelength , we may accept that the optical parameters of the films and the substrate are constant within the limit of these three measurements and correspond to the second (average) value of . Taking the thickness of TiONx film from additional measurements (e.g., by means of AFM) and choosing the refractive index of the silicon substrate according to the data shown in Figure 6, we may solve the direct and inverse problems of ellipsometry and determine film parameters.
Repeating this procedure many times with other groups of three points , we obtain all the optical parameters of TiOx films including dispersion, the spectral dependence of refractive index, and absorption of light on wavelength .
For comparison, at the bottom of Figure 5 we show curve 1 calculated without taking dispersion into account. It depicts the dependence of ellipsometric parameters , for curve 3a related to the film with thickness nm and refractive index corresponding to the red light with wavelength μm (see Figure 7). Comparing the two curves 3a and 3b in Figure 5, one can see that the account of the effect of dispersion in precise ellipsometric measurements does not bear only a character of correction but it holds the defining importance. Actually, the account of this fact does not hinder processing and interpretation of changes in film structures but even allows one, relying on the optical dispersion of the silicon substrate, to extract the unknown spectral characteristics of TiONx and similar films themselves. That is, it allows us to know the behaviour of both the refractive index and extinction coefficient versus the wavelength of light in any thin films on silicon. If the thickness and the optical parameters of the films are determined correctly, curves 3a and 3b in Figure 5 should coincide. We carried out this work, and the results for three samples when two curves 3a and 3b (in Figure 5) completely coincide are presented in Figure 7. Of course, due to the mentioned mathematical ambiguity of the solution of the inverse problem of ellipsometry, thickness should be determined also using an independent method.
(a)
(b)
(c)
The real thickness of TiONx films in all the three samples was to 22 nm. In order to demonstrate how film thickness may affect its optical characteristics, we processed seven successive thickness values. All the results with different thickness values shown in Figure 7 correspond to the experimental curves shown in Figure 5. It may be concluded from the behaviour of the absorption curve (Figure 7(c)) that the films with thickness nm and less cannot exist. We conclude this because starting from this thickness a “passive” film should possess amplifying properties because for nm the absorption coefficient of the film material changes its sign (see the lower dash line in Figure 7(c)).
One can see in the plot that the monotonous decrease in the refractive index of the TiONx film, so typical for quartz glass, is present only in the second (longer wavelength) half of the optical range, for radiation with the wavelengths longer than 500 nm. In the first (shorter wavelength) part of the optical range, the refractive index of the film material changes substantially, from 1.5 to 3.5. One can notice that light extinction attenuation in the film is very large. The absorption coefficient changes in the nonmonotonic manner and is large, within the range from 1 to 3. Because of this, the described procedure based on the data of spectral ellipsometry is applicable mostly to the thin films, that is, for TiONx films with thickness 2 to 50 nm (20–500 Å), and the films composed of any other material, if they are transparent enough (their thickness should be small enough to allow us to see the substrate). The procedure proposed by us for the determination of the behaviour of the optical characteristics of films uses the known dispersion in silicon as the calibration standard, the radiation must “sense” the material of the silicon substrate, and strongly absorbing films should not be too thick. In other words, in order to reduce the extent of ambiguity, we used the effect of dispersion of the refractive index, that is, its dependence on wavelength in Si. Usually even a small change of the refractive index versus wavelength within the visible light range is considered as a parasitic (undesirable) phenomenon, because it causes errors in the measurement of film parameters. However, it is this effect (dispersion of light in silicon) that turned out to be essential when thoroughly measured within a broadened wavelength range, and we used it as the basic information for calibration. In the case of thin films on silicon, either TiONx or other compositions, when the incident light wave partially penetrates and thus “senses” the Si substrate, this is possible.
So, we see that an increase in the number of ellipsometric measurements and broadening of the range of light angles for one, two, and more layered structures conserve ambiguity of the determination of film parameters. This ambiguity is not due to the physical limitations of the method but it bears the mathematical character of socalled illconditioned problem [23–25]. Because of this, in spite of the high sensitivity of the method, determination of refractive index and thickness of absorbing thin film requires additional data. For example, additional measurements of absorption or thickness of the film are to be carried out using other methods.
Because of this, the most promising approach is to combine ellipsometry with any other methods of surface studies, for example, Auger spectroscopy, UV and Xray spectroscopy, electron diffraction and ion scattering, and so forth. The advantages of ellipsometry are simplicity and rapidity of measurements (automatic ellipsometers are available), the possibility to carry out measurements in course of the process (in situ), in vacuum, at high temperatures, and in aggressive media; in addition, the surface is not contaminated or destroyed during measurements.
4. Conclusion
In spite of the high differential sensitivity of ellipsometric method, for titanium nitride films obtained by means of nitridation (metaltodielectric transition), determination of the absolute values of refractive index, attenuation, and film thickness requires independent additional information.
For onelayer, twolayer structures and the structures with larger number of layers, an increase in the number of ellipsometric measurements and broadening of the light angle range or wavelength range (for spectral ellipsometry) conserve the ambiguity of film parameter determination. Ambiguity is not due to the physical limitations of the method but has the mathematical character.
Solution of the inverse problem for two and more layered structures showed that the measurement of the same films within different angle ranges, as well as an increase in the number of measurements, does not allow one to obtain new information on film composition, and no adjustment of film parameters occurs in principle. Because of the mathematical specificity of equations, coincidence of solutions for different models is possible not only in discrete points and but also along the entire curve.
We considered one, two, three, and more layered structures that belong completely to the same ellipsometric curve, in spite of the substantial differences between their parameters.
This situation is conserved also for spectral ellipsometry. We demonstrated quite different structures that gave coincidence (within the measurement error) of the , curves within their entire continuous range under continuous change of the wavelength of light.
In order to eliminate ambiguity or decrease its extent in the ellipsometric control of growing TiONx or other films, we used dispersion of refractive index, that is, its dependence on the light wavelength in Si. Thoroughly measured within the whole wavelength range, this effect turned out to be essential. We used it as the major calibration information. This is possible in the case of absorbing thin films on silicon, when the incident light wave partially penetrates and “senses” the Si substrate.
Thus, confirming the necessity of additional information for the determination of film parameters using optical methods that are convenient, highly precise, and nondestructive, for TiONxSi films as example we demonstrated the whole set of the reasons of ambiguity inherent to ellipsometry, either discrete or continuous, either angular or spectroscopic.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This work was supported by the Ministry of Education and Science of the Russian Federation.
References
 T. Tanaka, T. Nakajima, and K. Yamashita, “Density functional study on the adsorption and surface reactions on SiO_{2} in TiNCVD using TiCl_{4} and NH_{3},” Thin Solid Films, vol. 409, no. 1, pp. 51–57, 2002. View at: Publisher Site  Google Scholar
 J. M. Wang, W. G. Liu, and T. Mei, “The effect of thermal treatment on the electrical properties of titanium nitride thin films by filtered arc plasma method,” Ceramics International, vol. 30, no. 7, pp. 1921–1924, 2004. View at: Publisher Site  Google Scholar
 N. Savvides and B. Window, “Electrical transport, optical properties, and structure of TiN films synthesized by lowenergy ion assisted deposition,” Journal of Applied Physics, vol. 64, no. 1, pp. 225–234, 1988. View at: Publisher Site  Google Scholar
 H. Z. Durusoy, O. Duyar, A. Aydinili, and F. Ay, “Influence of substrate temperature and bias voltage on the optical transmittance of TiN films,” Vacuum, vol. 70, no. 1, pp. 21–28, 2003. View at: Publisher Site  Google Scholar
 O. V. Naumova, Y. V. Nastaushev, S. N. Svitasheva et al., “Molecularbeam epitaxygrown Si whisker structures: morphological, optical and electrical properties,” Nanotechnology, vol. 19, no. 22, Article ID 225708, 2008. View at: Publisher Site  Google Scholar
 Y. Igasaki, H. Mitsuhashi, K. Azuma, and T. Moto, “Structure and electrical properties of titanium nitride films,” Japanese Journal of Applied Physics, vol. 17, no. 1, pp. 85–96, 1978. View at: Publisher Site  Google Scholar
 G. Zhao, T. Zhang, J. Wang, and G. Han, “Electrical and optical properties of titanium nitride coatings prepared by atmospheric pressure chemical vapor deposition,” Journal of NonCrystalline Solids, vol. 354, no. 1213, pp. 1272–1275, 2008. View at: Publisher Site  Google Scholar
 S. R. Kurtz and R. G. Gordon, “Chemical vapor deposition of titanium nitride at low temperatures,” Thin Solid Films, vol. 140, no. 2, pp. 277–290, 1986. View at: Publisher Site  Google Scholar
 L. Imhoff, A. Bouteville, and J. C. Remy, “Kinetics of the formation of titanium nitride layers by rapid thermal low pressure chemical vapor deposition from TiCl_{4}NH_{3}H_{2},” Journal of the Electrochemical Society, vol. 145, no. 5, pp. 1672–1677, 1998. View at: Publisher Site  Google Scholar
 L. Ryves, M. M. M. Bilek, T. W. H. Oates et al., “Synthesis and insitu ellipsometric monitoring of Ti/C nanostructured multilayers using a highcurrent, dual source pulsed cathodic arc,” Thin Solid Films, vol. 482, no. 12, pp. 133–137, 2005. View at: Publisher Site  Google Scholar
 U. I. Schmidt, B. Schröder, and H. Oechsner, “In situ ellipsometric study of the influence of powder formation on the initial growth of glow discharge aSi:H,” Thin Solid Films, vol. 233, no. 12, pp. 297–300, 1993. View at: Publisher Site  Google Scholar
 A. von Keudell and W. Jacob, “Growth and erosion of hydrocarbon films investigated by in situ ellipsometry,” Journal of Applied Physics, vol. 79, no. 2, pp. 1092–1098, 1996. View at: Publisher Site  Google Scholar
 A. Granier, G. Borvon, A. Bousquet, A. Goullet, C. Leteinturier, and A. van der Lee, “Mechanisms involved in the conversion of ppHMDSO films into SiO_{2}like by oxygen plasma treatment,” Plasma Processes and Polymers, vol. 3, no. 45, pp. 365–373, 2006. View at: Publisher Site  Google Scholar
 V. G. Polovinkin and S. N. Svitasheva, “Determination of the number of inverse eliipsometric problem solutions in the given domain of parameters,” Optoelectronics Instrumentation Data Processing, vol. 4, p. 79, 1999. View at: Google Scholar
 F. N. Dultsev, S. N. Svitasheva, Y. V. Nastaushev, and A. L. Aseev, “Ellipsometric investigation of the mechanism of the formation of titanium oxynitride nanolayers,” Thin Solid Films, vol. 519, no. 19, pp. 6344–6348, 2011. View at: Publisher Site  Google Scholar
 G. Hass and R. E. Thun, Eds., Physics of Thin Films, vol. 2, Academic Press, New York, NY, USA, 1964.
 R. M. A. Azzam and N. M. Bashara, Ellipsometry and Polarized Light, NorthHolland, Amsterdam, The Netherlands, 1977.
 R. Rosa, “The inverse problem of ellipsometry: a bootstrap approach,” Inverse Problems, vol. 4, no. 3, pp. 887–900, 1988. View at: Publisher Site  Google Scholar
 A. V. Rzhanov, Ed., Ellipsometry as a Method of Surface Investigation, Nauka, Novosibirsk, Russia, 1983, (Russian).
 F. Scandonne and L. Ballerini, “Teorie de la transmission et de la reflexion dans les systems de couches minces multiples,” Nuovo Cemento, vol. 3, pp. 81–91, 1946. View at: Google Scholar
 L. A. Ageev, K. S. Beloshenko, E. D. Makovetsky, and V. K. Miloslavsky, “Implantation of periodic structures formed by silver particles into quartz glass,” Optics and Spectroscopy, vol. 107, no. 5, pp. 796–802, 2009. View at: Publisher Site  Google Scholar
 D. Aspnes and A. Studna, “Dielectric functions and optical parameters of Si, Ge, GaP, GaAs, GaSb, InP, InAs, and InSb from 1.5 to 6.0 eV,” Physical Review B, vol. 27, no. 2, pp. 985–1009, 1983. View at: Publisher Site  Google Scholar
 A. Tikhonov and V. Arsenin, Solutions of IllPosed Problems, Winston and Sons, Washington, DC, USA, 1977.
 A. Akbalik, S. Soulan, J.H. Tortai et al., “An inverse ellipsometric problem for thin film characterization: comparison of different optimization methods,” in Metrology, Inspection, and Process Control for Microlithography 23, vol. 7272 of Proceedings of SPIE, San Jose, Calif, USA, March 2009. View at: Publisher Site  Google Scholar
 H. Arwin, “Is ellipsometry suitable for sensor applications?” Sensors and Actuators, A: Physical, vol. 92, no. 1, pp. 43–51, 2001. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 F. N. Dultsev and E. A. Kolosovsky. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.