Abstract

A numerically solvable engineering model has been proposed that predicts the sensitivity of metal oxide- (MOX-) based potentiometric pH sensors. The proposed model takes into account the microstructure and crystalline structure of the MOX material. The predicted pH sensitivities are consistent with experimental results with the difference below 6% across three MOX (RuO2, TiO2, and Ta2O5) analysed. The model distinguishes the performance of different MOX phases by the appropriate choice of surface hydroxyl site densities and dielectric constants, making it possible to estimate the performance of MOX electrodes fabricated through different high-temperature and low-temperature annealing methods. It further addresses the problem, cited by theoreticians, of independently determining the C1 inner Helmholtz capacitance parameter while applying the triple-layer model to pH sensors. This is done by varying the C1 capacitance parameter until an invariant pH sensitivity across different electrolyte ionic strengths is obtained. This invariance point identifies the C1 capacitance. The corresponding pH sensitivity is the characteristic sensitivity of MOX. The model has been applied across different types of metal oxides, namely, expensive platinum group oxides (RuO2) and cheaper nonplatinum group MOX (TiO2 and Ta2O5). High temperature annealed, RuO2 produced a high pH sensitivity of 59.1082 mV/pH, while TiO2 and Ta2O5 produced sub-Nernstian sensitivities of 30.0011 and 34.6144 mV/pH, respectively. Low temperature annealed, TiO2 and Ta2O5 produced Nernstian sensitivities of 59.1050 and 59.1081 mV/pH, respectively, illustrating the potential of using cheaper nonplatinum group MOx as alternative sensor electrode materials. Separately, the usefulness of relatively less investigated, cheap, and readily available MOX, viz. Al2O3, as the electrode material was analysed. Low-temperature-annealed Al2O3 with a Nernstian sensitivity of 59.1050 mV/pH can be considered as a potential electrode material. The proposed engineering model can be used as a preliminary prediction mechanism for choosing potentially cheaper alternative sensor electrode materials.

1. Introduction

Metal oxide (MOX) sensors are of great interest in air and water quality detection systems, with gas sensors being a prime application area [1]. In recent times, development of new sensors and improving the efficiency of existing sensors have received the attention of researchers [2]. The use of MOX enhances the selectivity and sensitivity of existing sensor devices and also allows for real-time sensing of multiple parameters concurrently [1]. This is especially relevant for MOX-based pH sensors which have been gaining popularity in application areas of daily importance such as food quality inspections [3] and in new-generation biomedical devices [4]. The use of WO3 nanoparticle-based conformable pH sensor in wearable biomedical devices [5] is one such example of active research taking place in this domain. Much of the development of electrochemical MOX-based pH sensors has focused on potentiometric pH sensors which involve a MOX-based sensitive electrode and a reference electrode [6]. Potentiometric pH sensors are based on the principle of measuring the potential difference between the reference electrode and the sensitive electrode which are immersed in an analyte solution of unknown pH. The sensitivity of such a sensor is a significant performance parameter since a greater sensitivity would allow for the detection of small pH changes in the analyte solution. The sensitivity performance of a sensor is determined by the degree to which the sensitivity factor is close to the ideal Nernstian response, as predicted by the Nernst equation (59.14 mV/pH at 25°C). A sensitivity factor close to this Nernstian value is thus desired.

The empirical determination of such sensitivity characteristics of MOX through experimental methods has been a matter of great interest to many researchers [711]. A theoretical model to predict the pH sensitivity of the MOX electrodes would be of a complementary value in optimizing the cost and time involved in experiment-based characterization. Furthermore, the sensitivity of a MOX-based pH sensor is found to be dependent not only on metal oxide used for the sensitive electrode but also on the fabrication methods which could influence the microstructure, porosity, surface homogeneity, and crystalline structure of the material as observed by Manjakkal et al. [6]. Thus, the need for a quantitative prediction model for sensitivity that involves both the electrode material and the fabrication method variations is abundantly clear. Furthermore, with current research works attempting to consider mixed metal oxides as possible cheaper alternatives as materials for the sensor electrodes [1214], a theoretical model that predicts the performance of pure metal oxides is of great value in deriving such models for mixed metal oxides.

1.1. The Faradaic and Nonfaradaic Approaches

There exist two broad approaches, one capturing the faradaic and the other nonfaradaic reactions, in modelling the sensing mechanism of MOX pH sensors. Among the five plausible mechanisms proposed by Fog and Buck [15] in their paper, the mechanism that involves ion exchange with the oxide surface layer imbibes both faradaic and nonfaradaic approaches. Both approaches involve the use of the Nernst equation, which captures the transduction of the pH variation as a measurable potential difference variation across the sensor electrodes. Faradaic processes involve the transfer of charged particles (electrons or ions) across the electrode without charge being progressively stored in the electrode [16]. These processes typically involve redox reactions. On the contrary, nonfaradaic processes involve the transfer of only ions, without electron transfer, but with progressive storage of charge. This allows nonfaradaic processes to be viewed as capacitive processes [16].

It is pertinent to discuss the limitations of applying the faradaic approach to model potentiometric sensor behaviour. This approach predicts the sensitivity based only on the number of electrons transferred through the redox reactions between the electrode and the electrolyte. This is because the Nernst equation serves as a linkage between the electrode potential and the concentration of participants in the electrode process [17]. Furthermore, it is an equation that is based on ideal surface conditions. One-electron electrode reaction would always result in an ideal slope factor of 59.14 mV/pH as the sensitivity of the sensor [6].

