**Abstract.**
Raman scattering is a nondestructive technique that is able to supply information on the crystal and electronic structures, strain, temperature, phonon–phonon, and electron–phonon interaction. In the particular case of semiconductor nanowires, Raman scattering provides additional information related to surfaces. Although correct, a theoretical approach to analyze the surface optical modes loses critical information when retardation is neglected. A comparison of the retarded and unretarded approaches clarifies the role of the electric and magnetic polarization in the Raman selection rules. Since most III–V compounds growing in the zincblende phase change their crystal structure to wurtzite when growing as nanowires, the polariton description will be particularized for these two important crystal phases. Confined phonons exist in cylindrical nanowires and couple with longitudinal and transverse modes due to the presence of the nanowire’s surface. This coupling vanishes in the case of rotational symmetry. The boundary conditions of the electromagnetic fields on small-size nanowires (antenna effect) have a dramatic effect on the polarization properties of a Raman spectrum.

# Review on Raman scattering in semiconductor nanowires: I. theory

**Andrés Cantarero**

University of Valencia, PO Box 22085, 46071 Valencia, Spain

*J. Nanophoton*. 7(1), 071598 (Nov 21, 2013). doi:10.1117/1.JNP.7.071598

#### Open Access

Nanowires (NWs) are quasi-one-dimensional (1-D) materials where one of the dimensions, usually in the 10- to 100-nm range, is much smaller than the other two (in the range of a few microns). They were first called nanowhiskers, but they have received other names in the literature such as nanorods or nanocolumns. They can be made of metals, polymers, or semiconductors, and the most common cross-sections are nearly cylindrical, hexagonal, triangular, or squared. The NW geometry has several advantages as compared to quantum wells or quantum dots, the most elemental one being the large surface-to-volume ratio as compared to a thin film (additional contribution of the lateral surface). Depending on the density and cross-section, an NW bunch can have 100 times more surface than a thin film of the same dimensions. This fact can be useful in the development of solar cells or chemical sensors.^{1} On the other hand, new concepts of heterostructures arise from the NW shape, like core/shell or axial heterostructures (quantum dots in a wire). Core/shell nanostructures can be particularly useful in the development of solar cells,^{2} while axial heterostructures are suitable for quantum computing.^{3}

Semiconductor NWs can be grown by chemical vapor deposition, chemical beam epitaxy, or molecular beam epitaxy.^{4}^{–}^{6} Usually, they are grown using gold as a catalyst,^{7} but depending on the semiconductor, they can also be grown without a catalyst.^{5} They are commonly grown on a convenient substrate like Si or GaAs, which can be prepatterned in order to develop ordered NW arrays.^{1} Due to the existence of a free surface, heterostructures can be grown with larger lattice mismatch than superlattices or quantum wells.^{8} The reason is the relaxation of the structure at the sample surface.^{9} In very thin NWs, the fact that the bond distances close to the surface are larger due to relaxation effects produces changes in the electronic band structure of the material. NWs have a huge amount of applications—in photonic as light emitting diodes or lasers, solar cells, as energy harvesting materials,^{1}^{,}^{2} in transport as single-electron transistors,^{3}^{,}^{10} or in phototransport devices.^{11}^{,}^{12}

One more particularity of semiconductor NWs is the fact that the balance between surface energy and bulk energy gives rise to new crystal structures, not stable when grown in bulk systems. In particular, most of the III–V semiconductors crystallizing in the zincblende (ZB) structure in bulk develop the wurtzite (WZ) structure when growing in the form of NWs;^{7} actually, very often, ZB staking faults appear within the WZ structure. A recent experiment on InAs NWs photoluminescence^{13} shows the existence of two bands in the emission at low temperatures, presumably one with origin in the ZB region and a second from the WZ region. From this variation, an energy gap of 2.45 eV for the $E1$ gap of InAs has been obtained. This value is in very good agreement with resonant Raman measurements^{14} and *ab initio* calculations.^{15} High-resolution transmission electron microscopy images clearly show the existence of these two regions. The same behavior has been observed in other semiconductors, like InP.^{6} A detailed band structure calculation in the WZ phase has been performed in these two compounds.^{15}^{,}^{16}

