## Origin of coherent G band phonon spectra in single wall carbon nanotubes

A. R. T. Nugraha , E. H. Hasdeo , G. D. Sanders , C. J. Stanton , R. Saito 1 1 2 2 1

1 Department of Physics, Tohoku University, Sendai 980-8578, Japan 2 Department of Physics, University of Florida, Box 118440, Gainesville, Florida 32611-8440, USA

(Dated: August 7, 2014)

## Abstract

Coherent phonons in single wall carbon nanotubes (SWNTs) are observed as oscillations of the differential absorption coefficient as a function of time by means of pump-probe spectroscopy. For the radial breathing mode (RBM) of a SWNT, the coherent phonon signal is understood to be a result of the modulated diameter-dependent energy gaps due to the coherent RBM phonon oscillations. However, this mechanism might not be the dominant contribution to other phonon modes in the SWNT. In particular, for the G band phonons, which correspond to bond-stretching motions, we find that the modulation of the interatomic dipole matrix element gives rise to a strong coherent G band phonon intensity comparable to the coherent RBM phonon intensity. We also further discuss the dependence of coherent G band and RBM phonon amplitudes on the laser excitation pulse width.

PACS numbers: 78.67.Ch,78.47.J-,73.22.-f,63.22.Gh,63.20.kd

## I. INTRODUCTION

Single wall carbon nanotubes (SWNTs), characterized by the chiral index ( n, m ), 1 have been important materials that provide us with a one-dimensional model system to study the dynamics and interactions between electrons, photons, and phonons. 2 In particular, rapid advances in ultrafast pump-probe spectroscopy have allowed researchers to observe lattice oscillations of SWNTs with the same phase in real time, known as coherent phonon spectroscopy. 3-7 The coherent phonon motions can be observed as oscillations of optical properties, such as the differential transmittance (∆ T/T ) or differential reflectivity (∆ R/R ) as a function of delay time between pump and probe pulses. By performing a Fourier transform of the oscillations of ∆ T/T or ∆ R/R with respect to time, we can obtain the coherent phonon spectra as a function of phonon frequency. Several peaks found in the coherent phonon spectra of a SWNT correspond to Raman active phonon modes, such as the radial breathing modes (RBMs), D bands, G bands, and G ′ bands. 8 Lim et al. showed that even the low-frequency acoustic phonon signals can be observed in purified (6 , 5) SWNTs by coherent phonon spectroscopy because of its ultrafine spectral resolution. 9 Moreover, ultrafast spectroscopy techniques allow us to directly measure phonon dynamics, including phase information, or life time of phonons, in the time domain. 3,4,6

It is known that oscillations of ∆ T/T or ∆ R/R as a function of delay time t between pump and probe pulses in coherent phonon spectroscopy are directly related to the modulations of the absorption coefficient α as a function of the probe energy E probe and t . 10 Therefore, in order to obtain the coherent phonon spectra theoretically, we need to calculate the absorption coefficient α E ( probe , t ) for a given coherent phonon amplitude. In the case of RBMs, the oscillations of α E ( probe , t ) have been understood as a result of energy gap modulations, which are inversely proportional to the nanotube diameter. 4,6 However, in the case of G bands, which are assigned to longitudinal-optical (LO) and in-plane transverse-optical (iTO) phonon modes, 2 it is known that these modes do not significantly modify the energy gaps because the SWNT diameters are not sensitive to the LO/iTO vibrations. While the coherent G band signals are experimentally observed to be on the same order of magnitude as the RBM signals, 9,11 our previous theoretical calculation predicted that the modulations of absorption coefficient due to the G band (LO) phonons are about 1000 times smaller than those caused by the RBM. 12 We expect that the reason for the discrepancy is because we

considered only the change of the energy gap as a main contribution for the coherent G band spectra and also the excitation pulse was too long (50 fs), whereas the G band oscillation period is about 20 fs. This fact indicates that a different mechanism is necessary to explain the coherent G band intensity and that the effects of laser pulse width on the coherent phonon intensity should be taken into account, both of which are the main subjects of this paper.

