InTechOpen uses cookies to offer you the best online experience. By continuing to use our site, you agree to our Privacy Policy.

Physics » Optics and Lasers » "Plasmonics - Principles and Applications", book edited by Ki Young Kim, ISBN 978-953-51-0797-2, Published: October 24, 2012 under CC BY 3.0 license. © The Author(s).

Chapter 1

Modelling at Nanoscale

By Paolo Di Sia
DOI: 10.5772/50755

Article top


R2 vs x=t/τ for some representative values ofτ, typical of doped Silicon [23] (ω0 = 0, T=300K) (Drude model). A complete description of R2 for Silicon requires the evaluation of the contribution of the Drude-Lorentz part.
Figure 1. R2 vs x=t/τ for some representative values ofτ, typical of doped Silicon [23] (ω0 = 0, T=300K) (Drude model). A complete description of R2 for Silicon requires the evaluation of the contribution of the Drude-Lorentz part.
R2 vs x=t/τ for 2 values of τ (ω0 = 1.12⋅1011 Hzdot-dashed line; ω0 = 2.24⋅1011 Hzdashed line) for TiO2 (m=6me, T=300K). Saturation values occur at sufficiently large t.
Figure 2. R2 vs x=t/τ for 2 values of τ (ω0 = 1.12⋅1011 Hzdot-dashed line; ω0 = 2.24⋅1011 Hzdashed line) for TiO2 (m=6me, T=300K). Saturation values occur at sufficiently large t.
R2 vs x=t/τ for 3 values of τ (ω0 = 0.5⋅1013 Hzsolid line; ω0 = 1013 Hzdot-dashed line; ω0 = 0.5⋅1014 Hzdashed line) for TiO2 (m=6me, T=300K).
Figure 3. R2 vs x=t/τ for 3 values of τ (ω0 = 0.5⋅1013 Hzsolid line; ω0 = 1013 Hzdot-dashed line; ω0 = 0.5⋅1014 Hzdashed line) for TiO2 (m=6me, T=300K).
Velocities correlation functionvs x=t/τ for some values of αI (m=6me, T=300K).
Figure 4. Velocities correlation functionvs x=t/τ for some values of αI (m=6me, T=300K).
Velocities correlation functionvs x=t/τ for two values of αR (m=6me, T=300K). Evident exponentially damped oscillations are displayed in this case.
Figure 5. Velocities correlation functionvs x=t/τ for two values of αR (m=6me, T=300K). Evident exponentially damped oscillations are displayed in this case.
Behaviour of the diffusion coefficient Dfor αI (Eq. (39) at classical level [1])varying in the interval [0.1-1] and τ = 10-12s. Note that D→0 for t→∞, indicating absence of diffusion at very long times.
Figure 6. Behaviour of the diffusion coefficient Dfor αI (Eq. (39) at classical level [1])varying in the interval [0.1-1] and τ = 10-12s. Note that D→0 for t→∞, indicating absence of diffusion at very long times.

Modelling at Nanoscale

Paolo Di Sia

1. Introduction

One of the most important aspects at nanoscale concerns the charge transport, which can be influenced by particles dimensions and assumesdifferent 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, we have a mesoscopic system, in which the transport depends on dimensions and in principle it is possible to correct the transport bulk theories by considering this phenomenon. A similar situation occurs also in a thin film, in which the smallest nanostructure dimension can be less than the free displacement and therefore requiresvariations to existing theoretical transport bulk models. Therefore a rigorous knowledge of transport properties is to be acquired. From a theoretical viewpoint, various techniques can be used for the comprehension of transport phenomena, in particular analytical descriptions based on transport equations. Many existing and used theories at today regard numerical approaches, not offering analytical results, which would be of great mathematical interest and in every case suitable to be implemented through the experimental data existing in literature and continuously found by the experimentalists.

To establish the applicability limit of a bulk model and to investigate the time response of systems at nanoscale, recently it has appeared a novel theoretical approach, based on correlation functions obtained by a complete Fourier transform of the frequency-dependent conductivity of the studied system [1]. With this model it is possible to calculate exactly the expressions of the most important connected functions, i.e. the velocities correlation function, the mean free displacement and the diffusion coefficient.