However, in reality, it is possible to approach this ideal Nernstian response across different fabrication methods only for certain “active” transition oxides such as RuO2 and IrO2 [6]. On the contrary, sensitivity studies on sensor electrodes fabricated using other MOX such as TiO2, Ta2O5, and CeO2, which involve the same one-electron electrode reactions, yield mixed results with some claiming high Nernstian sensitivities while the others sub-Nernstian sensitivities. The differences in the sensor fabrication methodology have been cited as a possible reason behind these observations [7, 18]. These variations however cannot be explained solely by considering faradaic processes.

Low-temperature annealing leads to large grain sizes that result in reduced leakage current due to a smaller number of grain boundaries [19]. On the contrary, high-temperature annealing leads to polycrystalline materials with smaller grain sizes and an increased leakage current due to increased number of grain boundaries [9]. The faradaic approach ignores this aspect of the material phase in polycrystalline materials as it assumes an ideal transfer of electrons across the electrolyte interface, resulting in an overestimation of the sensitivity of such electrodes. In contrast, a nonfaradaic approach, such as the site-binding theory, takes into account parameters such as the total number of surface hydroxyl groups per unit area, NS, and the dielectric constant of the electrode material, , which allow for the distinction of different MOX based on phase differences. The value of Ns depends on the crystalline structure of the MOX material [6].

Since the fabrication method or the annealing temperature of MOX materials affects the phase transformation of the materials [20], any theoretical approach to determine pH sensitivity should take into account Ns and as input parameters. Furthermore, the faradaic approach ignores the ionic interactions between the hydroxyl groups present on the MOX surface and the ions present in the electrolyte bulk. The actual amount of hydroxyl groups on the MOX surface available for sensing hydrogen ions would thus be overestimated in this approach resulting in a prediction of greater than expected sensitivities.

Only the faradaic approach has thus far been suggested for potentiometric sensors as seen in Al-Hilli and Willander’s work on ZnO nanorods [21], while the nonfaradaic approach (site-binding theory) has been suggested for ISFET sensors [22]. Kurzweill, in his work [23], mentioned that metal oxides, in general, can be viewed as mixed electronic and ionic conductors due to their oxygen defect stoichiometry. He claimed that modelling surface phenomena at ion exchanging surface sites is by itself sufficient to explain the mechanism of the pH response of metal oxides and that it is not necessary to include pH-dependent redox transitions. However, the scope of that work was limited to discussing platinum group oxides or “active” transition oxides such as RuO2 in this context.

1.2. A New Perspective on Potentiometric Sensing

A new perspective on the complex mechanism of potentiometric pH sensors has been developed in this work to unify various observations and theories made by experimenters and theoreticians alike in this field. Our perspective considers two groups of metal oxides; the first group is made of active transition oxides and metal oxides fabricated through low-temperature annealing methods, while the second group constitutes metal oxides fabricated through high-temperature annealing methods and nonstoichiometric MOX. The first group can be viewed as ideal pH sensor materials that produce sensitivities close to the predicted theoretical Nernstian sensitivity value. The potential difference developed when this group of metal oxides is in contact with the electrolyte solution can be modelled using a faradaic and a nonfaradaic voltage source connected in parallel, across the space charge and triple-layer capacitance represented as capacitors connected in series as seen in Figure 1.

The faradaic voltage source represents the potential difference developed across the sensitive MOX electrode and reference electrode through Faradaic processes. These processes refer to the redox reactions that take place on the interface of the MOX electrode surface and the electrolyte solution. The charge (electrons and ions) transfer that takes place in these reactions contributes to the voltage that is developed. Trasatti developed a mechanism that captures this process for RuO2 [24]. The following mechanism represents the general faradaic process:where MOx(OH)y is a higher valency metal oxide, while is a lower valency metal oxide.

The other voltage source in the model refers to the potential difference developed through nonfaradaic processes occurring across MOX and reference electrodes. These processes refer to the ion transfer that occurs during surface protonation/deprotonation reactions and surface complexation reactions. The protonation/deprotonation reactions occur between species such as SOH, SO, , and H+ ions in the electrolyte solution according to the following equilibrium [25], where S refers to a generic metal constituent of metal oxide in the study:

Surface complexation reactions occur between species such as SO and and electrolyte ions that are present. The equilibrium can be expressed in the following with K+ and Cl as the electrolyte ions as an example:

The developed voltage across the sensitive MOX electrode and the reference electrode drops across a series of capacitances in the MOX-electrolyte interface, represented as CSC and CSS, and the electrolyte itself, represented as CH and CG [26, 27].

CSC refers to the capacitance developed across the space-charge layer of the MOX electrode, while CSS refers to the surface-state capacitance. CH refers to the Helmholtz layer capacitance, while CG refers to the Gouy diffuse-layer capacitance.

The second group of metal oxides is more likely to produce sub-Nernstian sensitivities and thus could be termed as nonideal MOX materials. This includes nonstoichiometric metal oxides since reduction in oxygen sites on the electrode surface caused by oxygen out-diffusion would result in reduced surface hydroxyl sites and thus reduced pH sensitivity [9]. The potential difference developed when this group of metal oxides is in contact with the electrolyte solution can be modelled as seen in Figure 2. In this model, the dominance of the nonfaradaic process (ion transfer) in causing the potential difference is captured by the single voltage source, which is across the space charge and triple-layer capacitance represented as capacitors connected in series.

2. Proposed Engineering Model

