The Effect of Annealing Temperature on Statistical Properties of $WO_3$ Surface
G. R. Jafari, A. A. Saberi, R. Azimirad, A. Z. Moshfegh, and S., Rouhani

TL;DR
This study investigates how annealing temperature influences the surface statistical properties of $WO_3$ thin films, revealing an optimal temperature around 400°C where surface features and microstructure undergo significant changes.
Contribution
It introduces a combined use of level crossing and structure function methods to analyze surface properties and identifies a phase transition at 400°C related to microstructural changes.
Findings
Maximum effective surface area at 400°C
Surface transitions from amorphous to crystalline structure
Phase transition indicated by changes in roughness parameters
Abstract
We have studied the effect of annealing temperature on the statistical properties of surface using atomic force microscopy techniques (AFM). We have applied both level crossing and structure function methods. Level crossing analysis indicates an optimum annealing temperature of around 400 at which the effective area of the thin film is maximum, whereas composition of the surface remains stoichiometric. The complexity of the height fluctuation of surfaces was characterized by roughness, roughness exponent and lateral size of surface features. We have found that there is a phase transition at around 400 from one set to two sets of roughness parameters. This happens due to microstructural changes from amorphous to crystalline structure in the samples that has been already found experimentally.
Click any figure to enlarge with its caption.
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6
Figure 1
Figure 2
Figure 2
Figure 4
Figure 5
Figure 6| T | - | |||
|---|---|---|---|---|
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsSurface Roughness and Optical Measurements
The Effect of Annealing Temperature on Statistical Properties of
Surface
G. R. Jafari a,b, A. A. Saberi c, R. Azimirad c, A. Z. Moshfegh c, and S. Rouhani c
a* Department of Physics, Shahid Beheshti University, Evin, Tehran 19839, Iran
b Department of Nano-Science, IPM, P. O. Box 19395-5531, Tehran, Iran
c Department of Physics, Sharif University of Technology, P.O. Box 11365-9161, Tehran, Iran *
Abstract
We have studied the effect of annealing temperature on the statistical properties of surface using atomic force microscopy techniques (AFM). We have applied both level crossing and structure function methods. Level crossing analysis indicates an optimum annealing temperature of around 400 at which the effective area of the thin film is maximum, whereas composition of the surface remains stoichiometric. The complexity of the height fluctuation of surfaces was characterized by roughness, roughness exponent and lateral size of surface features. We have found that there is a phase transition at around 400 from one set to two sets of roughness parameters. This happens due to microstructural changes from amorphous to crystalline structure in the samples that has been already found experimentally.
I Introduction
Transition metal oxides represent a large family of materials possessing various interesting properties, such as superconductivity, colossal magneto-resistance and piezoelectricity. Among them, tungsten oxide is of intense interest and has been investigated extensively for its distinctive properties. With outstanding electrochromic Granqvist ; Bueno ; Azimirad ; Siokou ; Kuai ; Takeda , photochromic Avellaneda , gaschromic Lee , gas sensor Kim ; Gy rgy ; Kawasaki , photo-catalyst Gondal , and photoluminescence properties Feng , tungsten oxide has been used to construct ”smart-window”, anti-glare rear view mirrors of automobile, non-emissive displays, optical recording devices, solid-state gas sensors, humidity and temperature sensors, biosensors, photonic crystals, and so forth.
thin films can be prepared by various deposition techniques such as thermal evaporation Azimirad ; Lee , spray pyrolysis Hao , sputtering Takeda , pulsed laser ablation Gy rgy ; Kawasaki , sol-gel coating Bueno ; Kuai ; Garcia , and chemical vapor deposition Seman .
The gas sensitivity of heavily depends upon film parameters such as composition, morphology (e.g. grain size), nanostructure and microstructure (e.g. porosity, surface-to-volume ratio). Film parameters are related to the deposition technique used, the deposition conditions and the subsequent annealing process. Annealing, which is an essential process to obtain stable films with well-defined microstructure, causes stoichiometry and microstructural changes that have a high influence on the sensing characteristics of the films Stankova . Moreover, the surface structure and surface morphology of the metal oxides are also important for different applications. In fact, the electrochromic devices are made of amorphous oxides Granqvist , while crystalline phase plays a major role in catalysts and sensors Stankova . This is because, the minor change in their chemical composition and crystalline structure could modify different properties of the metal oxides.
In practice, one of the effective ways to modify the surface morphology is annealing process at various temperatures. So far, most of morphological analysis related to the surface were accessible through the experimental methods. Usually these analysis are rigorous and time consuming. Moreover, lack of the suitable analysis for AFM data to find the nano and microstructural properties of surfaces was feeling perfectly.
In this article, we introduce the methods: roughness analysis and level crossing as suitable candidates and show that we can get easily the structure and morphological properties of a surface in a fast manner, only using the AFM observation as an initial data.
The roughness of a surface has been studied as a simple growth model using analytical and numerical methods Barabasi ; Jafari10 ; Jafari11 ; Irajizad ; Halpin ; Krug ; Meakin ; Kardar ; Masoudi . These studies quite generally proposed that the height fluctuations have a self-similar character and their average correlations exhibit a dynamic scaling form. Also some authors recently use the average frequency of positive slope level crossing to provide further complete analysis on roughness of a surface Tabar . This stochastic approach has turned out to be a promising tool also for other systems with scale dependent complexity, such as in surface growth where one would like to measure the roughness Tabar1 . Some authors have applied this method to study the fluctuations of velocity fields in Burgers turbulence Movahed and the Kardar-Parisi-Zhang equation in (d+1)-dimensions Bahraminasab and analyzing the stock market Jafari2 .
In this work, we have used the scaling analysis to determine the roughness, roughness exponent and the lateral size of surface features. Moreover, level crossing analysis has been utilized to estimate the effective area of a surface.
This paper is organized as follows: In section II, we have discussed about the film preparation and experimental results obtained from AFM, XPS and UV-visible spectrophotometer for the annealed samples at the various temperatures. In section III, we have introduced the analytical methods briefly. Data description and data analysis based on the statistical parameters of surface as a function of annealing temperatures are given in section IV.Finally, section V concludes presented results.
II Experimentals
Thin films of were deposited on microscope slide glass using thermal evaporation method. The deposition system was evacuated to a base pressure of . Thickness of the deposited films was considered about 200 measured by the stylus and optical techniques. More details about the other deposition parameters of the films are recently reported elsewhere Moshfegh .
To study the effect of annealing temperature on surface structure and optical properties of the samples, they were annealed at 200, 300, 350, 400, 450, and 500 in air for a period of 60 . Optical transmission and reflection measurements of the deposited films were performed in a range of 300-1100 wavelength using a Jascow V530 ultraviolet (UV)-visible spectrophotometer with resolution of 1 .
X-ray photoelectron spectroscopy (XPS) using a Specs EA 10 Plus concentric hemispherical analyzer (CHA) with anode at energy of 1486.6 was employed to study the atomic composition and chemical state of the tungsten oxide thin films. The pressure in the ultra high vacuum surface analysis chamber was less than . All binding energy values were determined by calibration and fixing the line to 285.0 . The XPS data analysis and deconvolution were performed by SDP (version 4.0) software. The nanoscale Surface topography of the deposited films was investigated by Thermo Microscope Autoprobe CP-Research atomic force microscopy (AFM) in air with a silicon tip of 10 radius in contact method. The AFM images were recorded with resolution of about in a scale of .
II.1 XPS Characterization
The elemental and chemical characterizations of the films were performed by XPS. Figure shows the core level spectra recorded on the ”as deposited” sample, and the results of its fitting analysis. To reproduce the experimental data, one doublet function was used for the component. This contains at 35.6 and at 37.8 with a full-width at half-maximum (FWHM) of . The area ratio of these two peaks is 0.75 which is supported by the spin-orbit splitting theory of levels. Moreover, the structure was shifted by 5 toward higher energy relative to the metal state.
It is thus clear that the main peaks in our XPS spectrum attributed to the state on the surface Granqvist ; Bueno ; Crist . In stoichiometric , the six valence electrons of the tungsten atom are transferred into the oxygen p-like bands, which are thus completely filled. In this case, the tungsten valence electrons have no part of their wave function near the tungsten atom and the remaining electrons in the tungsten atom experience a stronger Coulomb interaction with the nucleus than in the case of tungsten atom in a metal, in which the screening of the nucleus has a component due to the valence electrons. Therefore, the binding energy of the level is larger in than in metallic tungsten. If an oxygen vacancy exists, the electronic density near its adjacent atom increases, the screening of its nucleus is higher and, thus, the level energy is expected to be at lower binding energy Granqvist .
By increasing the annealing temperature it was observed that the position of peak did not obviously change. But for thin film annealed at (Fig. 1b), the peak moved to a lower binding energy so that position was observed at 35.0 . This can be related to oxygen vacancy at this high annealing temperature and formation of .
II.2 Optical Characterization
The transmittance and reflectance spectra in the visible and infrared range recorded for the thin films before and after annealing at different temperatures (Fig. 2a). It is seen that, the transmittance of the ”as deposited” films in the visible range varies from about 80 up to nearly (without considering the substrate contribution). Correspondingly, maximum value of the reflectance for both the film and the substrate is about (the reflectance from the bare glass substrate was measured about ). The sharp reduction in the transmittance spectrum at the wavelength of is due to the fundamental absorption edge that was also reported previously Granqvist ; Bueno ; Azimirad .
The oscillations in the transmission and reflection spectra are caused by optical interference. The optical transmittance of films strongly depends on the oxygen content of the films. In fact, non-stoichiometric films with composition of show a blue tinge for Lee2 .
The ”as deposited” pure tungsten oxide films were highly transparent with no observable blue coloration, under our experimental conditions. As can be seen from Fig. 2a after annealing process at 200 to 400, the transmittance and reflectance of the films have not changed significantly. Only, the position of the oscillations altered due to thickness reduction and film condensation after the heat-treatment process Granqvist . At transmittance and reflectance of the annealed film is reduced about , therefore at this temperature, the film turn into non-stoichiometric composition, so that it could be seen from changing color of the film.
The optical gap () was evaluated from the absorption coefficient () using the standard relation: , in which depends on the kind of optical transition in semiconductors, and was determined near the absorption edge using the simple relation: ,where is thickness of the film. More useful explanation about the optical band gap calculation reported in Moshfegh . The relationship between the optical band gap energy and annealing temperature for thin films has been shown in Fig. 2b. As can be seen from it, the optical band gap for the ”as deposited” evaluated 3.4 . Amorphous structure of the ”as deposited” causes to is bigger than 2.7 . After annealing samples at 200 and , the optical band gap decreased slightly about 0.1 which can be related to condensation of the films. But the optical band gap of the annealed at reduced to 3.1 due to crystallization of the film. This reduction continues to 2.5 for the sample annealed at . Reason of the becomes smaller than 2.7 is oxygen vacancy at this temperature as was seen in Fig 1b. It is to note that for evaporated films one has found Granqvist .
II.3 AFM Analysis
To study the effect of the annealing process on the surface morphology of the films, we have shown AFM images of the surfaces annealed at the different temperatures : 200, 300, 350, 400, 450, 500 in Figure 1. As can be seen from Fig. 1, for the annealed film at 200, it seems that the surface morphology of the film is relatively the same with a smooth surface, amorphous structure and nanometric grain size, as also reported by other investigators for films Antonik ; Mohammad . We have also observed similar image for the ”as deposited” which is not shown here. For thin films, increasing annealing temperature to did not significantly effect on surface parameters because it is low temperature for crystallization of Granqvist . But at higher annealing temperatures 400, 450 and 500, surface grain size and roughness begin to increase. The more precise analysis of these surfaces are given in the next section.
III Statistical quantities
III.1 Roughness Analysis
It is also known that to derive the quantitative information of the surface morphology one may consider a sample of size and define the mean height of growing film and its roughness by:
[TABLE]
where is growing time and denotes an averaging over different samples, respectively. Moreover, growing time is a factor which can be applied to control the surface roughness of thin films.
Let us now calculate the roughness exponent of the growing surface. Starting from a flat interface (one of the possible initial conditions), it is conjectured that a scaling of lenght by factor and of time by factor ( is the dynamical scaling exponent), rescales the roughness by factor as follows Barabasi :
[TABLE]
which implies that
[TABLE]
For large and fixed saturate. However, for fixed and large and , one expects that correlations of the height fluctuations are set up only within a distance and thus must be independent of . This implies that for , with . Thus, dynamic scaling postulates that
[TABLE]
The roughness exponent and the dynamic exponent characterize the self-affine geometry of the surface and its dynamics, respectively. In the present work, we see the surfaces at the limit and so we will only obtain the exponent.
The common procedure to measure the roughness exponent of a rough surface is use of the surface structure function depending on the length scale which is defined by :
[TABLE]
It is equivalent to the statistics of height-height correlation function for stationary surfaces, i.e. . The second order structure function , scales with as .
III.2 Level Crossing Analysis
Let denotes the number of positive slope crossing of for interval L.
Since the process is homogeneous, if we take a second time interval of immediately following the first we shall obtain the same result, and for two intervals together we shall therefore obtain Tabar1 :
[TABLE]
from which it follows that, for a homogeneous process, the average number of crossing is proportional to the interval . Hence
[TABLE]
or
[TABLE]
where is the average frequency of positive slope crossing of the level . We now consider how the frequency parameter can be deduced from the underlying probability distributions for .
Consider a small length scale of a typical sample function. Since we are assuming that the process is a smooth function of , with no sudden ups and downs, if is small enough, the sample can only cross with positive slope if at the beginning of the interval . Furthermore, there is a minimum slope at if the level is to be crossed in interval depending on the value of at position . So there will be a positive crossing of in the next interval if, at ,
[TABLE]
As shown in Tabar1 , the frequency can be written in terms of joint PDF (probability distribution function ) of as follows:
[TABLE]
and then the quantity which is defined as:
[TABLE]
will measure the total number of crossing the surface with positive slope. So, the and square area of growing surface are in the same order. Concerning this, it can be utilized as another quantity to study further the roughness of a surface Tabar .
IV Results and Discussion
Thin films of were deposited by using thermal evaporation method and then surface micrographs of samples were obtained by AFM technique after annealed at different temperatures (Fig.3).
These micrographs were then analyzed using methods from stochastic data analysis have introduced in the last section. Figure 3 shows AFM images of thin films annealed at 200 ,300 ,350 ,400 ,450, and 500. The ”as deposited” and annealed sample at 200 (Fig. 3a) have columnar structure, indicating that up to 200 no significant changes in the microstructure occurs. However, at higher temperatures (figs. 3b-3f) we have observed increased grain size and rougher surface. Specifically at 500 (Fig. 3f) we observe stark changes in the micrograph which is accompanied by composition changes in the surface. This can be related to the phase transition to Magneli phase e.g. in the annealing process Mohammad . This is also confirmed by our XPS and UV-visible spectrophotometry analysis (Sec. II). These are shown the significant formation of state in the surface at 500.
Also our analysis shows that below 400 the surfaces are in amorphous phase with the same behavior for all scales, but as soon as the crystalline phase appears the system behaves differently which diagnostics at small and large scale for temperatures above 400. By using parameters of the analytical method given here, these transitions can be quantified.
Now, we will use the statistical parameters introduced in the last section and will obtain some quantitative information about the effect of annealing temperature on the surface topography of the samples.
The structure function as defined in Eq.(5) can be used to quantify the topology of a rough surface. The structure function is plotted against the length scale of the sample in Fig.4 . The saturated is an indication of the surface roughness, as . The most obvious observation indicates that roughness is raised with increasing annealing temperature. Roughness has a minimum of at 27 and 200 and a maximum of at 500. This is because higher temperatures create higher peaks (i.e. peaks with more deviations from the average) . All exponents which is derivable from have been summarized and given in Table I.
As depicted in Fig.4 , the structure function , has a different behavior in the various temperatures. So that, in the annealing temperature range 27-350 it has a typical behavior in all scales, but in the higher temperature range 400-500 its behavior is different in the small and large scales. In the other words, the phase transition is occurred at 400, because for higher temperatures, there are two sets of roughness parameters needed to simulate the surface morphology. It can be related to the phase transition in the structure of the surface from amorphous to crystalline phase has been yielded from the band gap energy (see the section II.B).
The slope of each curves at the small and large scales yields the roughness exponents and of the corresponding surface. Hence, it is seen that the mono roughness exponent increases with the addition of annealing temperature up to 400. In the higher temperatures, we have obtained two roughness exponents( -) equal to the 0.40-0.14, 0.71-0.20, and 0.69-0.24 for temperatures 400, 450 and 500, respectively. Difference in the values, in these temperatures, are in agreement with changes of correlation length. Where the correlation length, is the distance at which the structure function behaves differently.
The range of the scaling upon correlation length listed in the forth column in table I. The value of denotes the correlation length at small scales and for large scales. The higher value represents a smoother surface (as we expected from Fig.3). The correlation length obtained from the structure function is also a measure of minimum lateral size of surface features at each annealing temperature.
The another important film parameter is the effective area of the sample which has an important role in the gas sensitivity of surfaces. To obtain a measure for this, we utilize the level crossing analysis. As shown in Fig.5, the average frequency as a function of height , is plotted for the various annealing temperatures. The broad curves indicate the higher magnitude of height fluctuations around the average, and sharp curves show that the most of fluctuations are around the height average. This conclusion is in the correspondence with the results obtained from Fig.3.
According to the Eq.(11) i.e. The total number of the crossing surface with positive slope is proportional to the square of area of the growing surface. To obtain the optimum value of the effective area, we have calculated the ratio of effective areas with respect to the area of the ”as deposited” surface (27).
The values are given in the last column in Table I. It means, although the roughness increases by the annealing temperature but the effective square area of the rough surface has a maximum value of Tabar .
For more clarity, we have calculated the temperature dependence of normalized numerically (Fig.6) around 400, and we have obtained the three following functions for this quantity :
[TABLE]
[TABLE]
[TABLE]
According to this figure, the maximum value of the effective area is at 400 (with respect to its value at 27) with the relative value equal to 2.00. Thus, applying this analysis easily shows that if one follows the condition which the effective area as an important parameter in the gas sensitivity of surfaces is optimum and furthermore, the film composition has not been changed (e.g. The Magneli phase transition has not been occurred), should choose the annealed surface at 400 for better performance.
V Conclusions
We have investigated the role of annealing temperature, as an external parameter, to control the statistical properties of a rough surface. The AFM microstructure of the surfaces is just needed to apply in our analysis. We have computed the statistical quantities such as roughness exponent, roughness and lateral size of surface features of the ”as deposited” and annealed surfaces at 200, 300, 350, 400, 450, and 500, using the structure function. We have seen a phase transition at 400, because for higher temperatures there are two sets of roughness parameters, due to structural changes from amorphous to the crystalline phase. Moreover, using the level crossing analysis we have obtained an optimum annealing temperature, 400 in which the surface of the has maximum value about twice relative to the ”as deposited” film without any changes in the film composition that may increase surface reaction of the film as the gas sensor or photo-catalyst.
VI Acknowledgment
GRJ and AAS would like to thank S.M.Fazeli for his useful comments and especially M.R.Rahimitabar for his useful lectures on ”stochastic data analysis”. AZM would like to acknowledge research council of Sharif University of Technology for financial support of the work.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1(1) C. G. Granqvist, Handbook of Electrochromic Materials (Elsevier, Amsterdam, 1995).
- 2(2) P. R. Bueno, F. M. Pontes, E. R. Leite, L. O. S. Bulh es, P. S. Pizani, P. N. Lisboa-Filho, and W. H. Schreiner, J. Appl. Phys. 96, 2102 (2004).
- 3(3) R. Azimirad, O. Akhavan, and A. Z. Moshfegh, J. Electrochem. Soc.153, E 11(2006).
- 4(4) A. Siokou, S. Ntais, S. Papaefthimiou, G. Leftheriotis, and P. Yianoulis, Surface Science, 566/568, 1168 (2004).
- 5(5) S.-L. Kuai, G. Bader, and P. V. Ashrit, App. Phys. Lett., 86, 221110 (2005).
- 6(6) Y. Takeda, N. Kato, T. Fukano, A. Takeichi, and T. Motohiro, J. Appl. Phys., 96, 2417 (2004).
- 7(7) C. O. Avellaneda and L. O. S. Bulh es, Solid State Ionics, 165, 117 (2003).
- 8(8) S.-H. Lee, H. M. Cheong, P. Liu, D. Smith, C. Edwin Tracy, A. Mascanrenhas, J. R. Pitts, and S.K. Deb, J. Appl. Phys., 88, 3076 (2000).
