Open access peer-reviewed chapter

Modelling at Nanoscale

By Paolo Di Sia

Submitted: January 10th 2012Reviewed: June 18th 2012Published: October 24th 2012

DOI: 10.5772/50755

Downloaded: 2524

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 vm107 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 Jproportional 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:

ε//(ω)=1ωp2ω(ω+i/τ)+ωs2ω02ω2iγωE1

where the amplitudeωs, the resonant frequency ω0and 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]:

εm(ω)=εωp2ω2+iγω+jεstjωTOj2ωTOj2ω2iΓjωE2

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:

εm(ω)=εωp2ω2+iγωE3

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]:

εm(ω)=ε+εstωTO2ωTO2ω2iγωE4

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:

σ(ω)=0j(t)exp(iωt)dtE5

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:

j(0)=n*e2/mE6

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

0Reσ(ω)dω=π2j(0)=ωp28E7

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

j(t)/j(0)=exp(t/τ)E8

from which the standard Drude formula is obtainable:

σ(ω)=(n*e2τ/m)/(1iωτ)E9

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

j(t)/j(0)=exp(t/τ)[1+n=1cn(t/τ)n/n!]E10
The cnfactors hold into account of the original electrons speed remained after the n-th collision. The analytical form of the complex conductivity is therefore [9]:
σ(ω)=n*e2τm(1iωτ)[1+n=1cn(1iωτ)n]E11

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:

ωp2=σ0ε0τ=ne2ε0mE12

with σ0conductivity 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:

ε(ω)=1ωp2ω2E13

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 ωpand 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:

σext=24πNAa3εm3/2λln(10)[εi(εr+χεm)2+εi2]E14

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 εrimaginary 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:

σext=2πVεm3/23λj[(1/Pj)2εj(εr+1PjPjεm)2+εj2]E15

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 αidetermined by the dielectric function of nanoparticles. The induced dipole Piof every element results Pi=αiEloc,iin an applied plane wave field. The cross section is determined by:

σext=4πk|Einc|2j=1NIm(Einc,j*Pj)E16

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:

H=H0+H1E17

with H1having small effects with respect toH0, and of the form:

H1=λAeiωteηtE18

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:

H1=eErE19

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

E=E0eiωtE20

The time dependent corresponding current is:

J(t)=σ(ω)E(t)E21

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/Vin the β direction with the electric field Ein the αdirection, where V is the volume of the system. This permits to deduce the susceptivityχ(ω), which is correlated to σ(ω)through the relation:

1+4πχ(ω)=1+4πiσ(ω)ωE22

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

Reσβα(ω)=e2ωπVSβα(ω)(1eω/KT)E23

with:

Sβα(ω)=+dtrα(0)rβ(t)TeiωtE24

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

v=ddtr=i[H,r]E25
,

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:

Reσβα(ω)=e22VKT+dtvα(0)vβ(t)TeiωtE26

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

<vα(0)vβ(t)>T=KTVπe2+dωReσβα(ω)eiωtE27

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:

R2(t)=[R(t)R(0)]2E28

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:

R2(t)=20tdt'(tt')v(t')v(0)E29

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:

D(t)=dR2(t)2dt=0tdt'v(t')v(0)E30

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:

mz¨=kzλz˙+eEoeiωtE31

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:

z(t)=z0eiωtE32

The current density in the field direction is:

jz=nez˙=nez0(iω)eiωt=σEE33

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

σ=iωne2m(ω2ω02+iλω)E34

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

Reσ(ω)=τne2mω2/τ2(ω02ω2)2+ω2/τ2E35

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:

ωpi2=4πNe2mfiE36

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:

v(0)v(t)=(KTm)i(fiexp(t2τi)[cos(αiR2tτi)1αiRsin(αiR2tτi)])E37

with:αiR=4τi2ωi21;(37)

v(0)v(t)=12(KTm)i((fiαiI)[(1+αiI)exp((1+αiI)2tτi)(1αiI)exp((1αiI)2tτi)])E38

with:αiI=14τi2ωi2.(39)

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

Integrating these expressions through Eqs. (28) and (29), we obtain the mean squared displacement R2and 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 R2for 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 R2approaching 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.12⋅1011 Hzdot-dashed line; ω0 = 2.24⋅1011 Hzdashed 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 R2may become larger than the size of the nanoparticles composing the films, depending on the parameters ω0andτ, indicating that carriers created at time t=0will 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 R2decreases 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.5⋅1013 Hzsolid line; ω0 = 1013 Hzdot-dashed line; ω0 = 0.5⋅1014 Hzdashed 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=1is 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 αI0corresponds 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 τ194psandτ22.4ns. From Eq. (38) it is possible to deduce a ratio of two different relaxation times τ1andτ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.5THzhas been suggested to explain the frequency-dependent conductivity obtained from THz experiments. A long characteristic obtained time τc1.1÷51012sis in accordance with the new model, corresponding to τc=τ/(1αI)=5.1012sand leading to αI=0.1÷0.8for τ=1012sandω0τ=0.3÷0.49.

Similar considerations can be applied to single walled carbon nanotubes (SWCN), where the resonance state is at ωo=0.5THzwith 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 ω0varies 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 ω0the 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].

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Paolo Di Sia (October 24th 2012). Modelling at Nanoscale, Plasmonics - Principles and Applications, Ki Young Kim, IntechOpen, DOI: 10.5772/50755. Available from:

chapter statistics

2524total chapter downloads

1Crossref citations

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Computational Electromagnetics in Plasmonics

By Guy A. E. Vandenbosch

Related Book

First chapter

Dual-Periodic Photonic Crystal Structures

By Alexey Yamilov and Mark Herrera

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More about us