International
Tables for
Crystallography
Volume D
Physical properties of crystals
Edited by A. Authier

International Tables for Crystallography (2013). Vol. D, ch. 2.3, pp. 334-348
doi: 10.1107/97809553602060000913

Chapter 2.3. Raman scattering

I. Gregoraa*

aInstitute of Physics, Academy of Sciences of the Czech Republic, Na Slovance 2, CZ-18221 Prague 8, Czech Republic
Correspondence e-mail: gregora@fzu.cz

Raman scattering is due to the inelastic scattering of photons by excitations in molecules, liquids or solids. This chapter deals specifically with the inelastic scattering of photons by lattice vibrations (phonons) in crystals, but the treatment is also applicable to other types of collective excitations in a crystal. The basic notions are introduced in the first section. First-order scattering by phonons is then described. A table summarizes the symmetry properties of Raman tensors in the 32 crystallographic classes. Further sections discuss morphic effects in Raman scattering, spatial-dispersion effects and higher-order scattering.

2.3.1. Introduction

| top | pdf |

The term Raman scattering, traditionally used for light scattering by molecular vibrations or optical lattice vibrations in crystals, is often applied in a general sense to a vast variety of phenomena of inelastic scattering of photons by various excitations in molecules, solids or liquids. In crystals these excitations may be collective (phonons, plasmons, polaritons, magnons) or single-particle (electrons, electron–hole pairs, vibrational and electronic excitation of impurities). Raman scattering provides an important tool for the study of the properties of these excitations. In the present chapter, we shall briefly review the general features of Raman scattering in perfect crystals on a phenomenological basis, paying special attention to the consequences of the crystal symmetry. Our focus will be mainly on Raman scattering by vibrational excitations of the crystal lattice – phonons. Nevertheless, most of the conclusions have general validity and may be (with possible minor modifications) transferred also to inelastic scattering by other excitations.

2.3.2. Inelastic light scattering in crystals – basic notions

| top | pdf |

Although quantum concepts must be used in any complete theory of inelastic scattering, basic insight into the problem may be obtained from a semiclassical treatment. In classical terms, the origin of inelastically scattered light in solids should be seen in the modulation of the dielectric susceptibility of a solid by elementary excitations. The exciting light polarizes the solid and the polarization induced via the modulated part of the susceptibility is re-radiated at differently shifted frequencies. Thus inelastic scattering of light by the temporal and spatial fluctuations of the dielectric susceptibility that are induced by elementary excitations provides information about the symmetry and wavevector-dependent frequencies of the excitations themselves as well as about their interaction with electromagnetic waves.

2.3.2.1. Kinematics

| top | pdf |

Let us consider the incident electromagnetic radiation, the scattered electromagnetic radiation and the elementary excitation to be described by plane waves. The incident radiation is characterized by frequency [\omega{_I}], wavevector [{\bf k}{_I}] and polarization vector [{\bf e}{_I}]. Likewise, the scattered radiation is characterized by [{\omega}{_S}], [{\bf k}{_S}] and [{\bf e}{_S}]: [{\bf E}_{I,S}({\bf r},t) = E_{I,S}{\bf e}_{I,S}\exp(i{\bf k}_{I,S}{\bf r}-{\omega}t). \eqno (2.3.2.1)]

The scattering process involves the annihilation of the incident photon, the emission or annihilation of one or more quanta of elementary excitations and the emission of a scattered photon. The scattering is characterised by a scattering frequency [{\omega}] (also termed the Raman shift) corresponding to the energy transfer [\hbar \omega] from the radiation field to the crystal, and by a scattering wavevector q corresponding to the respective momentum transfer [\hbar {\bf q}]. Since the energy and momentum must be conserved in the scattering process, we have the conditions[\eqalignno{{\omega}_I - {\omega}_S &= {\omega}, &\cr {\bf k}{_I} - {\bf k}{_S} & = {\bf q}.&(2.3.2.2)}]Strictly speaking, the momentum conservation condition is valid only for sufficiently large, perfectly periodic crystals. It is further assumed that there is no significant absorption of the incident and scattered light beams, so that the wavevectors may be considered real quantities.

Since the photon wavevectors ([{\bf k}{_I}], [{\bf k}{_S}]) and frequencies ([{\omega}{_ I}], [{\omega}{_S}]) are related by the dispersion relation [{\omega}= ck/n], where c is the speed of light in free space and n is the refractive index of the medium at the respective frequency, the energy and wavevector conservation conditions imply for the magnitude of the scattering wavevector q[c^2 q^2 = n_I^2 \omega _I^2 + n_S^2 (\omega _I - \omega)^2 - 2n_I n_S \omega _I (\omega _I - \omega)\cos \varphi, \eqno (2.3.2.3)]where [{\varphi}] is the scattering angle (the angle between [{\bf k}{_I}] and [{\bf k}{_S}]). This relation defines in the ([\omega, q]) plane the region of wavevectors and frequencies accessible to the scattering. This relation is particularly important for scattering by excitations whose frequencies depend markedly on the scattering wavevector (e.g. acoustic phonons, polaritons etc.).

2.3.2.2. Cross section

| top | pdf |

In the absence of any excitations, the incident field [{\bf E}{_I}] at frequency [\omega_I] induces in the crystal the polarization P, related to the field by the linear dielectric susceptibility tensor [{\chi}] ([\varepsilon {_0}] is the permittivity of free space): [{\bf P} = \varepsilon {_0} {\bf \chi}(\omega {_I}) {\bf E}{_I}. \eqno (2.3.2.4)]The linear susceptibility [ \chi(\omega{_I})] is understood to be independent of position, depending on the crystal characteristics and on the frequency of the radiation field only. In the realm of nonlinear optics, additional terms of higher order in the fields may be considered; they are expressed through the respective nonlinear susceptibilities.

The effect of the excitations is to modulate the wavefunctions and the energy levels of the medium, and can be represented macroscopically as an additional contribution to the linear susceptibility. Treating this modulation as a perturbation, the resulting contribution to the susceptibility tensor, the so-called transition susceptibility [{\delta}{ \chi}] can be expressed as a Taylor expansion in terms of normal coordinates [Q{_j}] of the excitations:[\chi \rightarrow \chi + \delta \chi, \hbox{ where }\delta \chi = \textstyle\sum\limits_j {\chi ^{(j)}Q_j + \textstyle\sum\limits_{j,j'}{\chi ^{(j,j')}Q_j Q_{j'} + \ldots}}. \eqno (2.3.2.5)]The tensorial coefficients [\chi ^{(j)},\chi ^{(j,j')},\ldots] in this expansion are, in a sense, higher-order susceptibilities and are often referred to as Raman tensors (of the first, second and higher orders). They are obviously related to susceptibility derivatives with respect to the normal coordinates of the excitations. The time-dependent polarization induced by [{\delta}{\chi}] via time dependence of the normal coordinates can be regarded as the source of the inelastically scattered radiation.

The central quantity in the description of Raman scattering is the spectral differential cross section, defined as the relative rate of energy loss from the incident beam (frequency [\omega{_I}], polarization [{\bf e}{_I}]) as a result of its scattering (frequency [\omega{_S}], polarization [{\bf e}{_S}]) in volume V into a unit solid angle and unit frequency interval. The corresponding formula may be concisely written as (see e.g. Hayes & Loudon, 1978[link]) [{{{\rm d}^2 \sigma }\over {{\rm d}\Omega \,\,{\rm d}\omega _S }}= {{\omega _S^3 \omega _I V^2 n_S }\over {(4\pi)^2 c^4 n_I }}\left\langle {\left| {{\bf e}_I \delta \chi {\bf e}_S }\right|^2 }\right\rangle _\omega. \eqno (2.3.2.6)]The symbol [\left\langle {\ldots}\right\rangle _\omega ] stands for the power spectrum (correlation function) of the transition susceptibility fluctuations. The spectral differential cross section is the quantity that can be directly measured in a Raman scattering experiment by analysing the frequency spectrum of the light scattered into a certain direction. By integrating over frequencies [\omega{_S}] for a particular Raman band and, in addition, over the solid angle, one obtains, respectively, the differential cross section ([{\rm d}\sigma/{\rm d}\Omega]) and the total cross section ([\sigma_{\rm tot}]): [{{{\rm d}\sigma }\over {{\rm d}\Omega }} = \int {\left({{{{\rm d}^2 \sigma }\over {{\rm d}\Omega \,\,{\rm d}\omega _S }}}\right)}{\rm d}\omega _S, \quad \sigma _{\rm tot} = \int {\left({{{{\rm d}\sigma }\over {{\rm d}\Omega }}}\right){\rm d}\Omega }.]These quantities are useful in comparing the integrated scattered intensity by different excitations.

2.3.2.3. Experimental aspects

| top | pdf |

In a scattering experiment on crystals, the choice of the scattering geometry implies setting the propagation directions [{\bf k}{_I}] and [{\bf k}{_S}] and the polarization of the incident and scattered light with respect to the crystallographic axes and defining thus the direction of the scattering wavevector q as well as the particular component (or a combination of components) of the transition susceptibility tensor [{\delta\chi}]. In practice, the incident radiation is almost exclusively produced by a suitable laser source, which yields a monochromatic, polarized narrow beam, with a well defined wavevector [{\bf k}{_I}]. The light scattered in the direction of [{\bf k}{_S}] is collected over a certain finite solid angle [\Delta\Omega]. Its polarization is analysed with a suitable polarization analyser, and the scattered intensity as a function of frequency [\omega{_S}] (or Raman shift [{\omega}]) is analysed using a spectrometer.

To characterize the Raman scattering geometry in a particular experimental arrangement, standard notation for the scattering geometry is often used, giving the orientation of the wavevectors and polarization vectors with respect to a reference Cartesian coordinate system, namely: [{\bf k}{_I}({\bf e}{_I}, {\bf e}{_S}){\bf k}{_S}]. Thus, for example, the symbol [x(zy)z] means that right-angle scattering geometry is used, where the incident beam polarized in the [z=[001]] direction propagates along the [x=[100]] axis, while the scattered beam is collected in the z direction and the polarization analyser is set parallel with the [y = [010]] direction. The measured intensity, being proportional to [|\chi_{zy}|{^2}], gives information on this particular component of the transition susceptibility tensor. By virtue of the momentum conservation, the scattering wavevector q in this case is oriented along the [[10\bar 1]] direction.

In a typical Raman experiment with visible light ([\omega \ll \omega{_I} \approx\omega{_S}]), the magnitudes of the wavevectors [ k{_I}\approx k{_s}= k] are of the order of 105 cm−1, much lower than those of the reciprocal-lattice vectors K ([\approx 10^{8}] cm−1). Consequently, the range of the magnitudes of the scattering wavevectors q accessible by varying the scattering geometry from [{\varphi}= 0]° (forward scattering) to [{\varphi}=180]° (back scattering) is [0 \leq q \leq 2k], i.e. by about three orders of magnitude lower than the usual dimensions of the Brillouin zone. The use of back-scattering geometry is imperative in the case of opaque samples, which show stronger absorption for the exciting (or scattered) light.

It should be noted that the general formula for the spectral differential cross section (2.3.2.5)[link] applies to the situation inside the crystal. Since in real experiments the observer is always outside the crystal, several corrections have to be taken into account. These are in particular due to refraction, reflection and transmission of the incident and scattered light at the interfaces, as well as absorption of light in the crystal. Attention must be paid in the case of anisotropic or gyrotropic crystals, where birefringence or rotation of the polarization direction of both incident and scattered light may occur on their paths through the crystal, between the interfaces and the scattering volume.

We conclude this section by remarking that, owing to the obvious difficulties in taking all the properties of the experimental setup and the corrections into consideration, measurements of absolute Raman intensities tend to be extremely rare. There exist, however, several crystals for which absolute determination of the cross section for particular excitations has been made with reasonable reliability and which may serve as secondary standards.

2.3.3. First-order scattering by phonons

| top | pdf |

In what follows, we shall be more specific and by underlying excitations we shall explicitly understand lattice vibrationsphonons – although the treatment is also applicable to other types of collective excitations in a crystal.

Let us recall (see Chapter 2.1[link] ) that atomic displacements in the crystal lattice can be expressed as linear combinations of the normal modes of vibrations – eigenvectors of the dynamical matrix. [{\bf u}_{\kappa {\bf l}} = \sum\limits_{\bf q}\sum\limits_j {\bf u}_{\kappa {\bf l}}({\bf q},j) = {1 \over {\sqrt{Nm_{\kappa}}}}\sum\limits_{\bf q}\sum\limits_j Q_j ({\bf q}){\bf e}_\kappa ({\bf q},j)\exp (i{\bf qr}_{\bf l}), \eqno (2.3.3.1)]where l denotes the primitive unit cell, [{\kappa}] ([{\kappa}=1,\ldots,s]) is the index of the atom (mass [m{_\kappa}]) in the unit cell and N is the number of unit cells in the crystal. The eigenvectors [{\bf e}{_\kappa}({\bf q},j)], also called polarization vectors, are normalized vectors describing the displacement pattern of atoms in the reference unit cell in a jth normal mode of vibration ([j = 1, 2, 3,\ldots, 3s]) with a wavevector q. The normal modes can be classified according to irreducible representations of the crystal space group and labelled correspondingly. In the case of degeneracy, the branch index j stands, in fact, for two indices: one for the irreducible representation, and the other distinguishing the degenerate partners of the same frequency. The coefficients [Q{_j}({\bf q})] are called normal coordinates and represent the time-dependent amplitudes of the normal modes, vibrating with frequencies [\omega{_j}({\bf q})]. The reality of the displacements and the hermiticity of the dynamical matrix impose the following conditions:[{\bf e}{_\kappa}({\bf q}, j) = {\bf e}_{\kappa}^*(-{\bf q}, j), \quad \omega{_j}({\bf q}) = \omega{_j}(-{\bf q}), \quad Q{_j}({\bf q}) = Q_{j}^*(-{\bf q}).\eqno (2.3.3.2)]In the limit of [{\bf q} \rightarrow 0], three of the vibrational branches correspond to homogeneous displacements of the crystal with vanishing frequency along three independent directions. These branches are termed acoustic. The remaining [3s-3] branches (provided that [s\,\gt \,1]) are called optic; they correspond to relative displacement of sublattices with nonzero frequency.

