Thermally Stimulated Currents in Nanocrystalline Titania

A thorough study on the distribution of defect-related active energy levels has been performed on nanocrystalline TiO2. Films have been deposited on thick-alumina printed circuit boards equipped with electrical contacts, heater and temperature sensors, to carry out a detailed thermally stimulated currents analysis on a wide temperature range (5–630 K), in view to evidence contributions from shallow to deep energy levels within the gap. Data have been processed by numerically modelling electrical transport. The model considers both free and hopping contribution to conduction, a density of states characterized by an exponential tail of localized states below the conduction band and the convolution of standard Thermally Stimulated Currents (TSC) emissions with gaussian distributions to take into account the variability in energy due to local perturbations in the highly disordered network. Results show that in the low temperature range, up to 200 K, hopping within the exponential band tail represents the main contribution to electrical conduction. Above room temperature, electrical conduction is dominated by free carriers contribution and by emissions from deep energy levels, with a defect density ranging within 1014–1018 cm−3, associated with physio- and chemi-sorbed water vapour, OH groups and to oxygen vacancies.


Introduction
The study of the electronic transport in nanocrystalline Titanium dioxide (nc-TiO 2 ) is motivated by its wide range of application, from catalysis to green energy systems such as Dye Sensitized Solar Cells (DSSCs) [1] and toxic gas sensing devices [2]. In fact, performances of devices based on nc-TiO 2 strongly depend on electron transport mechanisms, which can be very different from those dominating in the bulk single-crystal semiconductor, due to the complex morphology structure and to the huge active surface of this porous material with respect to its volume. Investigation of the electrical transport properties of nc-TiO 2 and its relationship with surface and bulk defects is thus strategic in the perspective of increasing their performance. A model of the electrical conductivity in nc-TiO 2 that takes into account all the complexities of this material (disorder, fractional dimensionality of the nanoporous material and potential barriers between the constituent nanoparticles) is still lacking; its implementation is a cumbersome task. In order to progress in this direction, a systematic conductivity. As a result, mixed conduction mechanisms with hopping and free carriers contributions should be taken into account when discussing the electrical transport in such a material, especially at low-to-moderate temperatures and high carrier densities.
As it is usually done, we consider the nc-TiO 2 as a strongly disordered intrinsic semiconductor whose charge carriers are electrons: then the Density of States (DoS) is modelled in the simplest way as the sum of the extended state contribution (i.e., the conduction band), with the typical square root dependence on energy, and the localized state contribution characterized by an exponential tail, below the conduction band, as for amorphous materials [19]. In fact, this model neglects the DoS features in the intermediate region connecting the exponential Urbach tail and the conduction band, whose description has been object of several studies [20]. Such expressions are rather complex and their use, even at the cost of the computation time, do not provide substantial improvements to the fitting procedure. Moreover, in our measurements the Fermi level never attaints this region and the hopping, as well the multi-trapping, take place over an extended part of the DoS, then coming rather insensitive to the DoS details. Then, the density of states in the overall energy range, putting at ε = 0 the border between extended and localized states, is given by the function: with θ(ε) Heaviside function, E L tailing factor, N L = 0 −∞ g(ε)dε total density of localized states within the band gap. N C0 is related to the effective density of states at the conduction band minimum. The value N C0 = 4 × 10 19 cm −3 · meV −3/2 , used for all subsequent calculations, has been chosen in order to get at room temperature N C = 5.2 × 10 21 cm −3 , typical for TiO 2 . The energy parameter: is introduced so as to have the continuity of g(ε) at ε = 0. The DoS obtained with this procedure is shown in Figure 1. should be taken into account when discussing the electrical transport in such a material, especially at low-to-moderate temperatures and high carrier densities.
As it is usually done, we consider the nc-TiO2 as a strongly disordered intrinsic semiconductor whose charge carriers are electrons: then the Density of States (DoS) is modelled in the simplest way as the sum of the extended state contribution (i.e., the conduction band), with the typical square root dependence on energy, and the localized state contribution characterized by an exponential tail, below the conduction band, as for amorphous materials [19]. In fact, this model neglects the DoS features in the intermediate region connecting the exponential Urbach tail and the conduction band, whose description has been object of several studies [20]. Such expressions are rather complex and their use, even at the cost of the computation time, do not provide substantial improvements to the fitting procedure. Moreover, in our measurements the Fermi level never attaints this region and the hopping, as well the multi-trapping, take place over an extended part of the DoS, then coming rather insensitive to the DoS details. Then, the density of states in the overall energy range, putting at ε = 0 the border between extended and localized states, is given by the function: with ( ) Heaviside function, EL tailing factor, = ( ) total density of localized states within the band gap. is related to the effective density of states at the conduction band minimum. The value = 4 × 10 19 cm −3 meV −3/2 , used for all subsequent calculations, has been chosen in order to get at room temperature = 5.2 × 10 21 cm −3 , typical for TiO2. The energy parameter: is introduced so as to have the continuity of ( ) at = 0. The DoS obtained with this procedure is shown in Figure 1. Density of states (DoS) as obtained with an exponential decay tail in the forbidden gap and a square root dependence in the conduction band region. The boundary is kept at = 0, minimum of the conduction band. Conduction is due to both free electrons in the conduction band and to hopping at localized states.
Mixed electrical conductivity, in case of a n-type nc-TiO2 layer, is the result of transport processes of both hopping electrons within the band tail and free electrons in the conduction band, as depicted in Figure 1, then for the total conductivity we have: Figure 1. Density of states (DoS) as obtained with an exponential decay tail in the forbidden gap and a square root dependence in the conduction band region. The boundary is kept at ε C = 0, minimum of the conduction band. Conduction is due to both free electrons in the conduction band and to hopping at localized states.
Mixed electrical conductivity, in case of a n-type nc-TiO 2 layer, is the result of transport processes of both hopping electrons within the band tail and free electrons in the conduction band, as depicted in Figure 1, then for the total conductivity we have: where σ H and σ C , µ H and µ C , n H and n C , are conductivity, mobility and electron concentration of hopping carriers and free carriers respectively; q is the electron charge. In (3) µ H is an effective mobility that must be calculated taking into account that the hopping probability depends on the site. Assuming the validity of the Einstein relation, the hopping mobility is related to the carrier energy through the diffusivity D(ε): with ν hopping rate and r(ε) average distance between hopping sites available for hopping for carrier with energy ε. The rate for hopping from site i to site j, respectively characterized by energies ε i and ε j , is given by the Miller-Abrahams [21] model: where r ij is the distance between sites at ε i and at ε j , α is the localization radius of the electron and ν 0 is the hopping frequency coefficient. The localization radius can be assumed as a constant, of the order of a few angstrom, for states deeply localized, but when approaching the boundary ε = 0 with the extended states this length diverges. Then, as the hopping probability depends exponentially on this length, the description of the conduction for carriers localized in states near the mobility edge must take into account this behavior. This divergence, according to several models, behaves as ε −γ with γ ranging between 0.5 and 1.5 [19]. Here, we assume for the localization length α: with α 0 carrier localization length of the deep states and z the energy to which α starts to diverge. The evaluation of the hopping conductivity requires some kind of approximation in order to get an expression useful for fitting experimental data. A considerable simplification is usually obtained referring to the "transport energy" level [8,14,22], that plays for hopping the same role of the mobility edge for the free carriers, but for a quantitave description of the TSC is necessary to go beyond this approach [13]. Then, to this purpose, we will follow the approach proposed by Nagy [15] where the hopping conductivity can be calculated as: where r(ε) is an average hopping distance for a site of energy ε, defined in term of the concentration of the available (unoccupied) final states: , Fermi-Dirac distribution function, α given by Equation (6) and B the percolative limit factor, namely the average number of site links, a parameter dependent on the system dimensionality (B = 2.7 for 3D hopping [14]). The integral in Equation (8) is splitted in two terms the first one taking into account, for the level of energy ε , the number of sites reachable with a down conversion process, having probability 1, the second the number of the sites reachable with a thermally activated up conversion process whose probability is given by the Boltzmann factor. The Fermi energy E F is obtained by evaluating the total electron concentration n as the sum of the electron density in the band tail (ε < 0), n H , and in the conduction band (ε > 0), n C ; then n = n C + n H , where: Finally, the contribution of free electrons to the electrical conductivity is given by the typical expression of the conductivity for electrons in the conduction band: with n C as calculated from Equation (9). A typical value for the TiO 2 mobility constant is µ C0 = 1 cm 2 meV 3/2 /Vs, which gives a room temperature mobility µ C = 7.5 × 10 −3 cm 2 /Vs [23].