One possible dominant contribution to the coherent G band intensity is the modulation of electron-photon interaction. For example, Gr¨neis u et al. discussed the optical absorption of graphene from π to π ∗ bands, where the atomic dipole matrix elements for the nearest neighbor carbon-carbon atoms, m opt , are essential. 13 The optical matrix elements are thus sensitive to the change in the carbon-carbon bond length, which can be significantly modified by the G band phonons. In this work, in addition to the changes in electronic structure which arise from the coherent phonons, we now consider changes to the optical matrix element which arise from the coherent phonon oscilllations. We find that modulation of m opt is particularly relevant to the coherent G band intensity and that the changes to the optical matrix element for the G band are larger than for the RBM oscillations. We calculate the coherent G band spectra for a specified SWNT chirality and compare them with the other coherent phonon modes in the SWNT. By a simple analytical model, we also study how the variation of the laser pulse width affects the coherent phonon intensity.

This paper is organized as follows. In Section II, we explain coherent phonon simulation methods which include a general theory for the generation and detection of coherent phonons in SWNTs. In Section III, we present the main results and discuss how the coherent G band intensity could have a stronger signal by considering the modulation of optical interaction and shorter pulse width. And finally, we give a conclusion in Section IV.

## II. SIMULATION METHODS

To calculate coherent phonon spectra, we follow the methods described in our earlier papers, 12,14 except that we will also now treat the effects of the coherent phonon modulations of the electron-photon interaction which we previously neglected for simplicity. We define a coherent phonon mode with wavevector q = 0 (Γ point phonon) whose amplitude satisfies a

driven oscillator equation

$$\begin{array} { c } \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \  \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \,
 \quad \quad \quad$$

where m and ω m denote the phonon mode (e.g. RBM, oTO, LO, iTO) and its frequency, respectively. Equation (1) is solved subject to the initial conditions Q m (0) = 0 and ˙ Q m (0) = 0. The driving function S m ( ) in the right hand side of Eq. (1) is given by t

$$S _ { m } ( t ) = - \frac { 2 \omega _ { m } } { \hbar } \sum _ { n k } \mathcal { M } _ { n } ^ { m } ( k ) \left ( f _ { n } ( k, t ) - f _ { n } ^ { 0 } ( k ) \right ).$$

where f n ( k, t ) is the time-dependent electron distribution function and f 0 n ( k ) is the initial equilibrium electron distribution function. Here n labels an electronic state, while k gives the electron wavevector. The electronic states of a SWNT are calculated within the extended tight-binding (ETB) approximation. 15 The electron-phonon matrix element M m n ( k ) in Eq. 2 is a shorthand for M m, 0 nk nk ; , where M m,q n k ′ ′ ; nk is the deformation potential electron-phonon matrix element in the ETB model with phonon wavevector q = k -k ′ and with a transition from the state n to n ′ . 16 .

From Eq. (2), we see that the driving function S m ( ) depends on the photoexcited electron t distribution functions, which can be calculated generally by taking photogeneration and relaxation effects into account. In coherent phonon spectroscopy, an ultrafast laser pulse generates electron-hole pairs on a time scale short in comparison with the coherent phonon period. The observed coherent phonon intensity is then proportional to the power spectrum of the oscillations of optical properties. 10 Within the scope of this work, we ignore relaxation effects of the photoexcited carriers and consider only the rapidly varying photogeneration term which can be calculated directly from the Fermi's golden rule. Neglecting carrier relaxation has a negligible effect on the computed coherent phonon signal since the relaxation time is much greater than the laser pulse duration and the coherent phonon period. 12 Using the Fermi's golden rule, we obtain the photogeneration rate for the distribution functions, 17

$$\mu & \quad \partial f _ { n } ( k ) = \frac { 8 \pi ^ { 2 } e ^ { 2 } \ u ( t ) } { \hbar { \ } n _ { g } ^ { 2 } \left ( E _ { \text{pump} } \right ) ^ { 2 } } \left ( \frac { \hbar { ^ { 2 } } } { m _ { 0 } } \right ) \sum _ { n ^ { \prime } } | P _ { n n ^ { \prime } } ( k, t ) | ^ { 2 } \\ & \times \left ( f _ { n ^ { \prime } } ( k, t ) - f _ { n } ( k, t ) \right ) \, \delta \left ( E _ { n n ^ { \prime } } ( k, t ) - E _ { \text{pump} } \right ), \\ \nu _ { n } \, \mu \, \quad \, | \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \, \mu \, \nu _ { n } \, \nu _ { n } \, \mu \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _{ n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n } \, \nu _ { n }$$