Figure 1 shows a set of InN NWs grown by molecular beam epitaxy under different temperatures, nitrogen flux, and In vapor pressure conditions. We can observe how the shape of the NWs can be completely different depending on the growth conditions.

**F1 :**

Indium nitride NW morphologies obtained under varied metal flux and substrate temperature (thermocouple temperature using the $1\xd71\u2009\u2009cm2$ with $8.88\xd78.88\u2009\u2009mm2$ opening) and constant growth time of 300 min and nitrogen flux $FN=2\u2009sccm$ ($cc/min\u2061$ at 0°C at 1 atm of pressure) and a plasma excitation power $PN=500\u2009\u2009W$, which corresponds to a growth rate at stoichiometry of $\u223c12.5\u2009\u2009m/min$. Arrows indicate the modifications produced with the increase of the substrate temperature $Ts$. (The author thanks Dr. Christian Denker for providing some SEM images of his PhD dissertation.)

Finally, the electronic bands bend at the surface of the NWs as in any surface in a solid, but in the case of an NW the effect is more notable due to the large contribution of the surface compared to the inner part of the NW. Although the origin is not completely clear, some NW semiconductors have a depletion layer at the surface, some others an accumulation layer.^{5}^{,}^{17} The accumulation layer influences both the optical^{18} and transport properties.^{19}

In this first work, we will focus on the theoretical development of the Raman theory, first in bulk and then particularizing for the case of NWs. In a second work, we will pay more attention to the comparison of the existing theories with the available experimental results.

In a semiconductor, most of the incident radiation is either reflected, absorbed, or transmitted; only a small amount is scattered. The elastically scattered light is called Rayleigh scattering. Actually, the elastic scattering is possible due to the existence of defects in the sample; in a nearly perfectly grown material with translational symmetry, the elastic scattering must be very low.^{20}^{,}^{21} A tiny amount of the scattered light, in general several orders of magnitude smaller than that scattered elastically, corresponds to inelastic light scattering or Raman scattering.

The Raman effect or light scattering was predicted by Smekal in 1923 (^{22}) and measured simultaneously by Landsberg and Mandelstam,^{23} and by Raman and Krishnan^{24} in solids and liquids, respectively. Actually, the first set of authors discovered the effect previously, but the publication of the second set of authors appeared earlier in 1928. For this discovery, Raman was the Nobel laureate in physics in 1930. Raman was the only Third-World scientist who received the Nobel laureate for a research performed in an undeveloped country.

Although the Raman effect is understood usually as the scattering of light due to phonons, any kind of inelastic light scattering is called Raman scattering. The scattering by electronic excitations or magnetic excitations are also Raman scattering. In this review, we are going to limit the analysis to the Raman scattering by phonons. The classical theory of the Raman effect was developed by the same authors who discovered the effect, Mandelstam et al.,^{25} and the quantum theory by Tamm.^{26} A rigorous development of the classical theory can be found in the book of Born and Huang.^{27} We will follow here the development of Cantarero et al.^{28}^{,}^{29} and Trallero-Giner et al.^{30}

In the infrared, there is a direct coupling of the electric field of the light with the phonons, but in Raman scattering, the coupling is via virtual electronic states, as it was shown many years ago by Loudon.^{31} If the semiconductor gap is around the energy of the laser or scattered light, the transition between electronic states becomes real and the Raman signal increases orders of magnitude. In that case, we have to take into account the width of the electronic states, and the process is usually referred as resonant Raman scattering (RRS).^{28} Double and triple resonances (multiple resonances with several electronic levels) can also increase, additionally, the scattering signal.^{32} Electric fields,^{33} magnetic fields,^{34} or high pressure^{35} can also be tuned to produce RRS. Excitons play an important role in first-order RRS (^{28}) or Raman scattering by one phonon. In second-order RRS (scattering by two phonons at the $\Gamma $-point), excitons also play a fundamental role in the resonant behavior^{36} because the intermediate excitonic state can also resonate. Electron–phonon interaction is thus the crucial part of the resonant process, since there is a direct interaction between the electrons and the phonons, even out of resonance. Electron–phonon interaction is actually very important within the electronic transitions itself since it renormalizes the electronic gap of the materials^{37} and gives a broadening to the electronic states. Also, when the electrons have enough energy to emit optical phonons, the photoluminescence processes become more efficient since the emission probability increases. Electron–phonon interaction is important not only in optics but also in transport, through the different relaxation mechanisms.^{38}