Rate Equations for Priming and Thermally Stimulated Process
The rate equation for hopping and free carriers is: where S(t) is the generation rate during priming (e.g., light exposure, null during the TSC scan).
Here conduction carriers decay is considered via annihilation on recombination centers or trapping from deep levels, characterized by an active energy level within the gap. The free electron lifetime, τ C , is typically dependent on the defect capture cross section σ t and concentration N t : where m is the effective mass of the free carriers (in TiO 2 it is about 7 times the electron mass m 0 [24]) and v th their average thermal velocity. In general, in a disordered semiconductor, defects may have a spread in energy, so its concentration is calculated through a gaussian distribution: The decay of the hopping carriers, in turn, is given as [14]: 1 In the following, we will describe results of decayed/fractionated Thermally Stimulated Currents (TSC) experiments [25] analyzed with the mixed conduction model described above. In this method, after priming the sample only once at a low temperature, successive cycles of heating/cooling are applied to fractionally deplete levels of lowering energy. The sample is heated up to a first maximum temperature T stop and cooled down to a first minimum temperature T start then it is heated and cooled again to higher T values and so on, eventually up to the final temperature.
At first, we have used the decayed TSC method to get useful information about the DoS shape. For each scan of the decayed TSC, we extracted the couple (E act , Q) with E act = activation energy as determined from the initial rise of the TSC and Q = emitted charge calculated by integrating the TSC of the corresponding scan [25]. In this way, we obtained the DoS shape as a function of the energy by plotting the charge released at each step as a function of the activation energy E act .
Then, delayed TSC measurements have been best fitted using the mixed conduction model, considering a constant heating/cooling rate: β = dT dt , starting from the initial condition due to priming, stated that in between two scans there is a time delay in which the electrical state of the sample evolve very slowly. So, assuming the initial state got from priming, it is possible to fit the entire sequence of the delayed TSC at once. The rate equation is solved starting from an evaluation of the Fermi level, then the distance between hopping sites, the hopping carrier density and the average hopping rate are obtained and finally iterating for each temperature the calculation is performed for the entire experimental data set.
In the high temperature range, above room temperature, when conductivity is mainly due to free carriers and hopping can be neglected, in our model TSC is considered as dominated by deep centers with discrete energy levels E t in the forbidden gap, characterized by a capture cross section σ n (for electrons) and a trap N t concentration. So, we consider the standard TSC expression as [26]: with e n (T) emission constant, Σ surface normal to electric field F, T i initial temperature of the scan. Due to the disorder in the nanocrystalline material, the TSC peak usually results in a peak broader in temperature than the standard one. This is due to the fact that the energy E t of a defect varies within a certain range due to local morphological changes. We have taken into account that by convoluting the TSC peak with a gaussian distribution, as given in Equation (13).

