# Receptor Function and Response of Semiconductor Gas Sensor

Authors: Noboru Yamazoe; Kengo Shimanoe
Journal: Journal of Sensors
 (2009)

---

## Abstract

Theoretical approaches to receptor function and response of semiconductor gas sensor are described, following the illustrations of some relevant key issues such as tunneling transport. Depletion in small semiconductor crystals is characterized by the occurrence of new type depletion (volume depletion) after conventional one (regional depletion), and inclusion of both types makes it possible to formulate the receptor function and response to oxygen (air base), oxidizing gas (nitrogen dioxide), and reducing gas (hydrogen). The equations derived theoretically using physical parameters of the semiconductor side and chemical parameters of the gases side appear to reproduce satisfactorily the sensing behavior to the aforementioned gases as well as the influence of changes in physical parameters such as grain size and donor density. Extension to the semiconductor crystals dispersed with surface electron-traps shows that the traps act as a sensitizer to promote sensor response.


---

## Body

## 1. Introduction

A semiconductor gas sensor (called device hereafter) possesses an electrical resistance made with a porous assembly of tiny crystals of an n-type metal oxide semiconductor, typically SnO2, In2O3, or WO3. The crystals are often loaded with a small amount of foreign substance (noble metals or their oxides) called a sensitizer. When operated at adequate temperature in air, the resistor changes its resistance sharply on contact with a small concentration of reducing gas or oxidizing gas, enabling us to know the concentration from the resistance change. For its inauguration with a report by Seiyama et al. [1] and a patent by Taguchi [2], this group of sensors has been subjected to a tremendous amount of R&D efforts world wide aiming at improvements of sensing performances and extensions to new applications. Thanks to these researches, the group not only has grown to provide important tools to detect and/or control gases in places in modern society but also has pioneered to founding a new technology field where the devices are called chemical sensors. Speaking more exactly, semiconductor gas sensors have been classified into two subtypes, that is, surface-sensitive type operating at temperatures below 500°C and bulk-sensitive one operating at high temperature (typically at 800°C) [3]. This article is concerned with those of the former type only.Apart from such remarkable achievements in practical applications, basic understandings of this group of sensors have hardly been satisfactory, despite tremendous efforts of so many researchers as summarized in reviews [5–8]. This is partly because there are many complex factors which affect sensing properties. Not only the selection of a proper oxide semiconductor is important but also the methods and conditions for fabricating sensor devices exert profound influences on gas sensing properties through changes in donor density, crystallite size, contacting geometry between crystals, packing density (or porosity), packing thickness, and so on [7]. In addition, the sensing properties are often modified largely with loading with foreign substances such as sensitizers. Understandings of these phenomena indeed have required interdisciplinary knowledge among semiconductor physics, surface chemistry, solid-state chemistry, and so on. In order to facilitate the understandings, we have proposed to assume that the sensing properties are determined by three main factors, that is, receptor function, transducer function, and utility factor, as schematically shown in Figure 1 [4]. The first factor is concerned with how each constituent crystal responds to the surrounding atmosphere containing oxygen and target gases (intraparticle issue). It is unanimous that oxygen is adsorbed on the crystals as negatively charged species, accompanied by the formation of a depletion layer inside the crystals. The target gases disturb the equilibrium through being adsorbed competitively or reacting with the adsorbed oxygen. The foreign substances like sensitizers dispersed on the crystals are assumed to affect these processes anyhow. The second factor is concerned with how the response of each particle is transformed into that of the whole device, and apparently this is related with the mechanism of electron transport between adjacent crystals (inter-particle issue). For a long time a double Schottky barrier model [9], which assumes migration transport of electrons over the barrier as shown, has been advocated for this process without critical check. The third one is concerned with the attenuation of the response due to the effect of diffusion and reaction of reactive target gases through the pores of the assembly of crystals (assembly issue) [10–12]. The above scheme has explained rather well qualitative nature of semiconductor gas sensors in several respects. However, it has failed to give quantitative understandings and, most importantly, to give new insights leading to innovations of this group of sensors. There should have been some serious defaults included in the scheme, particularly regarding the receptor and transducer functions.Figure 1
Three factors determining the response of semiconductor gas sensors [4].Fortunately, we encountered an interesting finding several years ago that thin film devices fabricated from hydrothermally prepared colloidal suspensions of SnO2 by a spin-coating technique showed temperature—almost independent resistances in air in the temperature range 150–400°C, as shown in Figure 2 [13]. Such thermal behavior of resistances is hardly consistent with the double Schottky barrier model mentioned above. Instead, tunneling transport of electrons across the contacts (or gaps) between adjacent crystals is strongly suggested. In addition, this transport mechanism has made much easier the theoretical modeling of receptor and transducer functions recently carried out [14–16], because the constituent crystals can now be treated independently from each other. As revealed during this process, depletion in small crystals easily goes beyond conventional one (regional depletion) to enter new type one (volume depletion).  Obviously, it is a lack of such information that has delayed fundamental understandings of this group of sensors, for most of their valuable gas sensing properties show up in the stage of volume depletion or nearby.Figure 2
Resistance of device under selected partial pressures of oxygen as measured at various temperatures (spin-coated SnO2 thin film) [13].In this review, we try to describe our theoretical approaches to semiconductor gas sensors. After brief descriptions of some experimental facts theoretically important, how to formulate receptor function and response of large and small constituent crystals to oxygen, nitrogen dioxide, and hydrogen is described to sufficient details followed by the comparisons with some experimental data. Also included is a recent extension [17] carried out to understand the sensitizing effects exerted by electron-accepting substances.


## 2. Key Issues about SnO2 Gas Sensors

### 2.1. Oxygen Adsorption on SnO2

The adsorption of oxygen on single-metal oxides was investigated by using a TPD technique by Iwamoto about 40 years ago [18]. He classified the transition of metal oxides into three groups depending on whether the amount of oxygen adsorbed and desorbed reversibly in the temperature range between room temperature and 500°C was large (Group A), small (B), or nondetectable (C). Group A consisted of p-type oxides, whereas some of n-type oxides such as SnO2, ZnO, and In2O3 belonged to B, and the others such as WO3 and V2O5 belonged to C. It is suggestive that Group B consists of typical oxides used for semiconductor gas sensors.Oxygen adsorption on SnO2 was investigated in more detail by our group, also by using the TPD technique but for an extended temperature range up to 600°C. The results are shown in Figure 3 [19]. There are four types of oxygen, α1,α2,β, and γ, recognized to be desorbed from SnO2 altogether. Of these, only the last two types, β (presumably O-) and γ (lattice oxygen O2-), are desorbed in the temperature range above 400°C after oxygen has been adsorbed during cooling from 600°C or at a fixed temperature of 400°C. On the other hand, after oxygen has been adsorbed in the lower temperature range (below 155°C), only α1 (neutral molecular O2) and α2(superoxide ion O2-) are desorbed in the temperature range up to 250°C. Clearly these molecular adsorbates are formed only when the other dissociated species are absent. Also it is clear that at usual sensor operating temperatures (typically 300–400°C) β species prevails on the SnO2 surface. The amount of adsorption of this species has been shown to be fairly small (less than 1% in surface coverage), suggesting that its adsorption is limited by the supply of electrons. Later, a report was published, which identified β to O2- [20]. This identification cannot fit well to the theoretical analysis of response to oxygen, however. So it is assumed to be O- in the following treatments.Figure 3
TPD chromatograms of oxygen desorbed from SnO2 [19].


### 2.2. Grain Size Effects

About twenty years ago, we explored experimental methods to prepare small grains (crystallites) of SnO2. During this study, we found remarkable grain size effects [21]. As shown in Figure 4, the resistances under air (base) as well as under exposure to H2 or CO in air, Ra and Rg, respectively, were found to increase sharply as the grain size (dm) decreased beyond a critical value (dc), while the sensor response to H2 or CO, Ra/Rg, also increased, though more gradually, as dm decreased. Intuitively we felt that the phenomena were associated with the completion of depletion up to the whole region of constituent grains. Our interpretation (or assumption) was that dc would be coincident with twice the thickness of depletion depth and that decreasing dm beyond this point would be responsible for the sharp increases of Ra, Rg, and Ra/Rg for some reasons yet unknown. Later, similar size effects were also observed for the response of lamellar or granular WO3-based devices to NO2, as shown in Figure 5 [22].Grain size effects on resistance in air and sensor response to H2 or CO in air observed with SnO2 (brush-coated thick films) [21].
(a)(b)Sensor response to 1 ppm NO2in air as correlated with granular size or lamellar thickness of constituent WO3 crystals [22]. The symbols attached indicate the preparation methods used for respective WO3 samples. See literature for details.
(a)(b)Although such grain size effects seemingly appeared to match well to the double Schottky barrier model, their real meanings were left open for a long time. Some researchers tried to explain the grain size effects [23–26]. Among them, Rothschild and Komen pointed out that the size effects could be explained as reflecting changes in the surface area/volume ratio of the constituent SnO2 grains [23, 24]. Such situation has turned out to appear typically in the stage of volume depletion in the present theory, as mentioned later.


### 2.3. Feasibility of Tunneling Transport

In sensor devices, constituent crystals are connected to adjacent ones either by contacts or by necks, the proportion of which depends on the methods and conditions of device fabrication. In case the device is fabricated through wet processes, the contact connections usually prevail, as seen in Figure6 [13]. In such cases, it is easily understood that the device resistance is determined as a sum of the resistance of each contact. This idea has already been conceived in the double Schottky barrier model, which assumes that electrons migrate over the barrier at the contact. However, this model cannot be consistent with the temperature-almost-independent behavior shown in Figure 2, as already mentioned. From this consideration, we explored the possibility of tunneling transport of electrons [13]. The tunneling model used is shown in Figure 7. The height of the wall to tunnel through is set to be equal to the electron affinity of the crystals, that is, the potential energy difference (V) between the vacuum level and the conduction band edge at the surface, while the thickness is set to the gap (L) between adjacent crystals. If V is known, tunneling probability (P) of an electron with energy E above conduction band edge can be estimated as a function of L by using a well-established equation. Reportedly, work function of SnO2 is dispersed between 4.7 and 5.7 eV [27, 28], while donor levels are shallow from conduction band edge. Therefore, V takes a value of about 5.7 eV at the maximum. The tunneling probability thus estimated assuming V=6 eV is shown in Figure 8. When the gap is zero (contact), P is unity naturally regardless of E. P decreases sharply with increasing L and decreasing E, but notably it still keeps a significant value (.01 or above) at L=0.1 nm and E=0.03 eV. Considering that thermal energy is 51.7 meV at 600 K, it is understood that electrons can tunnel through such a gap with a significant probability. What this means is that electrons, exited thermally, can be transported by tunneling not only through the directly contacting region (L=0) but also through a small gap in its vicinity, as schematically shown in Figure 9. It follows that even in the case of point contact (contacting area zero) electrons can be transported by tunneling but this is not the case by migration. As shown in Figure 7, V is maintained at the same value in either of flat-band state (a) and bent-band state (b). Therefore, the tunneling probability is kept the same in both, only the number of electrons involved being decreased with a progress of band bending. In this way, the tunneling mechanism is confirmed to be feasible well to the electron transport between adjacent crystals. There is no reason to exclude it. Rather it can even eliminate difficulties encountered by the migration mechanism.Figure 6
State of constituent grains as observed on SEM (spin-coated SnO2 film) [13].Model for electron tunneling across a gap between SnO2 grains [13]. (a) Flat band condition, (b) Bent band condition. V: Wall height V.L.: Vacuum level, C.B.: Conduction band, F.B.: Forbidden band, V.B.: Valence band.
(a)(b)Figure 8
Tunneling probability calculated for various values ofE as a function of gap width (L) [13].Figure 9
Schematic drawing of tunneling conduction across a gap between contacting spheres.


### 2.4. Models of Electron Transport between Grains

The mechanism of electron transport between neighboring crystals is directly related with the transducer function of sensors. There are probably three representative models, as schematically drawn in Figure10. Spherical crystals (uniform in size), connected with neighbors through a contact or a neck, are assumed to be depleted of electrons in the outer region only. Double Schottky barrier model (a) is a traditional one. It assumes that electrons are transported from nondepleted (core) region of one sphere to that of another by migration beyond a barrier in between. The barrier height not only determines the conductance of the contacting part but also gives rise to the activation energy of conduction on changing temperature. Tunneling transport model (b) assumes that electrons located at the periphery of one sphere are transported to that of another by tunneling through a small gap (typically 0–0.1 nm) in between. The conductance is proportional to the density of electrons at the periphery, which is determined by the surface barrier height. On changing temperature, no activation energy is associated with the conductance provided that the barrier height is kept the same. Neck (or conduction channel) model (c) assumes that electrons migrate through a conduction channel which is formed by connecting the core regions of neighboring spheres. The channel width is narrowed at the neck parts so that the conductance is determined by the geometric relation between neck size and depletion depth. Of these models, (a) and (c), though looking likely at a glance, are not always free from various difficulties. Probably the most serious one appears after the depletion depth has reached the radius of spheres or necks. Such situation can happen easily when the size of spheres is reduced or gas adsorption is strengthened. Both models hardly seem to be applicable to such a situation. In our opinion, (b) seems to survive in various situations. Actually theoretical treatises in this article are based on this model.Models for electron transport between adjacent crystals.
(a)(b)(c)It is remarked that important information concerning the transducer function is available from the Hall effects measurement. As generally accepted, conductance,σ, for an n-type semiconductor is given by the density, [e], and mobility, μ, of charge carriers (electrons) as follows:
(1)σ=-q[e]μ.
Here q is the electric charge of proton. The values of [e] and μ can be estimated independently from the Hall effects measurement. The measurements on semiconductor gas sensor devices have been carried out by Japanese experts [29, 30]. According to their results, σ and [e] changed in parallel to each other with increasing concentration of reducing gases while μ was kept almost constant, as measured for sputtered SnO2 thick film devices. This means that the change of [e] is responsible exclusively for the change of σ, thus supporting the tunneling transport model (b) in Figure 10. An exceptional case has been reported with an SnO2 device fabricated from tin metal by an evaporation and oxidation method [29]. In this case, both [e] and μ increased with increasing concentration of target gas (ethanol), resulting in enhanced sensitivity of σ as compared with the sputtered devices. This difference in the behavior of μ probably arose from a change in the morphology of the constituent crystals. The authors explained such gas sensitivity of μ as resulting from very small crystals on the basis of a neck model like one shown as (c) in Figure 10, but it seems difficult to explain the gas sensitivity of [e] based on this model. The tunneling transport model, on the other hand, seems to account for the sensitivity of both [e] and μ better provided that dispersion in crystal size is admitted, though this is yet to be done.