For completeness, we note that in a phenomenological treatment the amplitude of the polarization set up by atomic displacements of the normal modes with wavevector q can be expressed as [P_\alpha ({\bf q}) = \textstyle\sum\limits_j {Z_{j\alpha } ({\bf q})Q_j ({\bf q})}, \eqno (2.3.3.3)]introducing the ([3s\times3]) effective charge matrix of the jth normal mode of vibration (per unit cell of volume [V{_c}]), [Z_{j\alpha } ({\bf q}) = \left({{{\partial P_\alpha ({\bf q})}\over {\partial Q_j ({\bf q})}}}\right) = {1\over {V_c}}\sum\limits_\kappa {z_{\alpha \beta }^{(\kappa)}{{e_{\kappa \beta } ({\bf q},j)}\over {\sqrt {Nm_\kappa }}}}. \eqno (2.3.3.4)]The quantity [z_{\alpha \beta }^{(\kappa)}] represents the microscopic effective charge of the [{\kappa}]th atom in the unit cell. Owing to dynamic contributions of ionic deformability, it is not a scalar but rather a ([3\times 3]) tensor in general. Hence, in general, the direction of the contribution of a mode to the electrical polarization (2.3.3.3)[link] need not correspond to its mechanical polarization, as given by the mode eigenvectors [{\bf e}{_\kappa}({\bf q},j)].

In some cases, the optic modes carrying a nonzero effective charge [{\bf Z}{ _j}] (so-called polar optic modes) may be classified as transverse (TO) or longitudinal (LO). Whenever applicable, this classification should be understood in the correct sense, i.e. according to the orientation of the associated electric polarization [{\bf P} {_j}({\bf q})={\bf Z}{_j}({\bf q})Q{_j}({\bf q})] relative to the wavevector q.

2.3.3.1. First-order scattering cross section and Raman spectral line shapes

| top | pdf |

In the first-order scattering by a single excitation of the normal mode ([{\bf q},j]) of frequency [{\omega}{_j}] and wavevector [{\bf q}{_j}], the energy and wavevector conservation conditions give [\eqalign{\omega &=\pm\omega{_j} \cr{\bf q} &= \pm {\bf q}{_j},}]with the [+] sign corresponding to a Stokes process (one excitation quantum is created) and the − sign to an anti-Stokes process (one quantum is annihilated).

Let us explicitly consider the Stokes component, described by the term linear in the normal coordinate [Q{_j}] of the excitation. Inserting the plane-wave expressions for the quantities into the first term of the expansion (2.3.2.5)[link] for [\delta{\chi}] and comparing the terms with common time dependence, we get for the Stokes polarization due to the normal mode ([{\bf q},j])[{\bf P}({\bf k}{_S}) = {\varepsilon}{_0}{\chi}^{(j)}({\bf q}, {\omega}{_I},-{\omega})Q_{j}^*({\bf q},{\omega}){\bf E}{_I}({\bf k}{_I}),]which corresponds to a plane wave at a frequency of [{\omega}{_S}={\omega}{_I}-{\omega}].

Since, in the harmonic approximation, the normal modes ([{\bf q},j]) are dynamically independent, i.e. uncorrelated, the cross section for the Stokes component of the first-order scattering by phonons can be written as a sum over contributions from individual excitations. A summation convention over repeated Cartesian indices is understood throughout this chapter: [\eqalignno{&{{{\rm d}^2 \sigma }\over {{\rm d}\Omega \,\,{\rm d}\omega _S }}&\cr&\quad= {{\omega _S^3 \omega _I V^2 n_S }\over {({4\pi })^2 c^4 n_I }}\sum\limits_j {\left| {e_{I\alpha }\chi _{\alpha \beta }^{(j)}({\bf q},\omega _I, - \omega)e_{S\beta }}\right|^2 \left\langle {Q_j ({\bf q})Q_j^* ({\bf q})}\right\rangle _\omega }. &\cr&&(2.3.3.5)}]

Let us first briefly review the last term in this expression, which – together with the universal frequency-dependent first factor – essentially determines the shape of the scattered light spectrum, i.e. the frequency and temperature dependence of the spectral differential cross section. It depends exclusively on the fluctuation properties of the excitations participating in the scattering process.

The power spectrum of the fluctuations can be calculated using the linear response theory. The normal coordinates, i.e. excitation amplitudes of the normal modes, satisfy the decoupled equations of motion (with phenomenological damping constant [{\gamma}_{j,{\bf q}}] added to take into account the finite lifetime of excitations in real crystals within the quasi-harmonic approximation): [\ddot Q_j ({\bf q}) + \gamma _{j,{\bf q}}\dot Q_j ({\bf q}) + \omega _j^2 ({\bf q}) = 0.]If a fictitious generalized force [F(t)] with Fourier components [F({\omega})] is applied to the system, the average of the Fourier components of the excitation amplitude and the force are proportional, [\bar Q_j ({\bf q},\omega) = T_j ({\bf q},\omega)F(\omega),]where the proportionality coefficient [T{_j}({\bf q},{\omega})] is called a linear response function.

The fluctuation–dissipation theorem (or, more exactly, its quantum version) relates the power spectrum [\langle\ldots\rangle{_\omega}] of a fluctuating quantity to the imaginary part of the corresponding response function. The results for the Stokes and anti-Stokes components of the scattering are, respectively: [\eqalignno{\left\langle {Q_j ({\bf q})Q_j^* ({\bf q})}\right\rangle _\omega &= (\hbar/ \pi )\left [{n(\omega) + 1}\right]{\rm Im}T_j ({\bf q},\omega)\hbox{ and}&\cr \left\langle {Q_j^* ({\bf q})Q_j ({\bf q})}\right\rangle _\omega &= (\hbar/ \pi )n(\omega) {\rm Im} T_j ({\bf q},\omega),& (2.3.3.6)}]where [n({\omega})] is the Bose–Einstein statistical factor [n(\omega) = [{\exp ({{{\hbar \omega }/{k_B T}}}) - 1}]^{- 1}, \eqno (2.3.3.7)]which gives the occupation number of phonon states.

The linear response function of the normal coordinates is readily obtained from the equation of motion (N being the number of primitive cells in the crystal): [\eqalign{T_j ({\bf q},\omega) &= {1 \over N}{{1}\over {\omega _j^2 ({\bf q}) - \omega ^2 - i\omega \gamma _{j,{\bf q}}}}, \hbox{ hence}\cr {\rm Im}T_j ({{\bf q},\omega } ) &= {1 \over N} {{\omega \,\gamma _{j,{\bf q}}}\over { [{\omega _j^2 ({\bf q}) - \omega ^2 } ]^2 + \omega ^2 \gamma _{j,{\bf q}}^2 }}.}]If the damping parameter [\gamma_{j, {\bf q}}] of the excitations is much smaller than the frequency [{\omega}{_j}({\bf q})], as is usually the case for phonons in perfect crystals, this function is closely approximated by a Lorentzian centred at [{\omega}{_j}({\bf q})], with a full width [\gamma_{ j,{\bf q}}] at half maximum. Hence, introducing a suitably normalized lineshape function [L{_j}({\omega})],[\textstyle\int\limits_{- \infty }^\infty {L_j (\omega)\,\,{\rm d}\omega } = 1,]we can write for the power spectra (2.3.3.5)[link] of the jth normal mode in the spectral differential cross section [\eqalignno{\left\langle Q_j ({\bf q})Q_j^* ({\bf q})\right\rangle _\omega &= {\hbar \over {2N\omega _j ({\bf q})}}\{n[\omega _j ({\bf q})] + 1\} L_j (\omega)\hbox{ and}&\cr\left\langle Q_j^* ({\bf q})Q_j ({\bf q})\right\rangle _\omega &= {\hbar \over {2N\omega _j ({\bf q})}} n [\omega _j ({\bf q}) ] L_j (\omega ) & (2.3.3.8)}]for the Stokes and anti-Stokes case, respectively.

Note that at low temperatures the differential cross section of the anti-Stokes component in the spectrum becomes vanishingly small, because the mean number of thermally excited phonons vanishes at [T = 0].

Information about the interaction of photons with individual excitations is contained in the central term in the expression for the cross section (2.3.3.5)[link], i.e. in the nonlinear susceptibility [\chi _{\alpha \beta }^{(j)}({\bf q},\omega _I, - \omega)].

2.3.3.2. Symmetry properties of the scattering cross section

| top | pdf |

The quantity that controls the symmetry properties of the scattering cross section due to excitation [Q{_j}({\bf q})] is the squared modulus of the corresponding second-order susceptibility (second-rank tensor), contracted with the two polarization vectors of the incident and scattered light: [\left| {{\bf e}_S \boldchi ^{(j)}({{\bf q},\omega _I, - \omega }){\bf e}_I }\right|^2 \equiv \left| {e_{S\alpha }\chi _{\alpha \beta }^{(j)}({\bf q},\omega _I, - \omega)e_{I\beta }}\right|^2. \eqno (2.3.3.9)]The nonlinear susceptibility tensor [{\boldchi}^{(j)}({\bf q},{\omega}{_I},-{\omega})] is usually referred to as the first-order Raman tensor (defined in the literature to within a factor). Before discussing the consequences of the crystal symmetry on the form of the Raman tensor, let us mention two important approximations on which conventional analysis of its symmetry properties is based.

In a general case, the second-order susceptibilities are not necessarily symmetric. However, they fulfil a general symmetry property which follows from the symmetry of the scattering with respect to time inversion. Since the anti-Stokes process can be regarded as a time-inverted Stokes process (exchanging the role of the incident and scattered photons), it can be shown that in non-magnetic materials the susceptibilities obey the relation [\chi _{\beta \alpha }^{(j)}({- {\bf q},\omega _S, \omega }) = \chi _{\alpha \beta }^{(j)}({{\bf q},\omega _I, - \omega }). \eqno (2.3.3.10)]

In the quasi-static limit, i.e. if the scattering frequency is negligibly small compared with the incident photon frequency ([0\approx {\omega}\ll{\omega}{_I}\approx{\omega}{_S}]), it follows that the susceptibilities of non-magnetic materials become symmetric in the Cartesian indices α, β. This symmetry is very well fulfilled in a great majority of cases. Appreciable antisymmetric contributions are known to occur, e.g. under resonant conditions, where the quasi-static approximation breaks down as the energy of the incident (or scattered) photon approaches those of electronic transitions.

Thus, in the first approximation, we set [{\omega}] equal to zero and remove the time dependence in the phonon amplitudes, treating the normal coordinates as static. Then the nonlinear susceptibilities correspond to susceptibility derivatives, [\chi _{\alpha \beta }^{(j)}({\bf q},\omega _I, - \omega) \to R_{\alpha \beta }^j ({\bf q}) \equiv \chi _{\alpha \beta }^{(j)}\left({{\bf q},\omega _I }\right) = {{\partial \chi _{\alpha \beta }(\omega _I)}\over {\partial Q_j^* ({\bf q})}}, \eqno (2.3.3.11)]where we suppressed the explicit dependence on [{\omega}{_I}] and introduced a simplified notation for the Raman tensor [{\bf R}{^j}({\bf q})], still keeping the dependence on the scattering wavevector.

In deriving the symmetry properties of the Raman tensor [{\bf R}{^j}({\bf q})] that follow from the crystal lattice symmetry, the main point is thus to determine its transformation properties under the symmetry operation of the crystal space group.

Since the magnitude of the scattering vector [q\equiv |{\bf q}|] is very small compared with the Brillouin-zone dimensions, another conventional approximation is to neglect the q dependence of the susceptibilities. Setting  [{\bf q}\rightarrow 0] enables us to analyse the symmetry of the Raman tensor in terms of the factor group [G{_0}], which is isomorphous to the point group of the crystal lattice. This approach is, again, appropriate for the vast majority of cases. An important exception is, for instance, the scattering by acoustic modes (Brillouin scattering) or scattering by longitudinal plasma waves in semiconductors (plasmons): in these cases the Raman tensor vanishes for [{\bf q} = 0], since this limit corresponds to a homogeneous displacement of the system. Possible q-dependent effects can be treated by expanding the Raman tensor in powers of q and using compatibility relations between the symmetries at [{\bf q} = 0] and at the full symmetries applicable in the [{\bf q} \neq 0] case.

Let us mention that another notation is sometimes used in the literature for the Raman tensor. Since the square modulus of a second-rank tensor contracted with two vectors can be written as a fourth-rank tensor contracted with four vectors, one can introduce a fourth-rank tensor I(j), [I_{\alpha \beta \mu \nu }^{(j)} = {{\partial \chi _{\alpha \beta }^* }\over {\partial Q_j^{}}}{{\partial \chi _{\mu \nu }^{}}\over {\partial Q_j^ * }},]so that the scattering cross section of the jth mode is [{{{\rm d}\sigma ^{(j)}}\over {{\rm d}\Omega \,\,{\rm d}\omega _S }} \approx e_{I\alpha }e_{I\mu }e_{S\beta }e_{S\nu }I_{\alpha \beta \mu \nu }^{(j)}.]If there are no antisymmetric components in the susceptibility derivatives, it can be shown that the fourth-rank tensor [{\bf I}=(I_{\alpha\beta\mu\nu})] has at most 21 independent components, as for the elastic constants tensor.

2.3.3.3. Raman tensor and selection rules at [{\bf q}\approx 0]

| top | pdf |

The scattering cross section, being a scalar quantity, must be invariant with respect to all symmetry elements of the space group of the crystal. This invariance has two important consequences: it determines which normal modes (j) can contribute to the scattering (Raman activity of the modes) and it also gives the restrictions on the number of independent components of the Raman tensor (polarization selection rules).

At [{\bf q}\approx 0], the transformation properties of the incident and scattered light are described by the three-dimensional polar vector representation [\Gamma_{\rm PV}] of the appropriate point group of the crystal, since the quantities that characterize the light ([{\bf E}{_I}], [{\bf E}{_S}], [{\bf P}\ldots]) are all polar vectors, i.e. first-rank polar tensors ([T{_\alpha}]). The transformation properties of a normal mode j must correspond to an irreducible matrix representation [\Gamma(j)] of the crystal point group. We recall that in cases of two- or three-dimensional representations (degeneracy), the index j represents two indices.

In order that a particular normal mode j in a given crystal be Raman active, i.e. symmetry-allowed to contribute to the (first-order) scattering cross section, the necessary condition is that the corresponding irreducible representation [\Gamma(j)] must be contained in the decomposition of the direct Kronecker product representation [\Gamma_{\rm PV}\otimes\Gamma_{\rm PV}] at least once: [{\Gamma}_{\rm PV}\otimes{\Gamma}_{\rm PV}\supset{\Gamma}(j). \eqno (2.3.3.12)]In this case, the Kronecker product [\Gamma_{\rm PV}\otimes{\Gamma}_{\rm PV}\otimes{\Gamma}(j)] contains the identity representation at least once, so the cross section remains invariant under the transformation of the crystal point group. In the phenomenological formulation, the susceptibility derivatives correspond to third derivatives of a particular potential energy [\Phi] (interaction Hamiltonian), [R_{\alpha \beta }^j = \left({{{\partial \chi _{\alpha \beta } }\over {\partial Q_j^ * }}}\right)\sim \left({{{\partial ^3 \Phi }\over {\partial E_{I\alpha }\partial E_{S\beta }^ * \partial Q_j^ * }}}\right), \eqno (2.3.3.13)]such that the product [\Gamma_{\rm PV}\otimes{\Gamma}_{\rm PV}\otimes{\Gamma}(j)] is the reducible representation of the Raman tensor. If condition (2.3.3.12)[link] holds, then the Raman tensor [{\bf R}{^j}] does not vanish identically and may have at least one independent nonzero component. As the representation [\Gamma_{\rm PV}\otimes\Gamma_{\rm PV}] is that of a second-rank polar tensor, equivalent formulation of the Raman activity of a normal mode j is that the corresponding normal coordinate [Q{_j}] must transform like one or more components of a polar tensor. The transition susceptibility [{\delta \chi}^{(j)}] transforms accordingly. The task of determining whether a given normal mode j is Raman active or not thus consists of simply decomposing the representation [\Gamma_{\rm PV}\otimes\Gamma_{\rm PV}] and identifying the irreducible components [\Gamma (j)].