The laser wavelength or laser energy is the other crucial part of the resonant process, since the variation of the laser energy approaching the electronic transition is the basis of the RRS. As an example, the Raman signal of strontium titanate/barium titanate (STO/BTO) quantum wells when illuminated with visible light was not good enough to even see the Raman signal. But ultraviolet light has been shown to be very efficient for analyzing STO/BTO quantum wells constituted even by a single layer^{39} due to resonance effects. In single wires, since there is a small amount of material, resonant effects are very important^{15} to make the Raman effect observable.

In nonpolar semiconductors, there is no splitting between the the transverse (TO) and longitudinal (LO) optical phonons, but in polar semiconductors like GaAs, there is a splitting related to the difference between the low- and high-frequency dielectric permittivities (or the effective charges). The splitting is a measure of the polarity of the semiconductor. In polar semiconductors, plasmons interact with LO phonons giving rise to an anticrossing when the plasma frequency coincides with the LO phonon frequency, which depends on the electron concentration.^{40} From the measurement of the two split frequencies, the electron concentration can be obtained.^{40} In heavily doped samples, the electronic background is of the order of the Raman signal; electrons interact strongly with the phonons giving rise to the Fano lineshape,^{40}^{,}^{41} produced by the sum of quantum probabilities of the two processes, which can have opposite signs. Since phonons depend on the mass of the ions, Raman scattering has been used in the last years to study the isotopic composition of isotopically engineered samples^{42}^{–}^{44} and phonon anharmonicities.^{45} Phonons are also sensitive to pressure or stress variations.^{46}^{,}^{47} In thin films and semiconductor nanostructures, Raman scattering is used very often to obtain the lattice mismatch through the phonon deformation potentials.^{48}

Recently, with the advances in the fabrication of holographic filters (notch or edge filters), a single spectrometer is enough to study the Raman scattering in most of the materials. A good holographic filter has a cutoff of $\u223c100\u2009\u2009cm2$, or even less. There are, in the market, several companies selling compact instruments, smaller than a shoebox, completely equipped connectable to a portable computer capable of obtaining reasonable Raman scattering signals at a competitive price. If we need to study phonon modes at lower frequencies, a double or triple spectrometer, in subtractive mode, is needed. Also, if we need to increase the resolution because the phonon lines are narrow, we can use a triple spectrometer in additive mode and lower the temperature to decrease the phonon anharmonicities. Most of the Raman systems are equipped with a charged coupled device, but if we really need resolution, we have to go back to the photomultipliers. In the book of Kuzmany,^{49} there is a complete description of the apparatus used in Raman scattering: spectrometers, gratings, and detectors.

During the last decades, part of the scientific community has been dedicated to characterize, by Raman scattering, thin films, quantum wells, quantum dots, quantum wires, or even more complex nanostructures.^{50} In the simple case of a thin film, if we use the typical collection, a photographic camera objective, producing a laser spot $\u223c10\u2009\u2009\mu m$, most of the Raman signal will come from the substrate and it can mask the really important signal. We need, in that case, the use of a microscope (micro-Raman or μ-Raman system), where the laser spot can be $<1\u2009\u2009\mu m$ (with a confocal microscope, we can go to the diffraction limit^{50}). However, the use of a microscope has a couple of disadvantages. First, all the light, including the reflected light, goes through the entrance slit of the spectrometer (unless we use a notch or edge filter). The second disadvantage is that the selection rules are not completely fulfilled. If we have an objective of $100\xd7$ and a numerical aperture^{50}$NA=0.9$, such as that shown in Fig. 2(a), the entrance angle with a refractive index of 3 (a typical value for semiconductors) will be $\u223c17.5\u2009\u2009deg$. In a typical. μ-Raman system, a confocal microscope is employed in order to perform the measurements on a specific point of the sample (below 1 *μ*m in lateral resolution and 100 nm in depth). A confocal microscope and a motorized xyz plate allow the mapping of small areas at different depths. In Fig. 2(b), there is a scheme of a confocal system where we observe how the light enters into a sample with a smaller angle due to the difference of refractive indices. Although we do not have a complete backscattering geometry, the selection rules are fulfilled to a great extent due to the difference in refractive index.