### 2.5. Sensitization by Foreign Substances

As mentioned in Section1, actual sensor devices are almost always loaded with small amounts of foreign substances so called sensitizers. Kinds of sensitizers used or confirmed are quite various, including noble metals (Pt and Au), noble metal oxides (PdO, Ag2O, Rh2O3, etc.), base metal oxides (Fe2O3, Co3O4, etc.), alkaline or rare earth metal oxides (CaO, La2O3, etc.), and alkali metal oxide (Rb2O). To make sensitization in effect, sensitizers should be dispersed well on the constituent crystals of the sensor devices. As shown in Figure 11 [31], sensitizers are usually deposited as small particles adhering well to the underlying crystals. However, the roles of sensitizers (sensitization mechanisms) have not always been made clear. Beside the sensitization through modifying the acid-base property of the crystal surface which is vitally important in acid-base catalyzed reactions like the oxidation of alcohols, we have proposed two types of sensitization mechanisms, that is, chemical sensitization and electronic sensitization, operative for reducing gases [32]. The chemical sensitization is brought about from the catalytic actions of sensitizers like Pt. When target gas is resistant to oxidation like methane, for example, it is activated (or transformed) into more reactive intermediates like aldehydes before being exposed to the crystal surface. The electronic sensitization is in effect when the foreign substances draw electrons from the crystals, as easily witnessed from an increase in device resistance. Many of the sensitizers known are p-type oxides so that formation of p-n junctions between the tiny deposits of sensitizers and underlying SnO2 (n-type) does not seem to be so unlikely. It has been shown that work function of SnO2 is increased significantly by the loading with PdO (about 0.1 eV), Ag2O (0.3 eV) and CuO (0.5 eV), as measured on XPS [33]. Here the figures in the parentheses are increments of work function in respective cases. In these cases, each sensitizer can undergo a redox change and, when it is reduced with the target gases, work function and device resistance have been confirmed to go back to their original values. It has been thus inferred that coupling of the formation of p-n junctions and the redox changes results in the sensitization. This mechanism, however, encounters difficulties in accounting for the sensitization exerted by the other sensitizers like Au, Fe2O3, Co3O4, and Rb2O, which increase work function of SnO2 more or less but are less facile or totally stable to redox changes. It seems that there are different mechanisms of electronic sensitization which cover the redox changes-free cases. Theoretical approaches to this issue are worth being carried out.Figure 11
State of deposits of PdO on SnO2 grain as observed on TEM [31].


## 2.1. Oxygen Adsorption on SnO2

The adsorption of oxygen on single-metal oxides was investigated by using a TPD technique by Iwamoto about 40 years ago [18]. He classified the transition of metal oxides into three groups depending on whether the amount of oxygen adsorbed and desorbed reversibly in the temperature range between room temperature and 500°C was large (Group A), small (B), or nondetectable (C). Group A consisted of p-type oxides, whereas some of n-type oxides such as SnO2, ZnO, and In2O3 belonged to B, and the others such as WO3 and V2O5 belonged to C. It is suggestive that Group B consists of typical oxides used for semiconductor gas sensors.Oxygen adsorption on SnO2 was investigated in more detail by our group, also by using the TPD technique but for an extended temperature range up to 600°C. The results are shown in Figure 3 [19]. There are four types of oxygen, α1,α2,β, and γ, recognized to be desorbed from SnO2 altogether. Of these, only the last two types, β (presumably O-) and γ (lattice oxygen O2-), are desorbed in the temperature range above 400°C after oxygen has been adsorbed during cooling from 600°C or at a fixed temperature of 400°C. On the other hand, after oxygen has been adsorbed in the lower temperature range (below 155°C), only α1 (neutral molecular O2) and α2(superoxide ion O2-) are desorbed in the temperature range up to 250°C. Clearly these molecular adsorbates are formed only when the other dissociated species are absent. Also it is clear that at usual sensor operating temperatures (typically 300–400°C) β species prevails on the SnO2 surface. The amount of adsorption of this species has been shown to be fairly small (less than 1% in surface coverage), suggesting that its adsorption is limited by the supply of electrons. Later, a report was published, which identified β to O2- [20]. This identification cannot fit well to the theoretical analysis of response to oxygen, however. So it is assumed to be O- in the following treatments.Figure 3
TPD chromatograms of oxygen desorbed from SnO2 [19].


## 2.2. Grain Size Effects

About twenty years ago, we explored experimental methods to prepare small grains (crystallites) of SnO2. During this study, we found remarkable grain size effects [21]. As shown in Figure 4, the resistances under air (base) as well as under exposure to H2 or CO in air, Ra and Rg, respectively, were found to increase sharply as the grain size (dm) decreased beyond a critical value (dc), while the sensor response to H2 or CO, Ra/Rg, also increased, though more gradually, as dm decreased. Intuitively we felt that the phenomena were associated with the completion of depletion up to the whole region of constituent grains. Our interpretation (or assumption) was that dc would be coincident with twice the thickness of depletion depth and that decreasing dm beyond this point would be responsible for the sharp increases of Ra, Rg, and Ra/Rg for some reasons yet unknown. Later, similar size effects were also observed for the response of lamellar or granular WO3-based devices to NO2, as shown in Figure 5 [22].Grain size effects on resistance in air and sensor response to H2 or CO in air observed with SnO2 (brush-coated thick films) [21].
(a)(b)Sensor response to 1 ppm NO2in air as correlated with granular size or lamellar thickness of constituent WO3 crystals [22]. The symbols attached indicate the preparation methods used for respective WO3 samples. See literature for details.
(a)(b)Although such grain size effects seemingly appeared to match well to the double Schottky barrier model, their real meanings were left open for a long time. Some researchers tried to explain the grain size effects [23–26]. Among them, Rothschild and Komen pointed out that the size effects could be explained as reflecting changes in the surface area/volume ratio of the constituent SnO2 grains [23, 24]. Such situation has turned out to appear typically in the stage of volume depletion in the present theory, as mentioned later.


## 2.3. Feasibility of Tunneling Transport

In sensor devices, constituent crystals are connected to adjacent ones either by contacts or by necks, the proportion of which depends on the methods and conditions of device fabrication. In case the device is fabricated through wet processes, the contact connections usually prevail, as seen in Figure6 [13]. In such cases, it is easily understood that the device resistance is determined as a sum of the resistance of each contact. This idea has already been conceived in the double Schottky barrier model, which assumes that electrons migrate over the barrier at the contact. However, this model cannot be consistent with the temperature-almost-independent behavior shown in Figure 2, as already mentioned. From this consideration, we explored the possibility of tunneling transport of electrons [13]. The tunneling model used is shown in Figure 7. The height of the wall to tunnel through is set to be equal to the electron affinity of the crystals, that is, the potential energy difference (V) between the vacuum level and the conduction band edge at the surface, while the thickness is set to the gap (L) between adjacent crystals. If V is known, tunneling probability (P) of an electron with energy E above conduction band edge can be estimated as a function of L by using a well-established equation. Reportedly, work function of SnO2 is dispersed between 4.7 and 5.7 eV [27, 28], while donor levels are shallow from conduction band edge. Therefore, V takes a value of about 5.7 eV at the maximum. The tunneling probability thus estimated assuming V=6 eV is shown in Figure 8. When the gap is zero (contact), P is unity naturally regardless of E. P decreases sharply with increasing L and decreasing E, but notably it still keeps a significant value (.01 or above) at L=0.1 nm and E=0.03 eV. Considering that thermal energy is 51.7 meV at 600 K, it is understood that electrons can tunnel through such a gap with a significant probability. What this means is that electrons, exited thermally, can be transported by tunneling not only through the directly contacting region (L=0) but also through a small gap in its vicinity, as schematically shown in Figure 9. It follows that even in the case of point contact (contacting area zero) electrons can be transported by tunneling but this is not the case by migration. As shown in Figure 7, V is maintained at the same value in either of flat-band state (a) and bent-band state (b). Therefore, the tunneling probability is kept the same in both, only the number of electrons involved being decreased with a progress of band bending. In this way, the tunneling mechanism is confirmed to be feasible well to the electron transport between adjacent crystals. There is no reason to exclude it. Rather it can even eliminate difficulties encountered by the migration mechanism.Figure 6
State of constituent grains as observed on SEM (spin-coated SnO2 film) [13].Model for electron tunneling across a gap between SnO2 grains [13]. (a) Flat band condition, (b) Bent band condition. V: Wall height V.L.: Vacuum level, C.B.: Conduction band, F.B.: Forbidden band, V.B.: Valence band.
(a)(b)Figure 8
Tunneling probability calculated for various values ofE as a function of gap width (L) [13].Figure 9
Schematic drawing of tunneling conduction across a gap between contacting spheres.


## 2.4. Models of Electron Transport between Grains

The mechanism of electron transport between neighboring crystals is directly related with the transducer function of sensors. There are probably three representative models, as schematically drawn in Figure10. Spherical crystals (uniform in size), connected with neighbors through a contact or a neck, are assumed to be depleted of electrons in the outer region only. Double Schottky barrier model (a) is a traditional one. It assumes that electrons are transported from nondepleted (core) region of one sphere to that of another by migration beyond a barrier in between. The barrier height not only determines the conductance of the contacting part but also gives rise to the activation energy of conduction on changing temperature. Tunneling transport model (b) assumes that electrons located at the periphery of one sphere are transported to that of another by tunneling through a small gap (typically 0–0.1 nm) in between. The conductance is proportional to the density of electrons at the periphery, which is determined by the surface barrier height. On changing temperature, no activation energy is associated with the conductance provided that the barrier height is kept the same. Neck (or conduction channel) model (c) assumes that electrons migrate through a conduction channel which is formed by connecting the core regions of neighboring spheres. The channel width is narrowed at the neck parts so that the conductance is determined by the geometric relation between neck size and depletion depth. Of these models, (a) and (c), though looking likely at a glance, are not always free from various difficulties. Probably the most serious one appears after the depletion depth has reached the radius of spheres or necks. Such situation can happen easily when the size of spheres is reduced or gas adsorption is strengthened. Both models hardly seem to be applicable to such a situation. In our opinion, (b) seems to survive in various situations. Actually theoretical treatises in this article are based on this model.Models for electron transport between adjacent crystals.
(a)(b)(c)It is remarked that important information concerning the transducer function is available from the Hall effects measurement. As generally accepted, conductance,σ, for an n-type semiconductor is given by the density, [e], and mobility, μ, of charge carriers (electrons) as follows:
(1)σ=-q[e]μ.
Here q is the electric charge of proton. The values of [e] and μ can be estimated independently from the Hall effects measurement. The measurements on semiconductor gas sensor devices have been carried out by Japanese experts [29, 30]. According to their results, σ and [e] changed in parallel to each other with increasing concentration of reducing gases while μ was kept almost constant, as measured for sputtered SnO2 thick film devices. This means that the change of [e] is responsible exclusively for the change of σ, thus supporting the tunneling transport model (b) in Figure 10. An exceptional case has been reported with an SnO2 device fabricated from tin metal by an evaporation and oxidation method [29]. In this case, both [e] and μ increased with increasing concentration of target gas (ethanol), resulting in enhanced sensitivity of σ as compared with the sputtered devices. This difference in the behavior of μ probably arose from a change in the morphology of the constituent crystals. The authors explained such gas sensitivity of μ as resulting from very small crystals on the basis of a neck model like one shown as (c) in Figure 10, but it seems difficult to explain the gas sensitivity of [e] based on this model. The tunneling transport model, on the other hand, seems to account for the sensitivity of both [e] and μ better provided that dispersion in crystal size is admitted, though this is yet to be done.


## 2.5. Sensitization by Foreign Substances

