Abstract

A novel procedure based on the use of the Boltzmann equation to model the parameter, the film deposition rate, and the optical band gap of BaxSr1−xTiO3 thin films is proposed. The BaxSr1−xTiO3 films were prepared by RF cosputtering from BaTiO3 and SrTiO3 targets changing the power applied to each magnetron to obtain different Ba/Sr contents. The method to calculate consisted of fitting the angular shift of (110), (111), and (211) diffraction peaks observed as the density of substitutional Ba2+ increases in the solid solution when the applied RF power increases, followed by a scale transformation from applied power to parameter using the Boltzmann equation. The Ba/Sr ratio was obtained from X-ray energy dispersive spectroscopy; the comparison with the X-ray diffraction derived composition shows a remarkable coincidence while the discrepancies offer a valuable diagnosis on the sputtering flux and phase composition. The proposed method allows a quick setup of the RF cosputtering system to control film composition providing a versatile tool to optimization of the process.

1. Introduction

The BaxSr1−xTiO3 solid solution (BST-SS) is a ceramic material with a great potential to be used in dynamic random access memories (DRAM) and integrated optic devices. It is known that ferroelectric and optical properties of BST-SS are stoichiometry dependent, for example, the electrical resistance, hysteresis cycle, and resistive switching [1], extinction coefficient and refraction index [2], the relative dielectric constant, breakdown strength, and energy density of multifunctional fibers [3], and the remanent polarization [4]. BST-SS usually crystallizes in a cubic lattice changing its lattice parameter between and associated with the extreme compositions: BaTiO3 (BTO) and SrTiO3 (STO) [ICCD, PDF International cards 00-075-0215 (BTO); 00-040-1500 (STO)]. The progressive Sr substitution in the BaxSr1−xTiO3 solid solution leads to a lattice contraction, reflected in the shift of Bragg angles, as reported in BST materials prepared by diverse techniques such as arc discharge [5], sol-gel [68] pulsed laser deposition [7], hydrothermal reaction [9], and chemical synthesis [10]. As much as both titanates are miscible in each other, the peak shift is linear and continuous between the extreme compositions; that is, it follows Vegard’s law. However, the presence of the metastable tetragonal phase is possible due to the presence of oxygen vacancies [11, 12]. On the other hand, reactive sputtering modelling has been widely studied since the early work of Berg et al. [13]. Berg’s and other models such as that developed by Strijckmans and Depla [14] take into account the nonlinearity inherent to the sputter deposition process involving reactive species. A general conclusion extracted from the models is that the hysteresis due to the change in the target (poisoned target) cannot be eliminated by the power fed to the RF system because the thickness of the compound film varies in proportion to the power and also depends on the reactive gas flux. The underlying mechanisms of the sputtering process using ceramic targets have been described since the early 1980s [15]. An accepted description of the process for ceramic targets is that the species are evaporated as individual atoms with different sputtering yields [16] causing an initial enrichment with the species of smaller yield, although a stationary point on the target composition compared with the sputtering flux is eventually reached. Usually oxygen is added to increase the control of film stoichiometry [11], although there is not yet an established procedure to control the sputtering process. Figure 1 depicts a scheme of the Boltzmann sigmoidal equation, where before and after of the exponential growth there are two zones of stabilization distinguished by a central turning point, , , and , respectively. This function has also been used to describe other phenomena such as phase transition in smart gels [17], the relationship between reaction kinetics and its dynamic influence on viscosity in the chemorheology behavior of photopolymerizable acrylates [18], and the stem and rosebud growth in roses [19], among others.

Boltzmann’s sigmoidal equation can be written as (1) when it is increasing and as (2) when it is decreasing.where and are final values of the dependent variable, the independent variable, is the value of at the inflection point, and is a coefficient that describes the behavior of the slope of the process during the transition.

In the present contribution, a further development of the proposed Boltzmann modelling for BST-SS cosputtered films is presented in order to predict film properties such as film thickness and optical band gap, using as initial input the information extracted from X-ray diffractograms. The Boltzmann model is followed closely when the lattice parameter is plotted as a function of applied power.