**F2 :**

(a) Cross-section of an objective (by courtesy of Carl Zeiss Microscopy, Germany). 1. Objective thread. 2. Stop face of the objective. 3. Spring system for the specimen-protection mechanism. 4. to 7. Lens groups for the correction of image errors. 8. Correction collar for adapting to deviating cover glass thicknesses or temperatures. 9. Front lens system. 10. Front lens holder. (b) Schematic representation of a confocal setup, where we can observe the angle change inside the sample due to the differences in refractive indices.

In order to compare the Raman signal or the Raman intensity of different materials, it is convenient to use the Raman scattering efficiency (RSE), instead of the scattering cross-section, since the first amount is volume independent. In a polar binary compound semiconductor, the RSE per unit solid angle $\Omega $ can be defined as^{28}

*μ*the reduced mass of the ions, $vc$ the volume of the unit cell, $\omega L$ ($\omega S$) the laser (scattering) angular frequency, $nL$ ($nS$) the refractive index at the laser (scattered) frequency, $\omega \Gamma $ the angular frequency of the corresponding phonon, $eL$ ($eS$) the polarization vector of the laser (scattered) light, and $R\Gamma $ the Raman tensor corresponding to the $\Gamma $ phonon mode. The factor $N(\omega \Gamma )+1$ has its origin in a quantum mechanical description of the phonon field,

^{28}and it corresponds to the phonon emission (Stokes shift). In the case of phonon absorption (anti-Stokes shift or anti-Stokes scattering), the factor is $N(\omega \Gamma )$. In practice, this factor is replaced by the Bose–Einstein distribution function, and this is the reason why at the anti-Stokes frequency, the Raman scattering signal is very weak as compared to that appearing at the Stokes shift, unless the temperature is very high. In a classical description, both Stokes and anti-Stokes components have the same intensity. An exception is stimulated Raman scattering, where phonons at the anti-Stokes frequency

^{50}are excited. For Stokes shift, even if $N(\omega \Gamma )\u223c0$, the factor 1 remains and the Stokes shift can be observed at any temperature. The relative intensity of the Stokes/anti-Stokes Raman shift can be used to monitor the temperature.

^{51}

The Raman tensor is a second-order tensor, sometimes written as a third-order tensor (or even as a fourth-order tensor^{52}) contracted with the phonon polarization ($Rijkq^k$, where $q^$ is a unitary vector in the direction of the phonon displacement). The form of the Raman tensor depends on the symmetry of the crystal and the electron–phonon interaction, since in the Raman effect, electrons play a crucial role as virtual intermediate states in the Raman process.^{31} For an LO phonon propagating along $x$ ([1,0,0]) direction, the Raman tensor for the optical modes of a ZB crystal can be written in terms of the Raman polarizability $\alpha $ as

^{53}Actually, we can write $Ryzx=Rzyx=\alpha $. The remaining nonzero Raman tensor components are $Rxzy=Rzxy=Rxyz=Ryxz=\alpha $. We will come back to this point in Sec. 4. The Raman polarizability has the expression

^{54}

^{31}are nonresonant, and they can be neglected. $E\mu $ ($E\nu $) and $\gamma \mu $ ($\gamma \nu $) are the energy and broadening of the electronic states (which can be accessed in a time-resolved Raman experiment

^{55}). If there are several initial and/or final states (as it is the case of several valence bands in a ZB or WZ compound), we have to sum over $I,F$ after squares $WFI$. The sum over intermediate states is performed before squaring the expression. $HeR$ and $HeL$ are the electron-radiation and electron(hole)-phonon Hamiltonians, which can be written in the second quantization formalism.

^{28}The important point is the electron–phonon coupling constant, which consists of two parts, one corresponding to the emission of the phonon by the electron (see Feynmann diagram in Fig. 3) and a second corresponding to the emission of the phonon by the hole. In ZB- (and WZ-) type materials, only the second option is nonzero, and the electron–phonon coupling constant is

^{56}In the case of the Fröhlich interaction,

^{30}the basic expressions are the same; the only change is in the electron–phonon interaction Hamiltonian. Its coupling constant, in that case, is

^{30}