As mentioned in Section1, actual sensor devices are almost always loaded with small amounts of foreign substances so called sensitizers. Kinds of sensitizers used or confirmed are quite various, including noble metals (Pt and Au), noble metal oxides (PdO, Ag2O, Rh2O3, etc.), base metal oxides (Fe2O3, Co3O4, etc.), alkaline or rare earth metal oxides (CaO, La2O3, etc.), and alkali metal oxide (Rb2O). To make sensitization in effect, sensitizers should be dispersed well on the constituent crystals of the sensor devices. As shown in Figure 11 [31], sensitizers are usually deposited as small particles adhering well to the underlying crystals. However, the roles of sensitizers (sensitization mechanisms) have not always been made clear. Beside the sensitization through modifying the acid-base property of the crystal surface which is vitally important in acid-base catalyzed reactions like the oxidation of alcohols, we have proposed two types of sensitization mechanisms, that is, chemical sensitization and electronic sensitization, operative for reducing gases [32]. The chemical sensitization is brought about from the catalytic actions of sensitizers like Pt. When target gas is resistant to oxidation like methane, for example, it is activated (or transformed) into more reactive intermediates like aldehydes before being exposed to the crystal surface. The electronic sensitization is in effect when the foreign substances draw electrons from the crystals, as easily witnessed from an increase in device resistance. Many of the sensitizers known are p-type oxides so that formation of p-n junctions between the tiny deposits of sensitizers and underlying SnO2 (n-type) does not seem to be so unlikely. It has been shown that work function of SnO2 is increased significantly by the loading with PdO (about 0.1 eV), Ag2O (0.3 eV) and CuO (0.5 eV), as measured on XPS [33]. Here the figures in the parentheses are increments of work function in respective cases. In these cases, each sensitizer can undergo a redox change and, when it is reduced with the target gases, work function and device resistance have been confirmed to go back to their original values. It has been thus inferred that coupling of the formation of p-n junctions and the redox changes results in the sensitization. This mechanism, however, encounters difficulties in accounting for the sensitization exerted by the other sensitizers like Au, Fe2O3, Co3O4, and Rb2O, which increase work function of SnO2 more or less but are less facile or totally stable to redox changes. It seems that there are different mechanisms of electronic sensitization which cover the redox changes-free cases. Theoretical approaches to this issue are worth being carried out.Figure 11
State of deposits of PdO on SnO2 grain as observed on TEM [31].


## 3. Theoretical Approach to Sensor Response to Oxygen

For a long time, the mechanism of gas sensing has been one of the most interested subjects in semiconductor gas sensors. Despite many efforts so far made, however, few have succeeded in revealing it from a fundamental point of view. As known well, chemisorption and catalysis of gas molecules on semiconductor surfaces were paid much attention in the age earlier than the inauguration of semiconductor gas sensors. An electronic theory was proposed to explain chemisorption and catalysis on semiconductors [34–36]. It explains well how the electron transfer between adsobates and semiconductor results in the formation of a depletion, accumulation, or inversion layer in the subsurface of semiconductors. When applied for semiconductor gas sensors, it has been useful in understanding their behavior qualitatively. Unfortunately, however, no quantitative understandings have been available with it. There are several reasons for this, and the most serious one is, in our opinion, that the theory cannot deal correctly with so small semiconductor crystals as used in gas sensors. As easily understood, response to oxygen is a starting base of semiconductor gas sensors. The sensors are active in air, being alert to the access of target gases. In the beginning of this section, this theory is applied to the response to oxygen to demonstrate what comes out. Then a new theory is introduced to deal with small crystals, which is indispensable for understanding the response behavior of semiconductor gas sensors. In the following treatments, the semiconductor crystals taken into consideration are shaped in plate unless otherwise noted, because mathematical expressions for other shapes, sphere and column, become far more complex and tedious due to shape effect.


### 3.1. Conventional Theory of Depletion (Large Crystals)

Depletion of large semiconductor crystals can be treated with the conventional theory. Let us consider the adsorption of electron-accepting molecules (typically O2) on the surface of a large n-type semiconductor crystal. This situation is well approximated by the adsorption on a plane of a large crystal cleaved in halves, as schematically drawn in Figure 12 [37]. As a result of electron transfer, a depletion layer is formed in the subsurface region, the layer thickness increasing with an increase in electron transfer. The electron transfer reaches equilibrium when Fermi level at the adsorbed species coincides with that in the semiconductor bulk.Figure 12
Schematic diagram of electrostatic equilibrium for large semiconductor crystals under exposure to oxygen. Ec: conduction band edge, Ec: valence band edge, EF: Fermi level, Qsc: surface charge density, w: depletion depth, qV: potential energy of electrons, qVs: surface potential energy barrier [37].The electrostatic equilibrium inside the crystal in such a case can be treated quantitatively under simplifying assumptions. Conventionally it is assumed that (1) donors are ionized completely, (2) there are no other electron-traps other than the adsorbates, and (3) all of conduction electrons up to a depth w are completely transferred to the surface (abrupt model). Since donors in SnO2 are reportedly present at shallow levels (0.05–0.15 eV) below the conduction band edge [38], the first assumption can be acceptable fairly well. The second one may be fit to an ionic semiconductor like SnO2 better than to a nonionic semiconductor like silicon. The third one, though useful in simplifying mathematical treatments, cannot always be rationalized as mentioned later. However, abrupt model is used here always unless otherwise noted.Under these assumptions, the density of surface charges,QSC, is nominally equal to -qNdw, where q is electric charge of proton and Nd the density of donors of semiconductor:
(2)QSC=-qNdw=-qNdLDm.
Here m is reduced depletion depth defined by m=w/LD. Under these conditions, electric potential, V, in the depletion region should satisfy the following Poisson’s equation, where x is the depth from the surface and ε is the permittivity of semiconductor times the free space permittivity:
(3)d2Vdx2=-qNdε.
By introducing the boundary conditions that dV/dx as well as V are zero at x=w, (3) is solved to give the following depth profile of electric potential:
(4)V(x)=-(qNd2ε)(x-w)2.
Potential energy of electrons, qV(x), as well as the height of surface potential barrier, qVS, are obtained from (4):
(5)qV(x)kT=(q2Nd2εkT)(x-w)2=(12){(x-w)LD}2,qV(X)kT=(12)(X-m)2,(6)qVSkT=(12)(wLD)2=(12)m2.
Here LD is the Debye length of semiconductor given by LD=(q2Nd/kT)-1/2, and X is reduced distance defined by X=x/LD. The density of electrons at the surface, [e]S, is given from (6) as follows:
(7)[e]S=Ndexp(-qVSkT)=Ndexp(-m22).
In this way, the surface density of adsorbates, -QSC/q, and electrons can be correlated uniquely through m once LD and Nd are given.


### 3.2. Response of Large Crystals

It is assumed that oxygen is adsorbed asO- exclusively on the semiconductor surface. Then the overall equation of oxygen adsorption and its equilibrium are expressed as follows:
R1O2+2e=2O-,(8)(KO2PO2)1/2[e]S=[O-].KO2 and PO2 are the adsorption constant and the partial pressure of oxygen, respectively, and [e]S and [O-] are the surface densities of free electrons and O-, respectively. Equation (8) indicates that (KO2PO2)1/2 determines the density ratio, [O-]/[e]S, not the respective quantities. These quantities are determined only after the chemical equilibrium is coupled with the electrostatic equilibrium just mentioned.WhenO- is solely responsible for the total charge density of the semiconductor as assumed,
(9)[O-]=-QSCq=NdLDm.
By inserting (2) and (8) into (7), we obtain
(10)(KO2PO2)1/2Ndexp(-m22)=NdLDm,X=mexp(m22).
Here X is reduced adsorption strength of oxygen defined by X=(KO2PO2)1/2/LD. Indication of (10) is that if X is given, m is determined uniquely. This in turn determines [e]S through (8). Thus we can correlate oxygen adsorption with the surface density of electrons.In order to obtain the correlation with sensor resistance (R), we need to know the relation between R and [e]S. The relation is simple if the tunneling transport model (Figure 10(b)) is assumed, as follows:
(11)RR0=[e]S0[e]S=exp(m22).
Here R0 and [e]S0 are the values of R and [e]S at the flat band state, respectively. The latter quantity coincides with Nd for the simplified mode used here. R and R0 are influenced by many factors other than oxygen adsorption, such as packing of crystals and geometry of sensor device, but those influences can be canceled out effectively by normalizing R by R0. The correlation between R/R0 and X thus obtained through (10) and (11) is shown on logarithmic scales in Figure 13 [14]. R/R0 increases with increasing X and its slope becomes equal to unity in the middle region of X, explaining that the power law-type dependence of R on PO2 with a power index of 1/2 is actually observed. However, no important information else can be obtained from such a treatise. Similar conclusions have been obtained for the responses to NO2 and H2 [14]. This is because we have used conventional theory of depletion developed for large crystals, thus missing important information about the size and shape of crystals.Figure 13
Correlation betweenR/R0 and X on logarithmic scales [14]. X=(KO2PO2)1/2/LD


### 3.3. Two Types of Depletion in Small Crystals

Let us consider how depletion develops in a plate crystal of thickness2a (a, half thickness), starting from the flat-band state, with increasing PO2 in the ambient. The electrostatic equilibrium of this system is schematically drawn in Figure 14(a). When PO2 is small, depletion takes place in the same way as previously seen in a large crystal, depletion depth increasing with increasing PO2. However, this way of depletion is finished at PO2(II) where depletion has extended to cover the whole region of the crystal. Obviously the value of PO2(II) decreases with decreasing half thickness (a). What would happen when PO2 is further increased to PO2(III)) in the figure? This is a main concern here. The above limit of depletion comes solely from the physics of the crystal. The chemistry of oxygen adsorption, on the other hand, does not care the physics. It requires (8) to hold at any stages, or, speaking more specifically, the ratio [O-]/[e]S to increase linearly with increasing X=(KO2PO2)1/2. In order for the physics side (electrostatic equilibrium) to fulfill this requirement, the only available way is to allow Fermi level to shift down by a designated amount, pkT, whereas the band bending profile is kept unchanged. With this shift, not only the whole region from the surface adsorbates (O-) to the bulk of the crystal can be placed at the same Fermi level but also the potential energy of electrons goes up to the level at which [e]S takes a value consistent with the chemical equilibrium condition (8). This way of establishing the electrostatic equilibrium is thus very different from one observed at lower PO2. It is noted that electron transfer from the bulk to the surface continues still in this stage, though it is only nominal when the abrupt model is used (Figure 14(b)). The electron transfer is made more apparent in case some electrons are left behind in the depletion layer (modified model) as illustrated in the same figure. In this case, the density of electrons inside the crystal decreases progressively as the Fermi level shift increases, showing that electrons are squeezed out from the whole crystal to be supplied to the surface. It is therefore convenient to distinguish the depletion in this stage from the conventional one in the initial stage up to PO2(II). We call them here volume depletion and regional depletion, respectively, while the border in between is called boundary depletion.Schematic diagrams of electronic equilibrium for semiconductor plate (half thicknessa) placed under various partial pressures of oxygen (PO2). (a) Depth profiles of potential energy. (b) Depth profiles of conduction electrons density assumed on abrupt model (broken lines) and modified model (full lines). See section 3.6 for modified model [15].
(a)(b)


### 3.4. Response at Regional Depletion (Plate)

As mentioned earlier, the occurrence of regional depletion followed by volume depletion is characteristic to small crystals. Receptor function and response of small plate crystals to oxygen in the stage of regional depletion are discussed here. The distribution of potential inside the plate is obtained by solving Poisson’s equation, which is the same as (3) but has an origin located at a center of the plate with coordinate x being taken as a distance normal to the surface from origin (see Figure 14). In the stage of regional depletion, potential energy of electrons is expressed as shown below:
(12)qV(x)=(q2Nd2ε){x-(a-w)}2for|x|≧a-w,=0for|x|<a-w.
These are rewritten in reduced form as
(13)qV(X)kT=(12){X-(n-m)}2for|x|≧n-m,=0for|x|<n-m.
Reduced quantities, X and m, are the same as defined before, and n is reduced size defined by n=a/LD. These expressions are essentially the same as those derived for large crystals, though depletion in this case starts from both surfaces. Surface potential energy (barrier) is expressed by the same equation (7), while (9) holds also as assumed. As a result, either of (10) and (11) derived for large crystals can be used to correlate between R/R0 and X also in this case. It is noted, however, that this can be valid only for plate crystals, which are free from the shape effect.


### 3.5. Response at Volume Depletion (Plate)

Potential energy in the stage of boundary depletion is obtained by insertingw=a in (12) or m=n in (13) and that in volume depletion is larger by pkT or p, as follows:
(14)qV(x)=(q2Nd2ε)x2+pkT.(15)qV(X)kT=(12)X2+p.
Here p=0 indicates boundary depletion. The surface potential energy is obtained by inserting X=n into (15):
(16)qVskT=(12)n2+p.
We can then express [e]S and R/R0 by using n and p. [O-] is obtained from (9) if m is replaced by n. By inserting these quantities into the mass action law (8), we can correlate X with n and p. Thus we obtain a couple of equations:
(17)X=nexp{(12)n2+p},RR0=exp{(12)n2+p}.
These equations determine the correlation between R/R0 and X. Unlike the case of regional depletion, however, R/R0 is an explicit function of X, because we obtain from the couple of equations
(18)RR0=Xn=(KO2PO2)1/2a.
Normalized resistance is thus linear to (KO2PO2)1/2 and the reciprocal of a.The correlations betweenR/R0 and X in the stages of regional depletion and volume depletion are drawn for smaller selected values of n in Figure 15. In principle, regional depletion always precedes another on increasing X, but actually when n is small it is not visible well in the drawing on this scale, allowing volume depletion to prevail. It is also clear that the slope of each correlation (sensitivity to oxygen) increases sharply with decreasing n and is inversely proportional to n or a in the stage of volume depletion. This indicates a size effect on the sensor response to oxygen. It is remarked that the derivative of R or normalized R with respect to the concentration (or partial pressure) of target gas in problem is usually defined as sensor sensitivity. In the case of oxygen adsorption, however, the sensitivity is better defined as the derivative of R/R0 with respect to PO21/2, which is equal to KO21/2/a in the stage of volume depletion.Figure 15
Correlations betweenR/R0 and X for plate crystals with reduced thickness n [17].


### 3.6. Effects of Shape and Size