At nanoscale we are in the middle between classical and quantum effects, macroscopic and microscopic properties of the matter [2]. Actually one of the most important experimental technique for the study of the frequency-dependent complex-valued far-infrared photoconductivity σ(ω) is the Time-resolved THz Spectroscopy (TRTS), an ultrafast non-contact optical probe; data are normally fitted via Drude-Lorentz, Drude-Smith or effective medium models [3]. The quantum effects in the nanoworld have opened new ways in a lot of old and new technological sectors.

In the following it is illustrated an overview of the fundamental models that, starting from the Drude model, have attempted to analyze and explain in increasing accuracy the transport phenomena of the matter at solid state level and in particular at nanoscale, until the recent developments concerning the variations of Drude-Lorentz-like models, which involve in particular the concept of plasmon.

2. The Drude model

The most meaningful characteristics of metals are their elevated properties of electric and thermal conduction. Over the years this fact has brought to think in terms of a model in which the electrons are relatively free and can move under the influence of electric fields. Historically two models of the elementary theory of metals are born:

The Drude model, published in 1900 and based on the kinetic theory of an electron gas in a solid. It is assumed that all the electrons have the same average kinetic energy Em;

The previous model integrated with the foundations of quantum mechanics, called Sommerfeld model.

In the Drude model the valence electrons of atoms are considered free inside the metal; all the electrons move as an electronic gas. It is assumed that all the electrons have the same energy Em and that it exists a mechanism of collisions among ions and electrons, allowing the thermal equilibrium for the electrons; this fact implies the application of the kinetic theory of gases to such electronic gas. Drude published the theory three years after the discovery of the electron from J. J. Thomson. The free electrons have only kinetic energy, not potential energy; the average energy Em is therefore (3/2)kBT. It can be correlated to an average quadratic speed vm from the relation Em = (3/2)kBT = mvm2/2, denoting m the mass of the free electron. At environment temperature it is vm 107 cm/s, representing the average thermal speed of the electrons. It was assumed also that the electrons have collisions as instant events, i.e. the time of the diffraction is very smaller with respect to every other considered time. Through such collisions, the electrons acquire a thermal equilibrium corresponding to the temperature T of the metal. If the electrons don´t collide, they move in linear way following the Newton laws. The presence of a constant electric field determines an extra average velocity (the drift velocity) given by vd = - (eE/m)t. The relaxation time τ is defined as the average time between two collisions; it is so possible to get a mean free path, defined by lmfp = vmτ. The current density isJ=σcondE, with σcond electric conductivity. This result, obtained by Drude, is an important goal of the classic theory for the conduction of metals (said “Drude theory”).

3. The Drude-Lorentz model

The Lorentz model (1905) is a refining of the Drude model, in which the statistical aspects are specified. The electrons are considered as free charges, with charge “-e”; they are described by a maxwellian velocity distribution. Considering an electron gas in a spatial region with a constant electric field, the drift velocity of the electrons is constant; this corresponds to a current density J proportional to the applied fieldJ=σ0E, with σ0=ne2τ/m (n is the electron density). Estimating the relaxation time τ, Drude and Lorentz have obtained values of conductivity in good accordance with the experiments. In presence of an electric field of the formE(t)=E0eiωt, the complex conductivityassumes the formσω=σ0/(1iωτ). Such model, said “Drude-Lorentz model”, has received great success, but has also underlined series difficulties.

4. The most utilized Drude-Lorentz-like models

One of the most utilized models for describing experimental transport data is the Drude-Lorentz model [4,5]; with such model the main transport parameters are obtainable. Starting from the Drude-Lorentz model, it is possible to obtain the velocities correlation function, from this the quadratic average distance crossed by the charges as a function of time and examine directly the possible compatibility with the Einstein relation.

Considerable variations of this model were made in this years; the most used are the following:

  1. the “Maxwell-Garnett (MG) model”: in this model the dielectric function is given by aDrude term with an additional “vibrational” contribution at a finite frequency ωo [6], leading to a dielectric function of the form:


where the amplitudeωs, the resonant frequency ω0 and the damping constant γ are material dependent constants. The MG model usually describesan isotropic matrix containing spherical inclusions isolated from each other, such as the metal particles dispersed in a surrounding host matrix.

  1. In the “effective medium theories (EMTs)” the electromagnetic interactions between pure materials and host matrices are approximately taken into account [7]. The commonly used EMTs include the Maxwell-Garnett (MG) model and the “Bruggeman (BR) model”, which is a variation of the MG model.