The new perspective discussed above for potentiometric pH sensing mechanisms paves the way for an engineering model to be suggested. The engineering model proposed in this work aims to provide a general mechanism to predict the pH response across different MOX materials and their fabrication methods using the “triple-layer model” (TLM), as developed by Davis et al. [28], as the basis. When the electrode surface comes into contact with the electrolyte solution, the electric charge is induced on the electrode surface. This induced charge disrupts the ions present in the electrolyte, resulting in ions of opposite charge to be adsorbed towards the surface [29]. This layer is known as the Helmholtz layer.

The TLM captures the reactions that take place in the diffuse, inner, and outer Helmholtz layers. The TLM provides the sum of potential drops across the Helmholtz capacitance, VH, and the Gouy diffuse-layer capacitance, VG, for a certain pH value. Some of the model parameters required to apply the TLM were estimated using existing theoretical methods. In our engineering model, the total electrode potential for specific pH is expressed as a summation of the potentials VH, VG and the potential drop across the space-charge layer capacitance, VSC. The contribution to the electrode potential by surface states (VSS) is ignored in this model. This model alleviates the problem of independently determining the C1 inner Helmholtz capacitance parameter required to apply the TLM for pH sensors cited by theoreticians. This has been done by varying the C1 capacitance parameter value until invariant pH sensitivity across different ionic strengths is obtained. It also provides a means to obtain a close estimate of the actual pH sensitivity empirically determined through experiments through the electrolyte ionic strength invariant pH sensitivity.

The annealing temperature of the MOX material to fabricate the electrode has a strong influence on the phase of MOX on the sensing electrode [20]. The model distinguishes the performance of different metal oxide phases by choosing appropriate surface hydroxyl site densities and dielectric constants. The electrode potentials are calculated through the above method for a defined range of pH values for a specific C1 capacitance. The calculated electrode potentials against pH are plotted, and a linear regression fit is used to model the dataset. The slope of the linear regression fit obtained is represented as the sensitivity of the MOX electrode material for the chosen C1 and electrolyte concentrations in this work. These slopes for various C1 capacitances and electrolyte concentrations are used to predict the sensitivity of the MOX material. This work analyses four different metal oxides, RuO2, TiO2, Ta2O5, and Al2O3, and their respective different phases determined by their fabrication method temperature.

There have been other previous attempts made to provide an estimate of C1. One such work focuses on independently evaluating C1 from the intrinsic properties of MOX and the electrolyte solution [30]. The work proposes a theoretical model predicting C1 for different metal oxide combinations through the use of Born solvation theory and the intrinsic properties of the MOX materials such as their dielectric constant. However, this model requires the use of surface charge data from titration to obtain C1 through experimental fits and does not provide a handle for electrolyte ionic strength variation. Since the TLM takes into account ionic strength as a parameter, a C1 value predicted using this method would result in different pH sensitivities when applied in the model. This is to be expected as electrolyte concentrations are claimed to influence the mechanism of adsorption on electrode surfaces [31]. Thus, this method is of limited use in a potentiometric pH sensor context where an invariant pH sensitivity across ionic strengths is required to reflect the experimental observations.

Other studies have also been conducted in a similar vein such as the one by Hwang and Lenhart [32], which evaluated the dependence of Ns on C1 capacitance values for haematite. This work relies on surface titration data obtained from experiments. The results and analysis are thus specific to the haematite material. Such experimental data would otherwise be unavailable for generic MOX to be input as parameters into the TLM, limiting the scope of the work to a specific material.

3. Theory

This engineering model would solely make use of the nonfaradaic approach to model the pH response mechanism of potentiometric pH sensors. The nonfaradaic processes have been demonstrated to account for the potential difference across both ideal and nonideal MOX types as seen from Figures 1 and 2.

3.1. Potential across Helmholtz Layer Capacitance

The nonfaradaic processes have been captured using the same TLM as the basis. It is to be noted that the TLM is a general thermodynamic model that is used to specifically capture the adsorption of electrolyte ions to MOX surfaces [33] and also the protonation/deprotonation reactions that occur on the electrode surface. The MOX surface is occupied by hydroxyl groups due to the dissociative adsorption of water followed by proton displacement [34]. This process is more prevalent for the ideal type of MOX such as RuO2, resulting in a larger amount of hydroxyl groups in these oxides. An equilibrium similar to the above equations (2) and (3) was developed as follows:where K1 and K2 represent the equilibrium constants for the deprotonation and protonation reactions, respectively. The adsorbed hydrogen ions, , can be expressed in the following using the Boltzmann relation:where represents the potential at 0-plane, while e, k, and T represent the elementary charge of a particle, Boltzmann constant, and temperature, respectively. Substituting equations (6) and (8) and equations (7) and (8), the expressions for K1 and K2 can be derived as follows:

Surface complexation reactions occur between species such as SO and and electrolyte ions that are present. The equilibrium is expressed in the following with and represented as the adsorbed electrolyte ions:where Kani and Kcat represent the equilibrium constants of the respective surface complexation reactions. Similar to the adsorbed hydrogen ions, the adsorbed electrolyte ions, and , can be represented using the Boltzmann relation as follows:where represents the potential at the -plane. Substituting equations (11) and (14) and equations (12) and (13), the expressions for Kani and Kcat can be derived as follows:

Combining the surface protonation/deprotonation reactions with the surface complexation reactions, expressions which relate the hydroxyl groups (SOH) with the adsorbed electrolyte ions and hydrogen ions can be derived as follows:where and represent the equilibrium constants for these reactions. These equilibrium constants can be expressed using the Boltzmann expansion:

The TLM is thus particularly applicable to model the pH response mechanism of metal oxide electrode surfaces. , , and represent the potentials in the 0-plane, -plane, and d-plane, respectively. represents the sum of the potential drops across the Helmholtz capacitance (CH) and Gouy diffuse-layer capacitance (CG).

and represent the charge densities in the 0-plane and -plane. They can be expressed as follows:where B can be expressed aswhere F represents Faraday’s constant, while A is the unit surface area of oxide available per unit volume of electrolyte solution. represents the charge density at the d-plane, and it can be expressed as follows:where C and z represent the bulk concentration and charge of the electrolyte counterion under consideration in the diffuse layer. For electroneutrality to hold true, the following expression needs to be true:

The inner Helmholtz layer, represented as the intervening gap between the 0-plane and plane, holds a capacitance of C1. The outer Helmholtz layer, represented as the intervening gap between the -plane and d-plane, holds a capacitance of C2. According to the TLM, these capacitances are assumed to be constant in the regions between the planes [28]. This allows for the potential in the planes to be related to the charge densities according to the following relationships:

Using the surface mass balance, the surface site densities (Ns) are related to the concentrations of various electrode surface species as follows:

Equations (9), (10), (19), and (20) can be rewritten in terms of the concentrations of surface species as shown in the following:

The TLM requires seven sets of parameters to be input, namely, the surface hydroxyl site density, Ns, the surface acidity equilibrium constants, K1 and K2, the complex ionization equilibrium constants, and , and the inner and outer Helmholtz capacitances, C1 and C2. The surface acidity equilibrium constants, K1 and K2, can be theoretically derived for each MOX material as per the theoretical equations derived by Sverjensky and Sahai [25] and can be expressed in a manner consistent with the use of symbols and notations in this work as follows:where represents a parameter which includes Pauling’s bond length of the material (s) and the metal-hydroxyl bond length (rMOH). This rMOH was computed as the sum of the metal-oxygen bond length (rMO) and the O-H bond length which was fixed at 1.01  [25].

The set of equations, equations (21)–(32), can be solved numerically for any pH value and concentration of electrolyte solution.

3.2. Potential across Space-Charge Layer Capacitance

The potential drop across the space-charge layer capacitance (VSC) could be estimated by first focusing on the width of the space-charge layer. The expression for the width of the space-charge layer (W) was derived by Al-Hilli and Willander [21] as seen in the following:where is the permittivity of vacuum, is the dielectric constant of the semiconductor, ND is the density of an electron donor, V is the electrode potential, and Vfb is the flat-band potential. To estimate the value of ND, the total volume of a unit cell, , of MOX in the study was obtained using the lattice parameters of the unit cell as expressed in the following:

The lattice parameters (, , and ) are tabulated in Table 1. ND can thus be expressed as follows:where , , and represent the oxygen vacancy percentage, number of oxygen atoms per unit cell, and number of carriers per vacancy, respectively. The -value was fixed at 2.5% as it was a physically realizable vacancy concentration [40] for MOX. The -value was fixed at a value of 2 as there were two electron carriers for every oxygen vacancy. The total number of oxygen atoms per unit cell for every MOX considered in this work is tabulated in Table 1. Al-Hilli and Willander [21] represented V(0), the potential at the electrode surface (x = 0), as follows:where represents the width of the Helmholtz layer. can be expressed according to Gongadze et al. [41] as follows:where represents the dielectric constant of the electrolyte solution, and this is tabulated in Table 1. The Helmholtz layer capacitance, CH, is represented as the sum of the inner and outer Helmholtz layer capacitances (C1 and C2). The potential at the 0-plane as captured in the TLM, , can be related to V(0) as follows:

Assuming V − Vfb = VSC, equation (40) can be substituted into equation (35) allowing VSC to be expressed as follows:

The sum of the potentials across the Helmholtz layer and the space-charge layer is taken as the total electrode potential, ξ:

The contribution to the total electrode potential by surface states has not been considered in this engineering model. The total electrode potential is calculated for each pH value across the pH range that is analysed. The calculated electrode potentials against pH are plotted, and a linear regression fit is used to model the dataset. The slope of the linear regression fit obtained is represented as the predicted sensitivity of the MOX electrode material in the study. This study analyses four different metal oxides, RuO2, TiO2, Ta2O5, and Al2O3, and their respective different phases which were controlled by their fabrication method temperature.

4. Methodology

4.1. Total Electrode Potential against pH Variation

To analyse the potentiometric sensitivity of the pH sensor, it is necessary to determine the total electrode potential, ξ, of the metal oxide pH sensor. This could be expressed as the sum of the potential differences across the Helmholtz layer and the interface as shown in equation (42).

The potential across the Helmholtz layer VH can be determined by numerically solving the set of equations (21)–(32) for each varying pH value, using a nonlinear solver. In this work, Python’s “fsolve,” which implements the dogleg numerical methods, has been used for the purpose. The potential difference across the space-charge layer can then be calculated from equation (41).

By plotting the total electrode potential against pH, the sensitivity of the MOX pH sensor could be determined for a particular set of parameters that is input into the TLM. Extreme acidic and alkaline pH regimes have been excluded. A pH range of 2 to 11 has been chosen for this analysis. With the exception of C1, the rest of the TLM parameters discussed in Section 4.1 were kept constant for each experimental run. These parameters are tabulated in Table 2. Furthermore, the intrinsic parameters that are required to determine the surface acidity equilibrium constants according to equations (33) and (34) are also provided in Table 2.