where E nn ′ ( k, t ) = | E n ( k, t ) -E n ′ ( k, t ) | are the k dependent transition energies at time t of a coherent phonon oscillation, E pump is the pump laser energy, u t ( ) is the time-dependent energy density of the pump pulse, e is the electron charge, m 0 is the free electron mass,

and n g is the refractive index of the surrounding medium. The pump energy density u t ( ) is related with the pump fluence F by a relation F = ( c/n g ) ∫ u t dt ( ) and u t ( ) is also assumed to be a Gaussian. Thus

$$u ( t ) = A _ { p } e ^ { - 4 t ^ { 2 } \ln 2 / 2 \tau _ { p } ^ { 2 } },$$

where A p = (2 n F g √ ln 2 /π / cτ ) ( p ), with c is the speed of light. In Eq. (4), τ p is defined as the pump duration or laser pulse width. Unless otherwise mentioned, we use parameters τ p = 10fs, F = 10 -5 Jcm -2 , and n g = 1. To also account for spectral broadening of the laser pulses, we replace the delta function in Eq. (3) with a Lorentzian lineshape

$$\delta ( E _ { n n ^ { \prime } } - E _ { \text{pump} } ) & \to \frac { \Gamma _ { p } / ( 2 \pi ) } { ( E _ { n n ^ { \prime } } - E _ { \text{pump} } ) ^ { 2 } + ( \Gamma _ { p } / 2 ) ^ { 2 } }, \\ \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \ dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots\dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \Delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \d$$

where Γ p = 0 15 eV is the spectral linewidth (FWHM) of the pump pulse. . 12

By considering light polarized parallel to the tube axis ( z axis) that contributed to the optical absorption, we can write the optical matrix element P nn ′ in Eq. (3) within the dipole approximation as 13

$$P _ { n n ^ { \prime } } ( k ) = \frac { \hbar } { \sqrt { 2 m _ { 0 } } } \sum _ { i, j N } C _ { i } ^ { * } ( n ^ { \prime }, k ) C _ { j } ( n, k ) e ^ { i \phi _ { N } ( k ) } m _ { \text{opt} } ( i, j N ),$$

where C n,k i ( ) and φ N ( k ) respectively denote the expansion coefficient and phase factor from the N th two-atom unit cell of the symmetry-adapted ETB wave functions. 18 The atomic dipole matrix element is given by

$$m _ { \text{opt} } = \int d { \mathbf r } \varphi _ { i 0 } ^ { * } ( { \mathbf r } - { \mathbf R } _ { i 0 } ) \frac { \partial } { \partial z } \varphi _ { j N } ( { \mathbf r } - { \mathbf R } _ { j N } ),$$

where φ i,N is the 2 p z orbital of the i th atom in the N th unit cell.

We should note that Eqs. (6) and (7) still do not have an explicit time dependence. The time-dependence of the optical matrix element comes from the coherent phonon amplitude Q m ( ) which allows the atomic matrix element t m opt to also vary as a function of time as the positions of the carbon atom change. Gr¨ uneis et. al calculated the integral in Eq. (7) for planar graphene analytically by expanding the orbital wavefunctions in terms of Gaussians and it was found that m opt explicitly depends on the bond length between two carbon atoms a CC . 13 If the bond length a CC is altered by coherent phonon oscillations, the atomic matrix element m opt is directly affected, as is the dipole optical matrix element P nn ′ . This is because

the deformation of the bond lengths alters the transfer integral and overlap matrix elements in the ETB model.

Based on above argument, the time-dependence of E nn ′ ( k, t ) and P nn ′ ( k, t ) can be obtained from the time-dependent lattice displacements due to the change in a CC by the coherent phonon oscillations, especially for the G band, which is the in-plane C-C bondstretching mode. From the coherent phonon amplitudes, the time-dependent macroscopic displacements of each carbon atom in an SWNT are given by

$$\mathbf U _ { r, N } ( t ) = \frac { \hbar } { \sqrt { 2 M } } \sum _ { m } \frac { \hat { \mathbf e } _ { r \vec { J } } ^ { m } } { \sqrt { \hbar { \omega } _ { m } } } \ Q _ { m } ( t )$$