Generally, in the THz regime, the dielectric function εm(ω) consists of contributions of the high-frequency dielectric constant, conduction free electrons, and lattice vibration [8]:


in which ε is the high-frequency dielectric constant, the second term describes the contribution of free electrons or plasmons and the last term stands for the optical phonons.

When the response originates mainly from the contribution offree electrons or plasmons, it is usually adopted the Drude model:


which described well the dielectric properties of metals and semiconductors [7].

If instead the interaction of a radiation field with the fundamental lattice vibration plays a dominant role and results in absorption of electromagnetic wave, due to the creation or annihilation of lattice vibration, the dielectric function εm(ω) mainly consists of the contributions of the lattice vibrations, expressed by the classical pseudo-harmonic phonon model in the first approximation [7]:


5. The Smith model

Smith has started from the response theory for the optical conductivity, considering an electric field impulse applied to a system, in order to examine the answer with respect to the current. He has utilized the following Fourier transform for the frequency-dependent complex conductivity:


A field impulse, which exceeds every other acting force on the system, permits to assume the hypothesis that the electrons initially can be considered totally free; therefore it holds:

with n* effective electron density. After calculation, the real part of σ(ω) results:


If the initial current decays exponentially to its initial value with relaxation time τ, it is possible to write:


from which the standard Drude formula is obtainable:


Eq. (9) can be considered the first term of a series of the form:

The cn factors hold into account of the original electrons speed remained after the n-th collision. The analytical form of the complex conductivity is therefore [9]:

6. Plasmonics

6.1. Introduction

A lot of nanostructures, in particular the nanostructured metals, show very complex and interesting optical properties. One of the most interesting and promising phenomena encountered in these structures are electromagnetic resonances, due to collective oscillations of the conduction electrons, said “plasmons”. Plasmon modes exist in different geometries and in various metals, in particular in noble metals such as gold and silver. Under determined circumstances, plasmons are excited by light; this leads to strong light scattering and absorption and to an increase of the local electromagnetic field. The interest in plasmon modes has started to the beginning of the 20th century, with Zenneck (1907), Mie (1908), Sommerfeld (1909) and other scientists; recent advances regarding the structure, the manipulation and the observation at nanoscale, have increased the study of this scientific topic. The theoretical efforts have encountered also the growing demand at technological level, in particular for semiconductor based integrated electronic components, optical applications and new nano-components. Actually it remains an important challenge for research and development processes, like the guide of light in integrated optical systems and the interface with electronic components. A lot of nanostructures are believed to be one of the key ingredients of the future optoelectronic devices [10].

6.2. Plasmons

The concept of plasma is very useful in the description of some aspects of the interaction radiation/conductive matter; the free electrons of the conductive material (for example a metal) are considered as an electron fluid with high density, of order of1023cm3. Such concept is the base of the classic Drude model, which assumes that the material contains stopped positive ions and a gas of classical not interacting electrons, whose motion results slowed by a force of viscous friction, due to the collisions with the ions, and characterized from a relaxation time τ. The motion of electrons results so casual, due to the continuous collisions with the lattice. The plasma frequency is defined as the proper frequency of the collective motion of electrons in the following way:


with σ0 conductivity of material and common meaning of the other mathematical symbols. The dielectric constant of metal can be written as a function to the plasma frequency; in first approximation it results:


From the comparison between ω andωp, we can deduce the behaviour of the electromagnetic waves arriving to the metal. If it holdsω<ωp, the dielectric constant is negative, therefore its square root is imaginary pure; this involves the reflection of the incident wave. In the contrary case, the square root of the dielectric constant is real and the incident wave can propagate in the medium with a small attenuation.

The plasma oscillations are the fluctuations of charge density, which happen to the frequency ωp and propagate in the metal. The quanta of such fluctuations inside the volume are the “volume plasmons”.

There is also a mechanism on the metal surface, characterized by the “surface plasmons”. Normally it appears in the wavelength range between the visible and the infrared for the interface air/metal. Localized surface electromagnetic excitations can exist at the surface of nanoparticles and metallic nanostructures: the “localized surface plasmons” (LSP) [11]. The frequency and the intensity of the radiation are very sensitive with respect to the dimension, form and morphology of the nanostructures [12]. Nanoparticles and nanostructures with smaller dimensions with respect to the wavelength of the exciting light are characterized by a wide absorption band, normally in the range of the visible and near infrared spectrum.

