Investigation of TiO2 thin film deposited by microwave plasma assisted sputtering and its application in 3D glasses

: TiO 2 deposition using separate regions for sputtering and oxidation is not well investigated. We optimized process parameter for such as oxygen ﬂow and microwave power to produce high quality TiO 2 ﬁlters for Stereo/3D imaging applications. This deposition technique was chosen for its unique advantages: high deposition rates while increasing the probability of obtaining stoichiometric oxides, reduces possibility of target poisoning and provides better stability of process. Various characterization methods, such as scanning electron microscopy (SEM), atomic force microscopy (AFM), Raman, X-ray diffraction (XRD), transmission spectroscopy, were used in compliment to simulations for detailed analysis of deposited TiO 2 thin ﬁlms. Process parameters were optimized to achieve TiO 2 ﬁlms with low surface scattering and absorption for fabricating multi-passbands interference ﬁlter for 3D glasses. From observations and quantitative analysis of surfaces, it was seen that surface roughness increases while oxygen ﬂow or microwave power increases. As the content of anatase phase also increases with higher microwave power and higher oxygen ﬂow, the formation of anatase grains can cause higher surface roughness. Optical analysis of samples validates these trends and provided additional information for absorption trends. Optimized parameters for deposition process are then obtained and the ﬁnal fabricated 3D glasses ﬁlters showed high match to design, within 0.5% range for thickness error.

TiO 2 is a n-type semiconductor with a wide energy band gap and high refractive index with structural forms of amorphous, rutile, anatase and brookite (or mixtures); rutile is the most stable phase among these structures. The electronic and optical properties of TiO 2 film has be widely investigated experimentally [7] and theoretically [8]. As reported, Rutile has a direct forbidden gap (3.03 eV), which is almost degenerate with an indirect allowed transition (3.05 eV). Due to the low probability of the direct forbidden transition, the indirect allowed transition dominates the optical absorption just above the absorption edge. The refractive index, n D , of rutile is 2.609. Anatase has an indirect allowed transition at 3.2 eV and n D of 2.488. Brookite is similar to anatase in which they both have indirect allowed transition-as there are only minor differences in their local crystal environment. The n D of Brookite is 2.583. The amorphous and anatase phases of TiO 2 are preferred for the applications of optical coatings in visible range as they have wider band gap compared with rutile (narrower band gap of 3.05 eV at 406 nm wavelength) and higher absorption around 400 nm. Properties such as UV absorption and high refractive index in visible range are important for the design of interference optical coatings.
Various techniques have been used for the deposition of TiO 2 films, such as hydrothermal synthesis [9], plasma enhanced chemical vapor deposition (PECVD) [10], sol-gel [11], atomic layer deposition (ALD) [12], pulsed laser deposition (PLD) [13] and sputtering [7]. TiO 2 properties depend on the deposition process employed and process condition parameters. Using Monto Carlo simulation, Xiao et al. [14] have found that an atomically smooth surface can only be obtained within a certain deposition parameter window, at a regime of intermediate temperatures and low impingement rate. Tang et al. [7] produced single phase anatase and rutile using a reactive triode sputtering by controlling oxygen concentration and substrate temperature. These indicate that deposition methods and process parameters need to be carefully selected and optimized for a specific application as device performance will be significantly affected by the phase of TiO 2 and coating surface roughness.
Stereo/3D imaging has a wide range of applications, particularly in entertainment and games. Mehrabi et al. have provided an excellent review of 3D display including visual depth cues, technologies and their applications [15]. The concept of interest in this paper is binocular parallax of physiological depth cue, which fuses two slightly different images in human brain and provides 3D perception. This concept is currently widely used in TVs and cinemas. Various techniques can be used to realize to provide 3D perception based on binocular parallax [15,16], such as active stereo technique or passive stereo technique. The passive stereo technique uses polarizers or interference filters to separate images for left and right eye. Interference filter technique utilizes optical filter to pass only one or more specific wavelength bands and reflect the others, image separation is realized by wavelength difference. The advantages of interference filter technique for 3D display comparing to polarization technique include: nondepolarizing silver screens are not required [15]; image ghosting can be avoided; image resolution is maintained by using time sequential method even it is passive method [17]. Additional advantage is lower power consumption compared with polarization technique (when filters of glasses are well designed along with displayer, wavelength alignment) as polarizer only allows certain polarization light going through therefore reducing image brightness. Low power consumption becomes a more important factor as portable laser illuminated projector emerges on market [18]. On the other hand, this technique requires trained personnel to adjust the wavelengths of colors on the projectors; which increases costs [19].
Plasma assisted oxidation for a drum-based sputtering has been widely used in reactive thin film deposition processes for increasing productivity [20]. In this technique, the target material can be sputtered under metal-like conditions and further oxidized in highly reactive (oxygen or nitrogen) plasma area in a separated region; this allows high deposition rates while increasing the probability of obtaining stoichiometric oxides, furthermore this technique reduces possibility of target poisoning and provides better stability of process. Though TiO 2 sputtering processes have been investigated, the effects of process parameters on TiO 2 phases, surface structure and properties for this particular process are not well investigated. This paper investigated the influences of oxygen partial pressure and microwave power on TiO 2 structure and surface properties. Process parameters were optimized to achieve TiO 2 films with low surface scattering and absorption; then multi-passbands interference filter for 3D glasses were produced.

