SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

This article is part of the series Online First articles in Volume 6 (2011). These articles were published as Online First on SpringerLink in 2010. They should be cited with a 2011 publication year.

Open Access Nano Express

A High Diffusive Model for Nanomaterials

P Di Sia12* and V Dallacasa13

Author Affiliations

1 Department of Computer Science, Faculty of Science, Verona University, Strada Le Grazie 15, 37134, Verona, Italy

2 Faculty of Computer Science, Free University of Bozen, Piazza Domenicani 3, 39100, Bozen-Bolzano, Italy

3 Accademia Nazionale dei Lincei, via della Lungara 10, Rome, Italy

For all author emails, please log on.

Nanoscale Res Lett 2011, 6:39  doi:10.1007/s11671-010-9783-y

The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/6/1/39

Received:29 May 2010
Accepted:9 September 2010
Published:28 September 2010

© 2010 Di Sia et al.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


Considerable attention is today devoted to the engineering of films widely used in photocatalytic, solar energy converters, photochemical and photoelectrochemical cells, dye-sensitized solar cells (DSSCs), to optimize electronic time response following photogeneration. However, the precise nature of transport processes in these systems has remained unresolved. To investigate such aspects of carrier dynamics, we have suggested a model for the calculation of correlation functions, expressed as the Fourier transform of the frequency-dependent complex conductivity σ(ω). Results are presented for the velocity correlation functions, the mean square deviation of position and the diffusion coefficient in systems, like TiO2 and doped Si, of large interest in present devices. Fast diffusion occurs in short time intervals of the order of few collision times. Consequences for efficiency of this fast response are discussed in relation to nanostructured devices.

Correlation Functions; Diffusion; Frequency-Dependent Complex Conductivity; Nanostructures; Semiconducting Oxides; Montecarlo Simulation

Nano Express

One of the most important aspects of nanostructures concerns charge transport, which can be influenced by particle dimensions and assume different characteristics with respect to those of bulk. In particular, if the mean free path of charges due to scattering phenomena is larger than the particle dimensions, one has a mesoscopic system, in which the transport depends on dimensions and one might correct the transport bulk theories by considering this phenomenon. These problems occur also in a thin film, in which the smallest dimension can be less than the free displacement and therefore require variations to existing theoretical transport bulk models. This situation occurs particularly in connection with metal oxide, like transparency, hardness, etc. Therefore a rigorous knowledge of transport properties is to be acquired. To establish the applicability limit of a bulk model and to investigate the time response of systems at nanoscale we have performed a new approach based on correlation functions obtained by a Fourier transform of the frequency-dependent complex conductivity of the system [1]. With this method it is possible to calculate these functions using experimental data obtained by various films, like TiO2 and ZnO also in the form of nanowires, which have increasing interest for their technological, chemical and biomedical applications and which are engineered to reach the desired technical features. Also, the mesoporous films play a very important role, for their applications in devices for energy generation, photocatalytic processes in environment remediations and for the useful electronic propertiestechniques, in particular the Time-resolved THz spectroscopy (TRTS) [2,3]. Starting from the Drude–Lorentz model [4,5] we have obtained directly the correlation function of velocities, the quadratic average distance crossed by the charges as a function of time and the diffusion coefficient D.

From a mathematical viewpoint, the Kubo relation of the linear response must be inverted. But, due to the presence of a half Fourier-transform, it is necessary to modify this relation in such a way that the whole time axis (-∞, +∞) occurs. This procedure is not trivial and not previously found in the literature.

This new formula can be obtained by relying on linear response theory; we have started considering a system with an hamiltonian of the form:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/39/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/39/mathml/M1">View MathML</a>


with H1 having small effects respect to H0, and negligible in the remote past (adiabatic representation). In the case of an electric field of frequency ω we have:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/39/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/39/mathml/M2">View MathML</a>


For an electric field constant in space and depending on time as:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/39/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/39/mathml/M3">View MathML</a>


the time dependent corresponding current is:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/39/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/39/mathml/M4">View MathML</a>