It is interesting to realize that the DP is a local interaction, and thus it is related to the periodic part of the Bloch function, while the Fröhlich interaction, only possible in polar materials, is a long-range interaction related to the polarization charges in the medium, and only the envelope of the Bloch function contributes to the matrix element.^{30} On the other hand, since both processes are possible in the case of LO phonons, the amplitude probability must be added before squaring it, and interference effects can appear^{29} when the Raman polarizabilities corresponding to deformation potential (DP) and Fröhlich interaction are of the same order of magnitude, i.e., close to an electronic transition, in resonance. This is the reason we usually say that the selection rules are not completely fulfilled in resonance. During many years, the huge enhancement of Raman scattering close to the gap could not be explained, and it was attributed to unknown values of one of the constants that appear in the expression of RSE or to the existence of impurity-induced Raman scattering.^{57} Finally, it was demonstrated that the consideration of excitons, instead of electrons,^{31} as intermediate states in the Raman process could explain the absolute value of Raman scattering in the case of deformation potential,^{28} Fröhlich interaction,^{30} and also the interference between both.^{29} The interference is produced due to the fact that we have to add the probabilities before squaring (they correspond to different intermediate steps).

In Sec. 3, the Raman tensor of an LO phonon propagating along the [100] direction of the ZB structure was defined in Eq. (2). We also defined the remaining components of the Raman tensor using the conventional unitary vectors along $x$, $y$, and $z$. If one is interested in a different crystal surface, a general method to proceed is to define the following three vectors:^{58}

It is also convenient to define the matrix formed by these three vectors written in rows.

Following Fig. 4, the $a1$ axis will be always along the $x$ direction, i.e., the incoming wave vector of the light $\kappa L$ and outgoing $\kappa S$ are along $a1$. The phonon wave vector $q$, in backscattering configuration, will also be directed along $a1$. Thus, the LO phonon will vibrate in the $a1$ direction, while the two TO phonons will vibrate along the $a2$ and $a3$ directions. The three vectors are unitary. As commented in Sec. 3, the Raman intensity will be proportional to

**F4 :**

Fix axes system, where we have chosen $z$ as the direction of the cylinder axis, $x$ as the direction of the light in backscattering configuration, and $y$ as the perpendicular direction, forming a right-handed system. The vectors $eL$ and $eS$ indicate the polarization of the incoming and scattered light, respectively, and $\theta $ and $\varphi $ the corresponding angles of rotation, with origin in the $z$ axis with counterclockwise rotation.

In the case of an LO phonon, the Raman intensity will be proportional to

The reason of summing after squaring is that these are two different final states, but since they appear at the same energy when looking at the Stokes shift, we have to add the intensities. Of course, depending on the selection rules for a given surface, only one of them can contribute. In the case of the ZB structure, it is very simple to calculate the resulting intensity for a $(11\xaf0)$ surface, for instance. In that case,

If the polarization vectors $eL$ and $eS$ are both along the (110) direction, the transformed vectors

The Raman matrices multiplied by the phonon polarization, $Ri\u2032ai$,

^{52}

Although we have given a general procedure, we will pay particular attention to the case where one of the axis is along the [111] direction, the NW growth direction or the cylinder axis. The reason is that this is the equivalent to the [0001] direction of the WZ structure and the results can somehow be compared. Figure 4 shows the coordinate system where the NW axis is along the laboratory $z$ direction. In that case, it is better to select $A$ as

Polar plots can be drawn by considering, in the case of $(11\xaf0)$ surfaces, in agreement with Fig. 4,

In the WZ structure, there are nine optical modes, but two of them are silent modes. The seven Raman active optical modes are drawn in Fig. 5. The first mode vibrates along the $z$ direction of the WZ structure. The Raman tensor is

The 1-D $A1$ mode is a polar mode, with a different polarizability when propagating in the $z$ direction or in the $xy$ plane. As stated before, group theory gives only the difference in anisotropy ($a$ and $b$ polarizabilities), but does not consider the $LO\u2212TO$ splitting, which is due to the polarity of the medium. The two-dimensional (2-D) $E1$ modes are also polar modes, with Raman tensor

When propagating along $x$, $E1(x)$ will be LO, while it will be TO if it propagates along $y$. There are actually two bidimensional $E1$ modes, but one of them corresponds to two acoustic modes, the third acoustical mode being a 1-D $A1$ mode. There are also two nonpolar bidimensional $E2$ modes, as can be observed in Fig. 5. Both modes actually mix, but the amplitude of the heavier (lighter) atom dominates in the $E2l$ ($E2h$) modes, vibrating at lower (higher) frequency. The corresponding Raman tensor can be written as