Experiments
Microwave plasma assisted magnetron pulse sputtering was chosen as the deposition process for this research. Two in-situ thickness monitoring methods are used in this system (broadband optical method and quartz crystal method), shown in Figure 1.

Experiments
Microwave plasma assisted magnetron pulse sputtering was chosen as the deposition process for this research. Two in-situ thickness monitoring methods are used in this system (broadband optical method and quartz crystal method), shown in Figure 1. In this system, process takes place in two separate areas for deposition and oxidation, substrate mounted on a rotating drum that passes through sputter area to form a metal-like thin layer and then passes through microwave oxygen plasma area for further oxidation to form a fully oxidized thin film layer. In order to optimize process for producing high quality TiO2 film for optical filters, a set of TiO2 samples are deposited under various process parameters: varying microwave (MW) power and oxygen flow as listed in Table 1, other parameters are kept same.  1  3  15  2  3  20  3  3  25  4  3  30  5  3  35  6  3  40  7  2.5  30  8  2  30 Samples are characterized using various techniques. Optical characterization data was acquired with a PerkinElmer Lambda 40 UV/Vis spectrometer (PerkinElmer, Waltham, MA, USA) measuring in transmission mode from 250-1100 nm. Raman and X-ray diffraction (XRD) were used for material and structure characterizations. Raman was done using Thermo Fisher Scientific-DXR Raman Microscope (Thermo Fisher Scientific, Waltham, MA, USA) equipped with lasers of 4 wavelengths (760 nm laser was used for this research); XRD was done using a Siemens D5100 Diffractometer (Siemens, Berlin, Germany) in locked-couple mode with 2θ range from 20°-30°. Scanning electron microscopy (SEM) and atomic force microscope (AFM) were also carried out for surface morphology analysis. Here, cross-sectional and top surface images were obtained using a Hitachi S4100 SEM (Hitachi, Tokyo, Japan) and the Dimension Icon AFM (Bruker, Billerica, MA, USA) in contact mode (with the tip OESTPA-R3 fitted) was used to analyses surface roughness. In this system, process takes place in two separate areas for deposition and oxidation, substrate mounted on a rotating drum that passes through sputter area to form a metal-like thin layer and then passes through microwave oxygen plasma area for further oxidation to form a fully oxidized thin film layer. In order to optimize process for producing high quality TiO 2 film for optical filters, a set of TiO 2 samples are deposited under various process parameters: varying microwave (MW) power and oxygen flow as listed in Table 1, other parameters are kept same. Samples are characterized using various techniques. Optical characterization data was acquired with a PerkinElmer Lambda 40 UV/Vis spectrometer (PerkinElmer, Waltham, MA, USA) measuring in transmission mode from 250-1100 nm. Raman and X-ray diffraction (XRD) were used for material and structure characterizations. Raman was done using Thermo Fisher Scientific-DXR Raman Microscope (Thermo Fisher Scientific, Waltham, MA, USA) equipped with lasers of 4 wavelengths (760 nm laser was used for this research); XRD was done using a Siemens D5100 Diffractometer (Siemens, Berlin, Germany) in locked-couple mode with 2θ range from 20 • -30 • . Scanning electron microscopy (SEM) and atomic force microscope (AFM) were also carried out for surface morphology analysis. Here, cross-sectional and top surface images were obtained using a Hitachi S4100 SEM (Hitachi, Tokyo, Japan) and the Dimension Icon AFM (Bruker, Billerica, MA, USA) in contact mode (with the tip OESTPA-R3 fitted) was used to analyses surface roughness.

