Characterization and Modeling of a Pt-In2O3 Resistive Sensor for Hydrogen Detection at Room Temperature

Sensitive H2 sensors at low concentrations and room temperature are desired for the early warning and control of hydrogen leakage. In this paper, a resistive sensor based on Pt-doped In2O3 nanoparticles was fabricated using inkjet printing process. The H2 sensing performance of the sensor was evaluated at low concentrations below 1% at room temperature. It exhibited a relative high response of 42.34% to 0.6% H2. As the relative humidity of 0.5% H2 decreased from 34% to 23%, the response decreased slightly from 34% to 23%. The sensing principle and the humidity effect were discussed. A dynamic current sensing model for dry H2 detection was proposed based on Wolkenstein theory and experimentally verified to be able to predict the sensing behavior of the sensor. The H2 concentration can be calculated within a short measurement time using the model without waiting for the saturation of the response, which significantly reduces the sensing and recovery time of the sensor. The sensor is expected to be a promising candidate for room-temperature H2 detection, and the proposed model could be very helpful in promoting the application of the sensor for real-time H2 leakage monitoring.


Introduction
Hydrogen fuel is regarded as a vital green energy to support the operation of human society in the future due to its non-polluting combustion and high energy conversion rate [1,2]. However, the flammable and explosive properties of H 2 and the broad explosive limit concentration range of 4-75% are serious safety concerns in the development of hydrogen energy [3]. Since hydrogen gas (H 2 ) is colorless and odorless, it is challenging to detect it in time, especially at low concentrations and room temperature [4,5]. Therefore, it is of great significance to develop sensors that can detect H 2 sensitively and reliably at room temperature for the early warning and control of H 2 leakage.
Resistive H 2 sensors based on metal oxide semiconductors (MOS) have been widely reported because of their simple fabrication, low cost, and good selectivity. At present, many MOS materials, such as WO 3 [6], In 2 O 3 [7], TiO 2 [8], ZnO [9], and SnO 2 [10], have been used in the detection of H 2 . Among them, In 2 O 3 , an n-type semiconductor with a wide-bandgap (~3.6 eV), shows great potential due to its low cost, high electrical conductivity, high stability, and excellent H 2 sensitivity [11]. Various forms of low-dimensional In 2 O 3 have been studied to meet the requirements of room temperature and fast and stable H 2 detection [12,13]. The nanostructure enables the MOS materials to have a high surface/volume ratio and an enhanced gas adsorption capacity [14][15][16]. Additionally, doping noble metals, such as Au, Ag, Pd, Pt, etc., was used as another effective method to improve the sensing performance of materials, since the catalysis of noble metals can promote the sensing reactions [17][18][19][20][21].
In this paper, In 2 O 3 nanoparticles with a diameter of about 100 nm doped by 10-weight-concentration (wt.%) Pt were used as the sensitive material for the detection of H 2 at room temperature. Platinum was adopted as the dopant because it has a lower hydrogen-solubility compared with other commonly used inert metals, which is more beneficial to improving the stability of the sensor [20,22]. A precursor solution containing In 2 O 3 nanopowders and H 2 PtCl 6 was deposited by inkjet printing technology on the integrated electrodes and was heated to form the sensing layer. The adopted inkjet printing technology makes the fabrication process of the sensor simple and convenient with good repeatability and mass production capability. Besides, the heating temperature is controlled at 300 • C, which will not affect other semiconductor devices when the sensor is integrated with other functional circuits on the same silicon wafer for forming a sensing system. Sensing behaviors of the proposed sensor for both dry and humid H 2 were characterized in this work, which presented a reasonable H 2 sensing properties at room temperature. Moreover, the resistive sensor with nano Pt-In 2 O 3 was also modeled and analyzed in order to deeply understand the H 2 sensing principle at room temperature and help to predict the sensing results of the sensor. The method mainly relies on analytical descriptions of the adsorption and desorption behavior (Wolkenstein and Langmuir theories) and energy band theory of semiconductors. A model of transient sensing characteristics of the proposed sensor was built to predict the dynamic change in the current response to H 2 of the sensor and validated experimentally. It also contributed to explain the humidity effect on H 2 sensitivity of the sensor. In this work, the Pt-In 2 O 3 sensor could be a promising candidate for room-temperature H 2 detection, and the proposed model could be very helpful in promoting the application of the sensor for real-time H 2 leakage monitoring and the establishment of complex gas sensing systems by integrating with circuits.