Experimental Set-Up and Procedure
To manufacture our samples, we used a colloidal system produced by Solaronix (Aubonne, Switzerland), containing about 11 wt %. nanocrystalline titanium dioxide mixed with optically dispersing anatase particles (13/400 nm, Ti-Nanoxide D/SP). This commercial product, specifically developed for prototypal electrodes in DSSCs is an anatase titania particle paste for the deposition of active opaque layers. A mixing of large, 400 nm average size, and small, 13 nm average size, nano-particles ensures both very high surface area and efficient light diffusion. We deposited the nc-TiO 2 paste on alumina substrates having two parallel gold contacts, 7 mm long and spaced 0.8 mm; thickness of the film is about 1 µm. A picture of the sample is shown in Figure 2. After deposition, the films have been syntherized in two steps, 30 min each, first at 280 • C and then at 450 • C. The current-voltage characteristics of the sample showed an ohmic behaviour in the overall investigated range (0-100 V) with room temperature resistivity of the order of 2 × 10 8 Ωm [27]. In a typical DSSCs, with a 2 µm-thick nc-TiO 2 film, a voltage of about 0.5 V is applied and an average electric field of about 2.5 × 10 5 V/m is settled. In our TSC measurements we therefore chose to apply a bias of 100 V across the sample, in view to get an electric field of the same order of magnitude, considering the increased distance between our planar electrodes. current-voltage characteristics of the sample showed an ohmic behaviour in the overall investigated range (0-100 V) with room temperature resistivity of the order of 2 × 10 8 Ωm [27]. In a typical DSSCs, with a 2 μm-thick nc-TiO2 film, a voltage of about 0.5 V is applied and an average electric field of about 2.5 × 10 5 V/m is settled. In our TSC measurements we therefore chose to apply a bias of 100 V across the sample, in view to get an electric field of the same order of magnitude, considering the increased distance between our planar electrodes.  To perform TSC measurements in the temperature range 5-300 K, the alumina substrates coated with the nc-TiO 2 have been placed in a sample holder equipped with a 4 Ω wounded-wire heating resistor and a silicon temperature sensor (Leybold GmbH, Köln, Germany). The sample-holder has been inserted into a dewar containing liquid He (LHe) and positioned over the LHe vapours to ensure stable temperatures down to 4.2 K, minimize thermal inertia and reduce possible mismatch between the sample and the thermometer. Details of the experimental setup are given in [28]. Polarization of the sample and current reading was performed by a Keithley 6517 electrometer (Tektronix Ltd., Berkshire, UK), the heater was biased by a TTi QL564P power supply (TTI, Inc., Maisach-Gernlinden, Germany) and temperature was read by a DRC91C temperature controller (Lake Shore Cryotronics, Inc., Westerville, OH, USA). Priming was performed by a Light Emitting Diode (LED) source placed in front of the sample inside the sample holder. We used two priming sources: a 400 nm Ultra Violet (UV) and a 355 nm UV LED (Roithner-Lasertechnik, Vienna, Austria), having 12 mW (typical) and 8.4 mW (maximum) output power, respectively. LEDs were driven by a Systron Donner 110D pulse generator (Systron Dr, Concord, CA, USA). The light spot on the sample during illumination has a diameter of about 2 mm.
Delayed TSC measurements have been performed as follows. We primed the sample at a low temperature T 0 with the LED source, biasing the sample at 100 V. Then, we waited a time interval to make fast transient effects relaxing and to get a constant temperature on the whole sample. Then, fractionated TSC analysis has been carried out performing different heating/cooling cycles up to 300 K. TSC has been also studied in the temperature range from 300 to 630 K using a different chamber where heating/cooling is performed by a system controlling temperature, pressure and gas composition. During each TSC measurement, both in the low and high temperature ranges, the scan rate was fixed at 0.1 K/s. Figure 3 shows a typical TSC spectrum observed in the overall range 5-630 K obtained with the procedure described in the previous section. A fractionated TSC is performed up to 300 K, after priming at 5 K. Then, a second priming is carried out at 300 K inside the high temperature TSC setup and a second TSC analysis is performed up to 630 K. The low temperature analysis is divided into 4 TSC fractions in the ranges: 5-20, 20-80, 80-180, 180-300 K, then, a unique TSC curve is measured after priming at 300 K up to 630 K. At last, a final cooling step from 300 to 250 K is measured to close the whole cycle.

