Brillouin scattering originates from the interaction of an incident radiation with thermal acoustic vibrations in matter. The phenomenon was predicted by Brillouin in 1922 (Brillouin, 1922
) and first observed in light scattering by Gross (Gross, 1930a
,b
). However, owing to specific spectrometric difficulties, precise experimental studies of Brillouin lines in crystals were not performed until the 1960s (Cecchi, 1964
; Benedek & Fritsch, 1966
; Gornall & Stoicheff, 1970
) and Brillouin scattering became commonly used for the investigation of elastic properties of condensed matter with the advent of laser sources and multipass Fabry–Perot interferometers (Hariharan & Sen, 1961
; Sandercock, 1971
). More recently, Brillouin scattering of neutrons (Egelstaff et al., 1989
) and X-rays (Sette et al., 1998
) has been observed.
Brillouin scattering of light probes long-wavelength acoustic phonons. Thus, the detailed atomic structure is irrelevant and the vibrations of the scattering medium are determined by macroscopic parameters, in particular the density ρ and the elastic coefficients
. For this reason, Brillouin scattering is observed in gases, in liquids and in crystals as well as in disordered solids.
Vacher & Boyer (1972
) and Cummins & Schoen (1972
) have performed a detailed investigation of the selection rules for Brillouin scattering in materials of various symmetries. In this chapter, calculations of the sound velocities and scattered intensities for the most commonly investigated vibrational modes in bulk condensed matter are presented. Brillouin scattering from surfaces will not be discussed. The current state of the art for Brillouin spectroscopy is also briefly summarized.
The fundamental equation of dynamics (see Section 1.3.4.2
), applied to the displacement u of an elementary volume at r in a homogeneous material is
Summation over repeated indices will always be implied, and T is the stress tensor. In non-piezoelectric media, the constitutive equation for small strains S is simply
The strain being the symmetrized spatial derivative of u, and c being symmetric upon interchange of k and
, the introduction of (2.4.2.2)
in (2.4.2.1)
gives (see also Section 1.3.4.2
)
One considers harmonic plane-wave solutions of wavevector Q and frequency ω,
For
small compared with the wavelength
, the total derivative
can be replaced by the partial
in (2.4.2.3)
. Introducing (2.4.2.4)
into (2.4.2.3)
, one obtains
where
is the unit vector in the propagation direction,
is the unit tensor and
, where
is the phase velocity of the wave. This shows that
is an eigenvector of the tensor
. For a given propagation direction
, the three eigenvalues
are obtained by solving
To each
there is an eigenvector
given by (2.4.2.5)
and an associated phase velocity ![[V^{(s)} = \sqrt {C^{(s)}/\rho }. \eqno (2.4.2.7)]](/teximages/dach2o4/dach2o4fd7.svg)
The tensor
is symmetric upon interchange of the indices (
) because
. Its eigenvalues are real positive, and the three directions of vibration
are mutually perpendicular. The notation
indicates a unit vector. The tensor
is also invariant upon a change of sign of the propagation direction. This implies that the solution of (2.4.2.5)
is the same for all symmetry classes belonging to the same Laue class.
For a general direction
, and for a symmetry lower than isotropic,
is neither parallel nor perpendicular to
, so that the modes are neither purely longitudinal nor purely transverse. In this case (2.4.2.6)
is also difficult to solve. The situation is much simpler when
is parallel to a symmetry axis of the Laue class. Then, one of the vibrations is purely longitudinal (LA), while the other two are purely transverse (TA). A pure mode also exists when
belongs to a symmetry plane of the Laue class, in which case there is a transverse vibration with
perpendicular to the symmetry plane. For all these pure mode directions, (2.4.2.6)
can be factorized to obtain simple analytical solutions. In this chapter, only pure mode directions are considered.
In piezoelectric crystals, a stress component is also produced by the internal electric field E, so that the constitutive equation (2.4.2.2)
has an additional term (see Section 1.1.5.2
),
where e is the piezoelectric tensor at constant strain.
The electrical displacement vector D, related to E by the dielectric tensor
, also contains a contribution from the strain,
where
is at the frequency of the elastic wave.
In the absence of free charges,
, and (2.4.2.9)
provides a relation between E and S,
For long waves, it can be shown that E and Q are parallel. (2.4.2.10)
can then be solved for E, and this value is replaced in (2.4.2.8)
to give
Comparing (2.4.2.11)
and (2.4.2.2)
, one sees that the effective elastic tensor
now depends on the propagation direction
. Otherwise, all considerations of the previous section, starting from (2.4.2.6)
, remain, with c simply replaced by
.
The change in the relative optical dielectric tensor
produced by an elastic wave is usually expressed in terms of the strain, using the Pockels piezo-optic tensor p, as
The elastic wave should, however, be characterized by both strain S and rotation A (Nelson & Lax, 1971
; see also Section 1.3.1.3
):
The square brackets on the left-hand side are there to emphasize that the component is antisymmetric upon interchange of the indices,
. For birefringent crystals, the rotations induce a change of the local
in the laboratory frame. In this case, (2.4.3.1)
must be replaced by
where
is the new piezo-optic tensor given by
One finds for the rotational part
If the principal axes of the dielectric tensor coincide with the crystallographic axes, this gives
This is the expression used in this chapter, as monoclinic and triclinic groups are not listed in the tables below.
For the calculation of the Brillouin scattering, it is more convenient to use
which is valid for small
.
Brillouin scattering occurs when an incident photon at frequency
interacts with the crystal to either produce or absorb an acoustic phonon at
, while a scattered photon at
is simultaneously emitted. Conservation of energy gives
where positive
corresponds to the anti-Stokes process. Conservation of momentum can be written
where Q is the wavevector of the emitted phonon, and
,
are those of the scattered and incident photons, respectively. One can define unit vectors q in the direction of the wavevectors k by
where n and
are the appropriate refractive indices, and
is the vacuum wavelength of the radiation. Equation (2.4.4.3b)
assumes that
so that
is not appreciably changed in the scattering. The incident and scattered waves have unit polarization vectors
and
, respectively, and corresponding indices n and
. The polarization vectors are the principal directions of vibration derived from the sections of the ellipsoid of indices by planes perpendicular to
and
, respectively. We assume that the electric vector of the light field Eopt is parallel to the displacement Dopt. This is exactly true for many cases listed in the tables below. In the other cases (such as skew directions in the orthorhombic group) this assumes that the birefringence is sufficiently small for the effect of the angle between
and
to be negligible. A full treatment, including this effect, has been given by Nelson et al. (1972
).
After substituting (2.4.4.3)
in (2.4.4.2)
, the unit vector in the direction of the phonon wavevector is given by
The Brillouin shift
is related to the phonon velocity V by
Since
, from (2.4.4.5)
and (2.4.4.3)
, (2.4.4.4)
one finds
where
is the angle between
and
.
The power
, scattered from the illuminated volume V in a solid angle
, where
and
are measured inside the sample, is given by
where
is the incident light intensity inside the material,
is the appropriate elastic constant for the observed phonon, and the factor
results from taking the fluctuation–dissipation theorem in the classical limit for
(Hayes & Loudon, 1978
). The coupling coefficient M is given by
In practice, the incident intensity is defined outside the scattering volume,
, and for normal incidence one can write
Similarly, the scattered power is observed outside as
, and
again for normal incidence. Finally, the approximative relation between the scattering solid angle
, outside the sample, and the solid angle
, in the sample, is
Substituting (2.4.4.9a,b,c)![[link]](/graphics/greenarr.gif)
![[link]](/graphics/greenarr.gif)
in (2.4.4.7)
, one obtains (Vacher & Boyer, 1972
)
where the coupling coefficient
is
In the cases of interest here, the tensor
is diagonal,
without summation on i, and (2.4.4.11)
can be written in the simpler form ![[\beta = {1 \over {n^4 ({n'})^4 }}{{| {e_i n_i^2 p'_{ijk\ell }\hat u_k \hat Q_\ell e'_j n_j^2 }|^2 }\over C}. \eqno (2.4.4.12)]](/teximages/dach2o4/dach2o4fd35.svg)
The tables in this chapter give information on modes and scattering geometries that are in most common use in the study of hypersound in single crystals. Just as in the case of X-rays, Brillouin scattering is not sensitive to the presence or absence of a centre of symmetry (Friedel, 1913
). Hence, the results are the same for all crystalline classes belonging to the same centric group, also called Laue class. The correspondence between the point groups and the Laue classes analysed here is shown in Table 2.4.5.1
. The monoclinic and triclinic cases, being too cumbersome, will not be treated here.
Crystal system
|
Laue class
|
Point groups
|
Cubic
|
|
|
|
|
Hexagonal
|
|
|
|
|
Tetragonal
|
|
|
|
|
Trigonal
|
|
|
|
|
Orthorhombic
|
O
|
|
|
For tensor components
and
, the tables make use of the usual contracted notation for index pairs running from 1 to 6. However, as the tensor
is not symmetric upon interchange of
, it is necessary to distinguish the order
and
. This is accomplished with the following correspondence:![[\matrix{1,1 \to 1\hfill & 2,2 \to 2\hfill &3,3 \to 3\hfill\cr 1,2 \to 6\hfill & 2,3 \to 4\hfill & 3,1 \to 5\hfill\cr 2,1 \to \bar 6\hfill & 3,2 \to \bar 4\hfill & 1,3 \to \bar 5.\hfill}]](/teximages/dach2o4/dach2o4fd36.svg)
Geometries for longitudinal modes (LA) are listed in Tables 2.4.5.2
to 2.4.5.8![[link]](/graphics/greenarr.gif)
![[link]](/graphics/greenarr.gif)
![[link]](/graphics/greenarr.gif)
![[link]](/graphics/greenarr.gif)
![[link]](/graphics/greenarr.gif)
. The first column gives the direction of the scattering vector
that is parallel to the displacement
. The second column gives the elastic coefficient according to (2.4.2.6)
. In piezoelectric materials, effective elastic coefficients defined in (2.4.2.11)
must be used in this column. The third column gives the direction of the light polarizations
and
, and the last column gives the corresponding coupling coefficient
[equation (2.5.5.11)
]. In general, the strongest scattering intensity is obtained for polarized scattering (
), which is the only situation listed in the tables. In this case, the coupling to light (
) is independent of the scattering angle
, and thus the tables apply to any
value.
Tables 2.4.5.9
to 2.4.5.15![[link]](/graphics/greenarr.gif)
![[link]](/graphics/greenarr.gif)
![[link]](/graphics/greenarr.gif)
![[link]](/graphics/greenarr.gif)
![[link]](/graphics/greenarr.gif)
list the geometries usually used for the observation of TA modes in backscattering (
). In this case,
is always perpendicular to
(pure transverse modes), and
is not necessarily parallel to
. Cases where pure TA modes with
in the plane perpendicular to
are degenerate are indicated by the symbol D in the column for
. For the Pockels tensor components, the notation is
if the rotational term vanishes by symmetry, and it is
otherwise.
Tables 2.4.5.16
to 2.4.5.22![[link]](/graphics/greenarr.gif)
![[link]](/graphics/greenarr.gif)
![[link]](/graphics/greenarr.gif)
![[link]](/graphics/greenarr.gif)
![[link]](/graphics/greenarr.gif)
list the common geometries used for the observation of TA modes in 90° scattering. In these tables, the polarization vector
is always perpendicular to the scattering plane and
is always parallel to the incident wavevector of light q. Owing to birefringence, the scattering vector
does not exactly bisect
and
[equation (2.4.4.4)
]. The tables are written for strict 90° scattering,
, and in the case of birefringence the values of
to be used are listed separately in Table 2.4.5.23
. The latter assumes that the birefringences are not large, so that the values of
are given only to first order in the birefringence.
Brillouin spectroscopy with visible laser light requires observing frequency shifts falling typically in the range ∼1 to ∼100 GHz, or ∼0.03 to ∼3 cm−1. To achieve this with good resolution one mostly employs interferometry. For experiments at very small angles (near forward scattering), photocorrelation spectroscopy can also be used. If the observed frequency shifts are
cm−1, rough measurements of spectra can sometimes be obtained with modern grating instruments. Recently, it has also become possible to perform Brillouin scattering using other excitations, in particular neutrons or X-rays. In these cases, the coupling does not occur via the Pockels effect, and the frequency shifts that are observed are much larger. The following discussion is restricted to optical interferometry.
The most common interferometer that has been used for this purpose is the single-pass planar Fabry–Perot (Born & Wolf, 1993
). Upon illumination with monochromatic light, the frequency response of this instrument is given by the Airy function, which consists of a regular comb of maxima obtained as the optical path separating the mirrors is increased. Successive maxima are separated by
. The ratio of the maxima separation to the width of a single peak is called the finesse F, which increases as the mirror reflectivity increases. The finesse is also limited by the planarity of the mirrors. A practical limit is
. The resolving power of such an instrument is
, where
is the optical thickness. Values of R around
to
can be achieved. It is impractical to increase
above ∼5 cm because the luminosity of the instrument is proportional to
. If higher resolutions are required, one uses a spherical interferometer as described below.
A major limitation of the Fabry–Perot interferometer is its poor contrast, namely the ratio between the maximum and the minimum of the Airy function, which is typically ∼1000. This limits the use of this instrument to samples of very high optical quality, as otherwise the generally weak Brillouin signals are masked by the elastically scattered light. To avert this effect, several passes are made through the same instrument, thus elevating the Airy function to the corresponding power (Hariharan & Sen, 1961
; Sandercock, 1971
). Multiple-pass instruments with three, four or five passes are common. Another limitation of the standard Fabry–Perot interferometer is that the interference pattern is repeated at each order. Hence, if the spectrum has a broad spectral spread, the overlap of adjacent orders can greatly complicate the interpretation of measurements. In this case, tandem instruments can be of considerable help. They consist of two Fabry–Perot interferometers with combs of different periods placed in series (Chantrel, 1959
; Mach et al., 1963
). These are operated around a position where the peak transmission of the first interferometer coincides with that of the second one. The two Fabry–Perot interferometers are scanned simultaneously. With this setup, the successive orders are reduced to small ghosts and overlap is not a problem. A convenient commercial instrument has been designed by Sandercock (1982
).
To achieve higher resolutions, one uses the spherical Fabry–Perot interferometer (Connes, 1958
; Hercher, 1968
). This consists of two spherical mirrors placed in a near-confocal configuration. Their spacing
is scanned over a distance of the order of
. The peculiarity of this instrument is that its luminosity increases with its resolution. One obvious drawback is that a change of resolving power, i.e. of
, requires other mirrors. Of course, the single spherical Fabry–Perot interferometer suffers the same limitations regarding contrast and order overlap that were discussed above for the planar case. Multipassing the spherical Fabry–Perot interferometer is possible but not very convenient. It is preferable to use tandem instruments that combine a multipass planar instrument of low resolution followed by a spherical instrument of high resolution (Pine, 1972
; Vacher, 1972
). To analyse the linewidth of narrow phonon lines, the planar standard is adjusted dynamically to transmit the Brillouin line and the spherical interferometer is scanned across the line. With such a device, resolving powers of ∼
have been achieved. For the dynamical adjustment of this instrument one can use a reference signal near the frequency of the phonon line, which is derived by electro-optic modulation of the exciting laser (Sussner & Vacher, 1979
). In this case, not only the width of the phonon, but also its absolute frequency shift, can be determined with an accuracy of ∼1 MHz. It is obvious that to achieve this kind of resolution, the laser source itself must be appropriately stabilized.
In closing, it should be stressed that the practice of interferometry is still an art that requires suitable skills and training in spite of the availability of commercial instruments. The experimenter must take care of a large number of aspects relating to the optical setup, the collection and acceptance angles of the instruments, spurious reflections and spurious interferences, etc. A full list is too long to be given here. However, when properly executed, interferometry is a fine tool, the performance of which is unequalled in its frequency range.