
Coherent acoustic and optical phonon oscillations in Bi nanowire samples were studied with a femtosecond pump-probe technique. Laser pulses of 50 fs excited simultaneously acoustic oscillations at a frequency of about 9.5 GHz and optical phonons in the terahertz range. The transmission signal of nanowires on a glass substrate and the signal of light scattered from freestanding nanowires were measured. The acoustic velocity in nanowires was found to be close to that of bulk polycrystalline material. The changes in the optical phonon frequency at different laser fluences were simulated taking into account excitation inhomogeneity, lattice anharmonicity, diffusion, and recombination of the carriers and gave good agreement with experimental results. ©2008 American Institute of Physics
In the past, probing of scattered light proved to be a powerful tool for studies of condensed matter,1 in particular Raman spectra of optical phonons in Bi were measured.2 More recently, excitation of acoustic3 and optical phonon4,5,6 modes in films of different materials, including Bi,7,8 with short laser pulses was performed. Such pulses excite all active vibrations with frequencies contained in the pulse envelope.9 The laser light produces electronic excitations, which are coupled to the lattice vibrations via deformation potential3 and stimulated Raman scattering.9 After action of a short laser pulse the equilibrium positions of the nuclei are displaced, and they start oscillating around these new equilibrium positions. Thus the lattice distortions give rise to coherent acoustic and optical phonons. The excitation of the lattice oscillations of optical phonons was explained via the displacive mechanism,10 which was shown to be related to a particular manifestation of the impulsive stimulated Raman scattering mechanism in absorbing materials.11,12 A significant interest presents the application of the diagnostics with short laser pulses to nanosystems, such as nanoparticles of various shapes.13,14,15,16 In nanomaterials high frequency acoustic modes appear due to confinement of coherent acoustic phonons. Acoustic oscillations reflect the elastic properties, shape, and dimensions of the nanostructure, while optical phonon frequencies depend on the rigidity and the effective mass of the atoms involved in the optical phonon oscillation.
We investigated optical transient responses from samples of Bi nanowires with a femtosecond pump-probe technique. Bi nanowires are of special interest because of their exceptional thermoelectric properties.17 In nanosystems the surface contribution to the free energy is essential, so that their properties can be size dependent and different from bulk materials.18 The structure of the electronic density of states and the carrier transport properties of a nanowire, which can be considered as a quasi-one-dimensional quantum system, is also modified compared to the bulk material due to a confinement effect in the lateral directions.18 The Bi crystalline lattice (A7, 3m point group symmetry) presents a distorted simple-cubic structure with two atoms per a primitive cell and is especially sensitive to electronic excitations, shifting the equilibrium lattice position.19 A short laser pulse excites Raman active modes near the center of the Brillouin zone; for Bi these are the totally symmetric breathing A1g and the doubly degenerate Eg mode.20 Note that in femtosecond pump-probe experiments with bulk Bi samples the magnitude of the A1g mode was dominant, so that for the observation of the Eg mode an increase in the fluence, a manipulation of the pump/probe polarizations, or lowering the temperature was necessary.8,21,22 For bulk Te and Bi softening of the lattice and reduction in the optical phonon frequency were observed at high excitation levels.6,8 We demonstrate detection of both acoustic and optical phonon modes in Bi nanowires excited by femtosecond laser pulses and also consider major factors affecting the dynamics of the optical phonons and changes of their frequency at high excitation levels.
For measurements we used femtosecond laser pulses at a central wavelength of 810 nm with a repetition rate of 1 kHz, a pulse energy of up to 0.75 mJ, and a pulse duration of about 50 fs. The initial laser beam was split into two parts. The 90% part was modulated by a mechanical chopper at a frequency of 500 Hz and used as a source of the pump pulses. The 10% beam was directed onto a beta barium borate crystal, and the generated second harmonic at 405 nm was used as a probe. The intensity of each beam was additionally adjusted by optical neutral density filters. The pump and probe beams were focused on the sample into spots with diameters of about 0.6 and 0.4 mm, respectively. The delay between the pump and probe beams was varied by a computer-controlled translation stage. For transmission measurements at low sample densities a photodiode was used. With high density samples the light scattered from the sample surface in a nonspecular direction was collected by a lens and registered by a cooled photomultiplier. The signal from the photodetector was then measured with a lock-in amplifier.
The samples of Bi nanowires were fabricated by electrochemically depositing metal into the pores of a membrane by a process described in more detail elsewhere.23 We studied nanowires deposited on a glass substrate (the covered surface fraction of ~10%) and also densely packed freestanding nanowires [see scanning electron microscope (SEM) and transmission electron microscope (TEM) images in Fig. 1]. In the latter case the nanowires were submerged in water to reduce heating. For the nanowires on a substrate, the transmission of the probe beam was measured and for freestanding nanowires the transmission through the sample was below 1%; therefore measurements with the scattered light of the probe beam were performed.
Figure 1. The resultant length of the nanowires was about 15 µm and the filling factor, defined as the volume fraction of the nanowires in the sample, was f~0.5. From TEM images of Bi nanowires [Fig. 1(c)] the diameters of the nanowires were measured to be d
200±20 nm. The electron diffraction analysis of the crystalline structure shows that nanowires consist of relatively large crystallites with a typical length of about 1 µm.
The transmission response of a nanowire sample on a glass substrate is shown in Fig. 2. When nanowires were present on the substrate surface, the tail of the response contained terahertz oscillations. The signal from the substrate quickly decreased on a time scale of about 100 fs. After this the signal reflects response of the nanowires, and it was fitted by an expression
The fit produced the following values of interest:
1=0.48 ps, f=1.7 THz,
=0.13,
=0.60 rad, and
2=2.1 ps. The fit is shown in Fig. 2(a) as a solid line and also for the oscillating contribution the first term of Eq. (1) is plotted as a dotted line in the inset.
Figure 2. The optical response of the probe beam scattered from a sample of freestranding Bi nanowires is shown in Fig. 3. The energy of the pump pulses at the sample was set to 6 µJ, which corresponded to a relatively low fluence F=1.5 mJ/cm2. The energy of the probe pulses was 100 nJ. The observed signal (curve A) exhibited an initial peak followed by oscillations and was fitted with the same expression of Eq. (1), providing fitting parameters:
1=283 ps, f=7.7 GHz,
=0.11,
=0.35 rad, and
2=278 ps.
Figure 3. The initial peak was also recorded with a higher temporal resolution, revealing terahertz frequency oscillations in the initial portion of the signal. These high frequency oscillations are shown in the inset of Fig. 3(a) (curve B), and they are characterized by the following fitting parameters:
1=2.3 ps, f=2.6 THz,
=0.01,
=0.2 rad, and
2=2.7 ps.
The spectrum of the scattered light signals reveals an intensity modulation with two spectral peaks at f1=9.5 GHz and f2=2.8 THz shown in Fig. 3(b). As the following analysis confirms, frequencies f1,2 are determined by the excited coherent acoustic and optical phonons, respectively.
To describe the observed signals, we note that coherent acoustic and optical phonon oscillations modulate the dielectric constant of the material,
where Q1 and Q2 are the lattice displacements corresponding to acoustic and optical phonons and
1 and
2 describe the coupling of these oscillations to the material polarization. The acoustic oscillations are related to radial longitudinal eigenmodes of a nanowire. The excitation of optical phonons is produced by the abrupt change in the electronic energy-state distribution function. Both processes can be described by the model of a damped harmonic oscillator,12,15,24
In the case of acoustic oscillations Q1 is the change in the nanowire diameter,
![<i>F</i><sub>1</sub> = <i>F</i><sub><i>l</i></sub> + <i>F</i><sub><i>e</i></sub>; <i>F</i><sub><i>l</i></sub> = <i>Omega</i><sub>0</sub><sup>2</sup><i>Q</i><sub>0</sub>[(1/3<i>B</i>)<i>gamma</i> <i>C</i><sub><i>l</i></sub>(<i>T</i><sub><i>l</i></sub> − <i>T</i><sub>0</sub>)]<i>g</i><sub><i>l</i></sub>(<i>t</i>),](103110_1m3.gif)
![<i>F</i><sub><i>e</i></sub> = <i>Omega</i><sub>0</sub><sup>2</sup><i>Q</i><sub>0</sub>[(2/9<i>B</i>)<i>C</i><sub><i>e</i></sub>(<i>T</i><sub><i>e</i></sub>)(<i>T</i><sub><i>e</i></sub> − <i>T</i><sub>0</sub>)]<i>g</i><sub><i>e</i></sub>(<i>t</i>),](103110_1m4.gif)
where Cl,e are thermal capacities of the lattice and the electrons, B is the bulk modulus, and the time dependences can be approximated by