The complex ionization equilibrium constants were calculated using equations (19) and (20) with corresponding surface acidity equilibrium constants and electrolyte adsorption equilibrium constants, Kcat and Kani. logKcat and logKani were fixed at 1.9 and 2.2, as derived from literature [28], for electrolyte KCl which was taken as the electrolyte in the study for all the experimental runs.

The outer layer Helmholtz capacitance C2 was kept constant at 20 F/cm2 for every run as it was considered an invariant across different metal oxide electrodes [57]. Since the fabrication temperature governs the number of surface hydroxyl sites present on a metal oxide electrode [20], it was imperative to vary the hydroxyl site density values, Ns, for metal oxides wherever applicable. RuO2, low band-gap metal oxide known for excellent pH sensitivity [58], has a tendency to have higher Ns values since water dissociation on RuO2 rutile surfaces tends to form greater amounts of hydroxyl groups than other rutile surfaces such as TiO2, for instance [59].

Furthermore, the hydrated and unhydrated forms of RuO2 tend to have different site densities. The site density for the unhydrated form of RuO2 of 12.642 sites/nm2 has been extracted directly from literature [42]. To estimate the site density of a hydrated form of RuO2, the contribution to the site density from water dissociation reactions was added to the site density of unhydrated metal oxide. To do this, the amount of water molecules that typically adsorb on rutile surfaces (0.05 water/) [60] and the percentage of water dissociation that occurs on the surface of rutile RuO2 (50%) [59] were extracted from the literature. The product of these values is multiplied by 2 to represent the fact that two hydroxyl groups are formed for every water molecule that gets dissociated. This leads to an addition of 5 sites/nm2 to the unhydrated site density, resulting in a total of 17.642 sites/nm2. These variations have been captured in the experimental runs. For Ta2O5 and Al2O3, two different regimes of phases, namely, low-temperature-annealed polycrystalline and high-temperature-annealed polycrystalline phases, were captured with higher site density and lower site density values, respectively. Since there are no concrete experimental data for the site density values of Ta2O5, two different site density values of 10 sites/nm2 and 12 sites/nm2 were tried in the low-temperature-annealed regime. Two different site density values of 4 and 6 sites/nm2 were tried in the high-temperature-annealed regime. A lack of consensus among experimenters with regard to site density values for TiO2 is evident with a wide range of values suggested by experimenters [52].

Furthermore, it is evident from the literature that surface hydroxyl groups decreased with annealing temperature which resulted in the three distinct phases of TiO2, rutile, anatase, and low-temperature-annealed polycrystalline phases, to be modelled with increasing site density values. Two different site densities, 2.7 and 5 sites/nm2, were tried for the rutile phase, while a single-site density of 7 sites/nm2 and 12.5 sites/nm2 was tried for the anatase and low-temperature-annealed polycrystalline phase, respectively.

Appropriate values of the static dielectric constant, , were extracted from the literature for the respective phases as seen in Table 2. For TiO2, the value for the low-temperature-annealed polycrystalline phase was estimated to be similar to that of the anatase phase [61]. In a separate analysis, it was found that the variation of the value has a very marginal impact on the value when applying the TLM and thus has limited impact on the sensitivity performance of the pH sensor.

4.2. Determination of C1 Helmholtz Inner Layer Capacitance

The Helmholtz inner layer capacitance, C1, is a parameter of the TLM that does not have any robust, predictive, theoretical estimation procedures that could be used. Preliminary attempts have been made by theoreticians such as Sverjensky and Sahai and Dimitri [25, 30] by using electrostatic and solvation theories to estimate C1 capacitances for different metal oxide and electrolyte combinations. Others, such as Yates et al. [62], assumed a value of C1 capacitance to be of a fixed value across all oxides and electrolytes to achieve agreement between calculations and experimentally derived values of surface charge, , and diffuse-layer potential, , data.

However, these methods to determine C1 capacitance do not incorporate the notion of invariant pH sensitivity across different ionic strengths of the electrolyte solutions. In experimental studies of metal oxide pH sensors, the sensitivities are expressed in terms of the Nernstian expression [10]. The Nernstian expression, being primarily redox-based, does not include ionic strength variation effects on pH sensitivity. On the contrary, the TLM is primarily based on ion transfer reactions occurring on the surface of the metal oxide electrode, resulting in varying pH sensitivities across different ionic strengths of electrolyte solutions.

It was thus hypothesized that there should exist an ionic strength invariant sensitivity for some particular C1 capacitance which overcomes the issue of different sensitivities being predicted by the TLM for different ionic strengths of the electrolyte solution. Therefore, the C1 capacitance parameter in the TLM was varied over a wide range of values, and the corresponding pH sensitivity obtained was plotted as a function of C1 capacitance. This procedure was repeated for multiple ionic strengths, namely, 0.1, 0.3, and 0.5 M, to cover the range for which the corresponding parameters such as surface acidity equilibrium constants K1 and K2 are valid [28].

This point of invariance could be observed at both low and high C1 capacitance values. It is therefore imperative to discern which C1 capacitance value provides the true ionic strength invariant pH sensitivity for a specific phase of a metal oxide. A metal oxide fabricated through high-temperature annealing would be expected to develop a greater Helmholtz capacitance in the electrolyte solution owing to its lower surface site density. A lower site density translates to lower active sites for protonation/deprotonation to take place. This allows for a buildup of charges along the inner and outer Helmholtz planes caused by the accumulation of ions from the bulk electrolyte solution. The buildup of these ions, which would have otherwise taken part in the ion transfer reactions, is a reason for the increase in C1 capacitance in the high-temperature-annealed phase. Therefore, the C1 capacitance range was varied according to the phase of the metal oxide.