where ˆ e m r /vector J ≡ ˆ e m sj ( q = 0), /planckover2pi1 ω m ≡ /planckover2pi1 ω m ( q = 0), and M is the mass of a carbon atom. The bond length a CC at each time t of a coherent phonon oscillation can then be calculated from the macroscopic carbon atom displacements. Therefore, the time-dependent optical matrix element can be evaluated by

$$P _ { n n ^ { \prime } } ( k, t ) = P _ { n n ^ { \prime } } ( k, 0 ) + \Delta P _ { n n ^ { \prime } } ( k, t ),$$

where ∆ P nn ′ ( k, t ) is directly proportional to the time-dependent m opt and we take an average of P nn ′ over three nearest neighbor atoms.

In coherent phonon spectroscopy, a laser probe pulse is used to measure the time-varying absorption coefficient of the SWNT. The time-dependent absorption coefficient α t ( ) at a probe energy E probe is given by the Fermi's golden rule

$$\alpha ( E _ { \text{probe} }, t ) \in & \sum _ { n n ^ { \prime } } \int d k \ | P _ { n n ^ { \prime } } ( k, t ) | ^ { 2 } \Big ( f _ { n } ( k, t ) - f _ { n ^ { \prime } } ( k, t ) \Big ) \\ & \times \delta \Big ( E _ { n n ^ { \prime } } ( k, t ) - E _ { \text{probe} } \Big ), \\. \quad. \quad. \quad. \quad. \quad. \quad \Gamma \quad \partial \quad \partial \quad. \quad. \quad. \quad. \quad. \quad. \quad.$$

We replace the delta function in Eq. (10) with a broadened Lorentzian spectral lineshape with a FWHM of γ = 0 15 . eV, 12 similar to that in Eq. (5). Excitation of coherent phonons by the laser pump modulates the optical properties of the SWNTs, which gives rise to a transient differential transmission signal, or the modulations of absorption coefficient. The time-resolved differential gain measured by the probe is then given by

$$\Delta \alpha ( E _ { \text{probe} }, t ) = - [ \alpha ( \hbar { \omega }, t ) - \alpha ( \hbar { \omega }, t \to - \infty ) ]$$

We take the theoretical coherent phonon signal (or intensity, I ) to be proportional to the

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0672v2-figure-cd7c2390-1184-49ca-8313-9341d7820535.jpg)

FIG. 1: (a) Atomic matrix element as a function of carbon-carbon bond length in SWNTs and (b) the change of bond length as a function of time due to a coherent LO phonon oscillation of in a (6 5) , SWNT at a laser (both pump and probe) energy of 1 15eV. . Inset in (a) shows an enlarged region between 1 40 ˚ . A and 1 45 ˚ . . A A dot in the inset corresponds to the bond length without coherent phonon oscillations, whereas the shaded region corresponds to the area in which the bond length oscillates as shown in (b). (c) Coherent phonon intensity as a function of phonon frequency, showing the RBM and LO peaks. Solid (dashed) line is the calculated result with (without) considering the modulation of optical matrix element.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0672v2-figure-9f496fff-c3d5-45d9-bd1e-3683a55f2bc6.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0672v2-figure-09ca7e08-541c-48a7-b16d-ff34028a8d13.jpg)

Fourier power spectrum of such absorption modulations at a given energy E probe ,

$$I ( \omega ) = \int e ^ { - i \omega t } \, | \Delta \alpha ( E _ { \text{probe} }, t ) | ^ { 2 } \, d t,$$

where ω represents the phonon frequency that contributes to the coherent phonon spectra.

## III. RESULTS AND DISCUSSION

## A. Modulation of optical interaction

First we discuss the effects of coherent phonon oscillations on the optical interaction. The changes in a CC modulate the atomic matrix element m opt because of the direct correspondence between these two quantities at time t . Fig. 1(a) shows the calculated m opt as a function of a CC based on the formula given by Gr¨neis u et al. 13,19 . It indicates that the strength of optical interaction monotonically decreases as a function of a CC . In the inset of Fig. 1(a), we show the atomic matrix element within an enlarged region around 1 40 ˚ . A and 1 45 ˚ . . A The shaded region corresponds to the possible values of a CC affected by the coherent LO phonon oscillation given in Fig. 1(b). From this figure, we can say that the