The same treatment can be extended to spherical crystals (radiusa) and columnar crystals (columnar radius a). For this purpose, two things should be taken into account. First, Poisson’s equation is transformed, for the convenience of solution, into one fitting to the symmetry of the crystal form. A polar spherical coordinates system and a columnar coordinates system are illustrated together with one-dimensional system so far used in Figure 16. In the respective systems, Poisson’s equation is expressed as a function of a single variable as follows:
(19)d2Vdx2=-qNdε,(Plate),(1r2)ddr(r2dVdr)=-qNdε,(Sphere),(1r)ddr(rdVdr)=-qNdε,(Column).
Here r is radial displacement from the center of crystals. Each equation can be solved by using the same boundary conditions as used before to give a solution shown below:
(20)V(x)=-(qNd2ε){x-(a-w)}2,(Plate),V(r)=-(qNd/2ε){r-(a-w)}2{1+2(a-w)/r}3,(Sphere),V(r)=-(qNd4ε){r2-(a-w)2[1+2ln(r(a–w))]},(Column).
Second, the relations between total surface charge density (QSC) and depletion depth (w or m) should be altered to meet crystal symmetry:
(21)QSCq=-Ndw=-NdLDm,(Plate),QSCq=-Nd{∫a-wa4πr2dr}(4πa2)=-(13)NdLDn{1-[(n-m)n]3},(Sphere),QSCq=-Nd{∫a-wa2πrdr}(2πa)=-(12)NdLDn{1-[(n-m)n]2},(Column).
Here m and n are reduced depletion depth (w/LD) and reduced size (a/LD), respectively.Figure 16
Coordinates systems adapted for plate, sphere, and column [15].Based on the abrupt model, receptor function and sensor response to oxygen can be formulated for spherical and columnar crystals in the same way as done for plate crystals. The correlations betweenR/R0 and X thus derived are summarized as follows.Regional depletion
Plate:(22)X=mexp(m22),RR0=exp(m22);
Sphere:(23)X=(n3){1-[(n-m)n]3}×exp{(m26)[1+2(n-m)n]},RR0=exp(m26){[1+2(n-m)n]};
Column:(24)X=(n2){1-[(n-m)n]2}×exp{(14)[n2-(n-m)2[1+2ln(n(n-m))]]},RR0=exp{(14)[n2-(n-m)2[1+2ln(n(n-m))]]}.Volume depletion
(25)RR0=(Sa)(KO2PO2)1/2,S=1(plate),2(column),or3(sphere).S is the shape factor, which coincides with the surface area to volume ratio for each crystal shape. In the stage of volume depletion, sensitivity to oxygen decreases in the order, sphere > column > plate, due to the shape factor when a is the same.It would be instructive to estimate roughly the values ofLD and n for SnO2 crystals. Martinelli et al. have reported that the permittivity and the donor density of SnO2 are 10-10 Fm-1 and 4.1×1018 cm-3, respectively, [39]. These values allow us to estimate LD=2.8 nm at 600 K. Actually, however, donor density values are scattered rather largely depending on reports and the highest one is 1.5×1019 cm-3 [27], which leads to LD=1.4 nm. When we assume the former LD value, SnO2 grains of about 20 nm in diameter usually adapted in actual devices have n=3.6. Then regional depletion in these grains is finished at the reduced depletion depth m=n=3.6, which leads to R/R0=8.7, X=10.4, and KO2PO2=1.5×103 nm2 atm-1 according to (23). It follows that this value of R/R0 is the maximum response to oxygen allowed by regional depletion, and therefore the maximum of conventional sensor response to reducing gases in air (Ra/Rg) as well. It is understood that larger responses in R/R0 and Ra/Rg can never be acquired without the help of volume depletion.


### 3.7. Corrections for the Tailing of Electron Distribution

The abrupt model used above cannot always be rationalized, as mentioned. The assumption of complete transfer of electrons from the depletion region to the surface can be too drastic. It is likely that some electrons are left behind in the vicinity of the depletion frontier where potential energy of electrons remains modest, resulting in a tailing of electron distribution. For usual semiconductor devices, a distribution tail due to thermal fluctuations of electrons has been taken into consideration for correction [37]. However, this correction having aimed at applications for thick semiconductor devices hardly seems to be effective in gas sensors. It is worth remembering that the electron concentration at the surface, [e]S, is estimated from VS by assuming the Boltzmann distribution law. Then it would not be so unreasonable to assume that the same distribution law operates over the whole depletion region. If this assumption is admitted, the amount of electrons transferred from the bulk to the surface, which contributes to -QSC/q, can be estimated exactly from the potential energy profiles as follows:
(26)-QSCq=Nd∫a-wa{1-exp(-qV(x)kT)}dx,(Plate),-4πa2QSCq=Nd∫a-wa{1-exp(-qV(r)kT)}4πr2dr,(Sphere),-2πaQSCq=Nd∫a-wa{1-exp(-qV(r)kT)}2πrdr,(Column).
Rearrangements result in
(27)-QSCqNdLD=m-∫n-mnexp(-qV(x)kT)dX,(Plate),(28)-QSCqNdLD=(n3){1-[(n-m)n]3}-(1n2)∫n-mnR2exp(-qV(R)kT)dR,(Sphere),(29)-QSCqNdLD=(n2){1-[(n-m)n]2}-(1n)∫n-mnRexp(-qV(R)kT)dR,(Column).
The second term in each equation represents the quantity to be corrected for the tailing effect. Once such a correction is available, the correlations between R/R0 and X can be formulated in the same way as done before. The results are summarized below.Regional depletion
Plate:(30)X=(m-AP(m))exp(m22),RR0=exp(m22);
Sphere:(31)X={(n3){1-[(n-m)n]3}-AS(n,m)}×exp{(m26)[1+2(n-m)n]},RR0=exp{(m26)[1+2(n-m)n]};
Column:(32)X={(n2){1-[(n-m)n]2}-AC(n,m)}×exp{(14)[n2-(n-m)2[1+2ln(n(n-m))]]},RR0=exp{(14)[n2-(n-m)2[1+2ln(n(n-m))]]}.A(m) or A(n,m) is the correction term appearing in (27) through (29), and suffices, P, S, and C standing for plate, sphere, and column.Volum depletion
Plate:(33)RR0=(1n)X+(1n)AP(n)exp(n22);
Sphere:(34)RR0=(3n)X+(3n)AS(n)exp(n26);
Column:(35)RR0=(2n)X+(2n)AC(n)exp(n24).A(n) is the value of A(m) or A(n,m) at the boundary depletion (m=n). In each equation, the second term (intercept) is close to unity for smaller values of n but it increases rather sharply with increasing n in a range of larger n.The correlations given by these equations are shown for the cases of planar and spherical crystals in (a) and (b) in Figure17, respectively. Compared with Figure 15, the correlations for plate crystals (a) are more convex upward in the stage of regional depletion and this tendency becomes more conspicuous as n increases. Those for volume depletion, on the other hand, are straight lines, which now have intercepts against ordinate. Such behavior of the correlations seems to fit better to actual observations. However, the assumption used for the correction needs to be checked thoroughly for verification.Correlations betweenR/R0 and X for plate crystals (a) and spherical crystals (b) with reduced size n [15].
(a)(b)


### 3.8. Comparison with Experiments

Although the measurement of device resistance as a function of PO2 is a very basic task, carrying out it satisfactorily has turned out to be surprisingly difficult. The resistance is sensitive to trace amounts of surface-active impurities, that is, contaminant reducing gases at ppm or sub-ppm levels in commercial tubes of oxygen or of oxygen diluted with nitrogen, and contaminant oxygen as well as contaminant reducing gases in high purity nitrogen. Since these tube gases are mixed together to control PO2, contaminants can be serious origins of disturbances. It is necessary to reduce contaminants as much as possible but complete reduction has been yet to achieve. With utmost precautions about contaminants, fairly reliable data acquisition is under way in our laboratory. For example, the data obtained with screen-printed thick films of SnO2 of two different grain sizes (12 and 16 nm in diameter) at 300 and 350°C are shown in Figure 18. To check fitting between the data to the present theory, we need to know R0 but it has been difficult to measure it precisely yet due to contaminant oxygen (about 100 ppm in volume) in nitrogen. Here the resistance in nitrogen, R(N2), is substituted for R0. The resistance data normalized by R(N2) are thus plotted as a function of square root of PO2 in the figure. In spite of the ambiguity just mentioned, all the data at 300°C fall on a straight line fairly well in agreement with (25) or (34) derived for volume depletion. The data at 350°C, however, falls on a line which is convex upward in the beginning and then becomes linear, a feature characteristic to a change from regional depletion to volume depletion. From the slope of each straight line, full or partial, and grain size, we can deduce KO2 using (25) or (34). The slopes for 12 nm and 16 nm grains are 157 and 118 nm atm-1/2 at 300°C, respectively, both giving KO2=9.9×104 nm2 atm-1. The corresponding data at 350°C are 89 and 83 nm atm-1/2, giving KO2=3.1 and 4.9 in unit of 104 nm2 atm-1, respectively. In principle, KO2 is independent of grain sizes and this is in fact confirmed with these experimental data. These results are fairly satisfactory and   seem to assure the soundness of the present theory.Figure 18
Normalized response to oxygen,R/R(N2), at two designated temperatures as correlated with PO21/2 for two screen-printed SnO2 thick film devices different in grain size. R(N2); resistance in N2.Let us consider next the grain size effect onRa already shown in Figure 4. The same data are plotted on logarithmic scales in Figure 19. When volume depletion is assumed to prevail in the whole grains size range tested, Ra/R0 is inversely proportional to a according to (23) or (34), while it can be shown that R0 is also inversely proportional to a if ideal packing of spheres is assumed. As a result, Ra should be essentially constant regardless of the variation of a. In fact, this is confirmed in the range of a above 3.5 nm (ac). The sharp increase of Ra observed with decreasing a below ac can be attributed to an increase in the number of grains without donors (insulating grains). This can happen when the grain size (diameter) is made smaller than an average separation between neighboring donors under the condition of a fixed donor density (Nd). From ac, Nd is estimated to be 5.6×1018 cm-3, which is fairly close to reported values mentioned before. This in turn leads to LD=2.4 nm at 600 K. In this way, the grain size effect on base resistance is understood in a way largely different from what was imagined in the past.Figure 19
Grain size effect on Ra shown in Figure4 as reproduced on the logarithmic scales [16].


## 3.1. Conventional Theory of Depletion (Large Crystals)

Depletion of large semiconductor crystals can be treated with the conventional theory. Let us consider the adsorption of electron-accepting molecules (typically O2) on the surface of a large n-type semiconductor crystal. This situation is well approximated by the adsorption on a plane of a large crystal cleaved in halves, as schematically drawn in Figure 12 [37]. As a result of electron transfer, a depletion layer is formed in the subsurface region, the layer thickness increasing with an increase in electron transfer. The electron transfer reaches equilibrium when Fermi level at the adsorbed species coincides with that in the semiconductor bulk.Figure 12
Schematic diagram of electrostatic equilibrium for large semiconductor crystals under exposure to oxygen. Ec: conduction band edge, Ec: valence band edge, EF: Fermi level, Qsc: surface charge density, w: depletion depth, qV: potential energy of electrons, qVs: surface potential energy barrier [37].The electrostatic equilibrium inside the crystal in such a case can be treated quantitatively under simplifying assumptions. Conventionally it is assumed that (1) donors are ionized completely, (2) there are no other electron-traps other than the adsorbates, and (3) all of conduction electrons up to a depth w are completely transferred to the surface (abrupt model). Since donors in SnO2 are reportedly present at shallow levels (0.05–0.15 eV) below the conduction band edge [38], the first assumption can be acceptable fairly well. The second one may be fit to an ionic semiconductor like SnO2 better than to a nonionic semiconductor like silicon. The third one, though useful in simplifying mathematical treatments, cannot always be rationalized as mentioned later. However, abrupt model is used here always unless otherwise noted.Under these assumptions, the density of surface charges,QSC, is nominally equal to -qNdw, where q is electric charge of proton and Nd the density of donors of semiconductor:
(2)QSC=-qNdw=-qNdLDm.
Here m is reduced depletion depth defined by m=w/LD. Under these conditions, electric potential, V, in the depletion region should satisfy the following Poisson’s equation, where x is the depth from the surface and ε is the permittivity of semiconductor times the free space permittivity:
(3)d2Vdx2=-qNdε.
By introducing the boundary conditions that dV/dx as well as V are zero at x=w, (3) is solved to give the following depth profile of electric potential:
(4)V(x)=-(qNd2ε)(x-w)2.
Potential energy of electrons, qV(x), as well as the height of surface potential barrier, qVS, are obtained from (4):
(5)qV(x)kT=(q2Nd2εkT)(x-w)2=(12){(x-w)LD}2,qV(X)kT=(12)(X-m)2,(6)qVSkT=(12)(wLD)2=(12)m2.
Here LD is the Debye length of semiconductor given by LD=(q2Nd/kT)-1/2, and X is reduced distance defined by X=x/LD. The density of electrons at the surface, [e]S, is given from (6) as follows:
(7)[e]S=Ndexp(-qVSkT)=Ndexp(-m22).
In this way, the surface density of adsorbates, -QSC/q, and electrons can be correlated uniquely through m once LD and Nd are given.


## 3.2. Response of Large Crystals