For high-temperature-annealed metal oxides that are in the polycrystalline phase, the C1 capacitance was varied over a wide range (from 10 F/cm2 to 800 F/cm2 in some cases) as these materials are known to be highly capacitive [63]. On the contrary, for low-temperature-annealed polycrystalline phase materials, the C1 capacitance was varied over a relatively lower range (from 0.01 F/cm2 to 150 F/cm2). The threshold of 150 F/cm2 was based on the maximum C1 capacitance used in literature [28] for a similar study that involves the application of the TLM in the low-temperature annealing regimen.

5. Results and Discussion

5.1. Potentiometric Analysis Using C1 Capacitances from the Literature

A theoretical potentiometric analysis was done by numerically solving the proposed engineering model for a TiO2 electrode based on parameters obtained from Dimitri [30], whose work suggested the use of C1 capacitance of 120 F/cm2 for that specific system.

Figure 3 shows the total electrode potential (ξ) plotted against pH across different ionic strengths (0.1 M, 0.3 M, and 0.5 M) for this system. Figure 3 shows excellent agreement of the data points to the linear regression fit line as seen from the high R2 value of greater than 0.999 for all the electrolyte ionic strengths. The slopes of the three different ionic strength regimes are also different from one another with the 0.3 M slope greater than the 0.1 M slope by 5.93% and the 0.5 M slope greater than 0.3 M and 0.1 M by 2.26% and 8.32%, respectively. The y-intercepts of the linear regression fits are also different from one another with the 0.3 M intercept greater than the 0.1 M intercept by 5.45% and the 0.5 M intercept greater than 0.3 M and 0.1 M by 2.07% and 7.63%, respectively.

The slopes of the regression lines reflect the sensitivities of the pH sensors. It can be seen that there is an increasing trend of pH sensitivity with increasing ionic strengths while using Sverjensky’s choice of C1 capacitance value. This is not in consonance with the experimental setups which do not involve ionic strength variations of the electrolyte for pH sensitivity measurements. This is largely due to the fact that these measurements are based on the ideal Nernstian equation which does not involve electrolyte ionic strength variations. These results are to be expected since Sverjensky’s work does not incorporate the notion of invariant pH sensitivity with differing electrolyte ionic strengths. It is also noteworthy that the predicted pH sensitivities are not in agreement with the experimentally observed value of 58.21 mV/pH [64] reported for the low-temperature-annealed TiO2 system.

5.2. C1 Capacitance Variation

A wide range of C1 capacitance values were applied into the TLM to obtain an invariant pH sensitivity across different ionic strengths. The range of C1 capacitances used for high-temperature- and low-temperature-annealed phases was different in that a smaller range was required for low-temperature-annealed phases and a larger range for high-temperature-annealed phases.

The absolute slopes of the linear regression fits obtained were plotted against their respective C1 capacitance values, for a range of C1 capacitance values from 0 to 150 F/cm2, across different ionic strengths (0.1 M, 0.3 M, and 0.5 M) for TiO2, Ta2O5, and Al2O3 as seen in Figure 4. A polynomial ninth-order fit was applied to each dataset with R2 values exceeding 0.999 for every fit.

The absolute sensitivity of the metal oxide electrode is observed to decrease with increasing C1 capacitance values. There is a rapid drop in sensitivity from 0 to 10 F/cm2 followed by a gradual drop with increasing capacitance. The drop in sensitivity is observed to be greater at lower ionic strengths across the three metal oxides. The regression fits for the three ionic strengths (0.1 M, 0.3 M, and 0.5 M) are observed to intersect at a low capacitance value with a corresponding pH sensitivity of 59.1080, 59.1081, 59.1075, and 59.1050 mV/pH for TiO2, Ta2O5 (Ns = 10), Ta2O5 (Ns = 12), and Al2O3, respectively. The exact C1 capacitances at which this invariant sensitivity is achieved are reported in Table 3. There is no variation in sensitivity across the metal oxides and marginal variation in C1 capacitance values. It can also be observed that, among the three metal oxides, Ta2O5 (Ns = 10) experiences the greatest drop in sensitivity across the capacitance range with a decrease of 13.869 mV/pH across the range of C1 capacitance as compared to Ta2O5 (Ns = 12) with a decrease of 12.712 mV/pH, TiO2 with a decrease of 11.021 mV/pH, and Al2O3 with a decrease of 7.041 mV/pH.

An analysis similar to the above was performed for the high-temperature-annealed phases in that the absolute slopes of the linear regression fits of ξ against pH analysis were plotted against their respective C1 capacitance values across different ionic strengths (0.1 M, 0.3 M, and 0.5 M) for RuO2, TiO2, Ta2O5, and Al2O3 as seen in Figures 5 and 6. The same polynomial ninth-order fit was applied to each dataset with R2 values exceeding 0.999 for every fit.