6.3. Related theoretical models

The models concerning plasmonics simulate the extinction of the “localized surface plasmons resonance” (LSPR) from nanoparticles and nanoarrays; the application regards the calculation of the light absorption and the scattering. The most used models are:

  • the Mie theory [13];

  • the Gans theory [14];

  • the discrete-dipole approximation method (DDA) [15];

  • the finite-difference time-domain method (FDTD) [16].

6.3.1. The Mie theory

The Mie theory is a theoretical approach concerning the optical properties of the nanoparticles. When the nanoparticle dimension is smaller than the wavelength of the incident light, such theory predicts that the extinction caused by a metallic nanosphere is estimated in the quasi-static and dipole limit. The relative cross section of the process results:


where NA is the surface density, a the radius of the metallic nanosphere, εmthe dielectric constant of the medium surrounding the nanosphere, λthe wavelength, εiand εr imaginary and real part of the dielectric function of the metallic nanosphere respectively and χ form factor concerning the nanoparticle. The localized resonance depends also on the interparticle space and on the dielectric constant of the substrate.

6.3.2. The Gans theory

The Gans theory extends the Mie theory to spheroidal particles case. In the case of polarization of incident light parallel to the symmetry axis of the spheroid, the cross section of the process is given by:


with V volume of nanoparticle and Pj depolarization factors along the three cartesian aces, which hold into account of the anisotropic form of particles.

6.3.3. The discrete-dipole approximation method (DDA)

With this method the nanoparticles are divided in a cubic array of N polarizable dipoles, with polarizabilities αi determined by the dielectric function of nanoparticles. The induced dipole Pi of every element results Pi=αiEloc,i in an applied plane wave field. The cross section is determined by:


with k=ω/c (in vacuum). Such method can be applied for the calculation of the absorption, scattering, extinction and other optical properties of nanoparticles of various forms and dimensions. It is considered one of the most important methods for the understanding of the structural characteristics and optical properties of nanomaterials and nanostructures.

6.3.4. The finite-difference time-domain method (FDTD)

This is a method of numerical calculation in order to resolve the Maxwell equations directly in the time domain. Being a time-domain method, the solutions can cover a wide frequency range with a single process of simulation. It is a versatile and useful technique in applications where the resonance frequencies are not exactly known. A great variety of magnetic and dielectric materials can be modelled in relatively simple way with such method [17].

7. Linear response theory: a new interesting idea

We start considering a system with an hamiltonian of the form:

with H1 having small effects with respect toH0, and of the form:


beingλ a real quantity and η positive, so that in remote past the perturbation is negligible (adiabaticrepresentation:limtH1=0). In the case of an electric field of frequency ω we have:


If the electric field is constant in space and its evolution depends on time, it is writable as follows:


The time dependent corresponding current is:


The conductivityσ(ω) is in general a complex function of the frequency ω and can be deduced from linear response theory. Following the standard time-dependent approach,it is possible to find a general formula for the linear response of a dipole moment density B=er/V in the β direction with the electric field E in the α direction, where V is the volume of the system. This permits to deduce the susceptivityχ(ω), which is correlated to σ(ω) through the relation:


Analytical calculations permit to write the real part of the complex conductivity σ(ω) as:




The part .....T in the integral (24) is the thermal average, and the exponential factor arises from equilibrium thermal weights for Fermi particles. By considering the identity


Eq. (23) can be written in a form containing the velocities correlation function instead of the position correlation function. Assuming the high temperature limit ω<<KT (valid in such contests), we obtain:


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


The new introduced key idea is the possibility to perform a complete inversion of Eq. (26) on temporal scale, i.e. considering the entire time axis (-∞, +∞), not the half time axis (0, +∞), as usually considered in literature [9,18]. This idea is viable if we consider the real part of the complex conductivityσ(ω). Via contour integration by the residue theorem in Eq. (26), the integral is determined by the poles ofReσ(ω). This leads to an exact formulation and gives a powerful method to describe the velocities correlation function (and consequently the mean square displacement and the diffusion coefficient) in analytical way.

8. The other important functions in the nano-bio-context