The WZ is an anisotropic structure. At this point, it is interesting to remark that there are two kinds of anisotropic structures. In the first one, the crystal anisotropy dominates;^{59}^{,}^{60} actually $\epsilon \u2225$ is very different from $\epsilon \u22a5$. This is the case of most of the layered compounds,^{46} or nonpolar compounds like graphite, but the most common case is that where the Coulomb interaction dominates.^{35} This is the case of III–V compounds grown in the WZ structure. In recent cases where typical ZB crystals grow in the WZ structure, the material anisotropy is unkonwn, but, on one side, we know the $LO\u2212TO$ splitting, i.e., the intensity of the Coulomb interaction, and on the other side, a small difference is expected in the crystal anisotropy since the structures are identical until second neighbors. This is very important when we want to draw polar plots in the case of WZ materials. When we rotate the polarization, the polar modes are mixed.

Instead in a strongly anisotropic material, $A1(TO)$ and $A1(LO)$ [or $E1(TO)$ and $E1(LO)$] mix. Since we are interested in crystals with a small anisotropy, the Raman polarizabilities $a$ and $b$ will be taken as equal to simplify the final expressions, i.e.,

With this approximation, we can easily obtain the Raman selection rules in backscattering from the lateral side of a single WZ-NW grown along the c-axis. Looking at the picture given by Sander et al.,^{61} the $x$ axis is perpendicular to the $m$-plane $(101\xaf0)$ [$(100)$ plane in three-index notation], while the $y$ axis corresponds to the $a$-plane $(112\xaf0)$ [(110) in three-index notation]. The $z$ axis is along $c$, the [0001] direction ([111] direction of the ZB).

The selection rules for a few standard configurations are summarized in Table 1. If one is interested in polar plots, we have to calculate the components of the Raman tensor. Let us do, as an example, the scattering configuration $x(\zeta (\varphi )\zeta (\varphi ))x\xaf$ where $\zeta $ is in the zy plane. The vectors $eL\u2032=(0,0,1)$ and $eS\u2032=(0,\u2212sin\u2009\varphi ,cos\u2009\varphi )$.

Scattering configuration | Raman component |
---|---|

$x(yy)x\xaf$ | $|aTO|2$, $|d|2$ |

$x(zz)x\xaf$ | $|aTO|2$ |

$x(yz)x\xaf$ | $|cTO|2$ |

$z(xx)z\xaf$ | $|aLO|2$, $|d|2$ |

$z(yy)z\xaf$ | $|aLO|2$, $|d|2$ |

$z(xy)z\xaf$ | $|d|2$ |

The Raman intensities will be, in the $x(\zeta \zeta )x\xaf$ configuration,

The number of phonon modes in a crystal depends on the number of atoms in the unit cell. As is well known from a general course of solid-state physics, if a unit cell (primitive cell) has $n$ atoms, there are $3n\u22123$ optical modes. Depending on the complexity of the crystal, we can have from a few to tens of optical modes. For instance, in silicon, there are three optical modes, as in ZB-type materials; in WZ materials, there are nine optical modes (seven are Raman active), while in a more complex material, like $Bi2S3$ and isostructural layered materials, there are 57 optical phonons, 27 of them Raman active.^{62} In the case of an NW, there is no translational symmetry in the plane perpendicular to the NW axis. If the NW is very thin (a few nanometers of radius), we have confined modes (see Sec. 7). With increasing NW radius increases, the confined phonon modes overlap due to the finite linewidth, giving rise to an asymmetry at the low frequency region of the LO or TO phonon peaks. NWs of the order of several tens of nanometers can be considered, for most of the purposes, as a bulk material. Just as an example, in transport phonons can be treated as bulk materials even for NW diameters down to 30 to 40 nm.^{63} The only difference will come from the existence of a surface or an interface (core/shell NWs), giving rise to new phonon modes, the so-called surface optical (SO) modes.