Materials and Methods
The resistive sensor platform was fabricated on a silicon wafer with an O/N/O passivation layer consisting of SiO 2 (10 nm)/Si 3 N 4 (20 nm)/SiO 2 (10 nm). The Cr (30 nm)/Au (50 nm) multilayer was deposited and patterned to serve as the interdigitated electrodes. The width of the electrodes and the space between them are 2 µm and 1.8 µm, respectively. Then, an SU-8 passivation layer formed by spin coating was patterned on top of the wafer by a lithography process to expose only the interdigitated-electrode area and the contact pads of the of resistor platform, which can protect other devices on the wafer.
The Pt-In 2 O 3 sensitive material was formed by an inkjet printing process. Firstly, In 2 O 3 nanopowders (about 100 nm in diameter) were dissolved in ethanol (99%) and stirred thoroughly. Next, an 8 wt.% H 2 PtCl 6 (in H 2 O) solution was further diluted by deionized water to the desired concentration and mixed with the In 2 O 3 solution together to serve as the precursor ink. The as-prepared ink was printed on the interdigitated-electrode area of the platform using Omni Jet 100 inkjet printer, followed by a 2-h annealing process at 300 • C in air to fully evaporate the solvent and decompose H 2 PtCl 6 into Pt [23]. The wt.% of Pt in the sensing layer was set to be 10 wt.% to focus principally on the analysis of hydrogen adsorption and effects of moisture on the sensing properties. The resistance of the obtained sensor is about 7.3 × 10 5 kΩ. All chemicals used in this paper were purchased from Sigma-Aldrich without further purification.
The experimental setup is shown in Figure 1. In this paper, the prepared sensor was characterized for both dry and humid H 2 detection. During the test, the sensor was placed in the test chamber, and the total flow rate of the gas flowing into the chamber was always maintained at 400 sccm. In the tests of dry H 2 detection, MFC 1 in Figure 1 was off. In the initial stage, dry air served as reference gas (400 sccm) and was controlled by MFC 2 to blow into the chamber, and the sensor was in a standby state. When starting to sense H 2 , the reference gas was switched to the sample gas with a fixed concentration of H 2 (400 sccm) injected into the test chamber by using the switcher in Figure 1. At the same time, the reference gas was converted to the outlet and discharged into the atmosphere. Hydrogen gas samples were prepared by diluting 1% H 2 with a dry carrier gas (air) to desired concentrations using MFC 3 and MFC 4.
In the tests of humid H 2 detection, MFC 1 was turned on to let another stream of dry air flow into the bubbler and then mixed with the reference dry air or the dry H 2 samples to bring water vapor into the test chamber. The total flow rate of the gas injected into the chamber was still fixed at 400 sccm, and the relative humidity (RH) was monitored by the humidity meter in Figure 1 and tuned by MFC 1 to a desired level. Figure 1 also presented SEM images of the sensor with a transparent Pt-In 2 O 3 layer and Pt-In 2 O 3 layer nanoparticles. A clear magnification of the SEM image of Pt-In 2 O 3 particles in Figure 1 and EDS analysis are provided in Figure S1 of the Supplementary Materials which confirms the formation of Pt-In 2 O 3 . During all sensing tests, the voltage between the two electrodes of the sensor was fixed to measure the changes in current that corresponded to the variation in the H 2 concentration. All sensing characteristics of all sensors were tested at 25 • C (room temperature). Electrical measurements were carried out by using an Agilent B1500A analyzer.