Another interesting quantity at nanolevel is the mean squared displacement of particles at equilibrium, defined as:


where R(t) is the position vector at time t. Through a coordinate transformation relative to the integration region it is possible to rewrite the relation (27) as follows:


The mean squared displacement can therefore be evaluated through the velocities correlation function. Through Eq. (26) it is possible to deduce also the diffusion coefficient D:


The three relations (26), (28), (29) are fundamental in deducing the most important characteristics concerning the transport phenomena.

9. About the complex conductivity σ(ω)

Let us consider a particle in a region in which there is an electric time-oscillating field, direct along z-axis, with an elastic-type and a friction-type force acting on system; the dynamic equation of the particle can be written as:


where m is the mass of the particle, k the strenght constant of the oscillators, λ=m/τand τ is the relaxation time. We consider solutions of the form:


The current density in the field direction is:


From Eq. (32), via analytical calculations, it is possible to obtain the frequency dependent complex conductivity:


The real part of the complex conducivity (33) results [1,19]:


10. About the Drude-Lorentz-like models: a new promising “plasmon model”

Recently it has been published a new formulation of the Drude-Lorentz model [1,20], based on linear response theory and resonant plasmonic mode; it is able to accommodate some of the observed departures and gives a detailed description of the dynamic response of the carriers at nano-level. One of the peculiarities of this new model is the inversion of relation (26), extending the integration on time to the entire time axis (-∞, +∞); this procedure is not trivial and it was introduced for the first time in the nano-bio-context. One of the main advantages is the possibility to obtain analytical relations for the description of the dynamic behaviour of nanosystems. Such relations are functions of parameters, experimentally obtainable through Time-resolved techniques.An important consequence of the application of such formulation turns out to be the mathematical justification of the unexpected experimental results of high initial mobility of charge carriers in devices based on semiconductor nanomaterials, as the dye sensitized solar cells (DSSCs). The increase of the diffusion coefficientimplies a raise of the efficiency of such electro-chemical cells. The new formulation is able to describe very adequately the properties of transport of nano-bio-materials, but it holds also for other objects, like ions [19].

Many experimental data have indicated that plasmon models describe nanostructured systems in particularly effective way. At quantum level, the key factors incorporating the quantum behaviour are the weightsfi, related to plasma frequencies by:


where N is the carrier density,m and e respectively the mass and the charge of the electron.

The velocities correlation function in quantum case assumes the form:





(K is the Boltzmann’s constant, T the temperature of the system,ωi and τifrequencies and decaying times of each mode respectively).

Integrating these expressions through Eqs. (28) and (29), we obtain the mean squared displacement R2 and the diffusion coefficient D[20,21].

11. Some interesting results of the new model

The model has demonstrated at today high generality and good accordance with experiments. It has been applied in relation to the most commonly used and studied materials at nano-level, i.e. zinc oxide (ZnO), titanium dioxide (TiO2), gallium arsenide (GaAs), silicon (Si) and carbon nanotubes (CN), also at nano-bio-sensoristic level [19,22].

In Figure 1 it is showed R2 for doped Silicon. For this systems the conductivity is the contribution of a Drude-Lorentz term and a Drude term. At large times the Drude-Lorentz term leads to an R2 approaching a constant value (Figure 2), while the Drude term alone is the dominant one. Therefore, for sufficiently large times, only the Drude term survives.


Figure 1.

R2 vs x=t/τ for some representative values ofτ, typical of doped Silicon [23] (ω0=0, T=300K) (Drude model). A complete description of R2 for Silicon requires the evaluation of the contribution of the Drude-Lorentz part.

We can 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.


Figure 2.

R2 vs x=t/τ for 2 values of τ (ω0=1.121011Hzdot-dashed line; ω0=2.241011Hzdashed line) for TiO2 (m=6me, T=300K). Saturation values occur at sufficiently large t.

The Drude behaviour is contrasted by the plasmon behaviour. It has reported representative cases at large times x=t/τ>>1 (Figure 2) for TiO2 nanoparticle films with small frequenciesω0, i.e. close to the Drude case. It is observable that all of the curves reach a plateau value at sufficiently long times and that the slope within a given time interval increases withτ. As consequences of this behaviour, the plateau of R2 may become larger than the size of the nanoparticles composing the films, depending on the parameters ω0 andτ, indicating that carriers created at time t=0 will have enhanced mobility in the nanoporous films at small times, on the order of fewτ, in contrast with a commonly expected low mobility in the disordered TiO2 network. Secondly, the time derivative of R2 decreases to zero at very large times, indicating localization by scattering and therefore leading to decreasing mobility and to the existence of a strictly insulating state for t→∞