Dispersive Model Used for Transmittance Fitting
Dispersive function, n/k or dielectric constants versus wavelength, is needed for transmittance data fitting to obtain the n and k values for TiO 2 thin film required for thin film optical filter design. TiO 2 crystal is a direct band-gap semiconductor for rutile and brookite, anatase is an indirect band-gap semiconductor as the bottom of the conduction band (CB) is at Γ and the top of the valence band (VB) is at M [8]. However, the energy at Γ is only 0.18 eV lower than the top of the VB, therefore, Mo et al. [8] approximated anatase as a direct band-gap semiconductor. The eventual films produced in this work are mainly amorphous whose optical absorption behaves differently from a perfect crystalline semiconductor. The absorption edge of crystal semiconductor terminates abruptly at the energy gap, whereas the absorption of amorphous semiconductor does not stop at the energy gap and extends into the energy gap region. Thus the O'Leary-Johnson-Lim Model (OJL) model, a model for amorphous semiconductor, was employed to describe optical dispersion of dielectric function and carry out transmittance fitting [21]. OJL model assumes that the distribution of states/Density of Sates (DOS), N(E), exhibits a square-root functional dependence in the band region and an exponential functional dependence in the tail region, described by the following equations: The diagram of energy band structure of amorphous semiconductor as described by these equations are also shown in Figure 2, where the DOS of a perfect crystal follows a parabolic curve and its electronic band gap in can be derived as: V c and V v are ground state energies for the CB and VB. In amorphous semiconductors, the parameters of γ v and γ c are used to describe the tail-states for the valence and conduction bands respectively. The physical meaning γ v and γ c is disorder degree of amorphous material. For example, while γ c → 0, Equation (1) of conduction band DOS function gives the expression Equation (4). Similar results are also observed for the valence-band DOS function.
Using OJL model, the absorption coefficient of amorphous semiconductor can be obtained and therefore derive the extinction coefficient (k). The results of k can then be used to obtain the refractive index (n) by using the Kramers-Kronig relations (KKR). n and k against wavelength or energy of photon, dispersive function of amorphous semiconductor, is finally obtained. This dispersive function is then used to fit measured transmittance data and obtain the fitting parameters used in OJL model. Code software (W Theiss, version 3.5) was used to carry out transmittance fitting.
One important parameter, band gap E g0 , can be obtained by fitting measured transmittance data using OJL model. However, Tauc Plot [22] is also a commonly used method to find band-gap for crystal semiconductors using Equation (5): Here α is absorption coefficient; hγ is photon energy; β is a constant; r is transition related parameter where r = 1/2 for direct allowed transitions, r = 3/2 for direct forbidden transitions, r = 2 for indirect allowed transitions and r = 3 for indirect forbidden transitions.
An extrapolation of a linear function of (αhγ) 1/2 against hγ [23] is used to define the energy gap of amorphous semiconductors. If using Equation (5), r = 2 can be used for amorphous transitions, which is same as indirect allowed transitions for crystal semiconductor. In this paper, energy band gaps will be obtained using both OJL model and Tauc plot as a comparison.
The diagram of energy band structure of amorphous semiconductor as described by these equations are also shown in Figure 2, where the DOS of a perfect crystal follows a parabolic curve and its electronic band gap in can be derived as: Vc and Vv are ground state energies for the CB and VB. In amorphous semiconductors, the parameters of γv and γc are used to describe the tail-states for the valence and conduction bands respectively. The physical meaning γv and γc is disorder degree of amorphous material. For example, while γc → 0, Equation (1) of conduction band DOS function gives the expression Equation (4). Similar results are also observed for the valence-band DOS function. Figure 2. The diagram of energy band structure of amorphous semiconductor as described by Equations (1) and (2), thin solid curve: for crystal band structure; thick dash curve: for amorphous band structure.  (1) and (2), thin solid curve: for crystal band structure; thick dash curve: for amorphous band structure.

Rough Surface and Light Scattering
The methods used for surface roughness treatment and the calculation of light scattering are based on the method from Rakels' paper [24], however, specular transmittance were derived instead of reflectance (Rakel et al. derived equations for reflectance). Light scattering can be considered as coherent scattering if surface roughness is much smaller than wavelength of measurement light. The electric field of light needs to be calculated by considering the phase modulation of electric field by rough surface profile. For simplification, normal incidence illumination is assumed in this paper as all transmittance measurements were carried out at normal incident in this research.
Gaussian distribution is used to describe surface roughness: In this function, y is surface height, the mean distribution is zero, R q is the standard deviation. It can be found that R q is also the root mean squared roughness of rough surface by integrating y 2 with Gaussian distribution.
Set y = 0 position as phase reference, where the initial phase is 0. Then the phase difference comparing to y = 0 is ∆ = 2πy(n − 1)/λ for the position with height y.
The specular transmittance measured at normal direction is the ratio of light intensity, which is the squared modulus of electric field, then: Here T s is the specular transmittance for normal incident for rough surface, T 0 is the transmittance of smooth surface.
For a rough surface, light scattering effect is not the only effects to be considered. The rough surface also needs to be considered as a new layer consisting of a mixture of air and material (film or substrate). For example, Chen and Sun reported anti-reflection effect of nanostructured surface [25]. Effective Medium Approximation (EMA) was used to model the mixture layer of a rough surface and the Bruggeman model of EMA [26] is shown, Equation (9). The volume ratio of air to material can be assume 50% as roughness distribution is symmetric (Gaussian distribution).
ε 1 and ε 2 are dielectric function for air and material; ε is effective dielectric function of the mixture; and f 1 and f 2 are volume ratio, which is 50% here.