Let us start our discussion analyzing bulk phonons (i.e., there is no surface) in a solid. From group theory (at $q=0$), there is no splitting of the optical modes; LO and TO modes have the same energy. Actually, the $LO\u2212TO$ splitting is produced at finite $q$ in polar semiconductors due to the coupling of the phonon field with the electromagnetic field of the light (phonon–polariton). In bulk materials, the polariton branches are derived from the solution of the Maxwell equations inside the material and the coupling of electric field produced by the phonon displacement with the electric field of the electromagnetic wave.^{64} From the LO to the TO frequencies, there is a forbidden region where the electromagnetic wave cannot propagate in bulk materials (the electromagnetic field is totally reflected). Actually, because of the losses (imaginary part of the dielectric function), the reflectivity is not exactly one in this region.

Although the coupling between the ionic displacement in a polar crystal and the electric field of the electromagnetic wave of the light is produced directly in the infrared region of the spectrum, it can be analyzed by Brillouin scattering,^{65} and this coupling is responsible for the existence of surface modes, as we will see later. In a Raman experiment in backscattering geometry, $q\u22484\pi nL/\lambda L$, where $nL$ and $\lambda L$ are the refractive index and wavelength of the laser light, respectively. For the typical green line of an Ar laser ($\lambda L=514.5308\u2009\u2009nm$), $q\u22480.08\u2009\u2009nm\u22121$. In a Brillouin scattering experiment, the complete measurement range is below this amount.

If we consider the harmonic movement corresponding to a phonon traveling along the crystal, the equation of motion will be that corresponding to the harmonic oscillator. On the other hand, if there is an electromagnetic wave in the polar crystal, there is a polarization vector proportional to the electric field, the proportionality constant being the electric susceptibility at a microscopic scale. But these fields are coupled through the ionic movement. Let $u=u+\u2212u\u2212$ be the vector displacement between the positive and negative ion in a simple polar semiconductor (as the ZB, for instance) and $w=\mu /vcu$, ($\mu /vc$ is basically the density $\rho M$, the mass density). The equations to be solved in order to obtain both the electric polarization and the vector displacement due to the existence of the external field $E$ are

Equation (28) must fulfill at the same time the Maxwell equations, which we write in the Gaussian system, for simplicity, and assume a plane wave solution $ei(k\xb7r\u2212\omega t)$,

Combining these equations, we arrive to the solution

At very high frequencies, the ionic movement cannot follow the electric field, only the electrons, and we arrive to the trivial solution

On the other hand, uncoupled transverse phonons do not carry electric field and thus vibrate at their resonant frequency $\omega TO$, thus $b11=\u2212\omega TO2$. Finally, in static conditions ($\u22022w/\u2202t2=0$), we obtain from the coupled Eq. (28)

The final expressions coupling the fields and atomic movement are

Taking the displacement from Eq. (33), substituting in Eq. (34), and using Eq. (30) in the particular case of transverse modes ($P\xb7k=0$), we arrive at the following equation:

Using Eq. (35), we arrive at the solution

The solution of Eq. (36) is given in Fig. 6 using the data of GaAs.^{42} At $k=0$, there are two solutions obtained from the last equation taking the limit $k\u21920$: $\omega =\omega LO$ and $\omega =(c/\epsilon 0)k$. For large values of $k$, we obtain, taking the limit $k\u2192\u221e$, $\omega =\omega TO$ and $\omega =(c/\epsilon \u221e)k$. The longitudinal solutions are obtained assuming $P\xb7k=Pk\u21d2E=\u22124\pi P\u21d2D=0$. From the coupled equations, $\omega =\omega LO$. In Sec. 6.1, we will see how these equations are modified due to the existence of a limiting surface in the solid. The value of the electric field can be recovered from Eq. (28). In the case of TO modes, Eq. (28) gives $E=0$, and for LO modes,

**F6 :**

Lower and upper polariton branch in bulk GaAs calculated from Eq. (36) using the data of GaAs.^{42}

We can follow the same strategy to analyze the polariton dispersion in a crystal with the WZ structure. There are two modes, which are Raman and infrared active, $A1$ (a 1-D mode vibrating along the $c$ axis of the WZ structure) and $E1$ (a 2-D mode vibrating in the plane perpendicular to the $c$ axis). In a uniaxial crystal, the dielectric function is a tensor, isotropic in the plane perpendicular to $c$. Let us write