Few cases corresponding to short time behaviour and large frequencies are showed in Figure 3, with the scattering time in the range 10-13-10-14s.


Figure 3.

R2 vs x=t/τ for 3 values of τ (ω0=0.51013Hzsolid line; ω0=1013Hzdot-dashed line; ω0=0.51014Hzdashed line) for TiO2 (m=6me, T=300K).

In relation to the behaviour of the velocities correlation function, excluding the Drude model (ω0=0), characterized by simple exponential decreasing functions with time, the correlation function of velocities corresponds to either a damped oscillatory behaviour (Eq. (36)), or to a superposition of two exponentials (Eq. (38)), with quite different decay times depending on the value of αI (Figures 4, 5).

Atα1, i.e. close to Drude behaviour, proper case of some systems [3], the current results the superposition of a short τ and long τ/(1αI) time decay modes. The current is a damped oscillating function of time (Eq. (36)) whenω0τ>>1, and will have double exponential behaviour for ω0τ<<1 (Eq. (38)).

From this analysis, it is viable the possibility that the previous results can give an explanation of the ultra-short times and high mobilities, with which the charges spread in mesoporous systems, of large interest in photocatalitic andphotovoltaic systems [24]; the relative short times (few τ), with which charges can reach much larger distances than typical dimensions of nanoparticles, indicate easy charges diffusion inside the nanoparticles. The unexplained fact, experimentally found, of ultrashort injection of charge carriers (particularly in Grätzel’s cells) can be related to this phenomenon.


Figure 4.

Velocities correlation functionvs x=t/τ for some values of αI (m=6me, T=300K).

In systems with charge localization, where the carrier mean free path is comparable to the characteristic dimension of the nanoparticles, the conductivity response becomes more complicated. The response at low frequencies is characterized by an increasing real part of the conductivity and by a negative imaginary part. Deviations from the Drude model become strong in nanostructured materials, such as photoexcited TiO2 nanoparticles, ZnO films, InP nanoparticles [25], semiconducting polymer molecules [26],and carbon nanotubes [27].

In isolated GaAs nanowires theelectronic response exhibits a pronounced surface plasmon mode, that forms within 300 fs, before decaying within 10 ps as a result of charge trapping at the nanowire surface. The conductivity in this case was fitted by using the Drude model for a plasmon and the mobility found to be remarkably high, being roughly one-third of that typical for bulk GaAs at room temperature.

The Smith model with c1=1 is obtained as a limit of the new introduced plasmon model whenαI0. On performing this limit in Eq. (38) one finds the expression obtained by Smith with a scattering time twice the plasmon scattering time τ. The situation αI0 corresponds toτ1/2ωo. From the other hand, both Smith’s and the new model reduce to the Drude model in the limitωo0. So, although the two models are analytically different, their predictions are expected to be quite similar. From Figure 4 it is possible to note that, forαI0, the lower curve is of the form of the Smith curve. The backscattering mechanism invoked by Smith arises in a natural way in this new model, without further assumptions on successive scattering events.


Figure 5.

Velocities correlation functionvs x=t/τ for two values of αR (m=6me, T=300K). Evident exponentially damped oscillations are displayed in this case.

Time evolution of charges can be systematically studied by means of time resolved techniques. The THz technique allows a detailed investigation of very short time behaviour, such as the photoinjection, as well as longer time behaviour such as thermalized motion. The studies have indicated a common mechanism of the short time domain in nanostructures, where carriers are close to Drude behaviour with a rather large diffusion coefficient, followed by a range with time decreasing mobility. This latter stage is characterized by decay of the response as the superposition of a short time and a longer time exponentials.