![<i>g</i><sub><i>l</i></sub>(<i>t</i>) = {0,<i>t</i> < 0;[1 − exp(−<i>t</i>/<i>tau</i><sub><i>e</i>-ph</sub>)],<i>t</i> >= 0},](103110_1m6.gif)
where
e-ph is the thermalization time of the electron-phonon system and
is the Grüneisen constant of the lattice.15,24 The two contributions to the force Fe,l are due to the impulsive action of the electronic pressure25,26 and the contribution of the lattice heating. The frequency
1 corresponds to the characteristic radial breathing acoustic mode of the nanowires,
where the eigen-number
=2.21 is determined by the boundary problem of radial oscillations of a long rod.27 Since the impedance of the metal is much larger than the impedance of the surrounding medium
cl(Bi)/
cl(water))
25, the free boundary conditions are valid. Taking the value of the longitudinal sound velocity for a polycrystalline Bi cl=2200 m/s,28 we obtain from Eq. (4) f1=7.75 GHz, which is within the range of frequencies of the first spectral peak in Fig. 3(b).
It follows from our previous study23 that the main contribution to the excitation of Bi nanowires comes from the lattice heating. Taking into account only lattice contribution, we obtain for the coherent acoustic vibration after the action of the laser pulse (t>
0, where
0 is the pulse duration)
where
(
) is the Fourier transform of the normalized temporal dependence of the laser pulse intensity, and
1,±=−i/
d,1±(
−1/
)1/2.
For description of optical phonon oscillations we use Eq. (3) with Q2 being the displacement of atoms in the optical phonon oscillation. The force F2 can be found as an average of the electron-phonon coupling operator F=