modulations of optical interaction is about 0 02 [a.u] . -1 for the change of vibration amplitude of about 0 02 ˚ . . A These modulations of optical interaction is thus approximately 10 % of m opt = 0 25[a.u.] . -1 , which is not negligible for calculating the absorption coefficient of a SWNT. The coherent phonon intensity is proportional to | ∆ α | 2 ∝ | ∆ m opt | 4 , which is the leading order of the spectra. In the previous study, however, this fact was not taken into account and the optical matrix element was considered constant as a function of time. 12

Next, from the time-dependent optical matrix elements, we proceed to the calculation of coherent phonon spectra by taking the Fourier transform of Eq. (11). The calculation is performed by allowing the probe energy in Eq. (10) to be varied independently while keeping the pump energy in Eq. (3) constant. We take a (6 5) SWNT chirality as a sample , for this calculation. This SWNT has the first and second optical transition energies (band gaps) of 1 27 eV and 2 42 eV, denoted by . . E 11 and E 22 , respectively. 20 In this calculation we neglect the exciton effects for simplicity. Basically, the exciton-photon matrix elements are about 100 times larger than the electron-photon matrix elements, 21 but such enhancement factors are common for all the phonon modes. Therefore, the exciton effects will not modify the relative intensity between the phonon modes. In Fig. 1(c), we show an example of the calculation for intensity as a function of phonon frequency by including or excluding the modulation of optical interaction. The coherent phonon spectra shows both the RBM and LO peaks for a particular laser excitation energy of 1 15 eV. . It can be seen that the LO intensity is enhanced significantly when taking the modulation of optical interaction into account, while the RBM intensity is just enhanced slightly.

To further understand the laser energy dependence of the spectra, we calculate the coherent phonon intensity for a given phonon frequency ω m by considering different laser probe energy from 1 0-3 0 eV with an interval of 0 1 eV. . . . In Fig. 2, we show absorption coefficients and coherent phonon spectra of the (6 5) SWNT as a function of probe energy. , For the coherent phonon spectra, we give the spectra both in the linear scale and logarithmic scale as shown in Figs. 2(b) and (c), respectively. The spectra are accompanied with the plot of absorption coefficient in Fig. 2(a) as a reference for showing the positions of the optical transition energy peaks. In Figs. 2(b) and (c), we compare the coherent G band phonon spectra (LO and iTO modes) with RBMs and also with oTO (out-of-plane TO) mode for the (6 , 5) tube. In Fig. 2(b), we can see that the coherent RBM intensity and LO intensity are on the same order, with the RBM intensity being slighly larger than that of the LO intensity

FIG. 2: (Color online) From top to bottom shows (a) absorption coefficient, (b) linearly scaled and (c) logarithmically scaled coherent phonon intensities as a function of probe energy for a (6,5) SWNT. The intensity is normalized to the maximum intensity of the RBM.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0672v2-figure-91f14791-caf0-4b97-931b-f942fb26c627.jpg)

by a ratio of about 2 5 and 2 1 at . . E 11 and E 22 , respectively. These results indicate that modulations of optical matrix elements become important in enhancing the coherent G band intensity. It should be noted that the coherent iTO intensity is hundred times smaller than the LO intensity. Therefore, the coherent G band phonon spectra are mainly dominated by the LO phonon modes.

It is also interesting to see in Fig. 2(c) that there is a dip at 2 eV for the RBM phonons, which might be related with the zero value of the coherent RBM phonon amplitude. 14

The dip of RBM coherent phonon spectra could give information of photon energy that would correspond to the transition from expansion to contraction (or vice versa) of the SWNT diameter. 14 Moreover, we obtain two peaks at each transition energy for all phonon modes, which are consistent with some earlier works that reported the excitation excitation energy dependence of coherent phonon intensity always shows a derivative-like behavior of the absorption coefficient. 4,11 The double-peak feature at each transition energy can be symmetric or asymmetric depending on whether or not the excitonic effects is taken into account. 22 In this work, the double-peak lineshapes are asymmetric because we neglect the exciton effects for simplicty. It is then worth comparing the ratio of the coherent RBM and LO intensity obtained in this study with that in the experiment. For example, a pump-probe measurement by Lim et al. gave the RBM intensity of about eight times larger than the LO intensity. 9 Although this discrepancy is not significant in our present discussion, we expect that it might come from the additional effect of the selection of laser pulse width τ p as we will discuss below.