The second consequence of the invariance condition is the imposition of restrictions on the Cartesian components of the Raman tensor for modes allowed to participate in the scattering. By virtue of the properties of the irreducible representations [\Gamma(j)], some components of the corresponding Raman tensor are required to vanish whereas others may have related values. This fact results in anisotropies in the observed cross section depending on the polarization directions of the incident and scattered light, and is usually referred to as polarization selection rules. As the scattering cross section of the excitation (j) is proportional to the scalar quantity [\left| {{\bf e}_S {\bf R}^j {\bf e}_I }\right|^2 \equiv \left| {e_{S\alpha }R_{\alpha \beta }^j e_{I\beta }}\right|^2,]one can generally `isolate' a given component of the Raman tensor by suitably arranging the scattering geometry in the experiment, i.e. by choosing the orientation of the wavevectors [{\bf k}{_I}] and [{\bf k}{_S}] and the polarization vectors [{\bf e}{_I}] and [{\bf e}{_S}] with respect to crystallographic axes.

For each normal mode (j) allowed in the scattering, the number of independent components of its Raman tensor is given by the multiplicity coefficients [c^{(j)}] of the irreducible representation [\Gamma(j)] in the decomposition [\Gamma_{\rm PV}\otimes\Gamma_{\rm PV} = c^{(1)}\Gamma(1)\oplus c^{(2)}\Gamma(2)\oplus \ldots, \eqno (2.3.3.14)]where the multiplicity coefficient [c^{(j)}] corresponds to the number of times the given irreducible representation [\Gamma(j)] enters the decomposition. If the representation [\Gamma(j)] is two- or three-dimensional, then for each occurrence of [\Gamma (j)] in (2.3.3.14)[link] there are two or three degenerate partners (of the same frequency) whose Raman tensors are symmetry-related.

The matrix form of the Raman tensor corresponding to a given irreducible representation – i.e. symmetry species[\Gamma(j)] can be readily constructed by finding the appropriate bilinear basis functions that transform according to the corresponding irreducible representation [\Gamma(j)]. The required number of such independent bases is given by the multiplicity coefficient [c^{(j)}]. Alternatively, one may construct invariant polynomials (transforming as scalars) of order four, i.e. of the same order as the product [E_{I\alpha}E_{S\beta}Q{_j}].

Making allowance for possible antisymmetric scattering, we have not explicitly supposed that the Raman tensor is symmetric. We recall that the derivative (2.3.3.13)[link] is not necessarily symmetric in the α and β indices as long as the fields [E_{I\alpha}] and [E_{S\beta}] correspond to different frequencies (inelastic scattering). However, each second-rank polar tensor [T_{\alpha\beta}] (nine components), transforming according to [\Gamma_{\rm PV}\otimes \Gamma_{\rm PV}], can be decomposed into a symmetric part [T'_{\alpha\beta}=T'_{\beta\alpha}] (six components), transforming like a symmetric polar tensor [[\Gamma_{\rm PV}\otimes \Gamma_{\rm PV}]{_S}], and an antisymmetric part [T''_{\alpha\beta}=-T''_{\beta\alpha}] (three components), transforming like an axial vector (for the definition of axial tensors, see Section 1.1.4.5.3[link] ) according to [\Gamma_{\rm AV}=[\Gamma_{\rm PV} \otimes \Gamma_{\rm PV}]{_A}].

The symmetry-restricted forms of the ([3\times 3]) Raman tensors corresponding to all Raman-active symmetry species are summarized in Table 2.3.3.1[link] (see e.g. Hayes & Loudon, 1978[link]) for each of the 32 crystal symmetry classes. Spectroscopic notation is used for the irreducible representations of the point groups. The symbols (x, y or z) for some Raman-active symmetry species in the noncentrosymmetric classes indicate that the respective components of polar vectors also transform according to these irreducible representations. Hence the normal coordinates of the phonons of these polar symmetry species (polar phonons) transform in the same way and, consequently, the corresponding component of the effective charge tensor [Z_{j\alpha}({\bf q}=0)], see (2.3.3.4)[link], is not required by symmetry to vanish. Polar phonons thus may carry a nonzero dipole moment and contribute to the polarization in the crystal, which manifests itself in infrared activity and also in the Raman scattering cross section (see Section 2.3.3.5[link]).

Table 2.3.3.1| top | pdf |
Symmetry of Raman tensors in the 32 crystal classes

The symbols a, b, c, d, e, f, g, h and i in the matrices stand for arbitrary parameters denoting possible independent nonzero components (in general complex) of the Raman tensors. Upper row: conventional symmetric Raman tensors; lower row: antisymmetric part. Alternative orientations of the point group are distinguished by subscripts at 2 or m in the class symbol indicating the direction of the twofold axis or of the normal to the mirror plane.

Triclinic

  [\pmatrix{ a & d & f \cr d & b & h \cr f & h & c \cr }]
1 [{\rm A}(x,y,z)]
[\bar 1] [{\rm A}{_g}]
  [\pmatrix{. & e & g \cr {- e}&. & i \cr {- g}& {- i}&. \cr }]

Monoclinic, unique axis z

  [{\pmatrix{ a & d &. \cr d & b &. \cr. &. & c \cr }}] [{\pmatrix{. &. & f \cr. &. & h \cr f & h &. \cr }}]
[2_z] [{\rm A}(z)] [{\rm B}(x,y)]
[m_z] [{\rm A'}(x,y)] [{\rm A}''(z)]
[2_z/m] [{\rm A}{_g}] [{\rm B}{_g}]
  [{\pmatrix{. & e &. \cr {- e}&. &. \cr. &. &. \cr }}] [{\pmatrix{. &. & g \cr. &. & i \cr {- g}& {- i}&. \cr }}]

Monoclinic, unique axis y

  [{\pmatrix{ a & . &d \cr . & b &. \cr d&. & c \cr }}] [{\pmatrix{. &f & . \cr f &. & h \cr . & h &. \cr }}]
[2_y] [{\rm A}(y)] [{\rm B}(x,z)]
[m_y] [{\rm A'}(x,z)] [{\rm A}''(y)]
[2_y/m] [{\rm A}{_g}] [{\rm B}{_g}]
  [{\pmatrix{. & . &e \cr .&. &. \cr{-e} &. &. \cr }}] [{\pmatrix{. &g & . \cr-g &. & i \cr .& {- i}&. \cr }}]

Orthorhombic

  [{\pmatrix{ a &. &. \cr. & b &. \cr. &. & c \cr }}] [{\pmatrix{. & d &. \cr d &. &. \cr. &. &. \cr }}] [{\pmatrix{. &. & f \cr. &. &. \cr f &. &. \cr }}] [{\pmatrix{. &. &. \cr. &. & h \cr. & h &. \cr }}]
222 [{\rm A}] [{\rm B}{_1}(z)] [{\rm B}{_2}(y)] [{\rm B}{_3}(x)]
[mm2] [{\rm A}{_1}(z)] [{\rm A}{_2}] [{\rm B}{_1}(x)] [{\rm B}{_2}(y)]
[mmm] [{\rm A}{_g}] [{\rm B}_{1g}] [{\rm B}_{2g}] [{\rm B}_{3g}]
    [{\pmatrix{. & e &. \cr {- e}&. &. \cr. &. &. \cr }}] [{\pmatrix{. &. & g \cr. &. &. \cr {- g}&. &. \cr }}] [{\pmatrix{. &. &. \cr. &. & i \cr. & {- i}&. \cr }}]

Tetragonal

  [{\pmatrix{ a &. &. \cr. & a &. \cr. &. & b \cr }}] [{\pmatrix{ d & e &. \cr e & {- d}&. \cr. &. &. \cr }}] [{\pmatrix{. &. & f \cr. &. & h \cr f & h &. \cr }}\quad{\pmatrix{. &. & {- h} \cr. &. & f \cr {- h}& f &. \cr }}]
[4] [{\rm A}(z)] [{\rm B}] [{\rm E}(x,y)]
[\bar 4] [{\rm A}] [{\rm B}(z)] [{\rm E}(x,-y)]
[4/m] [{\rm A}{_g}] [{\rm B}{_g}] [{\rm E}{_g}]
  [{\pmatrix{. & c &. \cr {- c}&. &. \cr. &. &. \cr }}]   [{\pmatrix{. &. & g \cr. &. & i \cr {- g}& {- i}&. \cr }}\quad{\pmatrix{. &. & {- i} \cr. &. & g \cr i & {- g}&. \cr }}]

  [{\pmatrix{ a &. &. \cr. & a &. \cr. &. & b \cr }}]   [{\pmatrix{ d &. &. \cr. & {- d}&. \cr. &. &. \cr }}] [{\pmatrix{. & e &. \cr e &. &. \cr. &. &. \cr }}] [{\pmatrix{. &. & f \cr. &. &. \cr f &. &. \cr }}\quad{\pmatrix{. &. &. \cr. &. & f \cr. & f &. \cr }}]
422 [{\rm A}{_1}] [{\rm A}{_2}(z)] [{\rm B}{_1}] [{\rm B}{_2}] [{\rm E}(-y,x)]
[4mm] [{\rm A}{_1}(z)] [{\rm A}{_2} ] [{\rm B}{_1}] [{\rm B}{_2}] [{\rm E}(x,y)]
[\bar 42_xm_{xy}] [{\rm A}{_1}] [{\rm A}{_2}] [{\rm B}{_1}] [{\rm B}{_2}(z) ] [{\rm E}(y,x)]
[\bar 4m_x2_{xy}] [{\rm A}{_1}] [{\rm A}{_2}] [{\rm B}{_1}] [{\rm B}{_2}(z) ] [{\rm E}(-x,y)]
[4/mmm] [{\rm A}_{1g}] [{\rm A}_{2g}] [{\rm B}_{1g}] [{\rm B}_{2g}] [{\rm E}{_g}]
    [{\pmatrix{. & c &. \cr {- c}&. &. \cr. &. &. \cr }}]     [{\pmatrix{. &. & g \cr. &. &. \cr {- g}&. &. \cr }}\quad{\pmatrix{. &. &. \cr. &. & g \cr. & {- g}&. \cr }}]

Trigonal

  [{\pmatrix{ a &. &. \cr. & a &. \cr. &. & b \cr }}] [{\pmatrix{ c & f & e \cr f & {- c}& d \cr e & d & . \cr }}\quad{\pmatrix{ f & {- c}& {- d} \cr {- c}& {- f}& e \cr {- d}& e &. \cr }}]
3 [{\rm A}(z)] [{\rm E}(x,y)]
[\bar 3] [{\rm A}{_g}] [{\rm E}{_g}]
  [{\pmatrix{. & h &. \cr {- h}&. &. \cr. &. &. \cr }}] [{\pmatrix{. &. & i \cr. &. & g \cr {- i}& {- g}&. \cr }}\quad{\pmatrix{. &. & {- g} \cr. &. & i \cr g & {- i}&. \cr }}]

  [{\pmatrix{ a &. &. \cr. & a &. \cr. &. & b \cr }}]   [{\pmatrix{f & . &. \cr . &-f & d \cr. & d &. \cr }}\quad{\pmatrix{ . &-f & {- d} \cr-f & .&. \cr {- d}&. &. \cr }}]
[32_x] [{\rm A}{_1}] [{\rm A}{_2}(z)] [{\rm E}(x,y)]
[3m_x] [{\rm A}{_1}(z)] [{\rm A}{_2}] [{\rm E}(y,-x)]
[\bar 3m_x] [{\rm A}_{1g}] [{\rm A}_{2g}] [{\rm E}_{g}]
    [{\pmatrix{. & h &. \cr {- h}&. &. \cr. &. &. \cr }}] [{\pmatrix{. &. &. \cr. &. & g \cr. & {- g}&. \cr }}\quad{\pmatrix{. &. & {- g} \cr. &. &. \cr g &. &. \cr }}]

  [{\pmatrix{ a &. &. \cr. & a &. \cr. &. & b \cr }}]   [{\pmatrix{. & f &. \cr f &. & d \cr. & d &. \cr }}\quad{\pmatrix{ f &.& {- d} \cr. & -f&. \cr {- d}&. &. \cr }}]
[32_y] [{\rm A}{_1}] [{\rm A}{_2}(z)] [{\rm E}(x,y)]
[3m_y] [{\rm A}{_1}(z)] [{\rm A}{_2}] [{\rm E}(y,-x)]
[\bar 3m_y] [{\rm A}_{1g}] [{\rm A}_{2g}] [{\rm E}_{g}]
    [{\pmatrix{. & h &. \cr {- h}&. &. \cr. &. &. \cr }}] [{\pmatrix{. &. &. \cr. &. & g \cr. & {- g}&. \cr }}\quad{\pmatrix{. &. & {- g} \cr. &. &. \cr g &. &. \cr }}]

Hexagonal

  [{\pmatrix{ a &. &. \cr. & a &. \cr. &. & b \cr }}] [{\pmatrix{. &. & e \cr. &. & d\cr e & d &. \cr }}\quad{\pmatrix{. &. & {- d} \cr. &. & e \cr {- d}& e &. \cr }}] [{\pmatrix{ c & f &. \cr f & {- c}&. \cr. &. &. \cr }}\quad{\pmatrix{ f & {- c}&. \cr {- c}& {- f}&. \cr. &. &. \cr }}]
6 [{\rm A}(z)] [{\rm E}{_1}(x,y)] [{\rm E}{_2}]
[\bar 6] [{\rm A}' ] [{\rm E}''] [{\rm E}'(x,y)]
[6/m] [{\rm A}{_g}] [{\rm E}_{1g}] [{\rm E}_{2g}]
  [{\pmatrix{. & h &. \cr {- h}&. &. \cr. &. &. \cr }}] [{\pmatrix{. &. & i \cr. &. & g \cr {- i}& {- g}&. \cr }}\quad{\pmatrix{. &. & {- g} \cr. &. & i \cr g & {- i}&. \cr }}]  

  [{\pmatrix{ a &. &. \cr. & a &. \cr. &. & b \cr }}]   [{\pmatrix{. &. &. \cr. &. & d \cr. & d &. \cr }}\quad{\pmatrix{. &. & {- d} \cr. &. &. \cr {- d}&. &. \cr }}] [{\pmatrix{. & f &. \cr f &. &. \cr. &. &. \cr }}\quad{\pmatrix{ f &. &. \cr. & {- f}&. \cr. &. &. \cr }}]
622 [{\rm A}{_1}] [{\rm A}{_2}(z) ] [{\rm E}{_1}(x,y)] [{\rm E}{_2}]
[6mm] [{\rm A}{_1}(z)] [{\rm A}{_2}] [{\rm E}{_1}(y,-x)] [{\rm E}{_2}]
[\bar 6m_x2_y] [{\rm A}{_1}'] [{\rm A}{_2}' ] [{\rm E}''] [{\rm E}'(x,y)]
[\bar 62_xm_y] [{\rm A}{_1}'] [{\rm A}{_2}' ] [{\rm E}''] [{\rm E}'(y,-x)]
[6/mmm] [{\rm A}_{1g}] [{\rm A}_{2g} ] [{\rm E}_{1g}] [{\rm E}_{2g}]
    [{\pmatrix{. & h &. \cr {- h}&. &. \cr. &. &. \cr }}] [{\pmatrix{. &. &. \cr. &. & g \cr. & {- g}&. \cr }}\quad{\pmatrix{. &. & {- g} \cr. &. &. \cr g &. &. \cr }}]  

Cubic

  [{\pmatrix{ a &. &. \cr. & a &. \cr. &. & a \cr }}] [{\pmatrix{ b &. &. \cr. & b &. \cr. &. & {- 2b} \cr }}\quad{\pmatrix{ {- \sqrt{3}b}&. &. \cr. & \sqrt{3}b&. \cr. &. &. \cr }}] [{\pmatrix{. &. &. \cr. &. & c \cr. & c &. \cr }}\quad{\pmatrix{. &. & c \cr. &. &. \cr c &. &. \cr }}\quad{\pmatrix{. & c &. \cr c &. &. \cr. &. &. \cr }}]
23 [{\rm A}] [{\rm E}] [{\rm F}(x,y,z)]
[m3] [{\rm A}{_g}] [{\rm E}{_g}] [{\rm F}{_g}]
      [{\pmatrix{. &. &. \cr. &. & d \cr. & {- d}&. \cr }}\quad{\pmatrix{. &. & d \cr. &. &. \cr {- d}&. &. \cr }}\quad{\pmatrix{. & d &. \cr {- d}&. &. \cr. &. &. \cr }}]

  [{\pmatrix{ a &. &. \cr. & a &. \cr. &. & a \cr }}] [{\pmatrix{ b &. &. \cr. & b &. \cr. &. & {- 2b} \cr }}\quad{\pmatrix{ {- \sqrt{3}b}&. &. \cr. & {\sqrt{3}}b &. \cr. &. &. \cr }}]   [{\pmatrix{. &. &. \cr. &. & c \cr. & c &. \cr }}\quad{\pmatrix{. &. & c \cr. &. &. \cr c &. &. \cr }}\quad{\pmatrix{. & c &. \cr c &. &. \cr. &. &. \cr }}]
432 [{\rm A}{_1}] [{\rm E}] [{\rm F}{_1}(x,y,z)] [{\rm F}{_2}]
[\bar 43m] [{\rm A}{_1}] [{\rm E}] [{\rm F}{_1}] [{\rm F}{_2}(x,y,z)]
[m3m] [{\rm A}_{1g}] [{\rm E}_{g}] [{\rm F}_{1g}] [{\rm F}_{2g}]
      [{\pmatrix{. &. &. \cr. &. & d \cr. & {- d}&. \cr }}\quad{\pmatrix{. &. & d \cr. &. &. \cr {- d}&. &. \cr }}\quad{\pmatrix{. & d &. \cr {- d}&. &. \cr. &. &. \cr }}]  

For convenience, the Raman tensors are explicitly split into a symmetric and possible antisymmetric part (upper and lower row of each part of the table, respectively, in each case). The conventional symmetric Raman tensors are appropriate for most cases of practical interest. Besides the resonant conditions mentioned above, there are other exceptions. For instance, there are optical phonons that transform like axial vectors, such as in the case of [{\rm A}{_2}] (or [{\rm A}_{2g}], [{\rm A}'_2]) modes in some uniaxial crystal classes, where the Raman tensor is purely antisymmetric. Antisymmetric scattering by these modes may become allowed at finite wavevector q. Antisymmetric Raman tensors are also needed for analysing the symmetry of scattering in magnetic materials (scattering by spin waves – magnons), or non-magnetic materials under a magnetic field, where the susceptibility itself is essentially nonsymmetric.

We note that the matrix form of the Raman tensors depends on the setting of the Cartesian axes with respect to the crystallographic axes. To avoid ambiguities and apparent disagreement with other sources, we give the results for alternative orientations of the point groups in several cases where different settings of the twofold axes or mirror planes with respect to the Cartesian axes are commonly used. This concerns all monoclinic classes (unique direction parallel to y or z), tetragonal class [\bar 4 2 m], trigonal classes 32, 3m and [\bar 3 m], as well as hexagonal class [\bar 6 2 m].

2.3.3.4. Centrosymmetric crystals

| top | pdf |

In those point groups that contain the inversion operation, i.e. in the eleven centrosymmetric (nonpolar) crystal classes [{\bar 1}, 2/m, mmm, 4/m, 4/mmm, {\bar 3}, {\bar 3m}, 6/m, 6/mmm, m3, m3m,]the irreducible representations are divided into two groups, odd and even, according to the parity. Since second-rank polar tensors must transform according to the even parity representations only, whereas polar vectors transform according to odd parity representations, the selection rules for electric dipole absorption (infrared activity) and for Raman scattering are incompatible. This is often expressed as the mutual exclusion rule or complementarity principle: The excitations in a crystal belonging to a centrosymmetric class cannot be simultaneously active in infrared absorption and in Raman scattering. Let us note, however, that even-parity excitations are not necessarily all Raman active, and that odd-parity excitations are not necessarily infrared active.

In the remaining noncentrosymmetric crystal classes, the excitations have no defined parity with respect to inversion and can be, in principle, both Raman and infrared active.

Example:  Consider a Raman scattering experiment on a crystal of tetragonal symmetry, class [4/mmm]. Raman-active phonons, allowed in conventional symmetric scattering, are of the symmetry species [{\rm A}_{1g}], [{\rm B}_{1g}], [{\rm B}_{2g}] and [{\rm E}{_g}]. (the [{\rm A}_{2g}] species admits purely antisymmetric scattering only). Straightforward application of Table 2.3.3.1[link] makes it possible to determine the polarization selection rules, i.e. to determine which symmetry species will contribute to the scattering cross section in various experimental configurations. Choosing the Cartesian axes [x \equiv [100]], [y \equiv [010]], [z \equiv [001]] consistent with the standard setting of the [4/mmm] point group, i.e. the fourfold rotation axis [4\parallel{\bf z}], let us further introduce the notation [x' \equiv [110]], [y' \equiv [\bar 110]]. Then the contributions to the cross section for different symmetry species can be distinguished by their dependence on the polarization vectors [{\bf e}{_I}] and [{\bf e}{_S}] of the incident and scattered light:[\eqalign{{\rm A}_{1g}&: |(e{_I}{_x}e{_S}{_x} + e{_I}{_y}e{_S}{_y})a + e{_I}{_z}e{_S}{_z}b|{^2}\cr {\rm B}_{1g}&:|(e{_I}{_x}e{_S}{_x}-e{_I}{_y}e{_S}{_y})d|{^2} \cr {\rm B}_{2g}&: |(e{_I}{_x}e{_S}{_y} + e{_I}{_y}e{_S}{_x}) e|{^2}\cr {\rm E}_{g}&: [(e{_I}{_x}e{_S}{_z} + e{_I}{_z}e{_S}{_x}){^2} + (e{_I}{_y}e{_S}{_z} + e{_I}{_z}e{_S}{_y}){^2}] |f|{^2}.\cr}]

Examples of some special scattering geometries that permit the separation of the contributions of different symmetry species are shown in Table 2.3.3.2[link] (five distinct configurations are sufficient to determine the five independent parameters a, b, d, e, f of the symmetric Raman tensors).

Table 2.3.3.2| top | pdf |
Raman selection rules in crystals of the [4/mmm] class

Scattering configurationCross section for symmetry species
Back scatteringRight-angle scattering[{\rm A}_{1g}][{\rm B}_{1g}][{\rm B}_{2g}][{\rm E}{_g}]
[\bar z(xx)z,\bar z(yy)z] [y(xx)z, x(yy)z] [\sim|a|^2] [\sim|d|^2]
[\bar x(zz)x, \bar y(zz)y] [x(zz)y] [\sim|b|^2]
[\bar z(xy)z] [y(xy)x, z(xy)x] [\sim|e|^2]
[\bar y(xz)y, \bar x(yz)x] [y(xz)x, x(yz)y] [\sim|f|^2]
[\bar z(x'x')z] [y'(x'x')z] [\sim|a|^2] [\sim|e|^2]
[\bar z(x'y')z] [y'(x'y')z] [\sim|d|^2]

If, for some reason, antisymmetric scattering is allowed, possible contribution of the [{\rm A}_{2g}] modes should be considered as well. The contribution to cross section from these modes is proportional to [|(e{_I}{_x}e_{Sy}-e{_I}{_y}e{_S}{_x}) c|{^2}], hence it can be distinguished from the contribution of the [{\rm B}_{2g}] symmetry species by a suitable choice of the scattering geometry.

2.3.3.5. Noncentrosymmetric crystals

| top | pdf |

Special care is required in treating the scattering by those optical phonons in the 21 noncentrosymmetric (polar) crystal classes (1, 2, m, 222, mm2, 4, [{\bar 4}], 422, 4mm, [{\bar 42m}], 3, 32, 3m, 6, [{\bar 6}], 622, 6mm, [{\bar 6m2}], 23, 432, [\bar 43m]) that are simultaneously infrared-active. Since these polar modes carry a nonzero macroscopic effective charge (2.3.3.4)[link], they contribute to the total polarization in the crystal, hence also to the macroscopic electric field, which in turn leads to a coupling between these modes. The polarization being a polar vector, the modes that contribute have the same symmetry character, i.e. they must also transform like the components of polar vectors.

An important consequence of the macroscopic field associated with polar modes in the crystal is the partial lifting of the degeneracies of the long-wavelength ([{\bf q}\approx 0]) mode frequencies (so-called TO–LO splitting). Since the macroscopic field in the crystal is longitudinal, it must be proportional to the longitudinal component of the polarization. Hence, the equations of motion for all polar modes carrying a nonzero longitudinal polarization (i.e. [{\bf P}_j\cdot{\bf q}\neq 0]) become coupled by the field and, consequently, their frequencies depend on the direction of q. This phenomenon is called directional dispersion and is connected with the fact that in the electrostatic approximation the dynamical matrix with long-range Coulomb forces shows non-analytic behaviour for [{\bf q}\to 0]. In lattice dynamics, the limit can be treated correctly by taking into account the retardation effects in the range where cq becomes comparable to [\omega{_j}({\bf q})], i.e. in the crossing region of free photon and optical phonon dispersion curves. As a result, one finds that for small q the true eigenmodes of the system – polaritons – have a mixed phonon–photon character and their frequencies show strong dispersion in the very close vicinity of [{\bf q}=0]. Experimentally, this polariton region is partially accessible only in near-forward Raman scattering [see (2.3.2.3)[link]]. For larger scattering wavevectors in the usual right-angle or back-scattering geometries, the electrostatic approximation, [cq\gg\omega], is well applicable and the excitations behave like phonons. Owing to the coupling via the longitudinal macroscopic electric field, however, the directional dispersion of these phonon branches remains.

Detailed analysis is complicated in the general case of a low-symmetry crystal with more polar modes (see e.g. Claus et al., 1975[link]). In crystals with at least orthorhombic symmetry, the principal axes of the susceptibility tensor are fixed by symmetry and for the wavevectors oriented along these principal axes the polar optic modes have purely transverse (TO) or longitudinal (LO) character with respect to the associated polarization. The character of a mode is usually mixed for a general direction of the wavevector.

Strictly speaking, conventional symmetry analysis in terms of irreducible representations of the factor group (point group) of the crystal, though giving a true description of polaritons at [{\bf q }= 0], cannot account for the lifting of degeneracies and for the directional dispersion of polar modes. A correct picture of the symmetries and degeneracies is, however, obtained by taking into account the finiteness of the wavevector q and classifying the vibrations according to the irreducible (multiplier) corepresentations of the point group of the wavevector [G{_0}({\bf q})], which is a subgroup of the factor group. Compatibility relations of the representations at [{\bf q}\to 0] can then be used to establish a correspondence between the two approaches.

The oscillating macroscopic field associated with long-wavelength LO polar modes acts as another source of modulation of the susceptibility. In addition to the standard atomic displacement contribution connected with the mechanical displacements of atoms, one also has to consider that the transition susceptibility also contains the electro-optic term arising from the distortion of electron shells of atoms in the accompanying macroscopic field E. Separating both contributions, we may write [\eqalignno{\delta \chi _{\alpha \beta }^{(j)}({\bf q}\approx 0,\omega _I) &= {{{\rm d}\chi _{\alpha \beta } }\over {{\rm d}Q_j }}Q_j &\cr&= {{\partial \chi _{\alpha \beta }}\over {\partial Q_j }}Q_j + {{\partial \chi _{\alpha \beta }}\over {\partial E_\gamma }}E_\gamma ^j &\cr&= \left({{{\partial \chi _{\alpha \beta } }\over {\partial Q_j }}+ {{\partial \chi _{\alpha \beta } }\over {\partial E_\gamma }}{{{\rm d}E_\gamma }\over {{\rm d}Q_j }}}\right)Q_j, &\cr&&(2.3.3.15)}]or, in terms of the Raman tensor, [\delta \chi = \textstyle\sum\limits_j {{\bf R}^j Q_j} = \textstyle\sum\limits_j {({\bf a}^j Q_j + {\bf bE}^j)} =\textstyle\sum\limits_j \left[{\bf a}^j + {\bf b}({\rm d}{\bf E}/{{\rm d}Q_j })\right]Q_j ,]where we introduce the notation [{\bf a}{^j}] and [{\bf b}({\rm d}{\bf E}/{\rm d}Q{_j})] for the atomic displacement and electro-optic contributions to the Raman tensor [{\bf R}{^j}]. As usual, [Q{_j}] stands for the normal coordinate of the jth mode and E for the total macroscopic electric field resulting from the longitudinal polarization of all optic modes. The modes that contribute to E are only LO polar modes; they transform as Cartesian components of polar vectors (x, y, z). Hence the electro-optic term contributes to the Raman cross section only if [{\bf E}{^j} = ({\rm d}{\bf E}/{\rm d}Q{_j})\neq 0], i.e. if the mode has at least partially longitudinal character. Hence, not only the frequencies but also the scattering cross sections of the TO and LO components of polar modes belonging to the same symmetry species are, in general, different.

Nevertheless, in view of the fact that the macroscopic electric field associated with LO polar phonons transforms in the same way as its polarization vector, the symmetry properties of both the atomic displacement and the electro-optic contributions to the Raman tensors of polar modes are identical. They correspond to third-rank polar tensors, which have nonzero components only in piezoelectric crystals. The symmetry-restricted form of these tensors can also be derived from Table 2.3.3.1[link] by combining the matrices corresponding to the x, y and z components. Note that these may belong to different irreducible representations in lower-symmetry classes (e.g. z cannot mix with x, y), and that in some uniaxial classes the z component is missing completely. Finally, in the noncentrosymmetric class 32 of the cubic system, the Raman tensors of the triply degenerate polar modes ([{\rm F}{_1}]) are purely antisymmetric; therefore all components of the piezoelectric tensor also vanish.

Example: To illustrate the salient features of polar-mode scattering let us consider a crystal of the [4mm] class, where of the Raman-active symmetry species the modes [{\rm A}{_1}(z)] and [{\rm E}(x,y)] are polar. According to Table 2.3.3.1[link], their ([{\bf q}=0]) Raman tensors are identical to those of the [{\rm A}_{1g}] and [{\rm E}{_g}] modes in the preceding example of a [4/mmm]-class crystal. Owing to the macroscopic electric field, however, here one has to expect directional dispersion of the frequencies of the long wavelength ([{\bf q}\approx 0]) [{\rm A}{_1}] and E optic phonon modes according to their longitudinal or transverse character. Consequently, in determining the polarization selection rules, account has to be taken of the direction of the phonon wavevector (i.e. the scattering wavevector) q with respect to the crystallographic axes. Since for a general direction of q the modes are coupled by the field, a suitable experimental arrangement permitting the efficient separation of their respective contributions should have the scattering wavevector q oriented along principal directions. At [{\bf q}\parallel{\bf z}], the [{\rm A}{_1}] phonons are longitudinal ([{\rm LO}{_\parallel}]) and both E modes ([2{\rm TO}_{\perp}]) are transverse, remaining degenerate, whereas at [{\bf q}\parallel{\bf x}] or [{\bf q}\parallel{\bf y}], the [{\rm A}{_1}] phonons become transverse ([{\rm TO}_{\perp}]) and the E phonons split into a pair of ([{\rm TO}{_\perp}], [{\rm LO}{_\perp}]) modes of different frequencies. The subscripts [\parallel] or [\perp] explicitly indicate the orientation of the electric dipole moment carried by the mode with respect to the fourfold axis ([4\parallel{\bf c}\equiv{\bf z}]).

Schematically, the situation (i.e. frequency shifts and splittings) at [{\bf q}\approx 0] can be represented by[\matrix{ & {\bf q}\parallel{\bf z}& {\bf q}\parallel{\bf x} &\cr &&-&{\rm A}_1({\rm TO}_\parallel)\cr {\rm A}_1({\rm LO}_\parallel) &- &&\cr &&-& {\rm E}_x({\rm LO}_\perp)\cr {\rm E}(2{\rm TO}_\perp)&-&-&{\rm E}_y({\rm TO}_\perp)\cr}]

For a general direction of q, the modes are of a mixed character and their frequencies show directional (angular) dispersion. The overall picture depends on the number of [{\rm A}{_1}] and E phonons present in the given crystal, as well as on their effective charges and on the ordering of their eigenfrequencies. In fact, only the [{\rm E}({\rm TO}_\perp)] modes remain unaffected by the directional dispersion.

Table 2.3.3.3[link] gives the corresponding contributions of these modes to the cross section for several representative scattering geometries, where subscripts TO and LO indicate that the components of the total Raman tensor may take on different values for TO and LO modes due to electro-optic contributions in the latter case.

Table 2.3.3.3| top | pdf |
Raman selection rules in crystals of the 4mm class

Scattering configurationCross section for symmetry species
A1E
[{\bf q}\parallel{\bf z}] [z(xx)z, z(yy)z] [\sim|a_{\rm LO}|^2]
[{\bf q}\perp{\bf z}] [x(zz)x, x(zz)y] [\sim|b_{\rm TO}|^2]
[\bar y(xz)y, \bar x(yz)x] [\sim|f_{\rm TO}|^2]
[x'(zx')y', x'(y'z)y'] [{\textstyle{1 \over 2}}|f_{\rm TO}|^2 + {\textstyle{1 \over 2}}|f_{\rm LO}|^2]

2.3.3.6. [\bf q]-dependent terms

| top | pdf |

So far, we have not explicitly considered the dependence of the Raman tensor on the magnitude of the scattering wavevector, assuming [{\bf q}\to 0] (the effects of directional dispersion in the case of scattering by polar modes were briefly mentioned in the preceding section). In some cases, however, the Raman tensors vanish in this limit, or q-dependent corrections to the scattering may appear. Formally, we may expand the susceptibility in a Taylor series in q. The coefficients in this expansion are higher-order susceptibility derivatives taken at [{\bf q}=0]. The symmetry-restricted form of these tensorial coefficients may be determined in the same way as that of the zero-order term, i.e. by decomposing the reducible representation of the third-, fourth- and higher-order polar Cartesian tensors into irreducible components [\Gamma(j)]. General properties of the q-dependent terms can be advantageously discussed in connection with the so-called morphic effects (see Sections 2.3.4[link] and 2.3.5[link]).

2.3.4. Morphic effects in Raman scattering

| top | pdf |

By morphic effects we understand the effects that arise from a reduction of the symmetry of a system caused by the application of external forces. The relevant consequences of morphic effects for Raman scattering are changes in the selection rules. Applications of external forces may, for instance, render it possible to observe scattering by excitations that are otherwise inactive. Again, group-theoretical arguments may be applied to obtain the symmetry-restricted component form of the Raman tensors under applied forces.

It should be noted that under external forces in this sense various `built-in' fields can be included, e.g. electric fields or elastic strains typically occurring near the crystal surfaces. Effects of `intrinsic' macroscopic electric fields associated with long-wavelength LO polar phonons can be treated on the same footing. Spatial-dispersion effects connected with the finiteness of the wavevectors, q or k, may also be included among morphic effects, since they may be regarded as being due to the gradients of the fields (displacement or electric) propagating in the crystal.