2. Experimental Details

BST-SS films were prepared by RF magnetron sputtering in a system equipped with two magnetrons: BaTiO3 (99.95%, SCI Engineered Materials Inc.) and SrTiO3 (99.9%, SCI Engineered Materials Inc.). Targets were 2′′ diameter and 0.125′′ thick. Before film deposition, the sputtering chamber was evacuated to a base pressure around  Pa; then, an Ar flushing was done filling the chamber to a pressure of 3.9 Pa during 10 minutes. For the film deposition, an Ar + O2 gas mixture was introduced into the chamber with an Ar/O2 = 90/10 ratio at an initial pressure of 6.6 Pa to ignite the plasma and perform a target presputtering during 15 minutes. Thereafter, the working pressure was set at 3.9 Pa to carry out the deposition. Quartz substrates  cm2 were successively rinsed with trichloroethylene, acetone, and ethanol before film deposit. A stainless-steel substrate holder was fixed at a distance of 8 cm from the magnetron in off-axis configuration. The substrate holder was rotated at 100 rpm for ensuring film uniformity; the substrate temperature was set at 549°C. The total applied RF-magnetron power was 120 watts, distributed between the two magnetrons as shown in Table 1, in order to produce BaxSr1−xTiO3 films with different compositions. The sputter deposit time for all of samples was 68 minutes.

Chemical composition was obtained by energy dispersive X-ray spectroscopy (EDS) employing a Jeol JSM-5300 electron microscope equipped with a Kevex Delta 1 spectrometer. The optical transmission spectra of the films were obtained in a Perkin-Elmer, Lambda 40 Spectrophotometer in the range of 250 to 800 nm. Film thicknesses were calculated from the transmittance spectra using the SCOUT software [20]. The X-ray diffractograms were acquired in a Phillips X’Pert diffractometer using CuK = 1.54060 Å, in a 2 theta range 20–95° with a step size of 0.02° and step time 1 s in the grazing incidence geometry with an incidence angle 0.5°.

3. Results and Discussion

Figure 2 presents the X-ray diffractograms of the BST-SS films described in Table 1. Diffraction peaks corresponding to (100), (111), (210), (211), (220), and (222) planes of the cubic BTO and STO phases were identified [21]. Diffraction peaks from (032) and (031) of the orthorhombic phase Ba1.91Sr0.09TiO4 (marked with ) [22], planes of (221) and (321), planes of the BTO tetragonal phase, and those of (211), (221), (310), (311), and (321) planes (marked with ) of the BST-SS tetragonal phase [8, 18, 23] were also identified. As STO sputtering power increases a shift of the BTO cubic phase reflections at higher diffraction angles is observed, indicating a lattice contraction due to the higher Sr composition in the films. Figure 3 depicts the values of the diffraction angles for the (110) plane of the measured films versus applied power, fitted employing the Boltzmann sigmoidal equation. A similar analysis was carried out for the shift observed for the positions associated with (111) and (211) planes.

The variables in (2) were assigned as follows: is representing the diffraction angle for the (110) plane; is corresponding to the diffraction angle when the STO magnetron power is 120 W, is , the diffraction angle with a BTO magnetron power = 120 W, is which represents the magnetron power values; is , the middle value indicated to be , and is which is the slope around ; thus (2) is rewritten as follows:The fitting values for , , , , and for (110), (111), and (211) planes of the cubic phase are presented in Table 2.

Composition values were obtained by considering that that peak shift follows Vegard’s law with the extreme values corresponding to those of SrTiO3 () and BaTiO3 ().