Following the standard time-dependent approach [4], we derived a general formula for the linear response of a dipole moment density <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/39/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/39/mathml/M5">View MathML</a> in the β direction with the electric field <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/39/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/39/mathml/M6">View MathML</a> directed in the α direction, where V is the volume of the system. This permits to deduce the susceptibility χ(ω), which is correlated to σ(ω) via the relation:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/39/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/39/mathml/M7">View MathML</a>


From Eq. 5 we have deduced the real part of σ(ω), denoted in the form σ'(ω) as:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/39/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/39/mathml/M8">View MathML</a>


where Sβα(ω)is the quantity:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/39/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/39/mathml/M9">View MathML</a>


The quantity 〈···〉T is the thermal average, and the exponential factor arises from equilibrium thermal weights for Fermi particles. By considering the identity <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/39/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/39/mathml/M10">View MathML</a>, Eq. 6 can be written in a form containing the velocity correlation function instead of the position correlation function. Assuming the high temperature limit ħω < < KT as usual in systems to be considered in this paper, we obtain:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/39/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/39/mathml/M11">View MathML</a>


The integral in Eq. 8 spans the entire t axis, so we can perform the complete inverse Fourier transform of this equation. It gives:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/39/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/39/mathml/M12">View MathML</a>


with V the volume of the system, K the Boltzmann's constant, T the temperature and σ'(ω) the real part of σ(ω), given by:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/39/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/39/mathml/M13">View MathML</a>


where n is the carrier density, ω0 the proper oscillator frequency, 1/τ the collision frequency [4,5].

The mean squared displacement in relation with the correlation function of velocities is given by:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/39/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/39/mathml/M14">View MathML</a>


By integration of Eq. 9 with Eq. 10, we deduced all the results for <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/39/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/39/mathml/M15">View MathML</a>, R2(t) and <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/39/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/39/mathml/M16">View MathML</a> with x = t/τ [6].

The main advantage of this new formulation is the disposal of exact results for describing the dynamic behaviour, as extracted by time-resolved techniques. In our analytical procedure we have distinguished the case ωo = 0 from the case ωo ≠ 0. For this latter, three cases occur in connection with the sign of the quantity <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/39/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/39/mathml/M17">View MathML</a>. After obtaining the respective σ'(ω), we have found the poles of these functions and then the residues for integration in the complex ω-plane via Cauchy theorem.

We have used our results for discussing transport in a conventional semiconductor such as doped Si, non conventional TiO2 and other systems where anomalous transport has been found.

The most important characteristics of the results are illustrated by concrete examples in Figures 1, 2, 3, 4, and 5.

thumbnailFigure 1. R2 vs. x = t/τ for some representative values of τ, typical of doped Si [7]0 = 0, T = 300 K) (Drude model). A complete description of R2 for Si requires the evaluation of the contribution of the Drude–Lorentz part (see text).

thumbnailFigure 2. R2 vs. x = t/τ for 2 values of τ (ω0 = 1.12 × 1011Hz dot-dashed; ω0 = 2.24 × 1011Hz dashed) for TiO2 (m = 6me, T = 300 K). Saturation values occur at sufficiently large t.

thumbnailFigure 3. R2 vs. x = t/τ at constant ω0τ, for 3 values of τ (ω0 = 0.5 × 1013Hz solid; ω0 = 1013Hz dot-dashed; ω0 = 0.5 × 1014Hz dashed) for TiO2 (m = 6me, T = 300 K).