## B. Effects of laser pulse width

To discuss the effects of laser pulse width ( τ p ) on the coherent phonon intensity, we can analytically model the driving function S m ( ) of Eq. (2) by using the laser pulse in the form t of Eq. (4) and then solve for Q m ( ). t By understanding the τ p dependence of Q m , we can qualitatively explain the trend of the coherent phonon intensity when τ p is varied. As we can see from Eq. (2), S m ( ) t is proportional to the carrier density, f n ( k, t ), which can be obtained by integrating Eq. (3) with respect to time. For simplicity, we can write S m ( ) to t be directly proportional to the integration of u t ( ),

$$S _ { m } ( t ) & \in \int _ { - \infty } ^ { t } A _ { p } e ^ { - 4 t ^ { \prime 2 } \ln 2 / 2 r _ { p } ^ { 2 } } d t ^ { \prime } \\ & \in \frac { n _ { g } F } { 2 c } \sqrt { \frac { \pi } { \ln 2 } } \left [ 1 + \text{erf} \left ( \frac { 2 t \ln 2 } { \tau _ { p } } \right ) \right ],$$

where erf( x ) = (2 / √ π ) ∫ x 0 e -x ′ dx ′ is the error function. To obtain the full equality between the left-hand and right-hand sides of Eq. (13), we can put an additional term of electronphonon matrix element as also indicated in Eq. (2). This additional term along with the prefactor in the right-hand side of Eq. (13) form a constant A m , which will change only

FIG. 3: (Color online) (a) Coherent phonon driving for the RBM and LO phonon modes of the (6 5) SWNT under , E 11 excitation, each calculated with three different values of τ p : 10fs, 20 fs, and 30fs. (b) Coherent phonon amplitude for the RBM and LO phonon modes obtained analytically as a function of pulse width. Inset shows the logarithmic plot of the same amplitude as a function of squared pulse width.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0672v2-figure-f859b19b-7fce-4ffc-a8fd-240115d20fc3.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0672v2-figure-1292be9a-2d6e-45fd-88c6-ada15678dc28.jpg)

when we have different phonon modes m . We can finally write the driving function as

$$S _ { m } ( t ) = \frac { A _ { m } } { 2 } \left [ 1 + \text{erf} \left ( \frac { 2 t \ln 2 } { \tau _ { p } } \right ) \right ],$$

and the corresponding solution for Q m ( ) with an initial condition of t Q m (0) = 0 and ˙ Q m (0) = 0 is

$$Q _ { m } ( t ) & = \frac { A _ { m } } { \omega _ { m } ^ { 2 } } \left [ 1 + e ^ { - \omega _ { m } ^ { 2 } \tau _ { p } ^ { 2 } / 1 6 \ln 2 } \cos ( \omega _ { m } t ) \right ]. & & ( 1 5 ) \\ \dots \cdots \cdots. & c \cap ( \omega _ { m } ) \cdots \cdots. & \ddots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots.$$

Having the solution of Q m ( ), we can now discuss its dependence on t τ p . First, the value of A m in Eqs. (14) and (15) can be obtained by fitting to the maximum value of the force S m ( ) simulated from the full microscopic treatment in Eq. (2). t In Fig. 3(a), we show the

simulated S m ( ) for the RBM and LO mode of (6 5) SWNT under t , E 11 excitation with three different values of τ p . We see that S m ( ) for all cases show a step-like behavior with a width t of τ p , and thus consistent with Eq. (14). The maximum values of S m ( ) only differ between t different phonon modes. The fitted values of A m for the RBM and LO mode are 1161 7 ps . -2 and 3516 5 ps . -2 , respectively. In Fig. 3(b), we show the magnitude of Q m , denoted as Q mag , as a function of pulse width τ p for the RBM and LO phonon modes. The definition for Q mag is