2.3.4.1. General remarks

| top | pdf |

Various types of applied forces – in a general sense – can be classified according to symmetry, i.e. according to their transformation properties. Thus a force is characterized as a polar force if it transforms under the symmetry operation of the crystal like a polar tensor of appropriate rank (rank 1: electric field E; rank 2: electric field gradient [\nabla{\bf E}], stress T or strain S). It is an axial force if it transforms like an axial tensor (rank 1: magnetic field H). Here we shall deal briefly with the most important cases within the macroscopic approach of the susceptibility derivatives. We shall treat explicitly the first-order scattering only and neglect, for the moment, q-dependent terms.

In a perturbation approach, the first-order transition susceptibility [\delta\chi] in the presence of an applied force F can be expressed in terms of Raman tensors [{\bf R}{^j}({\bf F})] expanded in powers of F: [\displaylines{\delta\chi ({\bf F}) = \textstyle\sum\limits_j {\bf R}^j ({\bf F})Q{_j },\cr \hbox{ where }{\bf R}{^j}({\bf F}) = {\bf R}{^j}{^0} + {\bf R}{^j}{^F}{\bf F }+ {\textstyle{1 \over 2}}{\bf R}{^j}{^F}{^F}{\bf F}{\bf F} + {\ldots}.\cr\hfill (2.3.4.1)}]Here, [{\bf R}{^j}{^0} = \chi ^{(j)}(0) = ({{{\partial \chi _{\alpha \beta }}/{\partial Q_j }}})] is the zero-field intrinsic Raman tensor, whereas the tensors [\eqalignno{{\bf R}{^j}{^F}{\bf F} &= \left({{\partial ^2 \chi _{\alpha \beta }}\over{\partial Q_j \partial F_\mu }}\right)F_\mu, &\cr {\bf R}^{jFF}{\bf FF} &= \left({{\partial ^3 \chi _{\alpha \beta }}\over {\partial Q_j \partial F_\mu \partial F_\nu }}\right)F_\mu F_\nu \,\,etc.&(2.3.4.2)}]are the force-induced Raman tensors of the respective order in the field, associated with the jth normal mode. The scattering cross section for the jth mode becomes proportional to [|{\bf e}{_S}({\bf R}^{j0} + {\bf R}^{jF}{\bf F} + {\textstyle{1 \over 2}}{\bf R}^{jFF}{\bf F}{\bf F} + {\ldots}){\bf e}{_I}|{^2}], which, in general, may modify the polarization selection rules. If, for example, the mode is intrinsically Raman inactive, i.e. [{\bf R}^{j0} = 0] whereas [{\bf R}^{jF} \neq 0], we deal with purely force-induced Raman scattering; its intensity is proportional to [F{^2}] in the first order. Higher-order terms must be investigated if, for symmetry reasons, the first-order terms vanish.