Similar to the low-temperature-annealed phase, the absolute sensitivity of the metal oxide electrode decreases with increasing capacitance value choice for the high-temperature-annealed phases. RuO2 displays a similar behaviour to its counterparts in the low-temperature-annealed phase in that there is a greater decrease in sensitivity with increasing capacitance value choice for a lower electrolyte ionic strength and a Nernstian pH sensitivity invariance point of 59.1094 mV/pH for a low C1 capacitance. This is expected due to its high site densities [42, 66] even in its high-temperature-annealed rutile phase as compared to other metal oxides such as TiO2 [59]. Its exceptional pH sensitivity, as compared to other metal oxides, has been cited by many experimenters [10, 11]. As expected, the hydrated RuO2 electrode was shown to have a marginal improvement in the sensitivity of 0.0012 mV/pH over the unhydrated RuO2 electrode. This improvement has also been reported in the literature although the improvement in sensitivity is observed to be greater in magnitude of around 3 mV/pH. There may be several factors that are influenced by the fabrication procedure such as point and surface defects present on the electrode that could have contributed to the increased sensitivity observed. The TLM and the determination of K1 and K2 involve only two parameters (Ns and ) that reflects the fabrication process, and thus, minor variations from literature values are expected. Nevertheless, both theoretically derived sensitivities and those reported in the literature are observed to be close to the Nernstian value.

For the other metal oxides, a contrasting behaviour is observed to that of the low-temperature-annealed phase in that there is a greater decrease in pH sensitivity with increasing capacitance value choice for a higher electrolyte ionic strength. Furthermore, the pH sensitivity invariance point is observed at much higher C1 capacitances. Therefore, the pH sensitivity is accordingly sub-Nernstian for these metal oxides in this phase.

For TiO2 and Ta2O5, it was observed that the variation of Ns had an impact on the pH sensitivity. Lower Ns resulted in a lower pH sensitivity and C1 value. The sensitivities for TiO2 ranged from 30.0 to 35.5 mV/pH for C1 ranging from 322 to 728 F/cm2 for increasing Ns values. The lower site density trial run for Ta2O5 had a pH sensitivity of 34.6 mV/pH and a C1 value of 299 F/cm2, while the higher site density run had a pH sensitivity of 36.5 mV/pH and a C1 value of 424 F/cm2. The single run conducted for Al2O3 yielded at a pH sensitivity of 29.7 mV/pH at a C1 value of 131 F/cm2. The results indicate that, among the high-temperature-annealed phases, Al2O3 is the least sensitive pH sensor, followed by TiO2 and Ta2O5. The results of these analyses along with their closest corresponding experimental sensitivities obtained from the literature are presented in Table 3. The sensitivities obtained from the model are broadly within the range of values reported in the literature. For Ta2O5, the sensitivities obtained from the model were within the lower range of the literature values. The differences between the theoretical and experimental sensitivities could be attributed to potential defects present in the electrode. These defects may act as potential sites for hydroxyl groups to be present, resulting in an increase in site density and thus an increase in sensitivity. There may be several factors that are influenced by the fabrication procedure which are not captured by this model. This may in turn cause minor variations from literature values to be observed.

With increased C1 capacitance effect, the ion transfer reactions (protonation/deprotonation) tend to be less effective in reflecting the pH change in the electrolyte solution. This is because a large C1 capacitance effectively acts as a barrier to free ion transfer between the surface of the electrode and the -layer. Since pH changes take place in the bulk electrolyte solution and the ion transfer reactions take place on the surface of the electrode, a large C1 capacitance in the interface would prevent much of the ions present in the bulk electrolyte from reaching the electrode surface. This would cause a lower concentration of ions in the bulk to take part in the protonation/deprotonation reactions that take place on the electrode surface, resulting in a lower sensitivity to be observed with the increase in the C1 capacitance effect. The behaviour of a greater decrease in pH sensitivity with increasing C1 capacitance for lower electrolyte ionic strengths observed for low-temperature-annealed phases contrasts with the same behaviour observed with higher electrolyte ionic strengths for high-temperature-annealed phases. This serves as useful criteria to determine whether the invariant pH sensitivity would occur at a lower C1 capacitance or at a higher C1 capacitance choice. When greater decrease in pH sensitivity is observed for lower electrolyte ionic strengths, this indicates that an invariant pH sensitivity is likely to be found at low C1 capacitance choice. On the contrary, greater decrease in pH sensitivity observed for higher electrolyte ionic strengths indicates that an invariant pH sensitivity is likely to be found at a higher C1 capacitance value choice. Hence, C1 capacitance value needs to be varied over a wider range to determine this sensitivity. It is noteworthy to mention that the sensitivities obtained from this analysis are in close agreement with published results from experiments conducted as seen in Table 3. Experimenters who used high-temperature annealing methods (over 850°C) such as screen printing for fabricating their metal oxide electrodes observed a sub-Nernstian pH sensitivity for nonideal MOX electrodes.

The observation of an invariant pH sensitivity value with regard to electrolyte ionic strength at a low capacitance value for the three different metal oxides in a low-temperature annealing phase indicates that this phase allows for low C1 capacitance buildup. This in turn allows for a high Nernstian sensitivity to be observed in these phases which is in line with the experimental results [9, 64]. The fact that Al2O3 manifests a similar behaviour to that of the other two metal oxides in terms of sensitivity performance in the low-temperature-annealed phase seems significant as this could position Al2O3 as another alternative material for metal oxide pH sensors.

The sub-Nernstian pH sensitivities seen in high-temperature-annealed phases, with the exception of RuO2, indicate that surface site density (Ns) has a significant impact on the sensitivity of the pH sensor. It could also be observed that, among the high-temperature-annealed metal oxides excluding RuO2, a greater Ns results in only a marginal improvement in pH sensitivity but a significant increase in C1 capacitance buildup. This is evident in the observation where an increase in Ns, for the Ta2O5 system, from 4 to 6 sites/nm2 results in an increase in the pH sensitivity of 5.5% from 34.6 to 36.5 mV/pH but an increase in the C1 capacitance of 42% from 298.5 to 423.9 F/cm2. This means that not only does Ns play a role in the sensitivity of the pH sensor but also C1 capacitance buildup. The exact role of Ns affecting C1 capacitance buildup and subsequently pH sensitivity merits further investigation.