$$Q _ { \max } = \frac { A _ { m } } { \omega _ { m } ^ { 2 } } e ^ { - \omega _ { m } ^ { 2 } \tau _ { p } ^ { 2 } / 1 6 \ln 2 },$$

which represents the difference between the maximum and minimum values of the coherent phonon oscillation amplitudes. We also have 2 π/ω RBM = 110 fs and 2 π/ω LO = 21 fs for the RBM and LO oscillation periods of the (6 , 5) SWNT, respectively. We can see from Fig. 3(b) that as the pulse width increases, the coherent phonon amplitude quickly decays following the Gaussian shape of the spectrum of the laser pulse. 23 However, the rate of the amplitude decay depends on the phonon mode oscillation frequency or period, as clearly shown in the inset of Fig. 3(b). If the pulse width is much smaller than the phonon oscillation period, the amplitude will be enhanced. In this case, the LO phonon mode is enhanced more quickly than the RBM mode after the pulse width becomes shorter than the LO oscillation period. Therefore, as we have used τ p = 10 fs in the simulation discussed earlier, the coherent LO intensity rapidly increase while at the same time the coherent RBM intensity increases more slowly. This could be the reason why we have a slightly different ratio of the RBM to the LO intensity since the coherent LO phonon amplitude is much more sensitive to the variation of the laser pulse width within sub-10 fs region compared to the coherent RBM phonon amplitude.

## IV. CONCLUSION

We have presented the mechanism by which a strong coherent G band signal could be generated in ultrafast pump-probe spectroscopy. Instead of the energy gap modulation mechanism which is dominant in the RBM case, we suggest that the modulations of electronphoton interaction as a function of time should be relevant to the coherent G band intensity. We also find an analytical formula that describes how a typical coherent phonon amplitude behaves as a function of laser pulse width. The formula indicates that the G band (LO

mode) intensity increases more rapidly than the RBM intensity, especially when the laser pulse width is much shorter than each of the phonon mode period.

## Acknowledgments

A.R.T.N acknowledges financial support from JSPS Fellowship for Young Scientists. R.S. acknowledges MEXT Grant No. 25286005. The University of Florida authors acknowledge NSF-DMR Grant No. 1105437 and OISE-0968405. We are all grateful to Prof. J. Kono (Rice University) and his co-workers for fruitful discussions which stimulated this work.

- 1 R. Saito, G. Dresselhaus, and M. S. Dresselhaus, Physical Properties of Carbon Nanotubes (Imperial College Press, London, 1998).
- 2 A. Jorio, M. S. Dresselhaus, and G. Dresselhaus, Carbon Nanotubes: Advanced Topics in Synthesis, Structure, Properties, and Applications (Springer-Verlag, Berlin, 2008).
- 3 A. Gambetta, C. Manzoni, E. Menna, M. Meneghetti, G. Cerullo, G. Lanzani, S. Tretiak, A. Piryatinski, A. Saxena, R. L. Martin, and A. R. Bishop, Nat. Phys. 2 , 515-520 (2006).
- 4 Y. S Lim, K. J. Yee, J. H. Kim, E. H. Haroz, J. Shaver, J. Kono, S. K. Doorn, R. H. Hauge, and R. E. Smalley, Nano Lett. 6 , 2696-2700 (2006).
- 5 K. Kato, K. Ishioka, M. Kitajima, J. Tang, R. Saito, and H. Petek, Nano Lett. 8 , 3102-3108 (2008).
- 6 J.-H. Kim, K.-J. Han, N.-J. Kim, K.-J. Yee, Y.-S. Lim, G. D. Sanders, C. J. Stanton, L. G. Booshehri, E. H. H´roz, and J. Kono, Phys. Rev. Lett. a 102 , 037402 (2009).
- 7 K. Makino, A. Hirano, K. Shiraki, Y. Maeda, and M. Hase, Phys. Rev. B 80 , 245428 (2009).
- 8 J.-H. Kim, A. R. T. Nugraha, L. G. Booshehri, E. H. Haroz, K. Sato, G. D. Sanders, K.-J. Yee, Y.-S. Lim, C. J. Stanton, R. Saito, and J. Kono, Chem. Phys. 413 , 55-80 (2013).
- 9 Y.-S. Lim, A. R. T. Nugraha, S.-J. Cho, M.-Y. Noh, E.-J. Yoon, H. Liu, J.-H. Kim, H. Telg, E. H. H@roz, G. D. Sanders, S.-H. Baik, H. Kataura, S. K. Doorn, C. J. Stanton, R. Saito, J. Kono, and T. Joo, Nano Lett. 14 , 1426-1432 (2014).
- 10 H. J. Zeiger, J. Vidal, T. K. Cheng, E. P. Ippen, G. Dresselhaus, and M. S. Dresselhaus, Phys. Rev. B 45 , 768-778 (1992).