Design of Triple Bandpass Filter for 3D Glasses
Though the filter for 3D glasses is triple passbands, the design here is started from a single bandpass filter then advances to multi bandpass filter design. One of the most basic design formulas for a single bandpass filter is a single-cavity design which uses only one Fabry-Perot cavity [27]: substrate |(HL) n 2H(LH) n | air, where H and L stand for high and low refractive index materials respectively. In this paper TiO 2 and SiO 2 were used as the high and low refractive index materials respectively (the index of TiO 2 is higher than the indices of ZrO 2 , Nb 2 O 5 and HfO 2 , therefore providing high difference between H and L layer which may benefit the interference filter design). To improve the performance of bandpass filter, multi-cavities were used, the transmittance band approaches a rectangle as the number of cavities increases. To achieve an approximate multi-passbands transmittance shape as a start design, an adapted structure of multi-cavities of (xHyL) n 2H (xLyH) n was used. One example here: high index = 2.5 (roughly sputtered TiO 2 index at 550 nm) and low index = 1.45 (roughly sputtered SiO 2 index at 550 nm) and structure is [(1.5H0.5L) 2 2H (0.5L1.5H) 2 L] 3 (1.5H0.5L) 2 2H (0.5L1.5H) 2 with reference wavelength of 650 nm for design of side 1 of 3D glasses. Using this design as an initial structure, an optimized design was found using TFCalc software (version 3.5) and needle method within a few minutes. Transmittances of this design before and after optimization are shown in Figure 3a. Similarly, the triple passband filter design for side 2 of 3D glasses is also achieved, with reference wavelength at 700 nm. The optical density (OD = −log(T)) performance are shown in Figure 3b for designs of both sides: ODs are higher than 3 for all necessary rejection bands.
ε and ε are dielectric function for air and material; ε is effective dielectric function of the mixture; and and are volume ratio, which is 50% here.

Design of Triple Bandpass Filter for 3D Glasses
Though the filter for 3D glasses is triple passbands, the design here is started from a single bandpass filter then advances to multi bandpass filter design. One of the most basic design formulas for a single bandpass filter is a single-cavity design which uses only one Fabry-Perot cavity [27]: substrate |(HL) n 2H(LH) n | air, where H and L stand for high and low refractive index materials respectively. In this paper TiO2 and SiO2 were used as the high and low refractive index materials respectively (the index of TiO2 is higher than the indices of ZrO2, Nb2O5 and HfO2, therefore providing high difference between H and L layer which may benefit the interference filter design). To improve the performance of bandpass filter, multi-cavities were used, the transmittance band approaches a rectangle as the number of cavities increases. To achieve an approximate multi-passbands transmittance shape as a start design, an adapted structure of multi-cavities of (xHyL) n 2H (xLyH) n was used. One example here: high index = 2.5 (roughly sputtered TiO2 index at 550 nm) and low index = 1.45 (roughly sputtered SiO2 index at 550 nm) and structure is [(1.5H0.5L) 2 2H (0.5L1.5H) 2 L] 3 (1.5H0.5L) 2 2H (0.5L1.5H) 2 with reference wavelength of 650 nm for design of side 1 of 3D glasses. Using this design as an initial structure, an optimized design was found using TFCalc software (version 3.5) and needle method within a few minutes. Transmittances of this design before and after optimization are shown in Figure 3a. Similarly, the triple passband filter design for side 2 of 3D glasses is also achieved, with reference wavelength at 700 nm. The optical density (OD = −log(T)) performance are shown in Figure 3b for designs of both sides: ODs are higher than 3 for all necessary rejection bands.

TiO2 Thin Film Process Optimization
Two groups of TiO2 samples were fabricated as described in the experimental session: fixed microwave power at 3 kW and varied oxygen flow; fixed oxygen flow at 30 sccm and varied

