Chapter 2.3. Raman scattering

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 , wavevector and polarization vector . Likewise, the scattered radiation is characterized by , and :

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 (also termed the Raman shift) corresponding to the energy transfer from the radiation field to the crystal, and by a scattering wavevector q corresponding to the respective momentum transfer . Since the energy and momentum must be conserved in the scattering process, we have the conditionsStrictly 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 (, ) and frequencies (, ) are related by the dispersion relation , 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 qwhere is the scattering angle (the angle between and ). This relation defines in the () 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.).

In the absence of any excitations, the incident field at frequency induces in the crystal the polarization P, related to the field by the linear dielectric susceptibility tensor ( is the permittivity of free space): The linear susceptibility 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 can be expressed as a Taylor expansion in terms of normal coordinates of the excitations:The tensorial coefficients 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 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 , polarization ) as a result of its scattering (frequency , polarization ) 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) The symbol 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 for a particular Raman band and, in addition, over the solid angle, one obtains, respectively, the differential cross section () and the total cross section (): These quantities are useful in comparing the integrated scattered intensity by different excitations.

In a scattering experiment on crystals, the choice of the scattering geometry implies setting the propagation directions and 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 . 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 . The light scattered in the direction of is collected over a certain finite solid angle . Its polarization is analysed with a suitable polarization analyser, and the scattered intensity as a function of frequency (or Raman shift ) 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: . Thus, for example, the symbol means that right-angle scattering geometry is used, where the incident beam polarized in the direction propagates along the axis, while the scattered beam is collected in the z direction and the polarization analyser is set parallel with the direction. The measured intensity, being proportional to , 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 direction.

In a typical Raman experiment with visible light (), the magnitudes of the wavevectors are of the order of 10⁵ cm⁻¹, much lower than those of the reciprocal-lattice vectors K ( cm⁻¹). Consequently, the range of the magnitudes of the scattering wavevectors q accessible by varying the scattering geometry from ° (forward scattering) to ° (back scattering) is , 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) 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.

In the first-order scattering by a single excitation of the normal mode () of frequency and wavevector , the energy and wavevector conservation conditions give 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 of the excitation. Inserting the plane-wave expressions for the quantities into the first term of the expansion (2.3.2.5) for and comparing the terms with common time dependence, we get for the Stokes polarization due to the normal mode ()which corresponds to a plane wave at a frequency of .

Since, in the harmonic approximation, the normal modes () 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:

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 added to take into account the finite lifetime of excitations in real crystals within the quasi-harmonic approximation): If a fictitious generalized force with Fourier components is applied to the system, the average of the Fourier components of the excitation amplitude and the force are proportional, where the proportionality coefficient is called a linear response function.

The fluctuation–dissipation theorem (or, more exactly, its quantum version) relates the power spectrum 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: where is the Bose–Einstein statistical factor 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): If the damping parameter of the excitations is much smaller than the frequency , as is usually the case for phonons in perfect crystals, this function is closely approximated by a Lorentzian centred at , with a full width at half maximum. Hence, introducing a suitably normalized lineshape function ,we can write for the power spectra (2.3.3.5) of the jth normal mode in the spectral differential cross section 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 .

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), i.e. in the nonlinear susceptibility .

The quantity that controls the symmetry properties of the scattering cross section due to excitation 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: The nonlinear susceptibility tensor 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

In the quasi-static limit, i.e. if the scattering frequency is negligibly small compared with the incident photon frequency (), 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 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, where we suppressed the explicit dependence on and introduced a simplified notation for the Raman tensor , still keeping the dependence on the scattering wavevector.