It is assumed that oxygen is adsorbed asO- exclusively on the semiconductor surface. Then the overall equation of oxygen adsorption and its equilibrium are expressed as follows:
R1O2+2e=2O-,(8)(KO2PO2)1/2[e]S=[O-].KO2 and PO2 are the adsorption constant and the partial pressure of oxygen, respectively, and [e]S and [O-] are the surface densities of free electrons and O-, respectively. Equation (8) indicates that (KO2PO2)1/2 determines the density ratio, [O-]/[e]S, not the respective quantities. These quantities are determined only after the chemical equilibrium is coupled with the electrostatic equilibrium just mentioned.WhenO- is solely responsible for the total charge density of the semiconductor as assumed,
(9)[O-]=-QSCq=NdLDm.
By inserting (2) and (8) into (7), we obtain
(10)(KO2PO2)1/2Ndexp(-m22)=NdLDm,X=mexp(m22).
Here X is reduced adsorption strength of oxygen defined by X=(KO2PO2)1/2/LD. Indication of (10) is that if X is given, m is determined uniquely. This in turn determines [e]S through (8). Thus we can correlate oxygen adsorption with the surface density of electrons.In order to obtain the correlation with sensor resistance (R), we need to know the relation between R and [e]S. The relation is simple if the tunneling transport model (Figure 10(b)) is assumed, as follows:
(11)RR0=[e]S0[e]S=exp(m22).
Here R0 and [e]S0 are the values of R and [e]S at the flat band state, respectively. The latter quantity coincides with Nd for the simplified mode used here. R and R0 are influenced by many factors other than oxygen adsorption, such as packing of crystals and geometry of sensor device, but those influences can be canceled out effectively by normalizing R by R0. The correlation between R/R0 and X thus obtained through (10) and (11) is shown on logarithmic scales in Figure 13 [14]. R/R0 increases with increasing X and its slope becomes equal to unity in the middle region of X, explaining that the power law-type dependence of R on PO2 with a power index of 1/2 is actually observed. However, no important information else can be obtained from such a treatise. Similar conclusions have been obtained for the responses to NO2 and H2 [14]. This is because we have used conventional theory of depletion developed for large crystals, thus missing important information about the size and shape of crystals.Figure 13
Correlation betweenR/R0 and X on logarithmic scales [14]. X=(KO2PO2)1/2/LD


## 3.3. Two Types of Depletion in Small Crystals

Let us consider how depletion develops in a plate crystal of thickness2a (a, half thickness), starting from the flat-band state, with increasing PO2 in the ambient. The electrostatic equilibrium of this system is schematically drawn in Figure 14(a). When PO2 is small, depletion takes place in the same way as previously seen in a large crystal, depletion depth increasing with increasing PO2. However, this way of depletion is finished at PO2(II) where depletion has extended to cover the whole region of the crystal. Obviously the value of PO2(II) decreases with decreasing half thickness (a). What would happen when PO2 is further increased to PO2(III)) in the figure? This is a main concern here. The above limit of depletion comes solely from the physics of the crystal. The chemistry of oxygen adsorption, on the other hand, does not care the physics. It requires (8) to hold at any stages, or, speaking more specifically, the ratio [O-]/[e]S to increase linearly with increasing X=(KO2PO2)1/2. In order for the physics side (electrostatic equilibrium) to fulfill this requirement, the only available way is to allow Fermi level to shift down by a designated amount, pkT, whereas the band bending profile is kept unchanged. With this shift, not only the whole region from the surface adsorbates (O-) to the bulk of the crystal can be placed at the same Fermi level but also the potential energy of electrons goes up to the level at which [e]S takes a value consistent with the chemical equilibrium condition (8). This way of establishing the electrostatic equilibrium is thus very different from one observed at lower PO2. It is noted that electron transfer from the bulk to the surface continues still in this stage, though it is only nominal when the abrupt model is used (Figure 14(b)). The electron transfer is made more apparent in case some electrons are left behind in the depletion layer (modified model) as illustrated in the same figure. In this case, the density of electrons inside the crystal decreases progressively as the Fermi level shift increases, showing that electrons are squeezed out from the whole crystal to be supplied to the surface. It is therefore convenient to distinguish the depletion in this stage from the conventional one in the initial stage up to PO2(II). We call them here volume depletion and regional depletion, respectively, while the border in between is called boundary depletion.Schematic diagrams of electronic equilibrium for semiconductor plate (half thicknessa) placed under various partial pressures of oxygen (PO2). (a) Depth profiles of potential energy. (b) Depth profiles of conduction electrons density assumed on abrupt model (broken lines) and modified model (full lines). See section 3.6 for modified model [15].
(a)(b)


## 3.4. Response at Regional Depletion (Plate)

As mentioned earlier, the occurrence of regional depletion followed by volume depletion is characteristic to small crystals. Receptor function and response of small plate crystals to oxygen in the stage of regional depletion are discussed here. The distribution of potential inside the plate is obtained by solving Poisson’s equation, which is the same as (3) but has an origin located at a center of the plate with coordinate x being taken as a distance normal to the surface from origin (see Figure 14). In the stage of regional depletion, potential energy of electrons is expressed as shown below:
(12)qV(x)=(q2Nd2ε){x-(a-w)}2for|x|≧a-w,=0for|x|<a-w.
These are rewritten in reduced form as
(13)qV(X)kT=(12){X-(n-m)}2for|x|≧n-m,=0for|x|<n-m.
Reduced quantities, X and m, are the same as defined before, and n is reduced size defined by n=a/LD. These expressions are essentially the same as those derived for large crystals, though depletion in this case starts from both surfaces. Surface potential energy (barrier) is expressed by the same equation (7), while (9) holds also as assumed. As a result, either of (10) and (11) derived for large crystals can be used to correlate between R/R0 and X also in this case. It is noted, however, that this can be valid only for plate crystals, which are free from the shape effect.


## 3.5. Response at Volume Depletion (Plate)

Potential energy in the stage of boundary depletion is obtained by insertingw=a in (12) or m=n in (13) and that in volume depletion is larger by pkT or p, as follows:
(14)qV(x)=(q2Nd2ε)x2+pkT.(15)qV(X)kT=(12)X2+p.
Here p=0 indicates boundary depletion. The surface potential energy is obtained by inserting X=n into (15):
(16)qVskT=(12)n2+p.
We can then express [e]S and R/R0 by using n and p. [O-] is obtained from (9) if m is replaced by n. By inserting these quantities into the mass action law (8), we can correlate X with n and p. Thus we obtain a couple of equations:
(17)X=nexp{(12)n2+p},RR0=exp{(12)n2+p}.
These equations determine the correlation between R/R0 and X. Unlike the case of regional depletion, however, R/R0 is an explicit function of X, because we obtain from the couple of equations
(18)RR0=Xn=(KO2PO2)1/2a.
Normalized resistance is thus linear to (KO2PO2)1/2 and the reciprocal of a.The correlations betweenR/R0 and X in the stages of regional depletion and volume depletion are drawn for smaller selected values of n in Figure 15. In principle, regional depletion always precedes another on increasing X, but actually when n is small it is not visible well in the drawing on this scale, allowing volume depletion to prevail. It is also clear that the slope of each correlation (sensitivity to oxygen) increases sharply with decreasing n and is inversely proportional to n or a in the stage of volume depletion. This indicates a size effect on the sensor response to oxygen. It is remarked that the derivative of R or normalized R with respect to the concentration (or partial pressure) of target gas in problem is usually defined as sensor sensitivity. In the case of oxygen adsorption, however, the sensitivity is better defined as the derivative of R/R0 with respect to PO21/2, which is equal to KO21/2/a in the stage of volume depletion.Figure 15
Correlations betweenR/R0 and X for plate crystals with reduced thickness n [17].


## 3.6. Effects of Shape and Size

The same treatment can be extended to spherical crystals (radiusa) and columnar crystals (columnar radius a). For this purpose, two things should be taken into account. First, Poisson’s equation is transformed, for the convenience of solution, into one fitting to the symmetry of the crystal form. A polar spherical coordinates system and a columnar coordinates system are illustrated together with one-dimensional system so far used in Figure 16. In the respective systems, Poisson’s equation is expressed as a function of a single variable as follows:
(19)d2Vdx2=-qNdε,(Plate),(1r2)ddr(r2dVdr)=-qNdε,(Sphere),(1r)ddr(rdVdr)=-qNdε,(Column).
Here r is radial displacement from the center of crystals. Each equation can be solved by using the same boundary conditions as used before to give a solution shown below:
(20)V(x)=-(qNd2ε){x-(a-w)}2,(Plate),V(r)=-(qNd/2ε){r-(a-w)}2{1+2(a-w)/r}3,(Sphere),V(r)=-(qNd4ε){r2-(a-w)2[1+2ln(r(a–w))]},(Column).
Second, the relations between total surface charge density (QSC) and depletion depth (w or m) should be altered to meet crystal symmetry:
(21)QSCq=-Ndw=-NdLDm,(Plate),QSCq=-Nd{∫a-wa4πr2dr}(4πa2)=-(13)NdLDn{1-[(n-m)n]3},(Sphere),QSCq=-Nd{∫a-wa2πrdr}(2πa)=-(12)NdLDn{1-[(n-m)n]2},(Column).
Here m and n are reduced depletion depth (w/LD) and reduced size (a/LD), respectively.Figure 16
Coordinates systems adapted for plate, sphere, and column [15].Based on the abrupt model, receptor function and sensor response to oxygen can be formulated for spherical and columnar crystals in the same way as done for plate crystals. The correlations betweenR/R0 and X thus derived are summarized as follows.Regional depletion
Plate:(22)X=mexp(m22),RR0=exp(m22);
Sphere:(23)X=(n3){1-[(n-m)n]3}×exp{(m26)[1+2(n-m)n]},RR0=exp(m26){[1+2(n-m)n]};
Column:(24)X=(n2){1-[(n-m)n]2}×exp{(14)[n2-(n-m)2[1+2ln(n(n-m))]]},RR0=exp{(14)[n2-(n-m)2[1+2ln(n(n-m))]]}.Volume depletion
(25)RR0=(Sa)(KO2PO2)1/2,S=1(plate),2(column),or3(sphere).S is the shape factor, which coincides with the surface area to volume ratio for each crystal shape. In the stage of volume depletion, sensitivity to oxygen decreases in the order, sphere > column > plate, due to the shape factor when a is the same.It would be instructive to estimate roughly the values ofLD and n for SnO2 crystals. Martinelli et al. have reported that the permittivity and the donor density of SnO2 are 10-10 Fm-1 and 4.1×1018 cm-3, respectively, [39]. These values allow us to estimate LD=2.8 nm at 600 K. Actually, however, donor density values are scattered rather largely depending on reports and the highest one is 1.5×1019 cm-3 [27], which leads to LD=1.4 nm. When we assume the former LD value, SnO2 grains of about 20 nm in diameter usually adapted in actual devices have n=3.6. Then regional depletion in these grains is finished at the reduced depletion depth m=n=3.6, which leads to R/R0=8.7, X=10.4, and KO2PO2=1.5×103 nm2 atm-1 according to (23). It follows that this value of R/R0 is the maximum response to oxygen allowed by regional depletion, and therefore the maximum of conventional sensor response to reducing gases in air (Ra/Rg) as well. It is understood that larger responses in R/R0 and Ra/Rg can never be acquired without the help of volume depletion.


## 3.7. Corrections for the Tailing of Electron Distribution

The abrupt model used above cannot always be rationalized, as mentioned. The assumption of complete transfer of electrons from the depletion region to the surface can be too drastic. It is likely that some electrons are left behind in the vicinity of the depletion frontier where potential energy of electrons remains modest, resulting in a tailing of electron distribution. For usual semiconductor devices, a distribution tail due to thermal fluctuations of electrons has been taken into consideration for correction [37]. However, this correction having aimed at applications for thick semiconductor devices hardly seems to be effective in gas sensors. It is worth remembering that the electron concentration at the surface, [e]S, is estimated from VS by assuming the Boltzmann distribution law. Then it would not be so unreasonable to assume that the same distribution law operates over the whole depletion region. If this assumption is admitted, the amount of electrons transferred from the bulk to the surface, which contributes to -QSC/q, can be estimated exactly from the potential energy profiles as follows:
(26)-QSCq=Nd∫a-wa{1-exp(-qV(x)kT)}dx,(Plate),-4πa2QSCq=Nd∫a-wa{1-exp(-qV(r)kT)}4πr2dr,(Sphere),-2πaQSCq=Nd∫a-wa{1-exp(-qV(r)kT)}2πrdr,(Column).
Rearrangements result in
(27)-QSCqNdLD=m-∫n-mnexp(-qV(x)kT)dX,(Plate),(28)-QSCqNdLD=(n3){1-[(n-m)n]3}-(1n2)∫n-mnR2exp(-qV(R)kT)dR,(Sphere),(29)-QSCqNdLD=(n2){1-[(n-m)n]2}-(1n)∫n-mnRexp(-qV(R)kT)dR,(Column).
The second term in each equation represents the quantity to be corrected for the tailing effect. Once such a correction is available, the correlations between R/R0 and X can be formulated in the same way as done before. The results are summarized below.Regional depletion
Plate:(30)X=(m-AP(m))exp(m22),RR0=exp(m22);
Sphere:(31)X={(n3){1-[(n-m)n]3}-AS(n,m)}×exp{(m26)[1+2(n-m)n]},RR0=exp{(m26)[1+2(n-m)n]};
Column:(32)X={(n2){1-[(n-m)n]2}-AC(n,m)}×exp{(14)[n2-(n-m)2[1+2ln(n(n-m))]]},RR0=exp{(14)[n2-(n-m)2[1+2ln(n(n-m))]]}.A(m) or A(n,m) is the correction term appearing in (27) through (29), and suffices, P, S, and C standing for plate, sphere, and column.Volum depletion
Plate:(33)RR0=(1n)X+(1n)AP(n)exp(n22);
Sphere:(34)RR0=(3n)X+(3n)AS(n)exp(n26);
Column:(35)RR0=(2n)X+(2n)AC(n)exp(n24).A(n) is the value of A(m) or A(n,m) at the boundary depletion (m=n). In each equation, the second term (intercept) is close to unity for smaller values of n but it increases rather sharply with increasing n in a range of larger n.The correlations given by these equations are shown for the cases of planar and spherical crystals in (a) and (b) in Figure17, respectively. Compared with Figure 15, the correlations for plate crystals (a) are more convex upward in the stage of regional depletion and this tendency becomes more conspicuous as n increases. Those for volume depletion, on the other hand, are straight lines, which now have intercepts against ordinate. Such behavior of the correlations seems to fit better to actual observations. However, the assumption used for the correction needs to be checked thoroughly for verification.Correlations betweenR/R0 and X for plate crystals (a) and spherical crystals (b) with reduced size n [15].
(a)(b)