Figure 4 compares the XRD-derived stoichiometry evolution curve versus the magnetron power with the experimental composition measured by EDS, correspondingly fitted with Boltzmann’s equation (Table 3). The Ba content was converted to with the following expression:Although the curves show a very close behavior, there is a separation in the middle part of the curve, that is, from BTO power = 30 W and from STO power = 15 W and a remarkable discontinuity between BTO magnetron power from 45 to 60 W. The inset in Figure 4 depicts the difference between the composition obtained from XRD () and that obtained from EDX measurements (). To evaluate the discrepancy, in the inset the difference versus the power is plotted. The difference curve displays a bimodal increase when the power in each magnetron increases with a minimum value at the point corresponding to the point in the Boltzmann fit. This point 59 W in the BTO magnetron corresponds to an value around 0.5. Above the tetragonal phase of BaxSr1−xTiO3 is reported [9, 14, 23], in accordance with the X-ray diffractograms in Figure 2.

The discrepancies observed between the EDS and XRD curves can be due to Ba-rich and Sr-rich composition volume zones in the film attributed to the variations in flux composition and Ba/Sr sputtering yields. When the power in a target, that is, the BTO target, begins to increase, decreasing power in STO target, an instability in the composition of the sputtered flux is expected because of the atomic masses in the system (Ba = 137.33 u, Sr = 87.62 u, Ti = 47.867 u, O2 = 15.999 u, and Ar = 39.948 u), affecting not only the sputtering yield but also scattering between atoms of similar masses and the impinging ions. The Ba-rich and Sr-rich conditions are similar to those obtained in direct current arc discharge plasma process [5]. Composition also may be affected by changes in target erosion and total gas pressure during the sputtering process [13]. Additional discrepancies may be due to stress originated by the difference in size of Ba and Sr cations and to the different phases present in addition to the main one, interplanar spacing from XRD data is sensitive to these EDS measurements.

The values for the Boltzmann fit of the EDS curve are presented in Table 3 and compared with those from the XRD curves for (110), (111), and (211) diffraction planes. The values in (1), for the      parameters, were set at 0 and 1, respectively.

Figures 5(a) and 5(b) present the lattice parameter calculated from diffraction peaks associated with (110), (111), and (211) planes of the prepared films versus the parameters derived from XRD and EDS, respectively. Points were fitted with a straight line corresponding to Vegard’s Law:where is the molar fraction in BaxSr1−xTiO3 and “” is the lattice parameter calculated from the XRD peaks (110), (111), and (211), respectively. Table 4 presents the , , and values, which indicate a linear relation of the composition of the cubic BST-SS with Vegard’s law as observed in Figure 5(a).

The inset in Figure 5(a) presents the thickness of the studied samples as a function of the Ba composition obtained from XRD. The thickness reduction of the film B15S105 compared with that of pure SrTiO3 (B0S120) can be explained due to the mass difference, the sputtering yield, and scattering events. As the cosputtering process starts Ti ions would suffer additional scattering events due to the presence of Ba atoms affecting the film deposition rate. Also it is known that the sputtering yield of SrTiO3 is smaller than that of BaTiO3 [24]. It is worth emphasizing that, in this work, due to the employed deposition characteristics, the variation in the lattice parameter is due to the substitutional incorporation of Ba as its ionic radium is greater than that of Sr [1, 7]. The effect of differences in the thermal expansion coefficient (TEC) between substrate and BST-SS films is very important for epitaxial growth giving place to values of critical film thickness from 1 to ten nanometers [25]. Reported values of the TEC of BST-SS are close to those of SrTiO3 and BaTiO3, and , respectively [2527]. In our case there is no an epitaxial growth not just because of the characteristics of the sputtering deposition technique but also due to the amorphous characteristics of the substrate. The reported value of the thermal expansion coefficient for quartz is [28, 29]. Even when the difference between the TEC values of the films and the substrate is important, the close values of the TEC for the different compositions of the BST-SS allow affirming that if an effect due to the TEC difference would exist, it should affect in the same way to the whole set of samples, ruling out an effect on the observed trend of XRD peak positions.

Figure 5(b) shows the dependence of lattice parameter against composition values obtained from EDS. The greater scatter of data is expected due to the limitations of the EDS technique, mainly its lack of sensitivity to distinguish between the contribution from different crystallographic phases. Although Rietveld phase refinement is beyond the scope of the present contribution, it is a future task to elucidate the contribution of each phase in the cubic BST-SS.