TiO2 Thin Film Process Optimization
Two groups of TiO 2 samples were fabricated as described in the experimental session: fixed microwave power at 3 kW and varied oxygen flow; fixed oxygen flow at 30 sccm and varied microwave power. Transmittance of samples were measured. Figure 4a shows the transmittance for the samples with varying oxygen flows. It is clearly divided into 2 groups (insert of Figure 4a): higher light loss (dotted lines, 30, 35 and 40 sccm oxygen) and lower light loss (solid lines, 15, 20 and 25 sccm oxygen). In the second group, the sample deposited at 15 sccm oxygen has less transmittance than others. Figure 4b shows the transmittance for the samples with varying microwave power. It shows that transmittance varies with microwave power, sample at 3 kW microwave power has highest light loss. There are two possible reasons for light loss in transmittance measurement: absorption and scattering. To identify the reasons, further characterizations are carried out: XRD, Raman, SEM and AFM.  Cross-section and top surface of TiO2 thin films are observed using SEM. The trends of surface morphology are: surface roughness increases while oxygen flow or microwave power increases. Figure 5 shows some SEM images of TiO2 thin films deposited under various oxygen flows and microwave power. However, SEM cannot give quantitatively surface roughness. AFM analysis ( Figure 6) was also employed and the AFM results confirm the trends seen using SEM: surface roughness increases while oxygen flow or microwave power increases. Mean square root of surface roughness Rq of samples are listed in Table 2 and the Rq of bare substrate is 0.31 nm. Cross-section and top surface of TiO 2 thin films are observed using SEM. The trends of surface morphology are: surface roughness increases while oxygen flow or microwave power increases. Figure 5 shows some SEM images of TiO 2 thin films deposited under various oxygen flows and microwave power.  Cross-section and top surface of TiO2 thin films are observed using SEM. The trends of surface morphology are: surface roughness increases while oxygen flow or microwave power increases. Figure 5 shows some SEM images of TiO2 thin films deposited under various oxygen flows and microwave power. However, SEM cannot give quantitatively surface roughness. AFM analysis (Figure 6) was also employed and the AFM results confirm the trends seen using SEM: surface roughness increases while oxygen flow or microwave power increases. Mean square root of surface roughness Rq of samples are listed in Table 2 and the Rq of bare substrate is 0.31 nm. However, SEM cannot give quantitatively surface roughness. AFM analysis ( Figure 6) was also employed and the AFM results confirm the trends seen using SEM: surface roughness increases while oxygen flow or microwave power increases. Mean square root of surface roughness R q of samples are listed in Table 2 and the R q of bare substrate is 0.31 nm.  While it is intuitive to understand how the roughness varies with microwave power as microwave energy will roughen the surface, effects of oxygen content on roughness needed further analysis: XRD and Raman analysis were done.
XRD can be used to analyze the phase structure of TiO2 thin films, such as anatase, rutile or amorphous. According XRD database, the strongest peaks of XRD (Cu Kα) are at 27.46° of 2θ for (110) rutile and at 25.35° of 2θ for (101) anatase [28]. Raman analysis can also demonstrate the phase structures of TiO2 and their Raman active modes are well characterized. The anatase phase has six Raman active modes: a strong peak at 144 cm −1 , followed by peaks located at 197, 399, 513, 519 and 641 cm −1 ; while the rutile phase has four active modes at 143, 236, 447 and 613 cm −1 [29].
XRD and Raman spectra for samples deposited under various oxygen flows at 3 kW microwave power are shown in Figure 7. In XRD, samples at 40 and 35 sccm oxygen flow showed anatase peak (25.35°), while the other samples did not show crystalline peaks. It can be seen that higher oxygen concentrations promote anatase formation, this agrees with the trend reported by Tang et al. [7] where pure anatase phase are obtained with oxygen concentration of 20%-30% and rutile phase obtained for oxygen concentration below 20%. In Raman spectra, samples at 40, 35 and 30 sccm oxygen flow showed anatase peak (144 cm −1 ). This confirms the XRD results. In addition, the higher sensitivity of Raman spectroscopy compared with XRD for thin film measurements, showed detection of anatase phase in the 30 sccm sample as well. The peak at 520 cm −1 seen in Figure 7b is due to crystalline silicon substrate.
Samples with varied microwave power and fixed oxygen flow (30 sccm) did not show any peak in XRD measurements. For Raman measurements, anatase peak (144 cm −1 ) was seen for samples at 3 and 2.5 kW but no anatase peak was seen for sample at 2 kW. This indicates that higher microwave  While it is intuitive to understand how the roughness varies with microwave power as microwave energy will roughen the surface, effects of oxygen content on roughness needed further analysis: XRD and Raman analysis were done.
XRD can be used to analyze the phase structure of TiO 2 thin films, such as anatase, rutile or amorphous. According XRD database, the strongest peaks of XRD (Cu Kα) are at 27.46 • of 2θ for (110) rutile and at 25.35 • of 2θ for (101) anatase [28]. Raman analysis can also demonstrate the phase structures of TiO 2 and their Raman active modes are well characterized. The anatase phase has six Raman active modes: a strong peak at 144 cm −1 , followed by peaks located at 197, 399, 513, 519 and 641 cm −1 ; while the rutile phase has four active modes at 143, 236, 447 and 613 cm −1 [29].
XRD and Raman spectra for samples deposited under various oxygen flows at 3 kW microwave power are shown in Figure 7. In XRD, samples at 40 and 35 sccm oxygen flow showed anatase peak (25.35 • ), while the other samples did not show crystalline peaks. It can be seen that higher oxygen concentrations promote anatase formation, this agrees with the trend reported by Tang et al. [7] where pure anatase phase are obtained with oxygen concentration of 20%-30% and rutile phase obtained for oxygen concentration below 20%. In Raman spectra, samples at 40, 35 and 30 sccm oxygen flow showed anatase peak (144 cm −1 ). This confirms the XRD results. In addition, the higher sensitivity of Raman spectroscopy compared with XRD for thin film measurements, showed detection of anatase phase in the 30 sccm sample as well. The peak at 520 cm −1 seen in Figure 7b is due to crystalline silicon substrate. Based the analysis results from XRD, Raman and AFM, it can be concluded that the films contain a mixture of anatase crystal grains and amorphous TiO2. Samples with higher anatase phase content also showed higher surface roughness, thus deducing that the formation of anatase grains can cause higher surface roughness. According to above analysis, higher oxygen flow and higher microwave power will lead higher surface roughness. This increased surface roughness will reduce optical performance thin film. However, low oxygen flow and microwave power may result in partially oxidized thin film, causing high optical absorption. A model was constructed by considering surface roughness and OJL model for analysis, using layer stack of substrate/ideal TiO2 layer/rough surface layer. For the rough surface layer, light scattering is taken into account while it is also treated as a thin film layer with complex refractive index (n and k) derived by using Bruggeman model (Equation (9)) with 50% volume fraction of air and TiO2 mixture. The dispersive dielectric function of TiO2 consists of dielectric background and one OJL model. Combining layer stack and dispersive function, measured transmittances of TiO2 thin film samples were fitted. Parameter of layers and dispersive function were obtained and used to calculate n and k versus wavelength. For conciseness, only the parameters (obtained through data fitting) which are used for later discussion are listed in Table 3. The TiO2 thicknesses, obtained through fitting, for all samples are around 490 nm. The optical fitted surface roughness Rq, obtained by using the Equation (8) in Section 3.2, are listed in Table 3, showing the same trend and similar level as the values obtained using AFM. Smoother surface is better for optical filter; however, the k value of film is also an important factor to consider when deciding process parameters; low k value give better optical filter performance. It is clear, therefore, that the best process parameters are 30 sccm oxygen flow and 2 kW microwave power for TiO2 thin film deposition for our specific drum-based system.  Samples with varied microwave power and fixed oxygen flow (30 sccm) did not show any peak in XRD measurements. For Raman measurements, anatase peak (144 cm −1 ) was seen for samples at 3 and 2.5 kW but no anatase peak was seen for sample at 2 kW. This indicates that higher microwave power promotes anatase formation.
Based the analysis results from XRD, Raman and AFM, it can be concluded that the films contain a mixture of anatase crystal grains and amorphous TiO 2 . Samples with higher anatase phase content also showed higher surface roughness, thus deducing that the formation of anatase grains can cause higher surface roughness.
According to above analysis, higher oxygen flow and higher microwave power will lead higher surface roughness. This increased surface roughness will reduce optical performance thin film. However, low oxygen flow and microwave power may result in partially oxidized thin film, causing high optical absorption. A model was constructed by considering surface roughness and OJL model for analysis, using layer stack of substrate/ideal TiO 2 layer/rough surface layer. For the rough surface layer, light scattering is taken into account while it is also treated as a thin film layer with complex refractive index (n and k) derived by using Bruggeman model (Equation (9)) with 50% volume fraction of air and TiO 2 mixture. The dispersive dielectric function of TiO 2 consists of dielectric background and one OJL model. Combining layer stack and dispersive function, measured transmittances of TiO 2 thin film samples were fitted. Parameter of layers and dispersive function were obtained and used to calculate n and k versus wavelength. For conciseness, only the parameters (obtained through data fitting) which are used for later discussion are listed in Table 3. The TiO 2 thicknesses, obtained through fitting, for all samples are around 490 nm. The optical fitted surface roughness R q , obtained by using the Equation (8) in Section 3.2, are listed in Table 3, showing the same trend and similar level as the values obtained using AFM. Smoother surface is better for optical filter; however, the k value of film is also an important factor to consider when deciding process parameters; low k value give better optical filter performance. It is clear, therefore, that the best process parameters are 30 sccm oxygen flow and 2 kW microwave power for TiO 2 thin film deposition for our specific drum-based system.
As discussed, gap energy, E g0 , is one of the important physical parameters for optical thin films. OJL gap energies for TiO 2 thin films are listed in Table 3 as well. Values are higher than reported data. For comparison, Tauc Plot is also used to extrapolate gap energies where r = 2 for Equation (5) where α values were calculated from fitted k values using α = 4πk/λ. One of extrapolation fitting is shown in Figure 8 as example and gap energies obtained using Tauc Plot are listed in Table 4. These values are much closer to reported values for anatase phase, though still slightly higher. As discussed, gap energy, Eg0, is one of the important physical parameters for optical thin films. OJL gap energies for TiO2 thin films are listed in Table 3 as well. Values are higher than reported data. For comparison, Tauc Plot is also used to extrapolate gap energies where r = 2 for Equation (5) where α values were calculated from fitted k values using α = 4πk/λ. One of extrapolation fitting is shown in Figure 8 as example and gap energies obtained using Tauc Plot are listed in Table 4. These values are much closer to reported values for anatase phase, though still slightly higher.