For force-induced Raman activity, in accordance with general rules, invariance again requires that a particular symmetry species [\Gamma(j)] can contribute to the first-order transition susceptibility by terms of order n in the force only if the identity representation is contained in the reducible representation of the nth-order Raman tensor.

An equivalent formulation is that the nth-order tensor-like coefficients in the corresponding force-induced Raman tensor, i.e. [R_{\alpha \beta \mu\ldots\nu }^{jF\ldots F} = \left({{{\partial ^{1 + \,n}\chi _{\alpha \beta }^{}}\over {\partial Q_j \partial F_\mu\ldots\partial F_\nu }}}\right) \hbox{ in the term }{\bf R}^{jF \ldots F}{\bf F}{\ldots}{\bf F},]vanish identically for symmetry reasons unless [[{\Gamma} _{\rm PV} \otimes {\Gamma}_{\rm PV}] \otimes [{\Gamma}({\bf F})]^{n}_{S}\supset{\Gamma}(j)]. Here [[\Gamma({\bf F})]^{n}_{S} =] [[\Gamma({\bf F})] [\otimes] [ \Gamma({\bf F})] [ \otimes] [\ldots \otimes\Gamma({\bf F})]{_S}] is the symmetrized nth power of the representation [\Gamma({\bf F})] according to which the generalized force F transforms under the operation of the point group. The requirement for the symmetrized part is dictated by the interchangeability of the higher-order derivatives with respect to the components of the force. We recall that the first factor representing the susceptibility, [[\Gamma_{\rm PV} \otimes \Gamma_{\rm PV}]], need not be symmetric in general. However, for most purposes (non-resonant conditions, non-magnetic crystals in the absence of a magnetic field) it can be replaced by its symmetrized part [[\Gamma_{\rm PV} \otimes \Gamma_{\rm PV}]{_S}].

Standard group-theoretical methods can be used to determine the force-induced Raman activity in a given order of the field and to derive the matrix form of the corresponding Raman tensors. Before treating several important cases of morphic effects in more detail in the following sections, let us make a few comments.

Beside the force-induced effects on the scattering tensors, there are also the direct morphic effects of the forces on the excitations themselves (possible frequency shifts, lifting of mode degeneracies etc.), which can be investigated by an analogous perturbation treatment, i.e. by expanding the dynamical matrix in powers of F and determining the corresponding force-induced corrections in the respective orders.

The lifting of degeneracies is a typical sign of the fact that the symmetry of the problem is reduced. The extended system crystal + applied force corresponds to a new symmetry group resulting from those symmetry operations that leave the extended system invariant. Consequently, the new normal modes (in the long-wavelength limit) can be formally classified according to the new point group appropriate for the extended system, which qualitatively accounts for the new reduced symmetries and degeneracies.

The force-induced Raman tensors referring to the original crystal symmetry should thus be equivalent to the Raman tensors of the corresponding modes in the new point group via the compatibility relations. The new point-group symmetry of the extended system is often used to investigate Raman-induced activity. It should be noted, however, that this approach generally fails to predict to what order in the force the induced changes in the Raman tensors appear. Such information is usually of prime importance for the scattering experiment, where appropriate setup and detection techniques can be applied to search for a force-induced effect of a particular order. Thus the perturbation method is usually preferable (Anastassakis, 1980[link]).

In the following sections, we shall briefly treat the most important cases in the conventional limit [{\bf q}\to 0] (neglecting for the moment the spatial dispersion).

2.3.4.2. Electric-field-induced scattering

| top | pdf |

Expanding the linear dielectric susceptibility into a Taylor series in the field, we write [\chi _{\alpha \beta }({\bf E}) = \chi _{\alpha \beta }(0) + {{\partial \chi _{\alpha \beta }}\over {\partial E_\gamma }}E_\gamma + {{\partial ^2 \chi _{\alpha \beta }}\over {\partial E_\gamma \partial E_\delta }}E_\gamma E_\delta +\ldots. \eqno (2.3.4.3)]The coefficients of the field-dependent terms in this expansion are, respectively, third-, fourth- and higher-rank polar tensors; they describe linear, quadratic and higher-order electro-optic effects. The corresponding expansion of the Raman tensor of the jth optic mode is written as [{\bf R}{^j}({\bf E}) =] [{\bf R}^{j0}] [+] [{\bf R}^{jE}{\bf E}] [+] [{\textstyle{1 \over 2}}{\bf R}^{jEE}{\bf EE}] [+] [\ldots].

Since the representation [\Gamma({\bf E}) = \Gamma_{\rm PV}], the coefficients of the linear term in the expansion for [\bf \chi], i.e. the third-rank tensor [b_{\alpha}{_\beta}{_\gamma} = (\partial {\chi}{_\alpha}{_\beta}/{\partial}E_{\gamma})], transform according to the reducible representation given by the direct product: [[{\Gamma}_{\rm PV} \otimes {\Gamma}_{\rm PV}]{_S} \otimes {\Gamma}_{\rm PV}.]First-order field-induced Raman activity (conventional symmetric scattering) is thus obtained by reducing this representation into irreducible components [{\Gamma}(j)]. Higher-order contributions are treated analogously.

It is clear that in centrosymmetric crystals the reduction of a third-rank polar tensor cannot contain even-parity representations; consequently, electric-field-induced scattering by even-parity modes is forbidden in the first order (and in all odd orders) in the field. The lowest non-vanishing contributions to the field-induced Raman tensors of even-parity modes in these crystals are thus quadratic in E; their form is obtained by reducing the representation of a fourth-rank symmetric polar tensor [[{\Gamma}_{\rm PV} \otimes {\Gamma}_{\rm PV}]{_S} \otimes [{\Gamma}_{\rm PV} \otimes {\Gamma}_{\rm PV}]{_S}] into irreducible components [{\Gamma}(j)]. On the other hand, since the electric field removes the centre of inversion, scattering by odd-parity modes becomes allowed in first order in the field but remains forbidden in all even orders. In noncentrosymmetric crystals, parity considerations do not apply.