thumbnailFigure 4. Velocity correlation function vs. x = t/τ for two values of αR (<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/39/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/39/mathml/M18">View MathML</a>) (m = 6me, T = 300 K). Clear exponentially damped oscillations are displayed in this case.

thumbnailFigure 5. Velocity correlation function vs. x = t/τ for some values of αI (<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/39/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/39/mathml/M19">View MathML</a>) (m = 6me, T = 300 K).

In Figure 1, we show R2 for doped Si. For this semiconductor, the conductivity is the contribution of two terms, a Drude–Lorentz term and a Drude term [7]. At large times the Drude–Lorentz term leads to an R2 approaching a constant value (see Figure 2), while the Drude term alone (Figure 1) is the dominant term at large times. Therefore for sufficiently large times, only the Drude term survives.

We observe that the linear relation at large times becomes quadratic at smaller times. The cross-over between the two regimes occurs at times comparable to the scattering time. This means that diffusion occurs after sufficient time has elapsed so that scattering events become significant, while at smaller times the motion is essentially ballistic.

In Figures 2 and 3R2 saturates at high t. The plateau value may assume high values so that R may be larger than the size of the nanoparticles composing the films. In general, these features indicate quite enhanced mobility of carriers in the nanoporous films at small times, in contrast with a commonly expected low mobility in a disordered network.

From these figures we can evaluate the diffusion coefficient <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/39/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/39/mathml/M20">View MathML</a>. It is remarkable that high D are obtained at t/τ of order unity. As an example, from Figure 3 the deduced that diffusion coefficient is D ~ 1 cm2/s for τ = 10-13s, i.e. comparable to the value ~1 cm2/s of the single crystal rutile [2].

From the other hand, much smaller D can develop at long times, with values D = 10-4–10-6 cm2/s typical of a disordered strong scattering system. So, our results indicate quite different behaviour in Si where normal diffusion occurs, and in TiO2 where the Drude–Lorentz model indicates anomalous diffusion.

The physical reason and mechanism of such increase can be traced back to ballistic-like motion of the carriers at early time when scattering is moderate yielding normal diffusion satisfying Eintein's rule and to strong localization due to the scattering at long times with anomalous diffusion with depression of D.

Figures 4 and 5 report the behaviour of the velocity correlation function.

We observe that, according to the equations of our model [8,9], the correlation function of velocities is never a single decreasing exponential of time, but it is in general a more complicated combination of exponentials, or an oscillating function of time.

When <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/39/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/39/mathml/M21">View MathML</a>, there is a change of sign of velocity with respect to initial velocity, a backscattering mechanism as indicated by Smith [10]: there are two regimes in the temporal response characterized by two different characteristic times, the inverted region being dominated by the longer decay time, and the positive velocity region being due to the shorter time; this region becomes the normal state diffusion region when ω0 = 0. These two regimes will give rise to small and large diffusion constants respectively, which will be discussed in connection with time-resolved techniques.

When <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/39/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/39/mathml/M22">View MathML</a>, we observe the presence of damped oscillations of the velocity in time with strong coupling leading to oscillating currents which average to zero in a sufficiently long interval along with the diffusion constant. This regime appearing at large frequencies for a given time constant does not seem to have been observed in real systems.

The results above give a precise indication on response times of a system subjected to charge motion. In the case of doped Si [7] we have verified the Einstein rule, giving rise to normal diffusion. In the case of TiO2 [2,3], anomalous diffusion is found with time-dependent diffusion coefficient vanishing at long times and oscillating behaviour in time of the transport parameters. We have compared our effective diffusion coefficients directly with experimental results [11-13] and with Monte Carlo simulations [14,15], which take into account the overall mechanisms of scattering, including phonons, imperfections and doping centres, and traps, finding that the diffusion coefficient reproduces the values of experimental or simulated coefficients.

We suggest the possibility that our results can give an explanation of the ultra-short times and of high mobilities with which the charges spread in mesoporous nanoparticle TiO2 systems, of deep interest in photocatalitic and photovoltaic systems [16,17]. In particular, the relative short times (few τ) with which charges can reach much larger distances than typical dimensions of nanoparticles indicate easy diffusion for charges photoproduced inside the nanoparticles towards the surface. The unexplained fact found experimentally of ultrashort injection of charge carriers (particularly in Graetzel's cells) can be related to this phenomenon [16,17].

Similar high diffusivity is found in a number of other devices, i.e. GaAs nanowires and ZnO nanoparticles on which terahertz time-resolved spectroscopy has revealed different time transport regimes with high diffusion processes at short times of the order of the scattering time and longer time localized motion due to the effects of scattering [18,19]. Interpretations of these results in terms of the model suggested here can be given.

Recently, an approach for converting nanoscale mechanical energy into electrical energy has been suggested by using piezoelectric zinc oxide (ZnO) nanowires and TiO2 [20]. Such devices have been shown to convert mechanical energy into electric energy with typical ∼1 nW output power per cm2 area. These unexpected efficiencies can be explained by anomalous high diffusion in the oxides of the type presented here.

In summary, we have evaluated the correlation functions for systems for which the Drude–Lorentz model is valid through the formulation of a new Drude–Lorentz-like model [8,9], in which such functions can be obtained as complete Fourier transform of the real part of the frequency-dependent complex conductivity σ(ω). From our results we deduce some important consequences connected with the nanometric film systems, in particular the possibility of a fast response of the transport of charge carriers with a direct consequence for the efficiencies of present devices based on such systems. Of particular interest for nanostructures is the fact that the limiting value of <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/39/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/39/mathml/M23">View MathML</a> reaches several nanometers in only few τ times, which means that R becomes comparable to dimensions of nanoparticles in few scattering events. This implies the possibility of having high mobility of carriers from and towards the surface of nanostructures. This result has possible and interesting implications in photocatalysis and in energy generators, i.e. in photochemical, photoelectrochemical cells and dye-sensitized solar cells (DSSCs) [3]. The principal consequence is the possibility to have high charge conversion efficiencies in particular time intervals. We can thus explain the rather unwaited experimental result that some film oxides as TiO2, in which the percolative layer structure would be expected to provide a low mobility, are in reality endowed with high response times of charge injection and with high mobility.


  1. Ventura PJ, Costa LC, Carmo MC, Roman HE, Pavesi L:

    J Porous Mater. 2000, 7:107. Publisher Full Text OpenURL

  2. Hendry E, Wang F, Shan J, Heinz TF, Bonn M:

    Phys Rev B. 2004, 69:081101(R). Publisher Full Text OpenURL

  3. Schmuttenmaer CA:

    Terahertz Sci Technol. 2008, 1(1):1. OpenURL

  4. Ziman M: Principles of the Theory of Solids. Cambridge University Press, New York; 1979.

  5. Kittel C: Introduction to Solid State Physics. Wiley, New York; 1995.

  6. Rudin W: Real and Complex Analysis. McGraw-Hill International Editions: Mathematics Series, McGraw-Hill Publishing Co.; 1987.

  7. Pirozhenko I, Lambrecht A:

    Phys Rev A. 2008, 77:013811. Publisher Full Text OpenURL

  8. Di Sia P:

    J Comput Theor Nanosci. 2011, 8:1. Publisher Full Text OpenURL

  9. Di Sia P, Dallacasa V, Dallacasa F:

    J Nano Res. 2010, 11:45. Publisher Full Text OpenURL

  10. Smith NV:

    Phys Rev B. 2001, 64:155106. Publisher Full Text OpenURL

  11. Vanmaekelberg D, de Jough PE:

    Phys Rev B. 2000, 61:4699. Publisher Full Text OpenURL

  12. Kambili A, Walker AB, Qiu FL, Fisher AC, Savin AD, Peter LM:

    Physica E Low Dimens Syst Nanostruct. 2002, 14:203. Publisher Full Text OpenURL

  13. Schlichtherle G, et al.:

    J Phys Chem B. 1999, 103:782. Publisher Full Text OpenURL

  14. Romanò L, Dallacasa V:

    J Phys Condens Matter. 1990, 2:4717. Publisher Full Text OpenURL

  15. Bettinelli M, Speghini A, Falcomer D, Daldosso M, Dallacasa V, Romanò L:

    J Phys Condens Matter. 2006, 18:S2149. Publisher Full Text OpenURL

  16. Graetzel M:

    Inorg Chem. 2005, 44(20):6841. PubMed Abstract | Publisher Full Text OpenURL

  17. Graetzel M:

    J Photochem Photobiol C: Photochem Rev. 2003, 4:145. Publisher Full Text OpenURL

  18. Dallacasa V, Dallacasa F, Di Sia P, Scavetta E, Tonelli D:

    J Nanosci Nanotechnol. 2010, 10:1043. PubMed Abstract | Publisher Full Text OpenURL

  19. Parkinson P, Lloyd-Hughes J, Gao Q, Tan HH, Jagadish C, Johnston MB, Herz LM:

    Nano Lett. 2007, 7(7):2162. Publisher Full Text OpenURL

  20. Wang ZL:

    Adv Funct Mater. 2008, 18:3553. Publisher Full Text OpenURL