The time response of ZnO films, nanowires, and nanoparticles to near-UV photoexcitation has been investigated in THz experiments. Films and nanoparticles show ultrafast injection, but with the addition of a second slower component. For ZnO nanoparticles, double exponential decay indicates characteristic times τ194ps andτ22.4ns. From Eq. (38) it is possible to deduce a ratio of two different relaxation times τ1 andτ2; it isτ1/τ2=(1αI)/(1+αI). Therefore, on using the experimental values, we findαI=0.92, from which we deduceω0τ0.2. For τ of the order of10-13s, we obtainω01.51012Hz, which is of the correct order of magnitude for the resonance in the infrared. The same procedure for injection times of films (500 nm grains) leads to ω02.51012Hz (forτ1013s).

In GaAs nanowires [27,28] a resonance at ωo=0.30.5THz has been suggested to explain the frequency-dependent conductivity obtained from THz experiments. A long characteristic obtained time τc1.1÷51012s is in accordance with the new model, corresponding to τc=τ/(1αI)=5.1012s and leading to αI=0.1÷0.8 for τ=1012s andω0τ=0.3÷0.49.

Similar considerations can be applied to single walled carbon nanotubes (SWCN), where the resonance state is at ωo=0.5THz with Drude-Lorentz scattering time τ ≈ 10 -13 s and to sensitized TiO2 nanostructured films. In the latter the conductivity of electrons injected from the excited state of the dye molecule into the conduction band of TiO2 is initially more Drude-like, and then evolves into a conduction dominated by strong backscattering as the carriers equilibrate with the lattice. These results conform to the backscattering mechanism of the new plasmon model.

In Figure 6 it is presented the fitting of data of GaAs photoconductivity [27,28] with the new model.


Figure 6.

Behaviour of the diffusion coefficient Dfor αI (Eq. (39) at classical level [1])varying in the interval [0.1-1] and τ = 10-12s. Note that D→0 for t→∞, indicating absence of diffusion at very long times.

In this figure the upper curve is the Drude result corresponding toω0=0. It is the highest value of the diffusion coefficient, as ω0 varies from zero to its maximum value. The dots in Figure 6 are experimental data derived from THz spectroscopy, the error bars on dots represent the distribution of the experimental data and the lines the predictions of the new plasmon model. D tends to a vanishing value as t→∞, the larger ω0 the faster this vanishing occurs. The results can be interpreted as follows: at early time, of the order of τ, the system behaves as Drude-like, irrespective ofω0, with carriers assuming large mobility values (the numerical evaluation gives D ≈ 7.5 cm2/s for the case of Figure 6); at increased times, charges become progressively localized as a result of scattering,significantly in agreement with the conclusions drawn by THz spectroscopy [29,30].

12. Conclusions

It is hoped that this new plasmon model will be useful to describe experimental data as an alternative to other generalizations of the Drude model. The principal findings of it are the reversal of the current in small systems like nanoparticles and the time dependence of the transport parameter D, which describes the dynamics of the system at short and long time, with increased mobilities at very short times followed by localization at longer times. This unusual behaviour, with respect to the predicted Drude-like, converges as a whole in the peculiar frequency dependence of the optical conductivity. The strength of the new model consist of its ability to accommodate this behaviour and include previous models, like the Smith model [19].