In deriving the symmetry properties of the Raman tensor 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 is very small compared with the Brillouin-zone dimensions, another conventional approximation is to neglect the q dependence of the susceptibilities. Setting   enables us to analyse the symmetry of the Raman tensor in terms of the factor group , 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 , 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 and at the full symmetries applicable in the 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), so that the scattering cross section of the jth mode is If there are no antisymmetric components in the susceptibility derivatives, it can be shown that the fourth-rank tensor has at most 21 independent components, as for the elastic constants tensor.

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 , the transformation properties of the incident and scattered light are described by the three-dimensional polar vector representation of the appropriate point group of the crystal, since the quantities that characterize the light (, , ) are all polar vectors, i.e. first-rank polar tensors (). The transformation properties of a normal mode j must correspond to an irreducible matrix representation 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 must be contained in the decomposition of the direct Kronecker product representation at least once: In this case, the Kronecker product 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 (interaction Hamiltonian), such that the product is the reducible representation of the Raman tensor. If condition (2.3.3.12) holds, then the Raman tensor does not vanish identically and may have at least one independent nonzero component. As the representation 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 must transform like one or more components of a polar tensor. The transition susceptibility transforms accordingly. The task of determining whether a given normal mode j is Raman active or not thus consists of simply decomposing the representation and identifying the irreducible components .

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 , 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 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 and and the polarization vectors and 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 of the irreducible representation in the decomposition where the multiplicity coefficient corresponds to the number of times the given irreducible representation enters the decomposition. If the representation is two- or three-dimensional, then for each occurrence of in (2.3.3.14) 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 – can be readily constructed by finding the appropriate bilinear basis functions that transform according to the corresponding irreducible representation . The required number of such independent bases is given by the multiplicity coefficient . Alternatively, one may construct invariant polynomials (transforming as scalars) of order four, i.e. of the same order as the product .

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) is not necessarily symmetric in the α and β indices as long as the fields and correspond to different frequencies (inelastic scattering). However, each second-rank polar tensor (nine components), transforming according to , can be decomposed into a symmetric part (six components), transforming like a symmetric polar tensor , and an antisymmetric part (three components), transforming like an axial vector (for the definition of axial tensors, see Section 1.1.4.5.3 ) according to .

The symmetry-restricted forms of the () Raman tensors corresponding to all Raman-active symmetry species are summarized in Table 2.3.3.1 (see e.g. Hayes & Loudon, 1978) 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 , see (2.3.3.4), 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).

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   1   Monoclinic, unique axis z Monoclinic, unique axis y Orthorhombic   222     Tetragonal     422         Trigonal   3   Hexagonal   6         622       Cubic   23           432

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 (or , ) 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 , trigonal classes 32, 3m and , as well as hexagonal class .

In those point groups that contain the inversion operation, i.e. in the eleven centrosymmetric (nonpolar) crystal classes 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.

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 (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 class Scattering configuration Cross section for symmetry species Back scattering Right-angle scattering — — — — — — — — — — — — — — — —

If, for some reason, antisymmetric scattering is allowed, possible contribution of the modes should be considered as well. The contribution to cross section from these modes is proportional to , hence it can be distinguished from the contribution of the symmetry species by a suitable choice of the scattering geometry.

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, , 422, 4mm, , 3, 32, 3m, 6, , 622, 6mm, , 23, 432, ) that are simultaneously infrared-active. Since these polar modes carry a nonzero macroscopic effective charge (2.3.3.4), 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 () 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. ) 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 . In lattice dynamics, the limit can be treated correctly by taking into account the retardation effects in the range where cq becomes comparable to , 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 . Experimentally, this polariton region is partially accessible only in near-forward Raman scattering [see (2.3.2.3)]. For larger scattering wavevectors in the usual right-angle or back-scattering geometries, the electrostatic approximation, , 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). 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 , 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 , which is a subgroup of the factor group. Compatibility relations of the representations at 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 or, in terms of the Raman tensor, where we introduce the notation and for the atomic displacement and electro-optic contributions to the Raman tensor . As usual, 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 , 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 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 () 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 class, where of the Raman-active symmetry species the modes and are polar. According to Table 2.3.3.1, their () Raman tensors are identical to those of the and modes in the preceding example of a -class crystal. Owing to the macroscopic electric field, however, here one has to expect directional dispersion of the frequencies of the long wavelength () 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 , the phonons are longitudinal () and both E modes () are transverse, remaining degenerate, whereas at or , the phonons become transverse () and the E phonons split into a pair of (, ) modes of different frequencies. The subscripts or explicitly indicate the orientation of the electric dipole moment carried by the mode with respect to the fourfold axis ().