Experimental Results and Discussion
priming at 5 K. Then, a second priming is carried out at 300 K inside the high temperature TSC setup and a second TSC analysis is performed up to 630 K. The low temperature analysis is divided into 4 TSC fractions in the ranges: 5-20, 20-80, 80-180, 180-300 K, then, a unique TSC curve is measured after priming at 300 K up to 630 K. At last, a final cooling step from 300 to 250 K is measured to close the whole cycle.  Main conductivity processes in the two temperature ranges 5-300 K and 300-630 K are different. In fact, in the low T range, hopping conduction is non-negligible against free conduction. Conversely, in the high T range conductivity is mainly due to free carriers. Moreover, in this latter case, the influence of the surrounding gas atmosphere to the charge state of deep discrete levels in the forbidden gap cannot be neglected. Thus, in the following, we will discuss separately results measured below and above room temperature. Figure 4 compares the TSC spectrum measured in the low T range together with TSC data reported in past for single crystal TiO 2 [29]. The experimental peaks for sc-TiO 2 have been multiplied respectively by 1700 (120 K) and 500 (230 K) to compare with those measured with the nc-TiO 2 ones. The lower TSC signal observed in single-crystal TiO 2 indicates that defect concentrations here is far below those encountered in the nanocrystalline morphology. Main conductivity processes in the two temperature ranges 5-300 K and 300-630 K are different. In fact, in the low T range, hopping conduction is non-negligible against free conduction. Conversely, in the high T range conductivity is mainly due to free carriers. Moreover, in this latter case, the influence of the surrounding gas atmosphere to the charge state of deep discrete levels in the forbidden gap cannot be neglected. Thus, in the following, we will discuss separately results measured below and above room temperature. Figure 4 compares the TSC spectrum measured in the low T range together with TSC data reported in past for single crystal TiO2 [29]. The experimental peaks for sc-TiO2 have been multiplied respectively by 1700 (120 K) and 500 (230 K) to compare with those measured with the nc-TiO2 ones. The lower TSC signal observed in single-crystal TiO2 indicates that defect concentrations here is far below those encountered in the nanocrystalline morphology.  [29] for single crystal TiO2 and with a best fit obtained considering a standard TSC emission from discrete energy levels. The experimental peaks for sc-TiO2 have been multiplied respectively by 1700 (120 K) and 500 (230 K) to compare with those measured with the nc-TiO2 ones.