1 - P. Di Sia, An Analytical Transport Model for NanomaterialsJournal of Computational and Theoretical Nanoscience201188489
2 - Wolf EL.Nanophysics and Nanotechnology:An Introduction to Modern Concepts in Nanoscience.John Wiley & Sons; 2006
3 - Schmuttenmaer CA.Using Terahertz Spectroscopy to Study Nanomaterials. Terahertz Science and Technology 20081118
4 - M. Ziman, Principles of the Theory of Solids.Cambridge University Press; 1979
5 - C. Kittel, Introduction to Solid State PhysicsWiley New York; 1995
6 - O. Levy, D. Stroud, Maxwell Garnett theory for mixtures of anisotropic inclusions: Application to conducting polymersPhysical ReviewB 1997561380358046
7 - J. Han, W. Zhang, W. Chen, S. Ray, J. Zhang, M. He, A. K. Azad, Z. Zhu, Terahertz Dielectric Properties and Low-Frequency Phonon Resonances of ZnO NanostructuresThe Journal of Physical Chemistry C2007111351300013006
8 - J. G. Han, F. Wan, Z. Y. Zhu, Y. Liao, T. Ji, M. Ge, Z. Y. Zhang, Shift in low-frequency vibrational spectra of transition-metal zirconium compounds.Applied Physics Letters 20058717172107172109
9 - Smith NV.Classical generalization of the Drude formula for the optical conductivityPhysical Review B20016415155106155111
10 - A. N. Shipway, E. Katz, I. Willner, Nanoparticle arrays on surfaces for electronic, optical and sensoric applications. ChemPhysChem200011852
11 - E. Hutter, J. H. Fendler, Exploitation of Localized Surface Resonance. Advanced Materials 2004; 161916851706 .
12 - Willets KA, Van Duyne RP.Localized Surface Plasmon Resonance Spectroscopy and Sensing.Annual Review of Physical Chemistry200758267297
13 - Stuart DA, Haes AJ, Yonzon CR, Hicks EM, Van Duyne RP.Biological applications of localised surface plasmonic Phenomenae.IEE Proceedings-Nanobiotechnology 200515211332
14 - A. Gulati, H. Liao, J. H. Hafner, Monitoring Gold Nanorod Synthesis by Localized Surface Plasmon Resonance.The Journal of Physical ChemistryB 2006110452232322327
15 - Zhang ZY, Zhao Y-P.Optical properties of helical Ag nanostructures calculated by discrete dipole approximation method.Applied Physics Letters200790221501221503
16 - Sherry LJ, Chang S-H, Schatz GC, Van Duyne RP.Localized Surface Plasmon Resonance Spectroscopy of Single Silver Nanocubes.Nano Letters200551020342038
17 - W. Zhang, Z. Yue, C. Wang, S. Yang, W. Niu, G. Liu, models. Theoretical, Fabrications, of. Applications, Surface. Localized, Resonance. Plasmon, sensors, IEEE 20104Project supported by the National Natural Science Foundation of China (Grant N. 60574091, 60871028).
18 - A. Sassella, A. Borghesi, B. Pivac, L. Pavesi, Characterization of porous silicon by microscopic Fourier transform infrared spectroscopy; 266267th International Conference on Fourier transform spectroscopy; 1994
19 - P. Di Sia, Classical and quantum transport processes in nano-bio-structures: a new theoretical model and applications. PhD Thesis. Verona University- Italy; 2011
20 - P. Di Sia, An Analytical Transport Model for Nanomaterials: The Quantum VersionJournal of Computational and Theoretical Nanoscience201293134
21 - P. Di Sia, Oscillating velocity and enhanced diffusivity of nanosystems from a new quantum transport model.Journal of Nano Research 2011164954
22 - P. Di Sia, New theoretical results for high diffusive nanosensors based on ZnO oxides. Sensors and Transducers Journal 2010122118
23 - I. Pirozhenko, A. Lambrecht, Influence of slab thickness on the Casimir forcePhysical Review A 200877013811013818
24 - M. Grätzel, Solar Energy Conversion by Dye-Sensitized Photovoltaic Cells.Inorganic Chemistry2005442068416851
25 - H. Nienhuys-K, V. Sundström, Influence of plasmons on terahertz conductivity measurements.Applied Physics Letters200587012101012103
26 - E. Hendry, M. Koeberg, J. M. Schins, H. K. Nienhuys, V. Sundström, L. D. A. Siebbeles, M. Bonn, Interchain effects in the ultrafast photophysics of a semiconducting polymer: THz time-domain spectroscopy of thin films and isolated chains in solutionPhysical Review B 200571125201125210
27 - P. Parkinson, J. Lloyd-Hughes, Q. Gao, H. H. Tan, C. Jagadish, M. B. Johnston, L. M. Herz, Transient Terahertz Conductivity of GaAs NanowiresNano Letters20077721622165
28 - P. Parkinson, H. J. Joyce, Q. Gao, H. H. Tan, X. Zhang, J. Zou, C. Jagadish, L. M. Herz, M. B. Johnston, Carrier Lifetime and Mobility Enhancement in Nearly Defect-Free Core−Shell Nanowires Measured Using Time-Resolved Terahertz Spectroscopy. Nano Letters 20099933493353
29 - P. Di Sia, V. Dallacasa, F. Dallacasa, conductivity. Transient, nanostructured. in, films, Journal of Nanoscience and Nanotechnology 20111116
30 - P. Di Sia, V. Dallacasa, charge. Anomalous, a. transport, “time. new, generalization. domain”, the. of, model. Drude, Plasmonics 20116199104