## 3.8. Comparison with Experiments

Although the measurement of device resistance as a function of PO2 is a very basic task, carrying out it satisfactorily has turned out to be surprisingly difficult. The resistance is sensitive to trace amounts of surface-active impurities, that is, contaminant reducing gases at ppm or sub-ppm levels in commercial tubes of oxygen or of oxygen diluted with nitrogen, and contaminant oxygen as well as contaminant reducing gases in high purity nitrogen. Since these tube gases are mixed together to control PO2, contaminants can be serious origins of disturbances. It is necessary to reduce contaminants as much as possible but complete reduction has been yet to achieve. With utmost precautions about contaminants, fairly reliable data acquisition is under way in our laboratory. For example, the data obtained with screen-printed thick films of SnO2 of two different grain sizes (12 and 16 nm in diameter) at 300 and 350°C are shown in Figure 18. To check fitting between the data to the present theory, we need to know R0 but it has been difficult to measure it precisely yet due to contaminant oxygen (about 100 ppm in volume) in nitrogen. Here the resistance in nitrogen, R(N2), is substituted for R0. The resistance data normalized by R(N2) are thus plotted as a function of square root of PO2 in the figure. In spite of the ambiguity just mentioned, all the data at 300°C fall on a straight line fairly well in agreement with (25) or (34) derived for volume depletion. The data at 350°C, however, falls on a line which is convex upward in the beginning and then becomes linear, a feature characteristic to a change from regional depletion to volume depletion. From the slope of each straight line, full or partial, and grain size, we can deduce KO2 using (25) or (34). The slopes for 12 nm and 16 nm grains are 157 and 118 nm atm-1/2 at 300°C, respectively, both giving KO2=9.9×104 nm2 atm-1. The corresponding data at 350°C are 89 and 83 nm atm-1/2, giving KO2=3.1 and 4.9 in unit of 104 nm2 atm-1, respectively. In principle, KO2 is independent of grain sizes and this is in fact confirmed with these experimental data. These results are fairly satisfactory and   seem to assure the soundness of the present theory.Figure 18
Normalized response to oxygen,R/R(N2), at two designated temperatures as correlated with PO21/2 for two screen-printed SnO2 thick film devices different in grain size. R(N2); resistance in N2.Let us consider next the grain size effect onRa already shown in Figure 4. The same data are plotted on logarithmic scales in Figure 19. When volume depletion is assumed to prevail in the whole grains size range tested, Ra/R0 is inversely proportional to a according to (23) or (34), while it can be shown that R0 is also inversely proportional to a if ideal packing of spheres is assumed. As a result, Ra should be essentially constant regardless of the variation of a. In fact, this is confirmed in the range of a above 3.5 nm (ac). The sharp increase of Ra observed with decreasing a below ac can be attributed to an increase in the number of grains without donors (insulating grains). This can happen when the grain size (diameter) is made smaller than an average separation between neighboring donors under the condition of a fixed donor density (Nd). From ac, Nd is estimated to be 5.6×1018 cm-3, which is fairly close to reported values mentioned before. This in turn leads to LD=2.4 nm at 600 K. In this way, the grain size effect on base resistance is understood in a way largely different from what was imagined in the past.Figure 19
Grain size effect on Ra shown in Figure4 as reproduced on the logarithmic scales [16].


## 4. Receptor Function and Response to NO2

### 4.1. Modeling

NO2 is selected here as a representative of oxidizing gases which are adsorbed as negatively charged ions on the semiconductor surface. As generally accepted, it is assumed to undergo adsorption as follows:
R2NO2+e=NO2-
The equilibrium of adsorption is written down:
(36)KNO2PNO2[e]S=[NO2-].KNO2 and PNO2 are the adsorption constant and the partial pressure of NO2, respectively, and [NO2-] is the surface density of NO2-. As noted, (R2) and (36) are quite similar to the corresponding equations for O2, (R1) and (8). This means that, in the absence of oxygen, receptor function and response to NO2 can be treated in exactly the same way as done to oxygen so far. That is, if we define the reduced adsorptive strength of NO2 as Z=KNO2PNO2/LD, all of the correlations between R/R0 and X(=(KO2PO2)1/2/LD) obtained to O2 above can be transferred to the correlations to NO2 by replacing X by Z.In the coexistence of oxygen, we are placed in a different situation. Now both O2 and NO2 are adsorbed on the surface of crystals to compete for electrons. The mass action laws for O2 and NO2, (8) and (36), should be satisfied simultaneously. That is,
(37)[O-]+[NO2-]={(KO2PO2)1/2+KNO2PNO2}[e]S.
This is a requirement from the chemistry side. In the physics side, it is required to have Fermi level kept the same at any sites among O-, NO2-, and bulk of crystal in the band energy diagram, and nothing else is changed except that the total surface charge density is now contributed by the two kinds of adsorbates:
(38)-QSCq=[O-]+[NO2-].
As a result, all the equations involving X so far derived to oxygen can survive if X is replaced by X+Z. That is, the response of plate crystals (half thickness a) is expressed as follows.Regional depletion
(39)X+Z=mexp(m22),RR0=exp(m22).Volume depletion
(40)RR0=(X+Z)n={(KO2PO2)1/2+KNO2PNO2}a.The meaning of (39) is the same as before; X+Z determines m unequely and then does R/R0 also unequely. Naturally, the correlations between R/R0 and X shown in Figure 15 are in effect in this case, too, if the abscissa is redefined to be X+Z. Correlations in the stage of volume depletion can be derived for other shapes of crystals and are summarized as follows:
(41)RR0=(X+Z)(Sn)={(KO2PO2)1/2+KNO2PNO2}(Sa),S=1(plate),2(column),and3(sphere).


### 4.2. Response to NO2 in Air

In the sensory measurement in air,PO2 is fixed at PO2 (air) and this determines the resistance level in air (air base), Ra/R0. Conventionally, sensor response has often been defined as the ratio of device resistance in gas (Rg) to that in air (Ra). The conventional response (Rg/Ra) thus defined corresponds to the ratio of the response to X+Z to that to X in the present treatment as follows:
(42)RgRa=(Rg/R0)(Ra/R0).
Let us consider the meaning of this using Figure 20, which depicts the correlations between R/R0 and X+Z for plates of two different thicknesses (n=1 and 3). If X is assumed to take a value of 100, the air base (Ra/R0) for n=1 is already on a linear correlation line of volume depletion, and addition of Z simply increases the response along with the same line. From (40) or (41), we can derive
(43)RgRa={KNO2PNO2(KO2PO2)1/2}+1.
In this case the conventional response reflects the ratio of adsorptive strengths of the two components. In the case of thicker plates (n=3), on the other hand, the air base point is still in the stage of regional depletion and, with an increase of Z, the response goes up toward volume depletion. The conventional response cannot be expressed in a simple equation like (43). The convex nature of the correlation lines tends to increase when the modified model is used (Figure 17), making the conventional response more nonlinear.Figure 20
Correlations betweenR/R0 and X+Z derived for plate crystals [17].


### 4.3. Comparison with Experiments

It has been shown that WO3 is very sensitive to NO2, though the sensitivity depends heavily on the methods of fabricating devices. Conventional sensor responses to NO2 observed with two devices different in fabrication method are shown in Figures 21 and 22. For the former device, WO3 was prepared from ammonium paratungstate through pyrolysis at 600°C for 5 hours and the resulting crystals were granular (grain size not measured) [40]. The response is seen to be linear to PNO2 in the tested range up to 80 ppm at 200 and 300°C but it changes to be apparently nonlinear at 400°C. On the basis of the present theory, the linear response at the lower temperatures reflects volume depletion, while that at the highest temperature does regional depletion. This change is easily understood from Figure 17(b) (after replacing X by X+Z) if we assume that Z is much larger than X, while KNO2 decreases exponentially with increasing temperature.Figure 21
Sensor response to NO2 at designated temperatures as correlated with NO2concentrations for granular WO3 crystals prepared by pyrolysis [16].Figure 22
Sensor response to NO2 at designated temperatures as correlated with NO2concentrations for lamellar WO3 crystals prepared by wet process [16].In the case of Figure22, WO3 was prepared through a sol-gel process into shape of thin lamellae (ca 30 nm thick) [41]. The sensor response is seen to be linear at both temperatures of 200 and 250°C reflecting the involvement of volume depletion. It is as high as being able to detect NO2 down to a level of 10 ppb at 200°C. Such high sensitivity comes from the smallness of lamellar size used, as understood from (40). From (40) and (42), the slope of each correlation, 9.4 or 7.4 in unit of 107 nm atm-1, respectively, gives rise to (KNO2/a)/(Ra/R0). Thus we can estimate KNO2 if Ra/R0 is known. Unfortunately, this value is not known yet, but even when Ra/R0 is assumed to take an ultimately small value of unity, KNO2 can amount to a value as large as 1.4 or 1.1 in unit of 109 nm atm-1, respectively. This confirms extremely strong nature of NO2 adsorption.


## 4.1. Modeling

NO2 is selected here as a representative of oxidizing gases which are adsorbed as negatively charged ions on the semiconductor surface. As generally accepted, it is assumed to undergo adsorption as follows:
R2NO2+e=NO2-
The equilibrium of adsorption is written down:
(36)KNO2PNO2[e]S=[NO2-].KNO2 and PNO2 are the adsorption constant and the partial pressure of NO2, respectively, and [NO2-] is the surface density of NO2-. As noted, (R2) and (36) are quite similar to the corresponding equations for O2, (R1) and (8). This means that, in the absence of oxygen, receptor function and response to NO2 can be treated in exactly the same way as done to oxygen so far. That is, if we define the reduced adsorptive strength of NO2 as Z=KNO2PNO2/LD, all of the correlations between R/R0 and X(=(KO2PO2)1/2/LD) obtained to O2 above can be transferred to the correlations to NO2 by replacing X by Z.In the coexistence of oxygen, we are placed in a different situation. Now both O2 and NO2 are adsorbed on the surface of crystals to compete for electrons. The mass action laws for O2 and NO2, (8) and (36), should be satisfied simultaneously. That is,
(37)[O-]+[NO2-]={(KO2PO2)1/2+KNO2PNO2}[e]S.
This is a requirement from the chemistry side. In the physics side, it is required to have Fermi level kept the same at any sites among O-, NO2-, and bulk of crystal in the band energy diagram, and nothing else is changed except that the total surface charge density is now contributed by the two kinds of adsorbates:
(38)-QSCq=[O-]+[NO2-].
As a result, all the equations involving X so far derived to oxygen can survive if X is replaced by X+Z. That is, the response of plate crystals (half thickness a) is expressed as follows.Regional depletion
(39)X+Z=mexp(m22),RR0=exp(m22).Volume depletion
(40)RR0=(X+Z)n={(KO2PO2)1/2+KNO2PNO2}a.The meaning of (39) is the same as before; X+Z determines m unequely and then does R/R0 also unequely. Naturally, the correlations between R/R0 and X shown in Figure 15 are in effect in this case, too, if the abscissa is redefined to be X+Z. Correlations in the stage of volume depletion can be derived for other shapes of crystals and are summarized as follows:
(41)RR0=(X+Z)(Sn)={(KO2PO2)1/2+KNO2PNO2}(Sa),S=1(plate),2(column),and3(sphere).


## 4.2. Response to NO2 in Air

In the sensory measurement in air,PO2 is fixed at PO2 (air) and this determines the resistance level in air (air base), Ra/R0. Conventionally, sensor response has often been defined as the ratio of device resistance in gas (Rg) to that in air (Ra). The conventional response (Rg/Ra) thus defined corresponds to the ratio of the response to X+Z to that to X in the present treatment as follows:
(42)RgRa=(Rg/R0)(Ra/R0).
Let us consider the meaning of this using Figure 20, which depicts the correlations between R/R0 and X+Z for plates of two different thicknesses (n=1 and 3). If X is assumed to take a value of 100, the air base (Ra/R0) for n=1 is already on a linear correlation line of volume depletion, and addition of Z simply increases the response along with the same line. From (40) or (41), we can derive
(43)RgRa={KNO2PNO2(KO2PO2)1/2}+1.
In this case the conventional response reflects the ratio of adsorptive strengths of the two components. In the case of thicker plates (n=3), on the other hand, the air base point is still in the stage of regional depletion and, with an increase of Z, the response goes up toward volume depletion. The conventional response cannot be expressed in a simple equation like (43). The convex nature of the correlation lines tends to increase when the modified model is used (Figure 17), making the conventional response more nonlinear.Figure 20
Correlations betweenR/R0 and X+Z derived for plate crystals [17].


## 4.3. Comparison with Experiments