TSC Analysis below Room Temperature
In single crystal TiO2, two peaks related to two discrete energy levels at 120 and 230 K are present, which can be described in terms of standard TSC emissions [26] as given in Equation (16). In nanostructured TiO2, we observe a much broader peak at 120 K, which cannot be described in terms of standard single-level TSC emission, and a peak at 220 K, rather similar to the one measured in single crystal TiO2. In the figure, a best fit of this latter peak is shown, obtained with a standard TSC analysis (Equation (16)), a very high capture cross section, σn ≈ 10 −10 m 2 and ≈ 0.8 eV, ≈ 10 14 cm −3 .
To investigate the origin of the broad band peaked at 120 K, we then performed a measurement in the same temperature range with more delayed heating steps (up to 10). Results are shown in the inset of Figure 5. We note that, in a standard thermally activated emission where carriers are emitted from discrete energy levels towards the corresponding extended band, the current measured at the foot of the TSC curve in each heating/cooling step has a dependence on temperature given by: ( ) ∝ , while if hopping conduction dominates, the dependence should be (Mott's expression):  [29] for single crystal TiO 2 and with a best fit obtained considering a standard TSC emission from discrete energy levels. The experimental peaks for sc-TiO 2 have been multiplied respectively by 1700 (120 K) and 500 (230 K) to compare with those measured with the nc-TiO 2 ones.
In single crystal TiO 2 , two peaks related to two discrete energy levels at 120 and 230 K are present, which can be described in terms of standard TSC emissions [26] as given in Equation (16). In nanostructured TiO 2 , we observe a much broader peak at 120 K, which cannot be described in terms of standard single-level TSC emission, and a peak at 220 K, rather similar to the one measured in single crystal TiO 2 . In the figure, a best fit of this latter peak is shown, obtained with a standard TSC analysis (Equation (16)), a very high capture cross section, σ n ≈ 10 −10 m 2 and E t ≈ 0.8 eV, N t ≈ 10 14 cm −3 .
To investigate the origin of the broad band peaked at 120 K, we then performed a measurement in the same temperature range with more delayed heating steps (up to 10). Results are shown in the inset of Figure 5. We note that, in a standard thermally activated emission where carriers are emitted from discrete energy levels towards the corresponding extended band, the current measured at the foot of the TSC curve in each heating/cooling step has a dependence on temperature given by:  [30]. Inset of Figure 5 shows current measured in the range 5-150 K in a Mott plot. Indeed, the foot of the logarithmic plots at each heating/cooling stage is linear with T −1/4 , in agreement with the fact that hopping conduction is prevailing in this temperature range and that, as suggested by [13], it should be related to a band tail deforming the DoS shape close to the mobility edge.
As described in the previous section, to investigate the DoS shape we evaluated, for each scan, the activation energy in the rising foot range of each peak and the corresponding total emitted charge. The result is shown in Figure 5, where the total charge is plotted as a function of the activation energy.
Results evidence a mono-exponential DoS: f (ε) = Q L e ε E L (E), in the energy range 0.1-0.6 eV with E L tailing factor. ( ) ∝ [30]. Inset of Figure 5 shows current measured in the range 5-150 K in a Mott plot. Indeed, the foot of the logarithmic plots at each heating/cooling stage is linear with T −1/4 , in agreement with the fact that hopping conduction is prevailing in this temperature range and that, as suggested by [13], it should be related to a band tail deforming the DoS shape close to the mobility edge.
As described in the previous section, to investigate the DoS shape we evaluated, for each scan, the activation energy in the rising foot range of each peak and the corresponding total emitted charge. The result is shown in Figure 5, where the total charge is plotted as a function of the activation energy.
Results evidence a mono-exponential DoS: ( ) = (E), in the energy range 0.1-0.6 eV with tailing factor. Best fit gives: = 47.5 meV, = 2.8 × 10 −5 C, values fairly in agreement with literature [31]. We note that in the low energy range of the plot, the charge is lower than the expected value indicated by the exponential trend. This can be explained considering that priming could not fill all the states in the highest part of the band tail, close to the boundary point at ε = 0. As a strong evidence of this, using a deeper UV LED (355 nm) the drop appears for shallower energies with respect to the Best fit gives: E L = 47.5 meV, Q L = 2.8 × 10 −5 C, values fairly in agreement with literature [31]. We note that in the low energy range of the plot, the charge is lower than the expected value indicated by the exponential trend. This can be explained considering that priming could not fill all the states in the highest part of the band tail, close to the boundary point at ε = 0. As a strong evidence of this, using a deeper UV LED (355 nm) the drop appears for shallower energies with respect to the shallower UV LED (400 nm). Considering a 2 mm diameter light spot and a sample thickness of about 1 µm, the effective volume involved in this process is: Vol ≈ π × 10 −12 m 3 and a rough estimate of the density of states is: N L = Q L qVol ≈ 5 × 10 19 cm −3 . To best-fit our TSC measurement, we considered as fixed a group of parameters, related to the crystalline quality of our sample, mainly N C0 , µ C0 , α and γ, already given in [13]. Then, we determined first the best-fit of a TSC peak from a unique scan through a χ-square procedure by opportunely changing variable parameters starting from values of E L and N L obtained from Figure 5. Significant parameters as frequency factor ν 0 , concentration of recombination centres, N rec , initial Fermi level position E F0 (this latter due to the filling procedure and calculated from the conduction band miminum) have been opportunely changed to improve our simulated peak. Figure 6 shows, as an example, how these variable parameters can affect the process. In the figure, the red curve is our best fit, obtained with E L = 60 meV, N L = 10 20 m −3 , N rec = 7.5 × 10 21 m −3 , frequency factor ν 0 = 4 × 10 12 s −1 , E F0 = 66.5 meV. Other curves in the plot are obtained with same parameters apart from one that has been intentionally changed to evidence its influence on the simulation. Green line, uses a higher recombination centre concentration, N rec = 6.0 × 10 22 m −3 , this fasten the decay of the TSC peak in the high temperature region, due to a reduced charge lifetime. Violet curve comes from a lower energy decay factor, E L = 50 meV, enhancing emission of charges at low temperature. Orange line is characterised by a lower frequency factor, ν 0 = 2 × 10 12 s −1 , bringing to a smoothed peak. Finally, the blue curve has been calculated using a deeper initial Fermi level, E F0 = 72.2 meV, corresponding to less initial charge within the hopping band tail after the priming process. In this case, a lower emission especially in the low temperature range is observed, as expected. So, each parameter is influencing a particular range of temperature and a peculiarity of the complex shape of the emission and, even if many parameters are inter-playing in the formation of the whole peak, each can be optimized almost individually to improve the overall simulation. 66.5 meV. Other curves in the plot are obtained with same parameters apart from one that has been intentionally changed to evidence its influence on the simulation. Green line, uses a higher recombination centre concentration, Nrec = 6.0 × 10 22 m −3 , this fasten the decay of the TSC peak in the high temperature region, due to a reduced charge lifetime. Violet curve comes from a lower energy decay factor, EL = 50 meV, enhancing emission of charges at low temperature. Orange line is characterised by a lower frequency factor, ν0 = 2 × 10 12 s −1 , bringing to a smoothed peak. Finally, the blue curve has been calculated using a deeper initial Fermi level, EF0 = 72.2 meV, corresponding to less initial charge within the hopping band tail after the priming process. In this case, a lower emission especially in the low temperature range is observed, as expected. So, each parameter is influencing a particular range of temperature and a peculiarity of the complex shape of the emission and, even if many parameters are inter-playing in the formation of the whole peak, each can be optimized almost individually to improve the overall simulation. Best-fit of TSC experimental data shown in Figure 4 in the low temperature range, obtained using our mixed conductivity model taking care of the band tail in the range down to 0.6 eV plus a discrete level at 0.8 eV, is shown in Figure 7. Of note the agreement between numerical and experimental data in the overall range, up to almost 4 orders of magnitude of the current. A disagreement is observed at high temperature, where the experimental current stabilizes itself on the pAs range, while in the numerical model it decreases to lower values, as a consequence of the Best-fit of TSC experimental data shown in Figure 4 in the low temperature range, obtained using our mixed conductivity model taking care of the band tail in the range down to 0.6 eV plus a discrete level at 0.8 eV, is shown in Figure 7. Of note the agreement between numerical and experimental data in the overall range, up to almost 4 orders of magnitude of the current. A disagreement is observed at high temperature, where the experimental current stabilizes itself on the pAs range, while in the numerical model it decreases to lower values, as a consequence of the decrease of the Fermi level towards midgap. The pAs contribution to the current could be due to the residual presence of water vapor physisorbed on the film surface, as discussed in the next section.

TSC Analysis above Room Temperature
In this temperature range, the effect of hopping conduction should become more and more negligible against free carrier one. Moreover, physisorption and chemisorption mechanisms at surface should also participate to conduction. In particular, dangling bonds at the nc-TiO2 surface are capturing and releasing oxygen depending on pressure and relative humidity, these effects should be possibly investigated separately. To this purpose, we performed different sets of measurements as follows.

Current vs. Temperature with No Priming
Measurements are performed without previous priming. Using a low heating/cooling rate (0.1 K/s) as a first approximation we can assume a quasi-stationary equilibrium. The sample is kept in dark with air at a pressure of 1100 mb, slightly higher than ambient pressure. A typical measurement is presented in Figure 8. The Arrhenius plot shows two distinct ranges: up to about 400 K the current decreases increasing the temperature T, then it increases with T. To explain this behavior we can consider the model proposed by [32], taking account of two dominant defects, one acting as a trap, the other as a Figure 7. Experimental TSC measured in the 5-250 K range and best fit obtained considering the mixed conductivity model taking care of both hopping and free carriers conduction plus emission from a discrete energy level at 0.8 eV.

TSC Analysis above Room Temperature
In this temperature range, the effect of hopping conduction should become more and more negligible against free carrier one. Moreover, physisorption and chemisorption mechanisms at surface should also participate to conduction. In particular, dangling bonds at the nc-TiO 2 surface are capturing and releasing oxygen depending on pressure and relative humidity, these effects should be possibly investigated separately. To this purpose, we performed different sets of measurements as follows.

Current vs. Temperature with No Priming
Measurements are performed without previous priming. Using a low heating/cooling rate (0.1 K/s) as a first approximation we can assume a quasi-stationary equilibrium. The sample is kept in dark with air at a pressure of 1100 mb, slightly higher than ambient pressure. A typical measurement is presented in Figure 8.

TSC Analysis above Room Temperature
In this temperature range, the effect of hopping conduction should become more and more negligible against free carrier one. Moreover, physisorption and chemisorption mechanisms at surface should also participate to conduction. In particular, dangling bonds at the nc-TiO2 surface are capturing and releasing oxygen depending on pressure and relative humidity, these effects should be possibly investigated separately. To this purpose, we performed different sets of measurements as follows.

Current vs. Temperature with No Priming
Measurements are performed without previous priming. Using a low heating/cooling rate (0.1 K/s) as a first approximation we can assume a quasi-stationary equilibrium. The sample is kept in dark with air at a pressure of 1100 mb, slightly higher than ambient pressure. A typical measurement is presented in Figure 8. The Arrhenius plot shows two distinct ranges: up to about 400 K the current decreases increasing the temperature T, then it increases with T. To explain this behavior we can consider the model proposed by [32], taking account of two dominant defects, one acting as a trap, the other as a recombination center (probably associated to dangling bonds at surface, releasing holes via a Figure 8. Current measured as a function of reciprocal temperature during a quasi-stationary heating process, in the range 300-600 K with no priming. Dark: data, red: best-fit. The Arrhenius plot shows two distinct ranges: up to about 400 K the current decreases increasing the temperature T, then it increases with T. To explain this behavior we can consider the model proposed by [32], taking account of two dominant defects, one acting as a trap, the other as a recombination center (probably associated to dangling bonds at surface, releasing holes via a thermally activating process). Neglecting the small contribution of the hopping, the rate equation for the charged carriers (free electrons and holes concentrations are denoted by n, p) is then given by: The first term of the right side is due to emission of electrons from the trap of energy E t1 , concentration N t1 and frequency factor c 1 . The second term describes the recombination of electrons with holes at the recombination center. We here assume that also hole capture is a thermally activated process: p = N t2 c 2 e −E t2 /K B T , with N t2 , E t2 , c 2 respectively concentration, energy and frequency factor of the recombination center and B a probability coefficient. The third term in Equation (17) takes account of other possible free electron removal mechanisms, with a coefficient γ, as trapping from deeper levels, acting on the excess concentration, n 0 being the equilibrium electron concentration.
During the current temperature measurements in dry fluxed air the system is actually in a quasi-stationary regime, so we can reasonably consider dn dt = 0, then: which gives a current dependence with temperature as: and where again Σ is the surface normal to electric field F. Best-fit of our data with Equation (19), shown in Figure 8, is obtained with energy values: E t1 = 1.30 ± 0.05 eV and E t2 = 0.40 ± 0.05 eV. To briefly comment on these values, we observe that e.g., in [33] shallower defect levels at ∼0.24−0.4 eV in anatase TiO 2 were attributed to Ti interstitials, while deeper ones at ∼0.9−1.1 eV to oxygen-vacancies.

TSC after Storage in Dark and Humid Environment
To analyze the effect of water vapor on TSC data, we first primed the sample by keeping it in a controlled humid environment (rh = 20%) in dark at room temperature (T = 300 K) for selected time intervals, up to 14 days. Then, a TSC heating/cooling cycle has been performed by fluxing dry air with a pressure slightly higher than atmosphere, in view to measure only emissions originated during charging in the storage period. As an example, Figure 9a,b show TSC experimental data obtained after 2 and 4 days storage respectively. Curves labeled "data" report experimental measurements, while those labeled TSC show the thermally stimulated currents obtained by subtracting the current measured during the cooling scan from that measured during the heating scan. TSC emissions are observed within two distinct ranges of temperatures: one from ambient temperature, up to 400 K, the other one above 400 K.  An evaluation of the main TSC components involved in these measurements has been carried out in order to identify the origin of the emissions. Measurements show statistically broadened emissions They can be fitted using TSC peaks convoluted with a gaussian as given in Equation (13). Best fits have been obtained by opportunely changing trap concentration N t for a same set of (E t , σ E t , σ n ), within errors, best fitting the two measurements. Up to seven energy levels are required to fit our data. Parameters are shown in Table 1, energy levels are peaked at E to = 0.7-1.14 eV and are characterized by an energy spread σ E t up to 70 meV. As a general trend, increasing the storage time, peaks at low temperatures decrease their concentration N t , while those at high temperatures increase N t . A source of uncertainty in the determination of concentration for the peak at ambient temperature is due to the increasing background current observed during the cooling stage, observed especially in Figure 9a. This effect can indicate reversible charging/discharging of the involved energy states, maybe due to adsorbing/desorbing from the porous alumina substrate. To comment on the origin of these peaks, we observe that our measurements are in agreement e.g., with temperature programmed desorption (TPD) analyses measured in past with TiO 2 after exposure to water. Four peaks at 155, 190, 295, and 490-540 K were observed by [34][35][36][37], the first three assigned to molecular desorption from multilayer, second layer, and first layer states, while the higher temperature feature was assigned to recombinative desorption. Effusion peaks from water were also observed in [34] by thermal desorption in the 150-350 K range from porous nanostructured TiO 2 and attributed to physiosorbed H 2 O, while in [35], two H 2 O effusion peaks were detected around 440 and 650 K. Moreover, physically adsorbed and dissociated H 2 O molecules in nanostructured anatase TiO 2 have been studied by Fourier Transform InfraRed (FTIR) emission spectroscopy at different temperatures in the range 100-300 • C [38]. A 3665 cm −1 band assigned to OH hydrogen bonded (adjacent) OH groups was observed to considerably decrease when the sample was heated from 373 to 573 K, while one at 3705 cm −1 , attributed to isolated OH groups, more difficult to remove from the surface than adjacent OH groups, only slightly changed. A 3250 cm −1 component attributed to the stretching vibration of water molecules that are hydrogen bonded was considerably weakened when heated up to 423 K, while the one at 3400 cm −1 , related to hydrogen-bound surface OH groups (Ti OH), became visible at this temperature. This latter and the 1625 cm −1 band (identified as the water bond-bending vibration mode) finally disappeared at 573 K.
Although a direct correlation between TPD, FTIR and our TSC measurements should be performed when same kind of samples and same experimental conditions are considered, we can qualitatively conclude that TSC peaks in the range up to 400 K should be related to adsorbed molecular water, while dissociated species, as hydrogen-bonded OH groups, should be involved in emissions in the range 400-600 K. In our measurements, increasing the storage time, peaks at low temperature are decreasing, while those at high temperature are increasing, and eventually saturating. This can be explained considering that molecular adsorbed water, in time, is slowly evolving into the formation of hydrogen bonded OH species, more stable at room temperature.

TSC after Illumination in He Atmosphere with Different Pressures
To evidence the effect on TSC of the oxygen-exchange at surface, we performed a set of measurements were the sample was primed in a dry He atmosphere at different pressure, from 10 −6 to 1 bar, at room temperature (T = 300 K), during illumination with a Xe lamp for a selected time interval. Then, TSC heating/cooling cycles were carried out by fluxing dry air with a pressure slightly higher than atmosphere, to measure only emissions originated during charging in the storage period. In fact, it is known that oxygen vacancies can be created by annealing TiO 2 at elevated temperatures in an oxygen-poor environment, such as a pure He gas atmosphere or vacuum condition [39].
Results of photocurrent measurements during priming are shown in Figure 10a. TSC curves after priming in these conditions are shown in Figure 10b.

TSC after Illumination in He Atmosphere with Different Pressures
To evidence the effect on TSC of the oxygen-exchange at surface, we performed a set of measurements were the sample was primed in a dry He atmosphere at different pressure, from 10 −6 to 1 bar, at room temperature (T = 300 K), during illumination with a Xe lamp for a selected time interval. Then, TSC heating/cooling cycles were carried out by fluxing dry air with a pressure slightly higher than atmosphere, to measure only emissions originated during charging in the storage period. In fact, it is known that oxygen vacancies can be created by annealing TiO2 at elevated temperatures in an oxygen-poor environment, such as a pure He gas atmosphere or vacuum condition [39].
Results of photocurrent measurements during priming are shown in Figure 10a. TSC curves after priming in these conditions are shown in Figure 10b. While at 1 bar the photocurrent is almost saturating during priming, at low pressure it increases superlinearly, a fact that can be explained considering the creation of extra oxygen-vacancies, which are releasing an ever increasing free carriers concentration, so favoring the passivation of deep traps during priming. To comment on TSC curves reported in Figure 10b (inset: logarithmic plot) we observe that, in the case of vacuum priming, the tail in the cooling stage of increasing current below 400 K observed in case of humid environment (Figures 8 and 9) is almost absent. Then, higher TSC emissions are observed in the high temperature range in vacuum, showing an increasing number of passivated deep traps in the priming stage. Best fits of the TSC measurements in He atmosphere have been performed starting with the same set of energy levels used in the previous section. They are shown in Figure 11a-c. Measured emissions have been calculated considering TSC peaks statistically While at 1 bar the photocurrent is almost saturating during priming, at low pressure it increases superlinearly, a fact that can be explained considering the creation of extra oxygen-vacancies, which are releasing an ever increasing free carriers concentration, so favoring the passivation of deep traps during priming. To comment on TSC curves reported in Figure 10b (inset: logarithmic plot) we observe that, in the case of vacuum priming, the tail in the cooling stage of increasing current below 400 K observed in case of humid environment (Figures 8 and 9) is almost absent. Then, higher TSC emissions are observed in the high temperature range in vacuum, showing an increasing number of passivated deep traps in the priming stage. Best fits of the TSC measurements in He atmosphere have been performed starting with the same set of energy levels used in the previous section. They are shown in Figure 11a-c. Measured emissions have been calculated considering TSC peaks statistically broadened as given in Equation (13). Results are shown in Figure 11a-c respectively for the cases of 1, 10 −3 , 10 −6 bar. Best-fit procedure turns out in a six-fold emission, with trap parameters listed in Table 2. The fit has been performed considering the same set of values (E t , σ E t , σ n ), within errors, for the three measurements, and best-fitting the TSC scans obtained by subtracting the response during cooling to the one measured during heating and by opportunely changing the trap concentration values N t . Trap parameters given in Table 2 are in agreement with the model in [40] indicating that localized donor states originating from oxygen vacancies are located at 0.75-1.18 eV below the conduction band of Titania. The six peaks used to best-fit our TSC curves are characterized by the same (E t , σ Et , σ n ), within errors, used to fit TSC measurements in the humid environment (apart of the shallowest two levels that here are present in only one component). Here, at every pressure analyzed, shallowest levels have negligible concentrations with respect to deepest levels. Looking to plots in Figure 11, we observe a good agreement between fit and data on a four-orders of magnitude scale. Logarithmic plots are shown as a function of 1/T: the observed linear trend is in favour with our previous observation, that hopping conduction in this high-temperature range is negligible.
To comment on these measurements, we observe that in low pressure/high temperature conditions the probability to form oxygen-vacancies on the surface of the nanostructured titania increases, an evidence widely discussed, e.g., in [41]. It is thus reasonable to hypothesize that TSC peaks shown in Figures 9-11 are related to such a phenomenon, as oxygen-atoms released from the surface leave behind electrons in the conduction band which are collected at electrodes and participate to the emission process. The fact that defects related to these TSC peaks are the same measured after exposure to humid environment, (apart from the shallowest levels which we attributed to adsorbed molecular water) evidences that also OH groups are more likely to be associated to oxygen vacancies, as e.g., suggested in [41].
errors, used to fit TSC measurements in the humid environment (apart of the shallowest two levels that here are present in only one component). Here, at every pressure analyzed, shallowest levels have negligible concentrations with respect to deepest levels. Looking to plots in Figure 11, we observe a good agreement between fit and data on a four-orders of magnitude scale. Logarithmic plots are shown as a function of 1/T: the observed linear trend is in favour with our previous observation, that hopping conduction in this high-temperature range is negligible. Figure 11. TSC response after priming in He atmosphere with a pressure in chamber (a) 1 bar; (b) 10 −3 bar (c) 10 −6 bar. TSC experimental data (black); Blue: TSC heating scan subtracted from background (cooling stage); dotted curves: TSC peaks convoluted with a gaussian, parameters are given in Table 2. Figure 11. TSC response after priming in He atmosphere with a pressure in chamber (a) 1 bar; (b) 10 −3 bar (c) 10 −6 bar. TSC experimental data (black); Blue: TSC heating scan subtracted from background (cooling stage); dotted curves: TSC peaks convoluted with a gaussian, parameters are given in Table 2.

Conclusions
Nanocrystalline Titanium dioxide is widely applied as a high gap semiconducting material in many optoelectronic devices, from solar cells to gas sensors, where its peculiar electronic properties and high chemical reactivity play a crucial role. To attain good performances in terms of efficiency and photoactivity, Titania is mostly used in form of porous nanocrystalline thin films, a material characterized by a high degree of microstructural disorder, which detrimental effect in transport properties should be taken in great care and possibly minimized.
The thermally stimulated current technique is one of the most effective tool for characterizing electrical defects in semiconductors. In this paper, we have used this method to get an overall picture of the defect distribution in nanocrystalline Titania used in state-of-art DSSC devices. The model of TSC used to interpret our experimental results, briefly described in this work, takes account of the heavily disordered microscopic nature of nc-TiO 2 . Mixed conductivity with non-negligible contributions from hopping between localized defects grouped in a band-tail below the conduction band is considered. Moreover, a broadening of the energy levels associated to discrete defects in the forbidden gap, has been accounted for, by convoluting the TSC standard emission with a gaussian distribution. Shallow-to deep energy levels ranging from 0.1 to 1.4 eV have been studied via a thermal spectroscopy spanning from 5 to 630 K and interpreted with this model: main results of our analysis can be summarized as follows.
An exponential DoS tail within the forbidden gap has been observed in the range 0.1-0.6 eV from the bottom of the conduction band, with energy tailing factor of 50-60 meV, characterized by density of states of the order of 1 × 10 20 cm −3 . This tail is responsible for a large TSC emission visible after priming with a UV source at 5 K for temperatures up to approximately 150 K. At higher temperatures, up to room temperature and above, the hopping contribution to conduction becomes more and more negligible against free carrier one and contributions from discrete energy levels emitting in the conduction band become visible. Similar to single crystal TiO 2 , a sharp TSC peak at 220 K is observed, with energy 0.8 eV, probably related to water adsorption. Above room temperature, dark current measured as a function of the temperature without any priming reveals to be non-negligible, as it should be due in pure intrinsic TiO 2 material. So, we studied it separately before any TSC analysis. Measurements as a function of temperature shows a double exponential trend, with a minimum at about 400 K. To explain this behavior a model taking account of two dominant defects, one acting as a recombination center, the other as a trap has been considered. The activation energy measured in this experiment in dark, for the recombination center, is about 0.4 eV a value compatible with dangling bonds at surface [22]. The trap energy, evaluated as about 1.3 eV, is in the highest energy range found for trap states observed in our further analyses.
Then, TSC response above background current has been studied as a function of the temperature, as a spontaneous emission observed after a prolonged storage of the sample in moderately (rh = 20%) humid ambient air. They are best-fitted using gaussian-broadened TSC emissions, with average energy 0.7-1.4 eV, characterised by uncertainties up to 70 meV. Two groups of peaks are measured, respectively below and above temperature. Peaks up to 350 K are probably related to molecular desorption from multilayer, second layer, and first layer states, while the higher temperature features should be assigned to desorption of OH groups [42]. Our measurements show that increasing the storage time molecolar adsorbed water slowly evolves into recombinative species. Finally, to evidence the relationship between TSC emissions and vacancy-oxygen defects, measurements have been carried out after priming in inert atmosphere (He) at different pressures. Main TSC emissions, observed in the high temperature range, 400-630 K, are best-fitted considering energy levels for localized energy states in the range 0.75-1.18 eV, in agreement with a model from [41] for donor states related to oxygen-vacancies below the conduction band.