5.3. Potentiometric Analysis Based on the Proposed Engineering Model

Figures 79 show the total electrode potential plotted against pH across different phases and different metal oxides based on derived C1 capacitances from the electrolyte ionic strength invariant engineering model proposed in this work. The equation of the linear regression fit has been shown. The slope of the equation represents the sensitivity of the corresponding MOX system.

The low-temperature-annealed metal oxide phases and RuO2 graphs of Figures 7(a), 7(b), and 9 show excellent agreement of the data points to the linear regression fit line as seen from the high R2 value of 1. On the contrary, the remaining high-temperature-annealed phases of Figures 7(c), 7(d), and 8 have marginally lower R2 value ranging from 0.9784 to 0.9947. This is especially prominent with numerical runs with lower Ns values such as in Figures 8(b) and 8(c). This deviation from linearity is more prominent at data points taken near the pHpzc of the respective metal oxides, which are 3.9 for Ta2O5, 5.3–5.7 for TiO2, and 8.8 for Al2O3. The pHpzc point is an important parameter that determines the affinity of the electrode surface to the ions in the electrolyte [67].

At the pHpzc point, the electrode surface has a net charge of zero. In low pH ranges (acidic regime) where solution pH is lower than the pHpzc, the electrode surface is positively charged. On the contrary, in higher pH ranges (basic regime), the electrode surface is negatively charged [68]. The acidic and basic regimes, away from the pHpzc point, are observed to have a more linear trend with notable exception for Al2O3 where the acidic regime tends to show a deviation from linearity.

These results show that a greater site density, apart from improving the magnitude of sensitivity, also ensures that the sensitivity is consistent across the pH range 2 to 11. On the contrary, a lower site density results in poorer sensitivity to occur in the immediate neighbourhood of the pHpzc than the acidic and alkaline regimes. Therefore, when using these metal oxides, it is imperative to find the pHpzc of the metal oxide material and determine if this lies within the pH range of interest for the end-use context of this pH sensor.

6. Conclusion and Future Work

One of the major impediments in applying the TLM in the context of estimating the sensitivity of pH sensors is the independent determination of the inner Helmholtz layer capacitance C1. In this work, we have proposed a new methodology to theoretically estimate the sensitivity of pH sensors which effectively addresses this problem by identifying an electrolyte concentration invariant point in the sensitivity response to C1 variation trends. The corresponding value of C1 and pH sensitivity at this invariance point can be considered as characteristics of MOX. The percentage difference between the calculated and experimentally observed values was below 6% across three MOX (RuO2, TiO2, and Ta2O5) and annealing methods where such experimental sensitivities were available. As observed for TiO2 and Ta2O5, the percentage difference between the values was lower (less than 2%) for low-temperature annealing methods as compared to high-temperature annealing methods. The consistent and close agreement of pH sensitivity thus determined, with experimentally observed values, across different MOX validates the methodology and its utility.

The proposed engineering model and the methodology distinguish the performance of different metal oxide phases by choosing appropriate surface hydroxyl site densities and dielectric constants. It is therefore possible to theoretically estimate the performance of metal oxide electrodes fabricated through different methods, namely, high-temperature and low-temperature annealing. Low-temperature annealing methods provide higher site densities, resulting in higher pH sensitivity values, while high-temperature annealing methods provide lower site densities, resulting in poorer pH sensitivity. This was evident from the fact that the sensitivities of nonplatinum group MOX such as TiO2 and Ta2O5 fabricated through low-temperature annealing methods were predicted to be close to Nernstian sensitivities of 59.1050 and 59.1081 mV/pH. With the exception of RuO2, an “active” platinum group MOX (sensitivity of 59.1082 mV/pH) and the other MOX analysed (TiO2 and Ta2O5), fabricated through high-temperature annealing methods, were predicted to have sub-Nernstian sensitivities of 30.0011 and 34.6144, respectively. These conclusions arrived by applying the model across three different metal oxides, namely, RuO2, TiO2, and Ta2O5, are consistent with experimental observations found in the literature.

The model was also applied to a relatively less investigated material Al2O3, in the context of potentiometric pH sensors. Low-temperature-annealed Al2O3 shows promise as a viable electrode material with a predicted Nernstian sensitivity of 59.105 mV/pH, which can be validated experimentally in the future. This is potentially a significant result since Al2O3 is widely available and is considerably cheaper as compared to other platinum group MOX traditionally used for fabrication sensor electrodes. On the contrary, high-temperature-annealed Al2O3, with a predicted sub-Nernstian sensitivity of 29.6951 mV/pH, displays characteristics similar to other nonplatinum group MOX such as TiO2 and Ta2O5.

The engineering model proposed helps in conducting focused experiments in this domain. Further investigation on the impact of site density by different fabrication methods would help provide a better estimate of the pH sensitivity of the electrode. This work can be further extended to investigate the sensitivity performance of mixed metal oxide mixtures. There can also be an investigation into the relationship between the surface site density and C1 inner Helmholtz capacitance as these two parameters are observed to be the major factors that influence pH sensitivity.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by Sreemagal Computers (Grant no. SC/2020-21/Sensors/284060626574) from Sreemagal Science Foundation, Chennai, India.