We will restrict ourselves to the case where the phonons propagate along the $z$ axis and the $x$ axis ($k$ is along $z$ or $x$—the $y$ direction is equivalent to $x$). If $k\u2225z^$, we have a $A1(LO)$ mode and two degenerate $E1(TO)$ modes (vibration along $x$ and $y$). If the phonons propagate along the $x$ (or $y$) direction, we will have the $A1(TO)$ mode vibrating along $z$ and propagating along $x$ and two $E1$ modes, a mode vibrating along $x$ and propagating along $x$, the $E1(LO)$, and a second mode vibrating along $y$ and propagating along $x$, the $E1(TO)$. Obviously, these two modes are not degenerate in that case.

For the case where $k\u2225z^$, the electric field can be along the $x$ axis (or $y$). In that case, Eq. (30) becomes

The equations coupling the fields with the atomic movement are

The same solution is obtained if the fields and ions vibrate in the $y$ plane [exchanging $x$ by $y$ in Eqs. (39) and (40)]. The solutions are the dispersion of the two $E1(TO)$ degenerate modes. By solving the three equations, we arrive at a solution equivalent to that of Eq. (36), replacing $\omega TO$ by $\omega \u2225$ and $\epsilon $ by $\epsilon \u2225$. The solution can be better written in an equivalent way as Eq. (30).

The limits $k\u21920$ and $k\u2192\u221e$ are easily extracted from these equations, $\omega =(c/\epsilon 0\u2225)k$ and $\omega =(c/\epsilon \u221e\u2225)k$, respectively (see Fig. 7).

**F7 :**

Lower and upper polariton branch in bulk ZnO calculated from Eq. (36) using the data of ZnO.^{66}

Here the electric field and the movement is along $z$; the solution is $Ez=\u22124\pi Pz$ ($D=0$), corresponding to the $A1(LO)$ with frequency $\omega =\epsilon 0\u2225/\epsilon \u221e\u2225\omega \u2225$. These solutions are plotted in Fig. 7 for ZnO. If $k$ is along the $x$ axis, the $E1$ mode propagating along the $x$ direction will be longitudinal, with frequency $\omega =\epsilon 0\u22a5/\epsilon \u221e\u22a5\omega \u22a5$. The dispersion of the modes propagating along $y$ and $z$ are obtained by solving the equations

In the case of Raman scattering, the value of $q$ is far away from the polariton region. In backscattering configuration and assuming that the selection rules are fulfilled (that depends on the cylinder radius), $q=kL\u2212kS$. For $\lambda L=514.53\u2009\u2009nm$, the laser frequency in $cm\u22121$ is $\omega L=1/\lambda L=19,435.213\u2009\u2009cm\u22121$. The scattered frequency, for instance, for Si is $\omega S=19,435.213\u2212519.5=19,383.263\u2009\u2009cm\u22121$ (515.909 nm). The value of $q$ will be $q=2\pi nL/\lambda L+2\pi nS/\lambda S\u22484\pi nL/\lambda L\u22480.0855\u2009\u2009nm\u22121=855,000\u2009\u2009cm\u22121$ (in the polariton region, $k\u226415,000\u2009\u2009cm\u22121$).

In finite-size samples, there are additional conditions on the electric and the magnetic fields since they must fulfill the boundary conditions at the sample surface. The existence of these boundary conditions generates additional solutions in the polariton equations. In general, there is a complete set of solutions ranging from below the TO to above the LO region, depending on the geometry.

Let us consider a single cylindrical semiconductor NW of radius $a$, with a dielectric permittivity $\epsilon (\omega )$, and an electromagnetic field with time dependence $e\u2212i\omega t$. If the external medium has a dielectric permittivity $\epsilon M$ (in case of air or vacuum, $\epsilon M=1$), the electric and magnetic fields outside and inside of the NW must fulfill the Helmholtz equations.

^{67}in particular, in cylindrical coordinates. In case of a hexagonal cross-section, the main conclusions from this section will not vary, but the calculation of the electromagnetic field must be done numerically due to the complexity of the boundary conditions.

Let us call $k2=\epsilon (\omega )\omega 2/c2$ and $k02=\epsilon M\omega 2/c2$. Equations (44), in cylindrical coordinates, become, for a scalar function $\psi $,

^{68}