Schematically, the situation (i.e. frequency shifts and splittings) at can be represented by

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 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 modes remain unaffected by the directional dispersion.

Table 2.3.3.3 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 configuration Cross section for symmetry species A1 E — — — —

So far, we have not explicitly considered the dependence of the Raman tensor on the magnitude of the scattering wavevector, assuming (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 . 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 . General properties of the q-dependent terms can be advantageously discussed in connection with the so-called morphic effects (see Sections 2.3.4 and 2.3.5).

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 , 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 in the presence of an applied force F can be expressed in terms of Raman tensors expanded in powers of F: Here, is the zero-field intrinsic Raman tensor, whereas the tensors 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 , which, in general, may modify the polarization selection rules. If, for example, the mode is intrinsically Raman inactive, i.e. whereas , we deal with purely force-induced Raman scattering; its intensity is proportional to 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 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. vanish identically for symmetry reasons unless . Here is the symmetrized nth power of the representation 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, , 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 .

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).

In the following sections, we shall briefly treat the most important cases in the conventional limit (neglecting for the moment the spatial dispersion).

Expanding the linear dielectric susceptibility into a Taylor series in the field, we write 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 .

Since the representation , the coefficients of the linear term in the expansion for , i.e. the third-rank tensor , transform according to the reducible representation given by the direct product: First-order field-induced Raman activity (conventional symmetric scattering) is thus obtained by reducing this representation into irreducible components . 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 into irreducible components . 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 and from field-induced elastic deformation (piezoelectric effect, d being the piezoelectric tensor). In order to separate these contributions, we write formally and get, to first order in the field,

The first term in these equations involves the susceptibility derivative at constant and S. The second term involves the second-order susceptibility derivatives with respect to the normal coordinates: .  Since , where the quantity is the effective charge tensor (2.3.3.4) of the normal mode p, its nonzero contributions are possible only if there are infrared-active optical phonons (for which, in principle, ) in the crystal. The third term is proportional to the field-induced elastic strain via the elasto-optic tensor 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 in this class reduces as follows: 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. Using these basis functions, one can readily construct the Cartesian form of the linear contributions to the electric-field-induced Raman tensors for all symmetry species of the -class crystals. The tensors are split into symmetric (conventional allowed scattering) and antisymmetric part.

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 Species Basis functions Symmetry A1 ; ; (s) (a) A2 − (s) ; − (a) B1 ; − (s) − (a) B2 ; (s) (a) E ; ; ; ]; , (s) ; , (a)

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)Further, in the absence of absorption, the susceptibility must be Hermitian, i.e. Hence, is neither symmetric nor real. Expanding in the powers of the field, 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 where the tensor f, referred to as the magneto-optic tensor, is real and purely antisymmetric in the first two indices: The representation of the magneto-optic tensor f may thus be symbolically written as 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 is equivalent to the representation of a general nonsymmetric second-rank tensor and reduces in exactly the same way (2.3.3.14).

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 must have the same symmetry properties as the product . The antisymmetric factor transforms, however, as , where is the fully antisymmetric third-rank pseudotensor (Levi–Civita tensor). Consequently, must transform in the same way as , which in turn transforms identically to . Therefore, comparison of the matrices corresponding to the irreducible components provides a simple mapping between the components of the Cartesian forms of the linear field-induced Raman tensors and the intrinsic Raman tensors . Explicitly, this mapping is given byFor 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.

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),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 so that the respective strain coefficients (conventional symmetric scattering) transform evidently as i.e. they have the same symmetry as the piezo-optic or elasto-optic tensor. Reducing this representation into irreducible components , 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.