Fabrication of 3D Glasses Filter
Following discussions form the above sections, the process parameters for TiO2 deposition were chosen to be 30 sccm oxygen flow and 2 kW microwave power. The n and k values for TiO2 and SiO2 used to produce 3D glasses filters are shown in Figure 9. These indices were used for replacing the constant index (2.5 for H layers and 1.45 for L layers in Section 3) in the initial designs. Then the initial designs of 3D glasses filters were re-optimized. The final design for side 1 has 49 layers with 4 um total thickness. The thinnest layer for side 1 is 24 nm. The final design for side 2 has 45 layers with 4.2 um total thickness. The thinnest layer for side 2 is 30 nm. The filters are deposited onto fused silica substrates, which are type JGS3 with index of 1.46 at wavelength of 550 nm and negligible absorption in the wavelength range of interest. Quartz crystals were used for in-situ thickness monitoring. Figure 10 shows the transmittance comparison between designs and measurements. It indicates thickness errors for both sides were controlled within the 0.5% range. It can also be seen that the transmittance discrepancies between design and measurement is larger for shorter wavelengths, this is due to more significant light scattering and absorption effects at shorter wavelengths. This effect can be more

Fabrication of 3D Glasses Filter
Following discussions form the above sections, the process parameters for TiO 2 deposition were chosen to be 30 sccm oxygen flow and 2 kW microwave power. The n and k values for TiO 2 and SiO 2 used to produce 3D glasses filters are shown in Figure 9. These indices were used for replacing the constant index (2.5 for H layers and 1.45 for L layers in Section 3) in the initial designs. Then the initial designs of 3D glasses filters were re-optimized. The final design for side 1 has 49 layers with 4 um total thickness. The thinnest layer for side 1 is 24 nm. The final design for side 2 has 45 layers with 4.2 um total thickness. The thinnest layer for side 2 is 30 nm. The filters are deposited onto fused silica substrates, which are type JGS3 with index of 1.46 at wavelength of 550 nm and negligible absorption in the wavelength range of interest. Quartz crystals were used for in-situ thickness monitoring. Figure 10 shows the transmittance comparison between designs and measurements. It indicates thickness errors for both sides were controlled within the 0.5% range. It can also be seen that the transmittance discrepancies between design and measurement is larger for shorter wavelengths, this is due to more significant light scattering and absorption effects at shorter wavelengths. This effect can be more dominated by scattering.  For numerical evaluation, the difference between a single layer and a multiple-layer structure need to be distinguished with regards to the characteristics of surface roughness. In Equation (8), only RMS (definition) of roughness parameter was considered, which would only be valid when the length between the lateral features of rough surface is very big, for example, the single layer surface roughness with reasonable thin thickness. However, for a multi-layer structure, particularly 45 layers or 49 layers structure in this paper, the feature of rough surface at lateral direction is much denser as roughness need to be considered for each interface and their effects accumulate. Thus, an additional parameter is used to define an arbitrary rough surface: autocovariance (ACV) length (l), describing lateral feature of a rough surface [30]. For normal incident of measurement light, when l is much smaller than the wavelength (λ) of incident light, Equation (8) can be re-written by considering the effect of l for transmission using similar as used by Harvey et al. [30]:  For numerical evaluation, the difference between a single layer and a multiple-layer structure need to be distinguished with regards to the characteristics of surface roughness. In Equation (8), only RMS (definition) of roughness parameter was considered, which would only be valid when the length between the lateral features of rough surface is very big, for example, the single layer surface roughness with reasonable thin thickness. However, for a multi-layer structure, particularly 45 layers or 49 layers structure in this paper, the feature of rough surface at lateral direction is much denser as roughness need to be considered for each interface and their effects accumulate. Thus, an additional parameter is used to define an arbitrary rough surface: autocovariance (ACV) length (l), describing lateral feature of a rough surface [30]. For normal incident of measurement light, when l is much smaller than the wavelength (λ) of incident light, Equation (8) can be re-written by considering the effect of l for transmission using similar as used by Harvey et al. [30]: For numerical evaluation, the difference between a single layer and a multiple-layer structure need to be distinguished with regards to the characteristics of surface roughness. In Equation (8), only RMS (definition) of roughness parameter was considered, which would only be valid when the length between the lateral features of rough surface is very big, for example, the single layer surface roughness with reasonable thin thickness. However, for a multi-layer structure, particularly 45 layers or 49 layers structure in this paper, the feature of rough surface at lateral direction is much denser as roughness need to be considered for each interface and their effects accumulate. Thus, an additional parameter is used to define an arbitrary rough surface: autocovariance (ACV) length (l), describing lateral feature of a rough surface [30]. For normal incident of measurement light, when l is much smaller than the wavelength (λ) of incident light, Equation (8) can be re-written by considering the effect of l for transmission using similar as used by Harvey et al. [30]: And, the transmission losses due to scattering when l << λ can be obtained to be: For a more generic format of the equation, it can be written with n as power for wavelength: To summarize, Equation (8) is valid for describing scattering effects when l >> λ, while Equations (13) and (14) are better for describing scattering effects when l << λ.
To compare transmission losses, data was extracted from Figure 10 where ∆T is the difference of transmission percentage between design and measurement. This data was plotted and compared with the estimations using Equations (13) and (14). The plotted graph, Figure 11, shows that when n = 4, the trends agree very well, while n = 6 gave the best fitting.
For a more generic format of the equation, it can be written with n as power for wavelength: To summarize, Equation (8) is valid for describing scattering effects when l >> λ, while Equations (13) and (14) are better for describing scattering effects when l << λ.
To compare transmission losses, data was extracted from Figure 10 where ΔT is the difference of transmission percentage between design and measurement. This data was plotted and compared with the estimations using Equations (13) and (14). The plotted graph, Figure 11, shows that when n = 4, the trends agree very well, while n = 6 gave the best fitting. Figure 11. Comparison of transmission losses extracted from Figure 10, to transmission losses due to scattering using Equation (14), with different n values.