Experimental Results
Firstly, a comparison between the resistive sensors with pure In 2 O 3 and 10 wt.% Pt-In 2 O 3 for the detection of 0.3% H 2 was conducted as shown in Figure 2a, which indicates that Pt was well-formed on the surface of In 2 O 3 and improved the H 2 sensitivity of the sensor. Then, dry H 2 sensing performance of the Pt-In 2 O 3 resistive sensor was measured. The results are shown in Figure 2b,c. In the initial stage of tests, 400 sccm of dry air was blown towards the sensor. From 20 s in Figure 2b, the air was replaced by a dry H 2 sample with a certain concentration and maintained for 100 s for the sensor to respond. After that (at 120 s in Figure 2b), the gas flowing into the chamber was changed back to the reference gas, and the sensor entered the recovery period. The above experimental procedure was repeated six times, and the H 2 concentration used each time was set from 0.1% to 0.6%, respectively, with a step of 0.1%. The responses, calculated by using Equation (1), are 3.04%, 9.25%, 17.60%, 26.07%, 33.69%, and 42.34% for 0.1-0.6% of H 2 , respectively, based on the data in Figure 2a.
where I H and I air denote the current of the sensor corresponding to H 2 samples and reference gas, respectively. The response and recovery times of the sensor were plotted as functions of H 2 concentration in Supplementary Figure S2. According to the results, the response of the sensor was very limited for the 0.1% H 2 . With the increasing H 2 concentration, the rising speed of the sensor current increased significantly. Figure 2c plots the response of the Pt-In 2 O 3 sensor as a function of H 2 concentration, showing a near-linear relationship. Each error bar corresponds to one standard deviation, considering the uncertainty of three response values. The detection limit (DL) was also calculated here, as it is a crucial indicator reflecting the sensing performance for practical application [24]. First, the standard deviation in the baseline of the response curve (rms noise ) was calculated from 20 data points, which is 0.003828, representing the level of the sensor noise. The slop of the linear fitting of gas response versus gas concentration at RT is 75.81454 according to Figure 2c. Therefore, the DL = 3rms noise /slop = 1.51 ppm according to literature [24,25]. Additionally, a comparison between the proposed sensor of this work and recent H 2 sensors based on various sensitive materials was conducted as shown in Table 1 [6,[26][27][28][29][30][31][32][33][34].
As the aim of this work is the detection of H 2 at room temperature, it is necessary to consider the effect of humidity on the sensor performance. Therefore, humid H 2 sensing tests were also conducted and the measurement results are shown in Figure 3. The experimental setup and gas controls can be found in Section 2. Figure 3a shows the dynamic response of the sensor towards 0.5% H 2 at 3.3% RH (corresponding to dry H 2 ), 9.4% RH, and 15.5% RH, respectively. The humid H 2 samples were maintained for 100s (from 120 s to 220 s in Figure 3a), and the change in current at 220 s was slightly reduced by the increasing RH. Figure 3b shows the variation of the response with RH for 0.5% H 2 . As the RH increased from 3.3% to 15.5%, the response decreases from 34% to 23%. The results illustrate that water vapor deteriorates the hydrogen sensitivity of the device. The rationale behind this was explored in-depth in the next section. Table 1. Comparison of sensing performances towards hydrogen gas of sensors with various materials [6,[26][27][28][29][30][31][32][33][34]. RT: room temperature, NW: nanowire, MOX: metal oxide, NF: nanofiber, R a : resistance in the air, R g : resistance in detection gas, -: not available. Adapted with permission from Ref. [27].

Sensor Materials T ( • C) Detection Range Response/Recovery Time (Response) at H 2 % Advantage/Disadvantage
Ref.

Current Model
In this section, by analyzing the sensing principle of the prepared sensor, a numerical sensing model was proposed based on adsorption kinetics and energy band theory and validated by experimental data. Finally, the proposed modeling approach was used to explain the mechanism of the humidity effect on the H 2 sensing performance of the inkjetprinted Pt-In 2 O 3 sensor.
For the sake of illustration, the Pt-In 2 O 3 resistive sensor was simplified to the structure shown in Figure 4a, i.e., two metal electrodes are connected by spherical In 2 O 3 nanoparticles embedded with Pt on the surface. In the figure, R c is the contact resistance between the electrode and the sensitive material, R s is the surface resistance at the interface between In 2 O 3 particles, R b and R d are the resistance of the bulk and surface depletion region of an In 2 O 3 particle, respectively; E F is the Fermi energy level, E b C and E s C are the conduction bands of the bulk and surface of the In 2 O 3 particle, respectively. At room temperature, the catalytic effect of Pt causes the adsorption of oxygen on the surface of In 2 O 3 particles in the form of O − x , as shown in Figure 4c (such as O − and O − 2 at low temperature) [35,36]. During this process, electrons in In 2 O 3 near the surface transfer to the adsorbed oxygen species to form a space charge region (depletion of electrons) according to a possible reaction shown as Equation (2) [37], which induces upward bending of the energy bands at the surface denoted as qV s (V s < 0) .
The average diameter of the In 2 O 3 particles used in this paper is 100 nm, which is larger than the Debye length, so it can be considered that the particles are in a partially depleted state. That is, there is a low-resistance electron accumulation region in the body of the particle, and only the surface is depleted by the chemisorption of oxygen. The thickness of the depletion region changes accordingly with the type and concentration of the gas molecules adsorbed on the surface and affect the overall resistance of the particle. Based on these analyses, the resistance of an In 2 O 3 particle can be simplified to the composition shown in Figure 4b. In the middle part of the particle, the resistances of the bulk and the surface depletion region are in parallel, and they connect with a surface resistance in series on each side of the particle, i.e., As In 2 O 3 is an n-type semiconductor, the total resistance of a particle can be written as Equation (4) according to literature [38].
Note that the contribution of R c to overall resistance of the sensor and the change of it by gas adsorption can be ignored [38,39]. Thus, the entire sensor can be approximately considered to be composed of numerous R p s in series and parallel, i.e., where γ is the coefficient related to the morphology of the sensitive material, the contact area between the particles, and the geometry chosen for simulating the particle (such as sphere, cylinder, and cube, etc.). Thus, the current can be written as where, β = U/(γ ·R b ), and U is the applied external DC voltage. Figure 4c shows the schematic diagram of the H 2 sensing mechanism of Pt-In 2 O 3 . When the sensor is exposed to H 2 , the H 2 molecules can be dissociated into H atoms by the catalytic effect of Pt even at room temperature according to a possible reaction shown as Equation (7) [22], which corresponds to process 1 in the figure.
The difference in H concentration between the surfaces of Pt and In 2 O 3 leads to the diffusion of H atoms from Pt to In 2 O 3 , which is so-called spillover corresponding to process 2 . Process 3 includes two reactions. On one hand, before diffusing to In 2 O 3 , hydrogen atoms can react with the adsorbed oxygen species which have already existed on Pt in the standby stage of the sensor, thereby generating water vapor leaving the surface of the sensitive material. On the other hand, as the H 2 gas sample is a mixture of 1 % H 2 and air, it contains both H 2 and O 2 , which can also react directly under the catalysis of Pt to produce water molecules when both gases meet with each other on the surface of Pt. When H atoms diffuse to In 2 O 3 , they react with the negatively charged oxygen species to generate water and release free electrons to In 2 O 3 , which is the reduction reaction that plays a major role in H 2 detection and denoted by 4 in Figure 4c. As a result, H 2 reduces the thickness of the depletion layer on the surface of In 2 O 3 , thereby reducing the |Vs| and increasing the current. where R c , R s , R b , and R d are the contact resistance between the metal electrode and the sensitive material, the surface resistance at the interface between In 2 O 3 particles, the bulk resistance inside the particle, and the resistance of the depletion region at the periphery of the particle, respectively; E F is the Fermi energy level, and E b C and E s C are the conduction bands of the bulk and surface of the In 2 O 3 particle, respectively. When the sensor is placed in an air atmosphere, O 2 naturally adsorbs on the surface of the In 2 O 3 , forming an upward bending of the surface energy band denoted as qV s (V s < 0); (b) Composition of the resistance of an In 2 O 3 particle after simplification. (c) Schematic diagram illustrating the H 2 sensing mechanism of Pt-In 2 O 3 . Here, 1 , 2 , 3 , and 4 represent the dissociative adsorption process of hydrogen at Pt, the spillover process, the direct reaction with adsorbed oxygen species on Pt to produce water vapor, and the reaction with adsorbed oxygen species on the surface of In 2 O 3 to produce water vapor and release electrons, respectively; O x and O − x represent the possible oxygen species adsorbed on Pt and In 2 O 3 .

Hydrogen Dynamic Sensing Model under Dry Conditions
There are two commonly used theories for the analytical descriptions of the adsorption and desorption behavior, which are Langmuir and Wolkenstein theories, respectively. The Langmuir model assumes that the bonding energy between the adsorbent and the adsorbed material is a fixed constant. However, both the oxygen adsorption and the H 2 detection involved in this paper are related to chemical reactions, and the processes include electron transfer between the adsorbent and the sensitive material. In this case, the binding energy or adsorption heat will vary with the coverage of the chemisorbed material, so that the Langmuir model cannot meet the requirements for establishing the dynamic adsorption model here. The Wolkenstein model, by contrast, takes into account both the electron interactions and their effects on the adsorption properties of the semiconductors. Therefore, the Wolkenstein model was adopted here. Figure 5 shows the energy band distributions and the corresponding internal and surface electron distributions of In 2 O 3 particles based on the Wolfenstein theory for four cases of pure surface (a), weak oxygen adsorption, (b) strong oxygen adsorption at equilibrium state (c), and hydrogen reduction (d). It is worth noting that Pt is not marked in the figure. The doped Pt and In 2 O 3 were studied as a whole in this section.
The core of Wolkenstein's model is to divide chemisorption of gases into two steps: weak chemisorption and strong chemisorption. In the ideal situation, the energy bands of the In 2 O 3 particle with an empty surface are assumed to be flat, as shown in Figure 5a. By regarding the In 2 O 3 particles as spheres, a coordinate axis x was established along the surface to the center of the sphere. Next, the energy band distribution was obtained at the bottom of Figure 5a, where r is the average radius of the particles. The energy bands of the particles are consistent from the inside to the outside, and the electrons are evenly distributed throughout the sphere. E C , E F , and E V represent the conduction band, Fermi energy level, and valence band, respectively.
When the sensor is exposed to a dry air, oxygen will be dissociated into O or adsorbed molecularly in the presence of Pt and diffuse to the surface of In 2 O 3 at low temperature [35,36,40]. As shown in Figure 4b, when O ∞ x reaches the indium oxide surface, a weak chemisorption without electron transfer is formed firstly, so that the energy bands are still flat and the concentration of electrons throughout the particle is still the same. During this process, the state of O species changes from O ∞ x to O 0 x , and the adsorption energy is q 0 . At the same time, an additional surface state (E ss ) is also introduced in the bandgap of In 2 O 3 .
Since the Fermi energy level is higher than E ss , the electrons in the conduction band bottom of In 2 O 3 oxide are taken away by O 0 x and the weakly adsorbed O 0 x is further strongly chemisorbed to form O − x ions under the room-temperature test condition of this paper [35]. Thus, at equilibrium, the energy bands bend upward, as shown in Figure 5c. At this point, the conduction band and valence band at the surface (E s C and E s V ) are no longer the same as the ones in the bulk (E b C and E b V ). The band bending is recorded as qV s0 (V s0 < 0), which leads a depletion of electrons near the surface of In 2 O 3 with a thickness of x d0 . The total adsorption energy of O dissociated by Pt on the In 2 O 3 is q 0 + E s C − E ss . When the sensor starts to detect H 2 , H 2 molecules will dissociate into H atoms under the catalysis of Pt, come to the surface of In 2 O 3 through the spillover process, and react with O − x species, thereby returning electrons to the conduction band of In 2 O 3 . Therefore, the upward bending degree of the energy bands is reduced from qV s0 to qV s , and the thickness of the depletion region is reduced from x d0 to x d , as shown in Figure 5d. If ∆H is denoted as the difference between the oxygen adsorption energy and the hydrogen reaction energy, the hydrogen reaction energy can be expressed as q 0 + E s C − E ss − ∆H. It is assumed that the O − species on the surface of In 2 O 3 plays the predominant role during the sensing process [35] and the concentration is N − . In the standby state of the sensor, the adsorption of O − is in equilibrium. When H 2 appears, N − is changed by the desorption of O − and the reduction reaction with hydrogen with time, which can be written as  The superscripts s and b represent the parameters of the surface and the core of an In2O3 particle, respectively. The binding energy of weak chemisorption is denoted as q0. O x ∞ , O x 0 , and O x are the oxygen species before adsorption, weakly adsorbed, and strongly adsorbed, respectively. ΔH represents the difference between the adsorption energy of oxygen species and the chemical energy of the reaction between the sensitive material and H2, r is the average radius of indium oxide particles, and xd0 and xd are the thickness of the depletion layer at oxygen adsorption equilibrium and after the reaction with hydrogen, respectively. after thermal equilibrium and electron depletion near the surface. (d) The particle undergoing a reduction reaction after the presence of hydrogen. E C , E F , and E V represent the conduction band, Fermi energy level, and valence band, respectively. The superscripts s and b represent the parameters of the surface and the core of an In 2 O 3 particle, respectively. The binding energy of weak chemisorption is denoted as q 0 . O ∞ x , O 0 x , and O − x are the oxygen species before adsorption, weakly adsorbed, and strongly adsorbed, respectively. ∆H represents the difference between the adsorption energy of oxygen species and the chemical energy of the reaction between the sensitive material and H 2 , r is the average radius of indium oxide particles, and x d0 and x d are the thickness of the depletion layer at oxygen adsorption equilibrium and after the reaction with hydrogen, respectively.
Assuming that the reduction reaction with hydrogen is the dominant factor effecting the concentration of O − on the surface of In 2 O 3 , then Equation (9) can be obtained as below according to Wolkenstein theory [41], where α 1 is the rate constant of the reaction, k is the Boltzmann constant, T is the thermodynamic temperature, and P H is the partial pressure of H 2 . As it is known that E b C = E s C + qV s (0 > V s > V s0 ), and P H ∝ C H , where C H is the concentration of H 2 , then Equation (9) can be rewritten as where α 1 is the coefficient related to P H and C H . Then, Equation (10) is rewritten as According to the full depletion approximation of the space charge region in semiconductor physics, i.e., there are no mobile carriers in the space charge region, the thickness of the depletion layer and the corresponding N − at a certain moment are obtained as follows where n b is the concentration of electrons inside the In 2 O 3 nanoparticle, and ε is the dielectric constant of the sensitive material. Substituting Equation (14) into 11 yields According to Equation (14), one can obtain: Thus, the change of N − with time becomes: As Equations (11) and (17) are equal, with u = qV s /kT, we obtain Define u 0 and u s to correspond to qV s /kT before H 2 detection and qV s0 /kT after detecting H 2 gas for time t, respectively. Integrate both sides of Equation (18) in the time range of 0-t to obtain Equation (19).
Let E i (x) = x −∞ e t /tdt, a = 2m and b = Ei(u 0 ), then Equation (19) and the expression of current become E i (u s ) = a·C H ·t + b (20) To simplify the calculation, the complex integral of Ei(x) was approximated as e x /x. Next, the current response data in Figure 2a from 60 s to 70 s were used for fitting to calculate a, b, and β, which can evade the errors caused by unstable experimental conditions, such as the changes of gases injected into the chamber at the beginning and the end of the response period. The final fitting results are shown in Figure 6, where a = −0.10094, b = −0.42957, and β = 1.55 × 10 −8 .
Thus, the final numerical model of hydrogen dynamic sensing under dry conditions is obtained as Equations (22) and (23).
To verify the model, the current values corresponding to the concentrations of each H 2 sample at 71 s and 72 s in Figure 2a were substituted into the model. The calculated H 2 concentrations were also plotted in Figure 6, as shown in the red hollow circles (71 s) and blue hollow triangles (72 s). The computation results using the data at 71 s and 72 s are in good agreement with the numerical model curve, except for a small error at 0.6% of H 2 .  (22) and (23)), respectively.

Analysis of Humidity Effects
When there are water molecules in the environment, they can also adsorb on the surface of In 2 O 3 , thus occupying the H adsorption site and reducing the concentration of the effective adsorption sites from N − to N − [42,43]. During the humid H 2 sensing measurement in this paper, the humidity was kept at a certain RH level in each round of testing, so it can be assumed that the adsorption of water molecules on the surface of In 2 O 3 was always in equilibrium in every humidity test. In addition to the detection of humid H 2 , the sensitivity to pure humidity of the Pt-In 2 O 3 sensor was also tested ( Figure S3 in the Supplementary Materials). There was no current response in the range of 0-15.5% RH, which indicates that there was no electron transfer between In 2 O 3 and water molecules. Therefore, according to the Langmuir theory, the equilibrium coverage of water molecule is expressed as where P is the partial pressure of water vapor, and α 2 is the Langmuir isotherm equilibrium constant. The relative humidity of the gas RH = P/P s , and P s is the saturation vapor pressure, then, with α 2 = α 2 ·P s With defining ∆N − max as the total change of the adsorption site concentration for water molecules caused by 100% RH at room temperature, the change of adsorption site concentration ∆N − at a certain RH level can be written as Thus, the effective adsorption site concentration N − for H 2 sensing is Then, when the sensor started to detect humid H 2 gas, by using the similar method in Section 3.2.2, according to Equations (9), (10), (14), (16) and (17), we can obtain It is found that |dV s /dt| decreases by comparing Equation (31) with the derivation in Section 3.2.2, indicating that the appearance of H 2 O slowed down the adsorption of H 2 . The above analysis on the effect of humidity is consistent with the test results in Figure 3.

Discussion
By using the dynamic sensing model obtain in Sections 3.2.2 and 3.2.3, the behavior of the sensor can be predicted in real time, and the H 2 concentration in the gas sample can be quickly calculated without waiting for the saturation of the sensor. For the errors in Figure 6, this section analyzed the main reasons and provided directions for further optimization of the sensing model as follows.

1.
Concentrations and types of the reactive oxygen ions. During the hydrogen reduction reaction, the concentration of O − adsorbed on the surface of In 2 O 3 is continuously reduced, which leads to the strong chemical adsorption of additional O ∞ x and O 0 x to generate O − x and contributes to dN − /dt. In addition, it is assumed that the predominant ionic species reacting with H atom is O − while modeling. However, O − 2 may also participate in the sensing reactions, which can influence the indexes of the parameters in Equation (9).

2.
Approximation of integral Ei. Integral Ei is one of the exponential integrals that gives solutions which cannot be evaluated using elementary functions. In order to obtain an analytical expression of dynamic current response, the function of e x /x was used to approximate Ei(x) in this article. The two functions are relatively close when |x| corresponding to |qV s /kT| is large, but as |x| gets closer to 0, the difference between the two becomes larger. Assuming that the initial value of qV s in air is −1 eV, the relative error of e x /x and Ei is 2.57%. Whereas, if |qV s | is reduced by half after H 2 sensing, the relative error of e x /x and Ei increases to 11.33%. Therefore, the more significant the reduction effect of hydrogen is (the higher H 2 concentration), the more obvious the error of the proposed model becomes.

3.
Experimental errors and noise effects. Firstly, the response of the sensor was weak for 0.1% H 2 , which was reflected in the deviation between the fitting calculation result and the real H 2 centration of 0.1% in Figure 6. In addition, the oxygen in the gas samples also consume H 2 under the catalysis of Pt, which may influence the response of the sensor and bring additional noise in the measurement results. Secondly, when the gas injected into the chamber is switched between the reference gas and the gas samples, it is inevitable that there are gas flow fluctuations at the beginning and end of the H 2 sensing period, which can introduce unwanted fluctuations and errors into the test results. This is also the reason that only the relatively stable data in 60 s-70 s in Figure 2 were selected to perform data fitting in Section 3.2.2.

4.
The effect of electric field. During the test, a fixed voltage was applied between the two electrodes of the device, thereby forming a certain electric field. This electric field makes the Fermi level of the sensitive material slightly different at different locations, which can also contribute to the inaccuracy of the model.

Conclusions
In summary, a Pt-In 2 O 3 resistive sensor was investigated for H 2 detection at room temperature. The sensor exhibits reliable sensitivities to variation of H 2 concentration in the range from 0.1 to 0.6%. The current response of the sensor slightly decreased with the increasing humidity of the tested H 2 samples because the water molecules partially occupied the sites for H 2 reaction at the surface of In 2 O 3 . The sensing principle of the sensor was discussed in-depth, and a dynamic current sensing model was proposed. The H 2 concentration of the tested gas sample can be figured out within relative short measurement time by using the dynamic model, which significantly reduces the sensing and recovery time of the sensor without waiting for the saturation of the response. Moreover, the model was also used to explain the effect of humidity on H 2 sensitivity. Finally, the limitations of the proposed model and the reasons behind it were discussed. Future work could involve improving the dynamic sensing model from the aspects mentioned in the discussion, optimizing the inkjet-printed Pt-In 2 O 3 sensor for better sensing properties, and establishing a robust sensing system on chip for real-time gas detection and identification.
Supplementary Materials: The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/s22197306/s1, Figure S1: SEM image and EDS analysis of the Pt-In 2 O 3 sensing material; Figure

Data Availability Statement:
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest:
The authors declare no conflict of interest.