For completeness, we note that, besides the direct electro-optic contribution to the Raman tensor due to field-induced distortion of the electronic states of the atoms in the unit cell, there are two additional mechanisms contributing to the total first-order change of the dielectric susceptibility in an external electric field E. They come, respectively, from field-induced relative displacements of atoms due to field-induced excitation of polar optical phonons [Q{_p}({\bf E}) \sim {\bf E}] and from field-induced elastic deformation [{\bf S}({\bf E})={\bf dE}] (piezoelectric effect, d being the piezoelectric tensor). In order to separate these contributions, we write formally [\chi({\bf E}) = \chi({\bf E}, Q{_p}({\bf E}),{\bf S}({\bf E}))] and get, to first order in the field,[\eqalignno{\delta \chi ({\bf E}) &= (\partial \chi /\partial {\bf E}){\bf E}+ \textstyle\sum\limits_p (\partial \chi /\partial Q_p)Q_p ({\bf E}) + (\partial \chi /\partial {\bf S}){\bf S}({\bf E}) &\cr&= \textstyle\sum\limits_j {\bf R}^{jE}{\bf E}Q_j \hbox{, where we define}&\cr {\bf R}^{jE} &= (\partial {\bf R}^j / \partial {\bf E}) +\textstyle\sum\limits_p (\partial {\bf R}^j /\partial Q_p)({\rm d}Q_p/{\rm d}{\bf E})+ (\partial{\bf R}^j/\partial{\bf S}){\bf d}. &\cr&&(2.3.4.4)}]

The first term in these equations involves the susceptibility derivative [{\bf b}=(\partial\chi/\partial{\bf E})] at constant [Q{_p}] and S. The second term involves the second-order susceptibility derivatives with respect to the normal coordinates: [{\chi}^{(j,p)}=(\partial^2\chi/\partial Q_j\partial Q_p)=(\partial R^j_{\alpha\beta}/\partial Q_p)].  Since [Q_p({\bf E})\sim Z_{p\nu}E_\nu], where the quantity [{\bf Z}{_p}=(Z_{p \nu})] is the effective charge tensor (2.3.3.4)[link] of the normal mode p, its nonzero contributions are possible only if there are infrared-active optical phonons (for which, in principle, [{\bf Z}{_p}\neq 0]) in the crystal. The third term is proportional to the field-induced elastic strain [{\bf S}({\bf E})={\bf dE}] via the elasto-optic tensor [{\bf p}=(\partial\chi/\partial{\bf S})] and can occur only in piezoelectric crystals.

Example: As an illustration, we derive the matrix form of linear electric-field-induced Raman tensors (including possible antisymmetric part) in a tetragonal crystal of the 4mm class. The corresponding representation [[{\Gamma}_{\rm PV} \otimes {\Gamma}_{\rm PV}] \otimes {\Gamma}_{\rm PV}] in this class reduces as follows: [\displaylines{[{\Gamma}_{\rm PV} \otimes {\Gamma}_{\rm PV}] {_S} \otimes {\Gamma}_{\rm PV}=3{\rm A}{_1} \oplus {\rm A}{_2} \oplus 2{\rm B}{_1} \oplus 2{\rm B}{_2} \oplus 5{\rm E},\cr [{\Gamma}_{\rm PV} \otimes {\Gamma}_{\rm PV}] {_A} \otimes {\Gamma}_{\rm PV}={\rm A}{_1} \oplus 2{\rm A}{_2} \oplus {\rm B}{_1} \oplus {\rm B}{_2} \oplus 2{\rm E}.}]Suitable sets of symmetrized (s) and antisymmetrized (a) basis functions (third-order polynomials) for the representations of the 4mm point group can be easily derived by inspection or using projection operators. The results are given in Table 2.3.4.1[link]. Using these basis functions, one can readily construct the Cartesian form of the linear contributions to the electric-field-induced Raman tensors [{\bf R}{^j}({\bf E})={\bf R}^{jE}{\bf E}] for all symmetry species of the [4mm]-class crystals. The tensors are split into symmetric (conventional allowed scattering) and antisymmetric part.[\matrix{&\hbox{Symmetric}\hfill&\hbox{Antisymmetric}\hfill\cr {\rm A}{_1}:\hfill &{\pmatrix{ {a_1^{}E_z }&. & {a_2^{}E_x } \cr. & {a_1^{}E_z }& {a_2^{}E_y } \cr {a_2^{}E_x }& {a_2^{}E_y }& {b_1^{}E_z } \cr }}\hfill & +{\pmatrix{. &. & {a_3^{}E_x } \cr. &. & {a_3^{}E_y } \cr {- a_3^{}E_x }& {- a_3^{}E_y }&. }}\hfill\cr \cr{\rm A}{_2}:\hfill &{\pmatrix{. &. & {c_2^{}E_y } \cr. &. & {- c_2^{}E_x } \cr {c_2^{}E_y }& {- c_2^{}E_x }&. \cr }}\hfill & +{\pmatrix{. & {c_1^{}E_z }& {c_3^{}E_y } \cr {- c_1^{}E_z }&. & {- c_3^{}E_x } \cr {- c_3^{}E_y }& {c_3^{}E_x }&. \cr }}\hfill\cr {\rm B}{_1}: \hfill&{\pmatrix{ {d_1^{}E_z }&. & {d_2^{}E_x } \cr. & {- d_1^{}E_z }& {- d_2^{}E_y } \cr {d_2^{}E_x }& {- d_2^{}E_y }&. \cr }}\hfill& + {\pmatrix{. &. & {d_3^{}E_x } \cr. &. & {- d_3^{}E_y } \cr {- d_3^{}E_x }& {d_3^{}E_y }&. \cr }}\hfill\cr {\rm B}{_2}:\hfill&{\pmatrix{. & {e_1^{}E_z }& {e_2^{}E_y } \cr {e_1^{}E_z }&. & {e_2^{}E_x } \cr {e_2^{}E_y }& {e_2^{}E_x }&. \cr }}\hfill &+ {\pmatrix{. &. & {e_3^{}E_y } \cr. &. & {e_3^{}E_x } \cr {- e_3^{}E_y }& {- e_3^{}E_x }&. \cr }}\hfill\cr {\rm E}:\hfill & {\pmatrix{ {(f_1 + f_2)E_x }& {f_4 E_y }& {f_5 E_z } \cr {f_4 E_y }& {(f_1 - f_2)E_x }&. \cr {f_5 E_z }&. & {f_3 E_x } \cr }}\hfill &+ {\pmatrix{. & {g_4^{}E_y }& {g_5^{}E_z } \cr {- g_4^{}E_y }&. &. \cr {- g_5^{}E_z }&. &. \cr }}\hfill\cr &{\pmatrix{ {(f_1 - f_2)E_y }& {f_4 E_x }&. \cr {f_4 E_x }& {(f_1 + f_2)E_y }& {f_5 E_z } \cr. & {f_5 E_z }& {f_3 E_y } \cr }}\hfill& + {\pmatrix{. & {- g_4^{}E_x }&. \cr {g_4^{}E_x }&. & {g_5^{}E_z } \cr. & {- g_5^{}E_z }&. \cr }}\hfill}]

Table 2.3.4.1| top | pdf |
Symmetrized (s) and antisymmetrized (a) sets of trilinear basis functions corresponding to symmetry species of the 4mm class

SpeciesBasis functionsSymmetry
A1 [(x{_1}x{_2}+y{_1}y{_2})z{_3}]; [z{_1}z{_2}z{_3}]; [(x{_1}z{_2}+z{_1}x{_2})x{_3}] [+] [(y{_1}z{_2} +z{_1}y{_2})y{_3}] (s)
[(x{_1}z{_2}-z{_1}x{_2})x{_3}] [+] [(y{_1}z{_2}-z{_1}y{_2})y{_3}] (a)
A2 [(x{_1}z{_2}+z{_1}x{_2})y{_3}][(y{_1}z{_2}+z{_1}y{_2})x{_3}] (s)
[(x{_1}y{_2}-y{_1}x{_2})z{_3}]; [(x{_1}z{_2}-z{_1}x{_2})y{_3}][(y{_1}z{_2}-z{_1}y{_2})x{_3}] (a)
B1 [(x{_1}x{_2}-y{_1}y{_2})z{_3}]; [(x{_1}z{_2}+z{_1}x{_2})x{_3}][(y{_1}z{_2}+z{_1}y{_2})y{_3}] (s)
[(x{_1}z{_2}-z{_1}x{_2})x{_3}][(y{_1}z{_2}-z{_1}y{_2})y{_3}] (a)
B2 [(x{_1}y{_2}+y{_1}x{_2})z{_3}]; [(x{_1}z{_2}+z{_1}x{_2})y{_3}] [+] [(y{_1}z{_2}+z{_1}y{_2})x{_3}] (s)
[(x{_1}z{_2}-z{_1}x{_2})y{_3}] [+] [(y{_1}z{_2}-z{_1}y{_2})x{_3}] (a)
E [[(x{_1}x{_2}+y{_1}y{_2})x{_3},(x {_1}x{_2}+y{_1}y{_2})y{_3}]]; [[z{_1}z{_2}x{_3}, z{_1}z{_2}y{_3}]]; [[(x{_1}z{_2}+z{_1}x{_2})z{_3}, (y{_1}z{_2}+z{_1}y{_2})z{_3}]]; [[(x{_1}x{_2}-y{_1}y{_2})x{_3},] [-(x{_1}x{_2}-y{_1}y{_2})y{_3}]]; [[(x{_1}y{_2}+y{_1}x{_2})y{_3}], [(x{_1}y{_2}+y{_1}x{_2})x{_3}]] (s)
[[(x{_1}z{_2}-z{_1}x{_2})z{_3}, (y{_1}z{_1}-z{_1}y{_2})z{_3}]]; [[(x{_1}y{_2}-y{_1}x{_2})y{_3}], [-(x{_1}y{_2}-y{_1}x{_2})x{_3}]] (a)

2.3.4.3. Raman scattering in a magnetic field

| top | pdf |

In a magnetic field, the dielectric susceptibility tensor of a crystal is known to obey the general relation (Onsager reciprocity theorem for generalized kinetic coefficients)[{\chi}_{\alpha\beta}({\bf H}) = {\chi}_{\beta\alpha}(-{\bf H}). \eqno (2.3.4.5)]Further, in the absence of absorption, the susceptibility must be Hermitian, i.e. [{\chi}_{\alpha\beta}({\bf H}) = {\chi}_{\beta\alpha}^*({\bf H}). \eqno (2.3.4.6)]Hence, [{\boldchi}({\bf H})] is neither symmetric nor real. Expanding [{\boldchi}({\bf H})] in the powers of the field, [\chi _{\alpha \beta }({\bf H}) = \chi _{\alpha \beta }(0) + {{\partial \chi _{\alpha \beta }}\over {\partial H_\mu }}H_\mu + {{\partial ^2 \chi _{\alpha \beta }}\over {\partial H_\mu \partial H_\nu }}H_\mu H_\nu + \ldots, \eqno (2.3.4.7)]it follows that all terms of the magnetic-field-induced Raman tensor that are of odd powers in H are purely imaginary and antisymmetric in α and β, whereas all terms of even powers in H are real and symmetric.

Let us discuss in more detail the symmetry properties of the first-order term, which can be written as [\Delta \chi _{\alpha \beta }({\bf H}) = i f_{\alpha \beta \mu }H_\mu, \eqno (2.3.4.8)]where the tensor f, referred to as the magneto-optic tensor, is real and purely antisymmetric in the first two indices: [f_{\alpha \beta \nu} \equiv - i(\partial{\chi}_{\alpha\beta}/\partial H {_\nu}) = - f_{\beta\alpha\nu}.]The representation [\Gamma({\bf f})] of the magneto-optic tensor f may thus be symbolically written as [\eqalignno{\Gamma({\bf f}) &= [\Gamma_{\rm PV} \otimes \Gamma_{\rm PV}]{_A} \otimes \Gamma_{\rm AV} = \Gamma_{\rm AV} \otimes \Gamma_{\rm AV} =\Gamma_{\rm PV} \otimes \Gamma_{\rm PV} &\cr &= \Gamma(T{_\alpha}T{_\beta}), &(2.3.4.9)}]since the antisymmetric part of the product of two polar vectors transforms like an axial vector, and the product of two axial vectors transforms exactly like the product of two polar vectors. Hence, the representation [\Gamma({\bf f})] is equivalent to the representation of a general nonsymmetric second-rank tensor and reduces in exactly the same way (2.3.3.14)[link]. [\Gamma({\bf f}) = \Gamma_{\rm PV} \otimes \Gamma_{\rm PV} = c^{(1)}\Gamma(1) \oplus c^{(2)}\Gamma(2) \oplus \ldots.]

We arrive thus at the important conclusion that, to first order in the field, only the modes that normally show intrinsic Raman activity (either symmetric and antisymmetric) can take part in magnetic-field-induced scattering. Moreover, the magnetic-field-induced Raman tensors for these symmetry species must have the same number of components as the general nonsymmetric Raman tensors at zero field.

In order to determine the symmetry-restricted matrix form of the corresponding field-induced Raman tensors (linear in H) in Cartesian coordinates, one can use the general method and construct the tensors from the respective (antisymmetric) basis functions. In this case, however, a simpler method can be adopted, which makes use of the transformation properties of the magneto-optic tensor as follows.

From the definition of the tensor f, it is clear that its Cartesian components [f_{\alpha\beta\nu}] must have the same symmetry properties as the product [[E{_\alpha}E{_\beta}]{_A}H{_\nu}]. The antisymmetric factor [[E{_\alpha}E{_\beta}]{_A}] transforms, however, as [{\varepsilon}_{\alpha \beta\mu}H{_\mu}], where [{\varepsilon}_{ \alpha \beta \mu}] is the fully antisymmetric third-rank pseudotensor (Levi–Civita tensor). Consequently, [f_{\alpha \beta \nu}] must transform in the same way as [{\varepsilon}_{\alpha \beta \mu}H{_\mu}H{_\nu}], which in turn transforms identically to [{\varepsilon}_{\alpha \beta \mu}E{_\mu}E{_\nu}]. Therefore, comparison of the matrices corresponding to the irreducible components [\Gamma(j)] provides a simple mapping between the components of the Cartesian forms of the linear field-induced Raman tensors [{\bf R}{^j}({\bf H})={\bf R}^{jh}{\bf H}] and the intrinsic Raman tensors [{\bf R}^{j0}]. Explicitly, this mapping is given by[R_{\alpha \beta \nu }^{jH} \equiv {{\partial ^2 \chi _{\alpha \beta }}\over {\partial Q_j \partial H_\nu }} = if_{\alpha \beta \nu }^{(j)}\,\, \leftarrow \,\, i\varepsilon _{\alpha \beta \mu } R_{\mu \nu }^{jo}. \eqno (2.3.4.10)]For any given symmetry species, this relation can be used to deduce the matrix form of the first-order field-induced Raman tensors from the tensors given in Table 2.3.3.1[link].