It has been shown that WO3 is very sensitive to NO2, though the sensitivity depends heavily on the methods of fabricating devices. Conventional sensor responses to NO2 observed with two devices different in fabrication method are shown in Figures 21 and 22. For the former device, WO3 was prepared from ammonium paratungstate through pyrolysis at 600°C for 5 hours and the resulting crystals were granular (grain size not measured) [40]. The response is seen to be linear to PNO2 in the tested range up to 80 ppm at 200 and 300°C but it changes to be apparently nonlinear at 400°C. On the basis of the present theory, the linear response at the lower temperatures reflects volume depletion, while that at the highest temperature does regional depletion. This change is easily understood from Figure 17(b) (after replacing X by X+Z) if we assume that Z is much larger than X, while KNO2 decreases exponentially with increasing temperature.Figure 21
Sensor response to NO2 at designated temperatures as correlated with NO2concentrations for granular WO3 crystals prepared by pyrolysis [16].Figure 22
Sensor response to NO2 at designated temperatures as correlated with NO2concentrations for lamellar WO3 crystals prepared by wet process [16].In the case of Figure22, WO3 was prepared through a sol-gel process into shape of thin lamellae (ca 30 nm thick) [41]. The sensor response is seen to be linear at both temperatures of 200 and 250°C reflecting the involvement of volume depletion. It is as high as being able to detect NO2 down to a level of 10 ppb at 200°C. Such high sensitivity comes from the smallness of lamellar size used, as understood from (40). From (40) and (42), the slope of each correlation, 9.4 or 7.4 in unit of 107 nm atm-1, respectively, gives rise to (KNO2/a)/(Ra/R0). Thus we can estimate KNO2 if Ra/R0 is known. Unfortunately, this value is not known yet, but even when Ra/R0 is assumed to take an ultimately small value of unity, KNO2 can amount to a value as large as 1.4 or 1.1 in unit of 109 nm atm-1, respectively. This confirms extremely strong nature of NO2 adsorption.


## 5. Response to Hydrogen

### 5.1. Modeling

Hydrogen is selected here as a representative of reducing gases, which decreases the device resistance. It is unanimous that H2 molecules react with O- adsorbates on the semiconductor crystals. We assume that H2 molecules collide directly with O- ions (Eley-Rideal mechanism) and H2O molecules formed are desorbed instantly from the surface:
R3O-+H2=H2O+e.
This reaction consumes O- ions, while they are resupplied from the gas phase through the reaction, (R1).The rate of accumulation of [O-] is expressed by
(44)d[O-]dt=k1PO2[e]S-k-1[O-]2-k3PH2[O-].
Here t is time, k1 and k-1 are the rate constants of forward and reverse reactions of (R1), respectively, k3 is that of forward reaction of (R3), and PH2 is the partial pressure of H2. At steady state, d[O-]/dt=0, so that we obtain
(45)(KO2PO2)1/2[e]S=[O-]{1+(cPH2[O-])}1/2.
Here KO2=k1/k-1 (equilibrium adsorption constant) and c=k3/k-1. In the absence of H2, this equation coincides with the mass action law of the oxygen adsorption equilibrium, (8). This means that the presence of H2modifies the value of [O-] simply according to (45). The electrostatic equilibrium of the crystal, on the other hand, does not see H2 but [O-], requiring the depletion depth to change correspondingly. Thus the equations relating VS and [O-] with m, n, and p in Section 3 can survive altogether in this new situation. Inserting those into (45) followed by rearrangement, we obtain an equation relating among X(=(KO2PO2)1/2/LD), Y(=(c/LDNd)PH2), and m, while R/R0 is always derived from VS. The sensor resistance measured to a mixture of H2 and oxygen is expressed here as Rg in order to distinguish it from that measured to oxygen only. The sought correlations for plates can be summarized as follows.Regional depletion
(46)Xexp(-m22)=m{1+(1m)Y}1/2,RgR0=exp(m22)(Plate).Volume depletion
(47)RgR0=S(X/n){1+(S/n)Y}1/2,S=1(plate),2(column),and3(sphere).These equations show that under fixedX (or PO2), Rg/R0 is an implicit or explicit function of Y. The equations in the stage of regional depletion for other shapes are excluded here for brevity.


### 5.2. Response to H2 in Air

In usual sensor operation,PO2 is fixed at that (0.21 atm) of the air. Sensor response is often defined as the ratio of the device resistance in air to that in the target gas (or normalized conductance), Ra/Rg. This corresponds to the response ratio in the present treatment as follows:(48)RaRg=(Ra/R0)(Rg/R0).
Each term in the right-hand side can be estimated as a function of Y for fixed X by (46) or (47) so that the conventional sensor response can also be done. For example, Figure 23 shows the sensor response as a function of Y1/2 for plate crystals, where n is selected between 0.1 and 3 and X is fixed at 1000. As easily understood, even when a large X is selected in order to secure deep stages of volume depletion in the absence of H2, the crystals are brought to regional depletion ultimately on increasing Y, and this tendency is enhanced with increasing n. Regional depletion is seen to begin at a relatively small value of Y1/2 for n=3. In the range of Y1/2 where volume depletion prevails, (47) holds so that (48) is transformed into an explicit function of Y: (49)RaRg={1+(1n)Y}1/2={1+(caNd)PH2}1/2.
When (c/aNd)PH2 is far larger than unity,(50)RaRg={Yn}1/2={1+(caNd)PH2}1/2.
This is why linear correlations appear in the range of volume depletion as shown. In the range of regional depletion, such linearity is not allowed, though detailed discussion is reserved here.Figure 23
Theoretical sensor response to H2as correlated with Y1/2 at X=1000 [17].


### 5.3. Comparison with Experiments

Sensor response (Ra/Rg) of two screen-coated SnO2 thick film devices different in grain size (12 and 16 nm in diameter) to H2 at 350°C is shown as a function of PH21/2 in Figure 24. The response data of each device fall on a correlation line which is slightly concave upward initially and then becomes almost straight. This behavior is easily understood from (46)–(50). For example, (49) and (50) show that the linear relation is obtained when PH2 is sufficiently large, and otherwise the concaved relation comes out. It is noted that the response as well as the slope of the correlation are lager with the smaller SnO2 grains. This is also consistent with what is expected from the derived equations. It is reserved, however, that really quantitive analysis of the response behavior is not mature yet at present because of a lack of knowledge about many  physicochemical constants needed. Probably we have to go up spirally repeating measurements and theoretical analyses.Figure 24
Observed sensor response to H2 in air at 350°C as correlated with PH21/2 for SnO2 devices different in grain size (screen-printed thick films).Finally let us consider the grain size effects previously shown in Figure4. The response data are reproduced in Figure 25, where the response is plotted against reciprocal of square root of a. Three data of the larger grains to either H2 (A) or CO (B) fall on a straight line passing through origin in agreement with (50), while the remaining data on the smaller grains deviate upward probably through improvements of utility factor due to the appearance of insulating grains. The ratio of the slopes of straight lines, A/B, gives the ratio of the rate constants (k2) of the oxidation reactions of H2 and CO. The ratio is analyzed to be 14. Similar response data to H2 (C) measured with Al-doped SnO2 [42] are also included in the figure. The slope ratio, C/A, which is 1/29 in this case, gives the ratio of donor density (Nd).Figure 25
Grain size effect on sensor response shown in Figure4 as reproduced in the relation between Ra/Rg and a-1/2.


## 5.1. Modeling

Hydrogen is selected here as a representative of reducing gases, which decreases the device resistance. It is unanimous that H2 molecules react with O- adsorbates on the semiconductor crystals. We assume that H2 molecules collide directly with O- ions (Eley-Rideal mechanism) and H2O molecules formed are desorbed instantly from the surface:
R3O-+H2=H2O+e.
This reaction consumes O- ions, while they are resupplied from the gas phase through the reaction, (R1).The rate of accumulation of [O-] is expressed by
(44)d[O-]dt=k1PO2[e]S-k-1[O-]2-k3PH2[O-].
Here t is time, k1 and k-1 are the rate constants of forward and reverse reactions of (R1), respectively, k3 is that of forward reaction of (R3), and PH2 is the partial pressure of H2. At steady state, d[O-]/dt=0, so that we obtain
(45)(KO2PO2)1/2[e]S=[O-]{1+(cPH2[O-])}1/2.
Here KO2=k1/k-1 (equilibrium adsorption constant) and c=k3/k-1. In the absence of H2, this equation coincides with the mass action law of the oxygen adsorption equilibrium, (8). This means that the presence of H2modifies the value of [O-] simply according to (45). The electrostatic equilibrium of the crystal, on the other hand, does not see H2 but [O-], requiring the depletion depth to change correspondingly. Thus the equations relating VS and [O-] with m, n, and p in Section 3 can survive altogether in this new situation. Inserting those into (45) followed by rearrangement, we obtain an equation relating among X(=(KO2PO2)1/2/LD), Y(=(c/LDNd)PH2), and m, while R/R0 is always derived from VS. The sensor resistance measured to a mixture of H2 and oxygen is expressed here as Rg in order to distinguish it from that measured to oxygen only. The sought correlations for plates can be summarized as follows.Regional depletion
(46)Xexp(-m22)=m{1+(1m)Y}1/2,RgR0=exp(m22)(Plate).Volume depletion
(47)RgR0=S(X/n){1+(S/n)Y}1/2,S=1(plate),2(column),and3(sphere).These equations show that under fixedX (or PO2), Rg/R0 is an implicit or explicit function of Y. The equations in the stage of regional depletion for other shapes are excluded here for brevity.


## 5.2. Response to H2 in Air

In usual sensor operation,PO2 is fixed at that (0.21 atm) of the air. Sensor response is often defined as the ratio of the device resistance in air to that in the target gas (or normalized conductance), Ra/Rg. This corresponds to the response ratio in the present treatment as follows:(48)RaRg=(Ra/R0)(Rg/R0).
Each term in the right-hand side can be estimated as a function of Y for fixed X by (46) or (47) so that the conventional sensor response can also be done. For example, Figure 23 shows the sensor response as a function of Y1/2 for plate crystals, where n is selected between 0.1 and 3 and X is fixed at 1000. As easily understood, even when a large X is selected in order to secure deep stages of volume depletion in the absence of H2, the crystals are brought to regional depletion ultimately on increasing Y, and this tendency is enhanced with increasing n. Regional depletion is seen to begin at a relatively small value of Y1/2 for n=3. In the range of Y1/2 where volume depletion prevails, (47) holds so that (48) is transformed into an explicit function of Y: (49)RaRg={1+(1n)Y}1/2={1+(caNd)PH2}1/2.
When (c/aNd)PH2 is far larger than unity,(50)RaRg={Yn}1/2={1+(caNd)PH2}1/2.
This is why linear correlations appear in the range of volume depletion as shown. In the range of regional depletion, such linearity is not allowed, though detailed discussion is reserved here.Figure 23
Theoretical sensor response to H2as correlated with Y1/2 at X=1000 [17].


## 5.3. Comparison with Experiments

Sensor response (Ra/Rg) of two screen-coated SnO2 thick film devices different in grain size (12 and 16 nm in diameter) to H2 at 350°C is shown as a function of PH21/2 in Figure 24. The response data of each device fall on a correlation line which is slightly concave upward initially and then becomes almost straight. This behavior is easily understood from (46)–(50). For example, (49) and (50) show that the linear relation is obtained when PH2 is sufficiently large, and otherwise the concaved relation comes out. It is noted that the response as well as the slope of the correlation are lager with the smaller SnO2 grains. This is also consistent with what is expected from the derived equations. It is reserved, however, that really quantitive analysis of the response behavior is not mature yet at present because of a lack of knowledge about many  physicochemical constants needed. Probably we have to go up spirally repeating measurements and theoretical analyses.Figure 24
Observed sensor response to H2 in air at 350°C as correlated with PH21/2 for SnO2 devices different in grain size (screen-printed thick films).Finally let us consider the grain size effects previously shown in Figure4. The response data are reproduced in Figure 25, where the response is plotted against reciprocal of square root of a. Three data of the larger grains to either H2 (A) or CO (B) fall on a straight line passing through origin in agreement with (50), while the remaining data on the smaller grains deviate upward probably through improvements of utility factor due to the appearance of insulating grains. The ratio of the slopes of straight lines, A/B, gives the ratio of the rate constants (k2) of the oxidation reactions of H2 and CO. The ratio is analyzed to be 14. Similar response data to H2 (C) measured with Al-doped SnO2 [42] are also included in the figure. The slope ratio, C/A, which is 1/29 in this case, gives the ratio of donor density (Nd).Figure 25
Grain size effect on sensor response shown in Figure4 as reproduced in the relation between Ra/Rg and a-1/2.


## 6. Influences of Electron-Traps Dispersed on the Crystals Surface (Plates)

### 6.1. Modeling of Receptor Function