The lattice parameter as a function of the RF-power is shown in Figure 6. The dependence follows a Boltzmann profile; a nonlinear behavior for a unique composition of Ba/Sr has been reported as a function of oxygen partial pressure and total gas pressure of the RF-magnetron sputtering process [11, 12].

The transmittance spectra of the deposited films are presented Figure 7(a) and analyzed according to Tauc model in Figure 7(b). In the transparent region (above the absorption edge) in the transmittance spectra, the number of interference maxima-minima reduces as the BTO sputtering power increases indicating a reduction in the film thickness, as will be discussed below. Correspondingly, the absorption edge shifts towards higher wavelength as is clearly observed in Figure 7(b).

By extrapolating the linear portion of the plot to , the band gap () values were calculated and presented in Figure 7(c) as a function of the composition parameter . For the sake of fitting, both axes are inverted. The sigmoidal dependence of the band gap with composition is evident. The tendency in the band gap values for BST-SS agrees with those reported in literature [2, 6, 30, 31].

Figure 8 presents the evolution of deposition rate, calculated from the film thickness and the sputter time, with respect to (a) the RF-magnetron power and (b) the parameter in the BaxSr1−xTiO3 solid solution. In Figure 8(a) a small decrement in the deposition rate is observed when the BTO target is turned on at the minimum power employed (105/15 power ratio STO/BTO). Thereafter the sputtering rate shows an exponential increase in the BTO magnetron power, fitted with the following expression:Correspondingly, when compared with the parameter, the deposition rate increases exponentially as follows: In Figure 8(a) the relationship between deposition rate Ba/Sr and 120 W in both targets is approximately 3/1.

The calculations of the area under the Boltzmann curve fitted to DRX and EDS in Figure 4 show a relationship between Ba/Sr roughly 3/2 and 1/1 respectively; these results indicate a better relationship between the DRX Boltzmann fit and deposition rate compared with the relationship between EDS Boltzmann fit and deposition rate due to reasons previously explained. Therefore, the growing of -parameter is directly related to the increment in the magnitude of power on BTO. These facts allow us to conclude that the calculations of -parameter based on DRX measurements are more accurate than those of EDS technique.

4. Conclusions

BaxSr1−xTiO3 (BST) films were grown on quartz substrates at 549°C by RF magnetron cosputtering from BTO and STO targets, using a total power of 120 W distributed between the two magnetrons. The films present a mixture of the cubic, orthorhombic, and tetragonal phases of BST. The Boltzmann fitting to the shift of diffraction peaks from the (110), (111), and (211) BST cubic phase planes as function of the applied RF-magnetron power allowed proposing a direct methodology to obtain the composition in BaxSr1−xTiO3 cubic solid solution. The “” parameter calculated with the Boltzmann expression was compared with the parameter obtained from the EDS measurements, and a close agreement was found. The lattice parameter variation versus power and the optical band gap versus the parameter were adjusted with similar Boltzmann expressions. The film deposition rate as a function of power and parameter follows an exponential growth behavior. The proposed method for estimation of the parameter from simple X-ray diffraction analysis has been found useful for presetting the initial sputtering conditions to adjust , film thicknesses, and optical band gap allowing a direct production of BST films with the adequate properties required for technological applications.

Disclosure

M. Meléndez-Lira is on sabbatical leave from Departamento de Fisica at Cinvestav-IPN.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

J. Reséndiz-Muñoz acknowledges financial support from the Ph.D. Program in Engineering and Materials Science at UASLP. F. Caballero-Briones contribution was supported by SIP-IPN Multidisciplinary Project 20161804-043, 20170215, and by EDI and SIBE-IPN grants. M. Meléndez-Lira acknowledges partial financial support from CONACyT during the sabbatical leave at IIT-UACJ, and thanks are due to L. López and M. Guerrero for their technical assistance at CINVESTAV-IPN.