Conclusions
Two groups of TiO2 samples were fabricated as described in the experimental session: fixed microwave power at 3 kW and varied oxygen flow; fixed oxygen flow at 30 sccm and varied microwave power. Their analysis, using SEM and AFM, showed surface morphology trend: surface roughness increases while oxygen flow or microwave power increases. Further analysis done using XRD and Raman showed presence of anatase phase. The results also indicated that higher microwave power and higher oxygen flow promotes anatase formation. Combining the analysis results from XRD, Raman and AFM, it can be concluded that the samples with higher anatase phase content also showed higher surface roughness: thus, the formation of anatase grains can cause higher surface roughness.
Optical analysis of samples validates that higher oxygen flow and higher microwave power will lead to higher surface roughness. A compromise between surface roughness and fully oxidizing thin film required a modelling of process parameters: oxygen flow rate and microwave power. A model constructed by considering surface roughness and OJL model was used for analysis, to eventually Figure 11. Comparison of transmission losses extracted from Figure 10, to transmission losses due to scattering using Equation (14), with different n values.

Conclusions
Two groups of TiO 2 samples were fabricated as described in the experimental session: fixed microwave power at 3 kW and varied oxygen flow; fixed oxygen flow at 30 sccm and varied microwave power. Their analysis, using SEM and AFM, showed surface morphology trend: surface roughness increases while oxygen flow or microwave power increases. Further analysis done using XRD and Raman showed presence of anatase phase. The results also indicated that higher microwave power and higher oxygen flow promotes anatase formation. Combining the analysis results from XRD, Raman and AFM, it can be concluded that the samples with higher anatase phase content also showed higher surface roughness: thus, the formation of anatase grains can cause higher surface roughness.
Optical analysis of samples validates that higher oxygen flow and higher microwave power will lead to higher surface roughness. A compromise between surface roughness and fully oxidizing thin film required a modelling of process parameters: oxygen flow rate and microwave power. A model constructed by considering surface roughness and OJL model was used for analysis, to eventually obtain n, k and surface roughness parameters. The roughness parameters calculated from transmittance measurements agrees with AFM results. Balancing n, k and surface roughness for high performance optical thin films, the parameters for our TiO 2 process was chosen to be 30 sccm oxygen flow and 2 kW microwave power. The modelling done using OJL also provided information on gap energy, values obtained are higher than reported data. Tauc Plot is also used to extrapolate gap energies, these values are much closer to reported values for anatase phase, though still slightly higher.
Following analysis, the process parameters for TiO 2 deposition were chosen to be 30 sccm oxygen flow and 2 kW microwave power for fabricating 3D glasses filters. The transmittance measurement of the filters matched the designs very well: transmittance comparison between designs and measurements showed that thickness errors for both sides was controlled within 0.5% range.