Semiconductor gas sensors are usually sensitized by loading the constituent crystals with sensitizers, as mentioned before. Among a variety of sensitizers known, we consider here only those which draw electrons from the constituent crystals, increasing work function of the crystals. Further, only the event of electron transfer between the sensitizers and the crystals is focused attention, excluding any other events such as redox changes of the sensitizers. It is assumed that the sensitizer in problem is deposited as tiny particles or clusters (denotedD), which form electron-trapping centers located in the forbidden band at the surface. Electrons would be transferred from the crystals to form ionized traps (D-). Unlike the case of gas adsorption, however, the electrons once trapped do not return back to the crystals no matter how gas atmosphere changes, if the trapping level is deep enough. In the treatment of electrostatic equilibrium, D- contributes to the total surface charge density (Qs) together with ionized adsorbates, typically O- and/or NO2-, as shown in Figure 26. For plate crystals placed in air, this is expressed as(51)Qs=-q([D-]+[O-])=-qmLDNd.
Here m is the reduced depletion depth in air. Obviously the presence of [D-] affects m depending on its magnitude. For ease of treatment, [D-] is expressed in terms of reduced depletion depth, mD, it develops in the absence of oxygen. That is,(52)[D-]=mDLDNd.
Two cases are distinguished here depending on whether the electron transfer between the traps and the bulk of crystals is reversible (reversible traps) or not (irreversible traps).Figure 26
Diagram of electrostatic equilibrium drawn for small crystals loaded with surface electron-traps.In the case of reversible traps, trapped electrons participate in the electrostatic equilibrium of the crystals; [D-], [O-], and the bulk of each crystal are kept at the same Fermi level. By using (8) and (52), (51) is rewritten as(53)mDLDNd+(KO2PO2)1/2[e]s=mLDNd.
It is reminded that m and mD are reduced depletion depths to [D-]+[O-] and [D-], respectively. We know that [e]s as well as R/R0 are expressed as a function of m, n, and p. The following equations result for correlating R/R0 with X (abrupt model).Regional depletion
(54)Xexp(-m22)=m-mD,RR0=exp(m22).Volume depletion
(55)RR0=exp{(n22)+p}=X(n-mD).The correlations derived formD=0,1,2, and 2.5, and n=3 are shown in Figure 27. For given X, R/R0 increases dramatically with increasing mD. In the stage of volume depletion, R/R0 of a trap-dispersed plate of reduced size n is the same as that exhibited by a clean plate of reduced size n-mD. The response in the stage of regional depletion is seen to be also enhanced dramatically on increasing mD.Figure 27
Correlations betweenR/R0 and X for plate crystals dispersed with reversible electron-traps [17].In the case of irreversible traps, trapped electrons are assumed to stay at the traps permanently. It is suspected that such case may happen when the foreign deposits react with the surface of the constituent crystals. For example, the component cations of the deposits are dissolved into the lattice of underlying crystals at the surface or nearby. The dissolved cations may act as strong acceptors of electrons (irreversible traps). Anyway those trapped electrons amounting to depthmD are kept off from the electrostatic equilibrium. As a result, oxygen adsorption cannot take place until PO2 exceeds a critical point (PO2(D)) at which depletion depth reaches mD in the absence of [D-]. Equation (51) is rewritten as(56)[D-]+{(KO2PO2)1/2-(KO2PO2(D))1/2}[e]s=mLDNd.
The resulting correlating equations are summarized as follows:Regional depletion
(57)(X-XD)exp(-m22)=m-mD,RR0=exp(m22).Volume depletion
(58)(X-XD)exp{-(n22)-p}=n-mD,RR0=exp{(n22)+p}=(X-XD)(n-mD).HereX(D)=(KO2PO2(D))1/2/LD. The correlations given by these equations are shown in Figure 28, where mD and n are set to the same values as those used in Figure 27. It is seen that the response does not appear up to X=XD (masking) and then rise up steeply with a slope increasing with increasing mD. The irreversible traps are thus effective also in enhancing the sensitivity to O2, though in a manner different from that with the reversible traps.Figure 28
Correlations betweenR/R0 and X for plate crystals dispersed with irreversible electron-traps [17].


### 6.2. Response to NO2 or H2 in Air

In this section, response of electron-trap dispersed plates to NO2 or H2 in air is discussed briefly. Under exposure to NO2 in air, total surface charge density is written as follows:(59)Qs=-q([D-]+[O-]+[NO2-])=-qmLDNd.
For either type of traps, reversible and irreversible, the sought correlations between R/R0 and X+Z are essentially the same as those obtained between R/R0 and X if X is replaced by X+Z. At fixed X, the sensitivity of R/R0 to Z increases sharply with increasing mD.The response to H2 in air can be developed in the same way as done for clean plate crystals. [O-] is given by (45), while [D-] is kept constant. Using these equations, it is possible to derive the sought correlations for reversible traps as follows.Regional depletion
(60)Xexp(-m22)=(m-mD){1+[1(m-mD)]Y}1/2,RgR0=exp(m22).Volume depletion
(61)RgR0={X(n-mD)}{1+[1(n-mD)]Y}1/2,(62)RaRg={1+[1(n-mD)]Y}1/2={1+[1(n-mD)][cLDNd]PH2}1/2.
Equation (62) formulates conventional sensor response. The response is drastically promoted by increasing mD, as illustrated in Figure 29, where correlations are sought between Ra/Rg and Y1/2 with n and X being fixed at 3 and 1000, respectively. However, the manner of promotion of the response is considerably different from that of the thickness effect observed with clean plates (Figure 23).Figure 29
Correlations betweenRa/Rg and Y1/2 for plate crystals dispersed with reversible electron-traps [17].In the case of irreversible traps, the correlations are given as follows.Regional depletion
(63)(X-XD)exp(-m22)=(m-mD){1+[1(m-mD)]Y}1/2,RgR0=exp(m22).Volume depletion
(64)RgR0={(X-XD)(n-mD)}{1+[1(n-mD)]Y}-1/2,(65)RaRg={1+[1(n-mD)][cLDNd]PH2}1/2.
Expression of conventional sensor response (65) is seen to be identical to (62). This is quite natural because it does not matter in the stage of volume depletion whether the traps are reversible or not.


## 6.1. Modeling of Receptor Function

Semiconductor gas sensors are usually sensitized by loading the constituent crystals with sensitizers, as mentioned before. Among a variety of sensitizers known, we consider here only those which draw electrons from the constituent crystals, increasing work function of the crystals. Further, only the event of electron transfer between the sensitizers and the crystals is focused attention, excluding any other events such as redox changes of the sensitizers. It is assumed that the sensitizer in problem is deposited as tiny particles or clusters (denotedD), which form electron-trapping centers located in the forbidden band at the surface. Electrons would be transferred from the crystals to form ionized traps (D-). Unlike the case of gas adsorption, however, the electrons once trapped do not return back to the crystals no matter how gas atmosphere changes, if the trapping level is deep enough. In the treatment of electrostatic equilibrium, D- contributes to the total surface charge density (Qs) together with ionized adsorbates, typically O- and/or NO2-, as shown in Figure 26. For plate crystals placed in air, this is expressed as(51)Qs=-q([D-]+[O-])=-qmLDNd.
Here m is the reduced depletion depth in air. Obviously the presence of [D-] affects m depending on its magnitude. For ease of treatment, [D-] is expressed in terms of reduced depletion depth, mD, it develops in the absence of oxygen. That is,(52)[D-]=mDLDNd.
Two cases are distinguished here depending on whether the electron transfer between the traps and the bulk of crystals is reversible (reversible traps) or not (irreversible traps).Figure 26
Diagram of electrostatic equilibrium drawn for small crystals loaded with surface electron-traps.In the case of reversible traps, trapped electrons participate in the electrostatic equilibrium of the crystals; [D-], [O-], and the bulk of each crystal are kept at the same Fermi level. By using (8) and (52), (51) is rewritten as(53)mDLDNd+(KO2PO2)1/2[e]s=mLDNd.
It is reminded that m and mD are reduced depletion depths to [D-]+[O-] and [D-], respectively. We know that [e]s as well as R/R0 are expressed as a function of m, n, and p. The following equations result for correlating R/R0 with X (abrupt model).Regional depletion
(54)Xexp(-m22)=m-mD,RR0=exp(m22).Volume depletion
(55)RR0=exp{(n22)+p}=X(n-mD).The correlations derived formD=0,1,2, and 2.5, and n=3 are shown in Figure 27. For given X, R/R0 increases dramatically with increasing mD. In the stage of volume depletion, R/R0 of a trap-dispersed plate of reduced size n is the same as that exhibited by a clean plate of reduced size n-mD. The response in the stage of regional depletion is seen to be also enhanced dramatically on increasing mD.Figure 27
Correlations betweenR/R0 and X for plate crystals dispersed with reversible electron-traps [17].In the case of irreversible traps, trapped electrons are assumed to stay at the traps permanently. It is suspected that such case may happen when the foreign deposits react with the surface of the constituent crystals. For example, the component cations of the deposits are dissolved into the lattice of underlying crystals at the surface or nearby. The dissolved cations may act as strong acceptors of electrons (irreversible traps). Anyway those trapped electrons amounting to depthmD are kept off from the electrostatic equilibrium. As a result, oxygen adsorption cannot take place until PO2 exceeds a critical point (PO2(D)) at which depletion depth reaches mD in the absence of [D-]. Equation (51) is rewritten as(56)[D-]+{(KO2PO2)1/2-(KO2PO2(D))1/2}[e]s=mLDNd.
The resulting correlating equations are summarized as follows:Regional depletion
(57)(X-XD)exp(-m22)=m-mD,RR0=exp(m22).Volume depletion
(58)(X-XD)exp{-(n22)-p}=n-mD,RR0=exp{(n22)+p}=(X-XD)(n-mD).HereX(D)=(KO2PO2(D))1/2/LD. The correlations given by these equations are shown in Figure 28, where mD and n are set to the same values as those used in Figure 27. It is seen that the response does not appear up to X=XD (masking) and then rise up steeply with a slope increasing with increasing mD. The irreversible traps are thus effective also in enhancing the sensitivity to O2, though in a manner different from that with the reversible traps.Figure 28
Correlations betweenR/R0 and X for plate crystals dispersed with irreversible electron-traps [17].


## 6.2. Response to NO2 or H2 in Air

In this section, response of electron-trap dispersed plates to NO2 or H2 in air is discussed briefly. Under exposure to NO2 in air, total surface charge density is written as follows:(59)Qs=-q([D-]+[O-]+[NO2-])=-qmLDNd.
For either type of traps, reversible and irreversible, the sought correlations between R/R0 and X+Z are essentially the same as those obtained between R/R0 and X if X is replaced by X+Z. At fixed X, the sensitivity of R/R0 to Z increases sharply with increasing mD.The response to H2 in air can be developed in the same way as done for clean plate crystals. [O-] is given by (45), while [D-] is kept constant. Using these equations, it is possible to derive the sought correlations for reversible traps as follows.Regional depletion
(60)Xexp(-m22)=(m-mD){1+[1(m-mD)]Y}1/2,RgR0=exp(m22).Volume depletion
(61)RgR0={X(n-mD)}{1+[1(n-mD)]Y}1/2,(62)RaRg={1+[1(n-mD)]Y}1/2={1+[1(n-mD)][cLDNd]PH2}1/2.
Equation (62) formulates conventional sensor response. The response is drastically promoted by increasing mD, as illustrated in Figure 29, where correlations are sought between Ra/Rg and Y1/2 with n and X being fixed at 3 and 1000, respectively. However, the manner of promotion of the response is considerably different from that of the thickness effect observed with clean plates (Figure 23).Figure 29
Correlations betweenRa/Rg and Y1/2 for plate crystals dispersed with reversible electron-traps [17].In the case of irreversible traps, the correlations are given as follows.Regional depletion
(63)(X-XD)exp(-m22)=(m-mD){1+[1(m-mD)]Y}1/2,RgR0=exp(m22).Volume depletion
(64)RgR0={(X-XD)(n-mD)}{1+[1(n-mD)]Y}-1/2,(65)RaRg={1+[1(n-mD)][cLDNd]PH2}1/2.
Expression of conventional sensor response (65) is seen to be identical to (62). This is quite natural because it does not matter in the stage of volume depletion whether the traps are reversible or not.


## 7. Discussion

To those who have been engaged in semiconductor gas sensor researches like us, it is a long-time dream to acquire a theoretical background of sensor researches. For a long time, many people have exerted efforts to upgrade capabilities of sensors or to extend applications of sensors empirically from various standpoints of sensing materials, material preparation and processing, sensitizers and stabilizers, sensing body structure and fabrication, device structure, electrodes, signal processing, transients based sensing, and so forth. These efforts have contributed to ever continuing progresses of gas sensors. However, it has not always been easy to understand well the meaning of a new progress achieved in a particular category, to those working in different categories. Even it has happened that the finders themselves do not understand their remarkable findings well, failing to connect them to another progress. All these have arisen mainly from a lack of a theoretical background commonly possessed by gas sensor researchers.The theoretical concept available to semiconductor gas sensors has so far been almost nothing more than the concept of semiconductor catalysis born a long time ago. It has been useful for qualitative interpretations of gas sensors but almost powerless for quantitative analyses. A main reason for this is that the concept has ignored the role or meaning of so small semiconductors as used in gas sensors. Behavior of small semiconductors is considerably different from that of large ones and this difference provides gas sensors with profound merits. This is a truth unrevealed so far. We have shown in this article that gas response behavior of small semiconductors can be formulated by using well-established concepts and well-known parameters in physics and chemistry. The equations derived seem to work well on understanding the response behavior experimentally observed, though available data are still limited at present. No doubt the equations need to be checked more rigorously by experiments. Nevertheless we believe we are at a new starting point for theoretical approaches to semiconductor gas sensors. It is anticipated that such approaches will contribute much to strengthening and expanding the background of gas sensors. The treatments described in this article have focused attention to receptor function of semiconductors only, setting transducer function to work in a simple mode. Elaboration of the transducer function is one of the most important subjects of theoretical approaches.


## 8. Conclusions

Through the theoretical approaches to semiconductor gas sensors described here, the following conclusions can be drawn.(1) Receptor function and response of small semiconductor crystals can be formulated in principle by using the chemical parameters of gases side, such as partial pressure, adsorption constant and rate constant, and the physical parameters of semiconductor side, such as shape and size, donor density, and Debye length.(2) Theoretical equations account for well-response behavior to oxygen, nitrogen dioxide, and hydrogen, as well as the grain size effects involved, though the tested examples are limited in number.(3) Electron-traps dispersed on the crystal surface influence strongly on receptor function and response to gases, thus exerting a kind of sensitization effect.



---

*Source: 875704-2009-07-12.xml*