Example: We consider again the 4mm class crystal. The representation [\Gamma({\bf f})] of the magneto-optic tensor f in the [4mm] class reduces as follows: [\Gamma({\bf f}) ={ \Gamma}_{\rm PV} \otimes { \Gamma}_{\rm PV} = 2{\rm A}{_1} \oplus {\rm A}{_2} \oplus {\rm B}{_1} \oplus {\rm B}{_2} \oplus 2{\rm E}.]Straightforward application of the mapping mentioned above then gives the following symmetry-restricted matrix forms of contributions to the magnetic-field-induced Raman tensors [{\bf R}^{jH}{\bf H}] for all symmetry species of the 4mm-class crystals. The number of independent parameters for each species is the same as in the intrinsic nonsymmetric zero-field Raman tensors: [\eqalign{A{_1}: \quad&\pmatrix{. & {ib'H_z }& {- ia'H_y } \cr {- ib'H_z }&. & {ia'H_x } \cr {ia'H_y }& {- ia'H_x }&. \cr } \cr A{_2}:\quad&\pmatrix{. &. & {ic'H_x } \cr. &. & {ic'H_y } \cr {- ic'H_x }& {- ic'H_y }&. \cr }\cr B{_1}:\quad &\pmatrix{. &. & {id'H_y } \cr. &. & {id'H_x } \cr {- id'H_y }& {- id'H_x }&. \cr }\cr B{_2}:\quad&\pmatrix{. &. & {- ie'H_x } \cr. &. & {ie'H_y } \cr {ie'H_x }& {- ieH_y }&. \cr }\cr E:\quad &\pmatrix{. & {ig'H_x }&. \cr {- ig'H_x }&. & {if'H_z } \cr. & {- if'H_z }&. \cr }\cr &\pmatrix{. & {ig'H_y }& {- if'H_z } \cr {- ig'H_y }&. &. \cr {if'H_z }&. &. \cr }.}]Let us note that the conclusions mentioned above apply, strictly speaking, to non-magnetic crystals. In magnetic materials in the presence of spontaneous ordering (ferro- or antiferromagnetic crystals) the analysis has to be based on magnetic point groups.

2.3.4.4. Stress- (strain-) induced Raman scattering

| top | pdf |

Stress-induced Raman scattering is an example of the case when the external `force' is a higher-rank tensor. In the case of stress, we deal with a symmetric second-rank tensor. Since symmetric stress (T) and strain (S) tensors have the same symmetry and are uniquely related via the fourth-rank elastic stiffness tensor (c),[T_{\alpha \beta} = c_{\alpha \beta \mu \nu} S_{\mu \nu},]it is immaterial for symmetry purposes whether stress- or strain-induced effects are considered. The linear strain-induced contribution to the susceptibility can be written as [\Delta \chi _{\alpha \beta } ({\bf S}) = \left({{{\partial \chi _{\alpha \beta } }\over {\partial S_{\mu \nu }}}}\right) S_{\mu \nu }]so that the respective strain coefficients (conventional symmetric scattering) transform evidently as [[{\Gamma}_{\rm PV} \otimes {\Gamma} _{\rm PV}]{_ S} \otimes [{ \Gamma}_{\rm PV} \otimes {\Gamma}_{\rm PV}]{_S}, ]i.e. they have the same symmetry as the piezo-optic or elasto-optic tensor. Reducing this representation into irreducible components [\Gamma(j)], we obtain the symmetry-restricted form of the linear strain-induced Raman tensors. Evidently, their matrix form is the same as for quadratic electric-field-induced Raman tensors. In centrosymmetric crystals, strain-induced Raman scattering (in any order in the strain) is thus allowed for even-parity modes only.

2.3.5. Spatial-dispersion effects

| top | pdf |

For [{\bf q}=0], the normal modes correspond to a homogeneous phonon displacement pattern (all cells vibrate in phase). Phenomenologically, the q-dependence of Raman tensors can be understood as a kind of morphic effect due to the gradients of the displacement field. Developing the contribution of the long-wavelength jth normal mode to the susceptibility in Cartesian components of the displacement of atoms in the primitive cell and their gradients, we obtain [\eqalignno{\delta \chi _{\alpha \beta }^{(j)}({\bf q}) &= \sum\limits_\kappa \left({{{\partial \chi _{\alpha \beta }}\over {\partial u_{\kappa, \gamma }^{(j)}}}}\right)_0 u_{\kappa, \gamma }^{(j)} ({\bf q}) + i \sum\limits_\kappa \left({{{\partial \chi _{\alpha \beta }}\over {\partial (\nabla u_{\kappa, \gamma }^{(j)})_\delta }}}\right)_0 q_\delta u_{\kappa, \gamma }^{(j)}({\bf q}), &\cr&& (2.3.5.1)}]where the derivatives are taken at [{\bf q}=0], and we use the obvious relation [\nabla u_{\kappa, \gamma }^{(j)} = i{\bf q}u_{\kappa, \gamma }^{(j)}].

Transforming to normal coordinates, using (2.3.3.1)[link], we identify the [{\bf q}=0] intrinsic Raman tensor [{\bf R}^{j0}] of the jth normal mode, explicitly expressed via Cartesian displacements of atoms, [R_{\alpha \beta }^{jo} \equiv \chi _{\alpha \beta }^{(j)}(0) \equiv \left({{{\partial \chi _{\alpha \beta }^{}}\over {\partial Q_j }}}\right) = \sum\limits_\kappa { \left({{{\partial \chi _{\alpha \beta } }\over {\partial u_{\kappa, \mu }^{(j)}}}}\right)}{{e_{\kappa, \mu } (0,j)}\over {\sqrt {Nm_\kappa }}}, \eqno (2.3.5.2)]and introduce the first-order q-induced atomic displacement Raman tensor coefficients [{\bf R}^{jq}]: [\eqalignno{R_{\alpha \beta \gamma }^{jq} &\equiv - i\left({{{\partial \chi _{\alpha \beta }^{(j)}}\over {\partial q_\gamma }}}\right) = - i\left({{{\partial ^2 \chi _{\alpha \beta }}\over {\partial Q_j \partial q_\gamma }}}\right) &\cr&= \sum\limits_\kappa { \left({{{\partial \chi _{\alpha \beta }^{(j)}}\over {\partial (\nabla u_{\kappa, \mu } )_\gamma }}}\right)}{{e_{\kappa, \mu } (0,j)}\over {\sqrt {Nm_\kappa }}}. &(2.3.5.3)}]Hence, to the lowest order in q, the transition susceptibility is expressed as [\delta \chi _{\alpha \beta }^{(j)}({\bf q}) \cong \left({R_{\alpha \beta }^{j0} + iR_{\alpha \beta \gamma }^{jq}q_\gamma }\right)Q_j(0). \eqno (2.3.5.4)]

In a more general case, spatial dispersion should be considered together with the electro-optic contributions due to the internal macroscopic field E and its gradients. Assuming the linear susceptibility to be modulated by the atomic displacements [Q{_j}] and the macroscopic electric field E as well as by their gradients [\nabla Q{_j}] and [\nabla{\bf E}], we can expand the transition susceptibility of the jth phonon mode [Q{_j}({\bf q})] to terms linear in q and formally separate the atomic displacement and electro-optic parts of the Raman tensor [see (2.3.3.15)[link]]: [\displaylines{\eqalign{\delta {\boldchi}^{(j)}({\bf q}) &= (\partial{\boldchi}/\partial Q_{j})Q_{j}({\bf q}) + i (\partial{\boldchi}/\partial \nabla Q_{j}){\bf q}Q_{j}({\bf q}) \cr &\quad+ (\partial{\boldchi}/\partial{\bf E}){\bf E}^{j}({\bf q}) + i (\partial{\boldchi}/\partial \nabla{\bf E}) {\bf qE}{^j}({\bf q}),}\cr\hbox{or concisely}\hfill\cr \delta {\boldchi}^{(j)}({\bf q}) ={\bf a}{^j}({\bf q})Q{_j}({\bf q}) + {\bf b}({\bf q}){\bf E}{^j}({\bf q}),\cr\hbox{with }\hfill\cr\hfill{\bf a}{^j}({\bf q})= ({\bf a}^{j0} + i {\bf a}^{jq}{\bf q}),\quad {\bf b}({\bf q}) = ({\bf b}{^0} + i {\bf b}{^q}{\bf q}). \hfill(2.3.5.5)}]Hence, setting [{\bf E}{^j}({\bf q}) = ({\rm d}{\bf E}/{\rm d}Q_j)Q_j({\bf q})], we write for the total Raman tensor [{\bf R}{^j}({\bf q}) = {\bf a}{^j}({\bf q}) + {\bf b}({\bf q}) ({\rm d}{\bf E}/{\rm d}Q{_j}). \eqno (2.3.5.6)]The definitions of the tensors [{\bf a}{^j}] and [{\bf a}^{jq}] correspond to (2.3.5.2)[link] and (2.3.5.3)[link]. Analogously, the tensors [{\bf b}{^0}] and [{\bf b}{^q}] are defined by [b_{\alpha \beta \gamma }^0 = \left({{{\partial \chi _{\alpha \beta } }\over {\partial E_\gamma }}}\right),\quad b_{\alpha \beta \gamma \delta }^q = - i\left({{{\partial ^2 \chi _{\alpha \beta } }\over {\partial q_\gamma \partial E_\delta }}}\right) = \left({{{\partial \chi _{\alpha \beta } }\over {\partial (\nabla E_\delta)_\gamma }}}\right). \eqno (2.3.5.7)]

The q-independent part [{\bf a}^{j0}] of the atomic displacement Raman tensor corresponds to the standard [{\bf q}= 0] Raman tensors [{\bf R}^{j0}], whose symmetry properties and matrix form were discussed in Section 2.3.3.3.[link]

Like [{\bf a}^{j0}], the form of the q-dependent contribution [{\bf a}^{jq}{\bf q}] also depends on the symmetry properties of the corresponding normal coordinate [Q{_j}({\bf q}\approx 0)]. Since q (or [\nabla Q{_j}]) and E are polar vectors, the symmetry properties of the [{\bf a}^{jq}] are identical to those of the coefficients [{\bf R}^{jE}] of the electric-field-induced Raman tensor discussed in Section 2.3.4.2[link]: they transform according to the [{\bf q} = 0] representation [[{\Gamma}_{\rm PV} \otimes {\Gamma}_{\rm PV}] {_S} \otimes {\Gamma}_{\rm PV} \otimes {\Gamma}(j)]. Hence, the symmetry-restricted matrix form of the q-dependent contribution to the atomic displacement Raman tensor [{\bf a}^{jq}{\bf q}] is exactly the same as that of the corresponding electric-field-induced Raman tensor [{\bf R}^{jE}{\bf E}]. In general, these linear terms must vanish for even-parity modes in centrosymmetric crystals (where the lowest-order non-vanishing contributions to the Raman tensor are quadratic). The third-rank tensor [{\bf b}{^0}] corresponds to the first-order susceptibility derivative b and the fourth-rank tensor [{\bf b}{^q}] to its first-order q-dependent part.

As mentioned above, the q-independent third-rank polar tensor [{\bf b}{^0}] is nonzero only in noncentrosymmetric (piezoelectric) crystals, where it contributes to Raman tensors for polar longitudinal optical (LO) phonons. The corresponding electro-optic terms in [{\bf R}{^j}], connected with the accompanying longitudinal electric field E, are given by [{\bf b}{^0}({\rm d}{\bf E}/{\rm d}Q_j)]. Symmetry arguments imply that the q-independent part of such terms must have the same form as the atomic displacement Raman tensor for polar LO phonons, since the corresponding normal coordinates transform as components of polar vectors.

The q-dependent part of the electro-optic contribution, the polar fourth-rank tensor [{\bf b}{^q}=(b^{q}_{\alpha\beta\gamma\delta})], transforms as [[{\Gamma}_{\rm PV} \otimes {\Gamma}_{\rm PV}]{^S} \otimes {\Gamma}_{\rm PV} \otimes {\Gamma}_{\rm PV}] and its symmetry properties are similar to those of the quadratic electro-optic tensor (however, as q and E are not interchangeable, there is no symmetry in the last two indices [\gamma,\delta]). Thus, for finite q, the term [{\bf b}{^q}{\bf q}] has nonzero components in all crystal classes. The corresponding contribution to the Raman tensor, however, is possible only in noncentrosymmetric crystals. Again, because of the ([{\rm d}{\bf E}/{\rm d}Q_j]) factor, the symmetry-restricted matrix form of this contribution for polar LO phonons will be equivalent to that of the [{\bf a}{^j}{\bf q}] term. As far as symmetry is concerned, the distinction between atomic displacement and electro-optic contributions is therefore immaterial.

The occurrence of q-dependent terms leads to polarization selection rules that are generally different from those of intrinsically ([{\bf q} = 0]) Raman-active modes. For this reason, this phenomenon is sometimes referred to as forbidden scattering. It is often observed under resonance conditions.

We recall that the terms linear in q, i.e. proportional to displacement gradients and elastic strains, are fundamental for the description of inelastic light scattering by those excitations for which the transition susceptibility identically vanishes in the limit [{\bf q} \to 0] regardless of lattice point symmetry. This is the case for scattering by acoustic phonons (Brillouin scattering) and also for scattering by plasma waves (plasmons in semiconductors).

We have explicitly considered only the q-dependent effects, due to gradients of phonon fields. In general, spatial dispersion may be also due to the gradients of the electric (or magnetic) field of the incident or scattered photons. The corresponding effects are often referred to as k-dependent effects. In view of the wavevector conservation condition (2.3.2.2)[link], the three wavevectors q, [{\bf k}{_I}] and [{\bf k}{_S}] are always related, which simplifies the symmetry analysis. Without going into details, we note that, microscopically, the k-dependent effects come from photon–electron interactions beyond the usual dipole approximation, i.e. from multipolar effects. The symmetry-allowed matrix form of the k-induced contributions, depending on the nature of the leading microscopic mechanism, can be obtained by standard group-theoretical techniques.

Example: As an example we give the symmetry-restricted form for the linear q-dependent contribution to the Raman tensors for Raman-inactive triply degenerate [{\rm F}_{1u}] modes in the m3m ([O{_h}]) class:[\eqalign{{\bf R}{^j} &= 0\semi\cr {\bf R}^{jq}{\bf q} &={\pmatrix{ {cq_x }& {bq_y }& {bq_z } \cr {bq_y }& {aq_x }&. \cr {bq_z }&. & {aq_x } \cr }}, {\pmatrix{ {aq_y }& {bq_x }&. \cr {bq_x }& {cq_y }& {bq_z } \cr. & {bq_z }& {aq_y } \cr }},\cr&\quad \,\,{\pmatrix{ {aq_z }&. & {bq_x } \cr. & {aq_z }& {bq_y } \cr {bq_x }& {bq_y }& {cq_z } \cr }}.}]

The same form of the q-dependent contribution holds for Raman-active optic F2 modes in the [\bar 43m] ([T{_d}]) class. The conventional ([{\bf q }= 0]) intrinsic Raman tensor of F2 is nonzero, but has off-diagonal components only. Since these modes are also infrared-active, there is a concomitant splitting of LO and TO frequencies as well as a possible electro-optic contribution to the Raman tensor due to the accompanying longitudinal macroscopic field. If one chooses [{\bf q}\parallel{\bf z}], for instance, the first two matrices correspond to two degenerate TO modes and the third one to the LO mode. Combining the q-independent and q-dependent contributions, we get for each triplet of (2TO [+] LO) F2 modes [\eqalign{{\bf R}^{j0} + {\bf R}^{jq}{\bf q} = &{\pmatrix{. &. & {bq_z } \cr. &. & {c_{\rm TO}} \cr {bq_z }& {c_{\rm TO}}&. \cr }}\quad j{_1 }={\rm F}{_2 }({\rm TO}{_x }),\cr & {\pmatrix{. &. & {c_{\rm TO}} \cr. &. & {bq_z } \cr {c_{\rm TO}}& {bq_z }&. \cr }}\quad j{_2 } = {\rm F}{_2 }({\rm TO}{_y }), \cr&{\pmatrix{ {aq_z }& {c_{\rm LO}}&. \cr {c_{\rm LO}}& {aq_z }&. \cr. &. & {dq_z } \cr }}\quad j{_3 } = {\rm F}{_2 }({\rm LO}{_z }).}]

The difference in the parameters [c_{\rm TO}] and [c_{\rm LO}] is due to the electro-optic contribution to the Raman tensor for the LO mode. In the back-scattering geometry for scattering from the (001) face of the crystal, only the LOz modes can take part. Intrinsic allowed components [|c_{\rm LO}|{^2}] are observable in the crossed (xy) polarization geometry, whereas the q-dependent terms [|aq{_z}|^2] appear as `forbidden' scattering in parallel (xx) or (yy) geometries.

2.3.6. Higher-order scattering

| top | pdf |

In higher-order processes, the scattering involves participation of two or more quanta (j and [j']) of the elementary excitations. Let us discuss briefly the second-order scattering by phonons, where the energy and wavevector conservation conditions read [\eqalignno{{\omega}{_I}-{\omega}{_S} & ={\omega}=\pm {\omega}{_j} \pm {\omega}_{j'},&\cr {\bf k}{_I}-{\bf k}{_S}&={\bf q}=\pm {\bf q}{_j}\pm{\bf q}_{j'}.&(2.3.6.1)}]The combinations of signs in these equations correspond to four possibilities, in which either both phonons, j and [j'], are created (Stokes process: [+ +]), both annihilated (anti-Stokes process: −−), or one is created and the other annihilated (difference process: [+ -], [- +]). If in the Stokes or anti-Stokes case both excitations are of the same type, [j = j'], one speaks of overtones. The corresponding terms in the transition susceptibility are the coefficients of a bilinear combination of normal coordinates in the expansion of [\delta\boldchi].

In the quasi-static limit, the transition susceptibilities for the second-order scattering correspond, again, to the susceptibility derivatives. Thus, the spectral differential cross section for the second-order scattering (Stokes component) can be formally written as [\eqalign{{{\rm d}^2 \sigma }\over {{\rm d}\omega \,\,{\rm d}\Omega }&\approx \sum\limits_{{j,j'} \atop {{\bf q}_j + {\bf q}_{j'} = {\bf q}\approx 0}}\left| e_{S\alpha}{{{\partial ^2 \chi _{\alpha\beta}}\over {\partial Q_j^* ({\bf q}_j)\partial Q_{j'}^* ({\bf q}_{j'})}}}e_{I\beta}\right|^2 \cr&\cr &\quad\times\delta [{\omega _j ({\bf q}_j) + \omega _{j'}({\bf q}_{j'}) - \omega }],}]with [{\omega} = {\omega}{_I} - {\omega}{_S} \,\gt \, 0]. In this formula, we have suppressed the universal factors [see (2.3.3.5)[link]] and the explicit expression for the response function (thermal factors). Instead, the delta function (response function in the limit of zero damping) expresses the energy-conservation condition.

The wavevector selection rules in the long-wavelength limit, with [{\bf q }= 0], imply that [{\bf q}{_j} = -{\bf q}{_{j'}}] (the same holds for anti-Stokes components, while [{\bf q}{_j} = {\bf q}_{j'}] for difference scattering), so the wavevectors themselves need not be small and, in principle, scattering by phonons with all wavevectors from the Brillouin zone can be observed.

Without invoking any symmetry arguments for the Raman activity, such as the restrictions imposed by crystal symmetry on the susceptibility derivatives, it is clear that the intensity of second-order scattering at a frequency [\omega] is controlled by the number of those combinations of phonons whose frequencies obey [\omega=\omega{_j}({\bf q})+\omega_{j'}(-{\bf q})]. The quantity determining this number is the combined density of states of phonon pairs, i.e. [\rho _2 (\omega) = \textstyle\sum\limits_{j,j'}{\textstyle\sum\limits_{\bf q}{\delta [{\omega _j ({\bf q}) + \omega _{j'}({\bf q}) - \omega }]}}. \eqno (2.3.6.2)]This function can be calculated provided the dispersion curves [\omega{_j}({\bf q})] of the excitations are known. The density of states is a continuous function and shows features known as the van Hove singularities corresponding to the critical points, where one or more components of the gradient [\nabla_{\bf q}[\omega_j({\bf q})+\omega_{j'}({\bf q})]] vanish. Most of the critical points occur for wavevectors on the boundary, where the vanishing gradients of the individual dispersion curves are often dictated by the crystal symmetry, but they also occur in those regions of the reciprocal space where both dispersion curves have opposite or equal slopes at the same wavevector q. To a first approximation, the second-order spectrum is thus essentially continuous, reflecting the two-phonon density of states, with peaks and sharp features at frequencies close to the positions of the van Hove singularities. This is to be contrasted with the first-order scattering, where (in perfect crystals) only single peaks corresponding to long-wavelength ([{\bf q}\approx 0]) phonons occur.

Group-theoretical arguments may again be invoked in deriving the selection rules that determine the Raman activity of a particular combination of excitations (Birman, 1974[link]). The susceptibility derivative again transforms as a tensor. For a given pair of excitations ([j,{\bf q}]) and ([j',-{\bf q}]) responsible for the modulation, the combined excitation symmetry is obtained by taking the direct product of the irreducible representations of the space group corresponding to the participating excitations, [\Gamma (j,j') = {\rm D}_{j,{\bf q}} \otimes {\rm D}_{j', - {\bf q}}. \eqno (2.3.6.3)]The representation [\Gamma (j,j')], unlike [{\rm D}_{j, {\bf q}}], corresponds to a zero-wavevector representation of the crystal space group and is therefore equivalent to a (reducible) representation of the crystal point group. It can be decomposed into irreducible components. Raman scattering of the pair is allowed if a Raman-active [{\bf q}=0] representation is contained in this decomposition of [\Gamma (j,j')] or, alternatively, if the product [[{\Gamma}_{\rm PV} \otimes {\Gamma}_{\rm PV} ]{_S} \otimes {\Gamma}(j,j')]  contains the totally symmetric representation [\Gamma(1)].

The selection rules for the second-order scattering are, in general, far less restrictive than in the first-order case. For example, it can be shown that for a general wavevector q in the Brillouin zone there are no selection rules on the participation of phonons in the second-order scattering, since the representations [\Gamma (j,j')] contain all Raman-active symmetries. In specific crystal structures, however, restrictions occur for the wavevectors corresponding to special symmetry positions (points, lines or planes) in the Brillouin zone. This implies that the selection rules may suppress some of the van Hove singularities in the second-order spectra.

Morphic effects in second-order scattering, due to an applied external force F (see Section 2.3.4.1[link]), may be investigated using the same criteria as in first-order scattering, i.e. decomposing the [{\bf q}=0] representation [{\Gamma}({\bf F}) \otimes [{\Gamma}_{\rm PV} \otimes {\Gamma}_{\rm PV} ]{_S} \otimes {\Gamma}(j,j')] and searching for the matrix form of the corresponding second-order Raman tensors.

Generalization to third- and higher-order processes is obvious.

Concluding this section, we note that in a Raman-scattering experiment, higher-order features in the spectra can in principle be distinguished from first-order features by different behaviour of the differential scattering cross section with temperature. For example, the respective thermal factors entering the expression for the second-order scattering cross section are given in Table 2.3.6.1[link].

Table 2.3.6.1| top | pdf |
Thermal factors for second-order Raman scattering

[n(\omega)] is given by (2.3.3.7)[link].

FactorProcessRaman shift
[[n({\omega}{_j}) + 1][n({\omega}_{j'}) + 1] ] Stokes [({\omega} = {\omega}{_j} + {\omega}_{j'})]
[[n({\omega}{_j}) + 1](n{\omega}_{j'})] Difference [({\omega} = {\omega}{_j} - {\omega}_{j'})]
[n({\omega}{_j})n({\omega}_{j'})] Anti-Stokes [({\omega} = -{\omega}{_j} -{\omega}_{j'})]

2.3.7. Conclusions

| top | pdf |

In this overview of Raman scattering in crystals, we have almost exclusively based our considerations on a phenomenological, semi-classical viewpoint without going into details of the underlying microscopic theory. This is surely an appropriate approach to a discussion of the fundamental consequences of crystal symmetry on the selection rules governing the varied phenomena of inelastic light scattering and on the symmetry-restricted form of the corresponding tensorial quantities encountered in this vast and fruitful field. We have attempted to treat the most important symmetry aspects of the inelastic scattering of light by collective excitations in perfect crystals, concentrating on scattering by optical phonons – in the traditional sense of Raman scattering studies. Within a limited scope, we tried to give some insight into the nature of the phenomena relevant in connection with this topic. Our coverage is certainly not exhaustive (nor original); we have also deliberately omitted all scattering phenomena connected with purely electronic excitations, although the corresponding symmetry aspects can be analysed on the same footing. The essence of the truth is rather simple: As long as the excitations are characterized by irreducible representations of the symmetry group of the system, the well proven tools of the theory of representations are at hand to work out the consequences (if common sense does not readily provide the answer).

For an experimental physicist, symmetry analysis is invaluable in designing the experimental arrangement expected to give evidence of some particular phenomenon. The type of answer one may expect from such analysis is characteristic of the group-theoretical predictions. One can learn whether a particular effect is allowed or forbidden, and – if it is allowed – one can learn how many independent coefficients may be needed to describe it. Symmetry analysis alone cannot predict how large the effect is, or why it occurs, though it may well indicate the way to make it more easily observable. In order to understand or justify its actual magnitude, one has to analyse in greater detail the underlying microphysical mechanisms; however, the analysis must be consistent with restrictions dictated by the symmetry of the problem and, here again, symmetry arguments with powerful group-theoretical tools provide reliable guidance.

For the sake of brevity, citations have been largely suppressed in the present overview. Instead, the relevant sources are included in the list of references.

For detailed information about all the varied aspects of light scattering in solids and recent advances in this vast and fruitful field, the reader is referred to specialized monographs, e.g. Turrell (1972[link]), Hayes & Loudon (1978[link]), and to a comprehensive series edited by Cardona & Güntherodt (1975–[link]).

2.3.8. Glossary

| top | pdf |

[\alpha, \beta, \gamma,\delta,\ldots,\mu,\nu]: Greek indices are used for Cartesian components of vectors and tensors (summation over repeated indices is understood).

[{\bf a}_j] atomic displacement Raman tensor of the jth phonon mode, defined as [\partial \boldchi/\partial Q_j]
[\bf b] susceptibility derivative [\partial\boldchi/\partial{\bf E}] (closely related to the electro-optic tensor [\bf r)]
c velocity of light
[\bf c] elastic stiffness tensor
[\bf d] piezoelectric tensor
[\boldchi] susceptibility tensor
[{\bf e}_{I,S}] polarization vector of incident and scattered light
[{\bf e}_{\kappa}({\bf q},j)] phonon eigenvector of the jth phonon mode of wavevector [\bf q]
[\bf E] electric field intensity
[\bf f] magneto-optic tensor
[\gamma_{j,{\bf q}}] damping constant of mode [(j,{\bf q})]
[\bf H] magnetic field intensity
[\varepsilon_0] permittivity of free space
[\boldvarepsilon] permittivity tensor
[{\bf k}_{I,S}] wavevector of incident and scattered light
[n_{I,S}] refractive index for incident and scattered light
[\bf P] polarization
[\bf q] phonon wavevector
[Q_j({\bf q})] normal coordinate (amplitude of the jth phonon mode of wavevector [\bf q)]
[{\bf R}^j] Raman tensor of jth phonon mode
[\bf S] strain tensor
[\bf T] stress tensor
[\sigma] scattering cross section
[\omega] angular frequency
[{\bf Z}({\bf q})] effective charge matrix

References

Agranovich, V. M. & Ginzburg, V. L. (1984). Crystal optics with spatial dispersion, and excitons. Berlin, Heidelberg, New York, Tokyo: Springer.
Anastassakis, E. M. (1980). In Dynamical properties of solids, edited by G. K. Horton & A. A. Maradudin, Vol. 4, pp. 159–375. Amsterdam, New York, Oxford: North-Holland.
Birman, J. L. (1974). Theory of crystal space groups and infrared and Raman lattice processes of insulating crystals. Handbuch der Physik, Vol. 25/2b. Berlin, Heidelberg, New York: Springer.
Burstein, E. & Pinczuk, A. (1971). The physics of opto-electronic materials, edited by W. A. Albers, pp. 33–79. New York: Plenum Press.
Callen, H. (1968). Crystal symmetry and macroscopic laws. Am. J. Phys. 36, 735–748.
Cardona, M. & Güntherodt, G. (1975–). Editors. Light scattering in solids, Vols. I–IV. Berlin, Heidelberg, New York, Tokyo: Springer.
Claus, R., Merten, L. & Brandmüller, J. (1975). Light scattering by phonon-polaritons. Berlin, Heidelberg, New York: Springer.
Hayes, W. & Loudon, R. (1978). Scattering of light by crystals. New York: John Wiley & Sons.
Lax, M. J. (1974). Symmetry principles in solid state and molecular physics. New York: John Wiley & Sons.
Nye, J. F. (1964). Physical properties of crystals. Oxford: Clarendon Press.
Poulet, H. & Mathieu, J. P. (1976). Vibration spectra and symmetry of crystals. Paris: Gordon & Breach.
Rousseau, D. L., Bauman, R. P. & Porto, S. P. S. (1981). Normal mode determination in crystals. J. Raman Spectrosc. 10, 253–290.
Turrell, G. (1972). Infrared and Raman spectra of crystals. London, New York: Academic Press.








































to end of page
to top of page