=
0.11,12 This approach takes into account also the resonant interaction of light with the energy levels, causing strong absorption, and expresses the force in terms of a modified Raman tensor.12 Consequently, the solution of Eq. (2) for t>
0 is
where Im(
) is the imaginary part of the dielectric permittivity of the material (Bi),
is the optical wavelength, V is the irradiated volume, Nph is the flux of incident photons,
2 is the optical phonon frequency, and
2,±=−i/
d,2±(
−1/
)1/2.
For relatively weak damping,
1,2

, the amplitude of the oscillation is proportional to the Fourier component of the intensity envelope at the phonon frequency.9,12 For a short pulse (as in our experiments)
(
1,2)![[approximate]](/stockgif3/ap.gif)
(0), and the oscillating part of the solution of Eq. (3) for t
0 can be presented as a damped cosine function, which corresponds to the oscillating contribution in Eq. (1).
There are several factors that also affect the observed optical response. Electronic excitations produce a contribution to the dielectric function that is changing in time as a result of such processes as thermalization of the electron-phonon system and the recombination of carriers promoted during the excitation to antibonding states. The contribution of these processes is reflected by the second and third terms in the fitting function of Eq. (1).
To clarify the relative roles of different factors on the formation of the optical response related to optical phonons we used the following model. Electronic excitations lead to a softening of the lattice and a redshift of the optical phonon frequency with increasing density of the excitations.6 The effect of the lattice anharmonicity was considered in Refs. 8,29. While in Ref. 8 the optical photon frequency change with the increasing excitation level was mainly attributed to the lattice anharmonicity, the calculations and experiments of Ref. 29 led to conclusion that below n
1.25%, where n is the density of the photoexcited electron-hole plasma expressed as a percentage of the excited valence electrons [for Bi n=1% corresponds to a carrier density of 1.43×1021 cm−3 (Ref. 29)], the role of the lattice anharmonicity is small. It was also found that with n approaching nc
2.04% the frequency of the optical phonon oscillations, calculated with the model of the displacive mechanism, approaches zero, and the effect of the anharmonicity is essential. At n
nc the energy of the oscillating motion reaches the middle peak (the Peierls distortion30) between two valleys in the Bi lattice potential. In the vicinity of this peak, which also somewhat decreases with n, the velocity of the motion decreases, resulting also in the increase of the period of oscillations. The frequency of the small-amplitude oscillations, reflecting the curvature of the potential near the equilibrium position becomes zero at a higher value, ncs
2.71%. To take into account these conclusions we approximated the frequency by the following expression:
where
0=2
f0 and f0 is the optical phonon frequency at the initial temperature [in our experiments this is the room temperature, and we take the frequency value corresponding to the bulk Bi, f0=2.92 THz (Refs. 8,20)]; A is the amplitude of the optical phonon oscillation and Amax(n) is the maximum value of A for a given n, which is realized just after the excitation, at t
0. If A=Amax, then
(n,Amax)=
0g1(n) and if A
0, then
(n,A)
0g2(n). The functions g1(n) and g2(n) were calculated in Ref. 29; however, our calculations have shown that in the region n
0.5% the decrease in frequency with the functions of Ref. 29 is much smaller compared to experimentally observed one (see further discussion). To account for this, we used the following approximations:



where the power indices p1(n)=p2(n)=
+
n change linearly with the excitation density. The expressions of Eq. (8) are a generalization of the square root dependence31 used for description of the red frequency shift of optical phonons in Te. The constants
=0.5 and
=0.3 were chosen to provide the average slope of the calculated dependence of
on F close to the value following from our observations and those of Ref. 32.
The excitation is strongly inhomogeneous and decays inside the material due to absorption. Since the diameter of the nanowires is much larger than the penetration depth of the optical radiation, we assumed a plane geometry in finding the solution for the diffusion of carriers. Taking into account diffusion and relaxation of the carriers and assuming that the surface recombination is relatively small, the solution for n can be presented in the form
where
r is the carrier relaxation time, D is the ambipolar diffusion constant, and l1 is the penetration depth of the pump beam [l1
22 nm at 810 nm (Refs. 28,33)]. We assumed that n is proportional to the laser energy flux, n(%)=F/Fn, where the value of the normalization constant Fn=4.0 mJ/cm2, which we estimated from the data of Refs. 29,34. We note that estimates based on the results of Ref. 8 gave a close value Fn=3.6 mJ/cm2. The decay of the optical phonon amplitude is dominated by the relaxation process with the time constant
d; consequently in Eq. (7) we use A/Amax
exp(−t/
d), and then the signal S(t) is proportional to the integral of the polarization modulation induced by the optical phonon [see Eq. (2)] in the material,
where n(x,t
0) is the excitation density just after the action of the laser pulse and l2 is the penetration depth of the probe beam [l2
24 nm at 400 nm (Refs. 33,35)]. We assume the response to be proportional to n(x,t
0) since the excitation of the coherent optical phonon occurs as a result of an abrupt change in the density of the excitations, and the following relaxation takes place with characteristic times substantially exceeding the period of the optical phonon. For a purely displacive excitation it holds
=0 10; however since the parameters 
0 and 1/(
r) are not infinitely small, also a phase shift can appear. With the frequency of optical phonon depending on the excitation level, as in our experiments, the accurate determination of the phase from experimental data becomes difficult, and therefore in calculations the initial phase shift
was disregarded.
The temporal behavior of the calculated optical phonon frequency is shown in Fig. 4. In the calculations we assumed the carrier relaxation time equal to the decay time of the average component of the signal, determined by the fit to the experimental data [time
2 in Eq. (1)]. Thus, we use
r=2.1 ps for F=10 mJ/cm2 and
r
2.7 ps for F=1.5 mJ/cm2, which are close to the values presented in Refs. 34,36. We assume the decay time of the optical phonon equal to the time
1, determined by the fit of Eq. (1), thus
d
0.48 ps for F=10 mJ/cm2 and
d
2.3 ps for F=1.5 mJ/cm2. These values are consistent with measurements of Ref. 36. The diffusion coefficient D=1.3 cm2/s was used for F=10 mJ/cm2 following the estimate37 based on the values of mobility measured in Ref. 38 and the Einstein relation as well as the discussion in Refs. 37,39. Although the diffusion coefficient tends to decrease with increasing density of excitations, for F=1.5 mJ/cm2 the photogenerated carriers diffused perpendicular to the length of the nanowires deposited on a substrate, in which case the mobility is expected to be reduced.40 Therefore we used the same relatively low value D=1.3 cm2/s also for F=1.5 mJ/cm2.
Figure 4. The observation of optical phonons was performed in nanowires deposited on a substrate and with freestanding nanowires. In the first case the transmission signal was measured at a relatively high laser fluence (required in view of the low level of the signal). For freestanding densely packed nanowires the signal of the scattered light was measured, and it contained oscillations from two frequency ranges: around 9.5 GHz and in the interval of 2.1–2.9 THz, which corresponded to acoustic and optical phonons, respectively.
For acoustic phonons the experimentally determined damping coefficient
=1/
d,1 includes homogeneous part
h due to inherent damping of the material and radiation losses into surrounding medium as well as an inhomogeneous part
i. The latter is determined by the size distribution of the nanowires, resulting in dephasing of the oscillations.24 Estimates show that in our experiments the damping of acoustic phonons was mainly due to the dephasing. Indeed, the estimate of the corresponding time constant
i
(1/
i)=d/(2
f
d),24 where
d
20 nm is the dispersion of the nanowire diameter distribution, gives the decay time
i
290 ps, which is close to the experimental value 280 ps.
For optical phonons some inhomogeneous damping can also be expected due to the inhomogeneous excitation of the material; however our simulations have shown that its effect is small for both a relatively small, F=1.5 mJ/cm2, and high, F=10 mJ/cm2, laser excitation fluences. We observed a significant reduction in the optical phonon frequency (down to 2.1 THz) and a faster decay at the high excitation levels. The reduction in the frequency was the largest just after the action of the laser pulse. Such redshifts in Bi films down to 2.45 THz (Ref. 8) and even 2.12 THz (Ref. 34) were observed at high excitation levels. The nature of this shift was the subject of discussion37,39 and a subsequent study.29 It was suggested8,39 that the anharmonicity of the lattice potential is mainly responsible for the observed shifts. According to Refs. 29,37 the effect of the anharmonicity in Bi is much smaller and practically negligible for n
1.25%; however for higher excitation levels [as in our case for F=10 mJ/cm2 (n=2.5%)] the influence of the anharmonicity can be noticeable.
We developed a model taking into account the softening of the lattice and its anharmonicity, inhomogeneity of the excitation, and relaxation and diffusion of carriers. This model reasonably well described the observed dynamics of the optical phonon frequency. As it follows from the results of Ref. 29, the anharmonicity effect plays a minor role at almost all excitation densities of interest, except for a relatively narrow interval in the vicinity of the critical density of excitations nc, at which the frequency becomes zero. Therefore except for this interval of the excitation densities close to nc (this value also somewhat changes due to anharmonicity) the inhomogeneity of the excitation and the carrier relaxation and diffusion mostly account for the observed dynamics of photoexcited optical phonons. We used approximation for the dependence of the optical phonon frequency on the density of the excitations somewhat different from the one calculated in Ref. 29 to better account for the slope of the dependence
(n) for relatively small n. For n
0.22% (F
0.9 mJ/cm2) the measured average slope 
/
F
−0.17 THz/(mJ/cm2) (Ref. 32) is close to the calculation with our model, 
/
F
−0.14 THz/(mJ/cm2), while
(n) of Ref. 29 (which account for inhomogeneous excitation and diffusion) gives a significantly lower value, 
/
F
−0.014 THz/(mJ/cm2). Note that the inhomogeneity of the excitation effectively reduces the observed frequency shift because the signal is observed from the penetration depth of the laser radiation including layers with a smaller excitation density.
In conclusion, with the femtosecond two-color pump-probe technique the simultaneous excitation of acoustic and optical phonons in Bi nanowires was registered. The presence of nanowires on a glass substrate resulted in a modulation of the transmitted probe beam at the optical phonon frequency. The transient responses were also measured for freestanding Bi nanowires, in which case both coherent acoustic and optical phonons were detected. The frequencies of the latter decreased with increasing laser fluence. The frequency of the acoustic breathing mode was close to the value calculated from the equation for eigenvalues of the oscillations with a longitudinal velocity corresponding to polycrystalline Bi. The developed model of the optical phonon dynamics accounted for softening and anharmonicity of the lattice, inhomogeneity of the excitation, as well as for relaxation and diffusion of the carriers. It reasonably well described the observed optical phonon frequency downshift for relatively low as well as for high excitation levels. Note that the parallel observation of the acoustic and optical phonons can be also useful in establishing quantitative relations between the main parameters responsible for these processes.
This work was supported by the Robert A. Welch Foundation under Grant Nos. A1546 and A1562 and by NSF under Grant Nos. DMR-0551813, DMR-0606529 and MRI-0722800 and by the Qatar National Research Fund, Grant NPRP30-6-7-35. The FE-SEM acquisition was supported by the NSF under Grant No. DBI-0116835.
Full figure (26 kB)Fig. 1. Investigated Bi nanowires: (a) SEM image of freestanding nanowires, (magnification of 2.6×103), (b) SEM image (view of nanowires from top, 1.4×104), and (c) TEM image (2×104). First citation in article
Full figure (28 kB)Fig. 2. (Color online) (a) The transmission signal from the Bi nanowires on a glass substrate showing oscillations (curve A) and the signal from a bare substrate (curve B, exhibiting no such oscillations); the laser fluence was F=10 mJ/cm2. The inset shows the oscillating contribution to signal A with a fit shown by a dotted line. (b) The spectrum of the oscillating part of signal A, showing a broad peak around 2.2 THz. First citation in article
Full figure (30 kB)Fig. 3. (Color online) Transient signals of scattered light from Bi nanowires at a relatively low pump fluence of F=1.5 mJ/cm2: (a) the temporal response is shown. For curve A the resolution, determined by the translational speed of the delay line, was set to ~1 ps, and in this case acoustic phonons were detected. For curve B (inset) the resolution was set to ~100 fs; this higher temporal resolution revealed the coherent optical phonon oscillations. (b) Spectral analysis of the oscillating portions of signals A and B: two spectral peaks are observed with maxima at f1=9.5 GHz and f2=2.8 THz. First citation in article
Full figure (25 kB)Fig. 4. (Color online) The variations in the optical phonon frequency from experimental data and simulations for F=10 and 1.5 mJ/cm2. The lines connecting experimental points are a guide to the eye. First citation in article
*Electronic mail: a-kolomenski@physics.tamu.edu.
†Also at Microscopy and Imaging Center, Texas A&M University, College Station, Texas 77843-2257.
Up: Issue Table of Contents
Go to: Previous Article | Next Article
Other formats: HTML (smaller files) | PDF ( kB)