- 11 J.-H. Kim, K.-J. Yee, Y.-S. Lim, L. G. Booshehri, E. H. H´roz, and J. Kono, Phys. Rev. B a 86 , 161415 (2012).
- 12 G. D. Sanders, C. J. Stanton, J.-H. Kim, K.-J. Yee, Y.-S. Lim, E. H. H´roz, L. G. Booshehri, a J. Kono, and R. Saito, Phys. Rev. B 79 , 205434 (2009).
- 13 A. Gr¨neis, R. Saito, Ge. G. Samsonidze, T. Kimura, M. A. Pimenta, A. Jorio, A. G. Souza u Filho, G. Dresselhaus, and M. S. Dresselhaus, Phys. Rev. B 67 , 165402 (2003).
- 14 A. R. T. Nugraha, G. D. Sanders, K. Sato, C. J. Stanton, M. S. Dresselhaus, and R. Saito, Phys. Rev. B 84 , 174302 (2011).
- 15 Ge. G. Samsonidze, R. Saito, N. Kobayashi, A. Gr¨neis, J. Jiang, A. Jorio, S. G. Chou, G. Dresu selhaus, and M. S. Dresselhaus, Appl. Phys. Lett. 85 , 5703 (2004).
- 16 J. Jiang, R. Saito, Ge. G. Samsonidze, S. G. Chou, A. Jorio, G. Dresselhaus, and M. S. Dresselhaus, Phys. Rev. B 72 , 235408 (2005).
- 17 S. L. Chuang, Physics of optoelectronic devices (Wiley, New York, 1995).
- 18 V. N. Popov and L. Henrard, Phys. Rev. B 70 , 115407 (2004).
- 19 A. Gr¨neis. u Resonance Raman spectroscopy of single wall carbon nanotubes . PhD thesis, Tohoku University, Sendai, Japan, September 2004.
- 20 A. R. T. Nugraha, R. Saito, K. Sato, P. T. Araujo, A. Jorio, and M. S. Dresselhaus, Appl. Phys. Lett. 97 , 091905 (2010).
- 21 J. Jiang, R. Saito, K. Sato, J. S. Park, Ge. G. Samsonidze, A. Jorio, G. Dresselhaus, and M. S. Dresselhaus, Phys. Rev. B 75 , 035405 (2007).
- 22 A. R. T. Nugraha, E. I. Rosenthal, E. H. Hasdeo, G. D. Sanders, C. J. Stanton, M. S. Dresselhaus, and R. Saito, Phys. Rev. B 88 , 075440 (2013).
- 23 It is apparent that the coefficient in front of the cosine term in Eq. (15) is related to the Fourier transform of u t ( ) given in Eq. (4). For a pulse shape in terms of Gaussian, the Fourier transform is also a Gaussian.
A. R. T. Nugraha E. H. Hasdeo G. D. Sanders C. J. Stanton R. Saito
Origin of coherent G band phonon spectra in single wall carbon nanotubes
Coherent phonons in single wall carbon nanotubes (SWNTs) are observed as
oscillations of the differential absorption coefficient as a function of time
by means of pump-probe spectroscopy. For the radial breathing mode (RBM) of a
SWNT, the coherent phonon signal is understood to be a result of the modulated
diameter-dependent energy gaps due to the coherent RBM phonon oscillations.
However, this mechanism might not be the dominant contribution to other phonon
modes in the SWNT. In particular, for the G band phonons, which correspond to
bond-stretching motions, we find that the modulation of the interatomic dipole
matrix element gives rise to a strong coherent G band phonon intensity
comparable to the coherent RBM phonon intensity. We also further discuss the
dependence of coherent G band and RBM phonon amplitudes on the laser excitation
pulse width.