ISRN Condensed Matter PhysicsVolume 2012 (2012), Article ID 382939, 9 pagesdoi:10.5402/2012/382939

Research Article

Phonon Conductance of Potassium- and Sodium-Doped Transpolyacetylene Chain

M. Belhadi1,2 and S. Kheffache¹

¹Laboratoire de Physique et Chimie des Matériaux (LPCM), Université Mouloud Mammeri, BP17 RP, 15000 Tizi-Ouzou, Algeria²Faculté des Sciences et de la Technologie, Université d’Adrar, Rue du 11 Décembre 1960, 01000 Adrar, Algeria

Received 4 August 2012; Accepted 30 September 2012

Academic Editors: M. Naito and A. Oyamada

Copyright © 2012 M. Belhadi and S. Kheffache. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

A theoretical approach for the study of phonon dynamics and scattering properties of doped transpolyacetylene chain is presented. The coherent reflection and transmission scattering cross-sections for phonons incident on the doped unit cell boundary are calculated in accordance with the Landauer-Buttiker electron scattering description, using the matching procedure with the nearest and next nearest neighbor elastic force constants. This is done for two different dopants, namely, the potassium and sodium atoms. Our numerical results yield an understanding of the transpolyacetylene chain dynamical properties and the effects on phonon conductance due to phonon incident on the doped unit cell boundary. The coherent reflection and transmission coefficients show characteristic spectral features, depending on the cutoff frequencies for the propagating phonons and on the nature of the dopants. They illustrate the occurrence of Fano resonances in the scattering spectra that result from the interactions of propagating elastic waves of the undoped transpolyacetylene chain with the localized modes due to the breakdown of the translation symmetry in the x direction.

1. Introduction

Conjugated polymers have been studied intensively because of their fascinating properties. These materials are of great fundamental interest in that they are exotic, low-dimensional semiconductors [1–3]. More recently, the use of such conjugated polymers has been explored in a variety of technological applications such as field effect transistors [4, 5], light emitting diodes [6, 7], and solar cells [8–10]. They are especially interesting in their doped state, when exhibiting the conductivity ranging from the conductivity typical for insulators to the conductivity typical for metals [11, 12]. For instance, as pointed out by Heeger [13], the nature of the elementary excitations, mobile charge carriers, or bound neutral excitants remains controversial.

Among these conjugated polymers, transpolyacetylene a quasi-one dimensional conductor, has been the subject of many investigations [14, 15]. Most of the theoretical works on elementary excitations in these materials have focused on a single-polyacetylene chain. In its doped state, polyacetylene chain forms a new class of conducting polymers, in which the electrical conductivity can be varied over a range of eleven orders of magnitude [16, 17]. It is well known that the existence of dopants in conjugated polymers can greatly modify the transport properties of the system. Previous studies have reported the doping dependence of the conductivity of conjugated polymers [18–20], and showed that the effects are crucial for the improvements of device performance. Understanding phonon scattering properties due to dopants in such systems will be essential to extend elementary excitations physics into the mesoscopic regime, where a variety of size and geometry effects are important for scattering and transport properties. Phonon conductance plays a fundamental role in controlling the performance of molecular devices.

Transport properties of low-dimensional systems have been subjects of interest for a long time. Most of the recent research in this area has been oriented toward the study of electronic scattering. The understanding of coherent electronic transport in the mesoscopic regime has been provided by the formalism of Landauer [21] and Büttiker [22] for electron transport in mesoscopic systems, who related the conductance of the system to its scattering matrix. Multiple scattering and quantum interference phenomena are not limited, however, to the electronic transport, but in a variety of problems they can become important for the coherent transport of other types of excitations such as phonons [23, 24] and magnons [25, 26] in low-dimensional systems.

The purpose of this work is to perform theoretical investigations on phonon dynamics and phonon conductance of potassium- and sodium-doped transpolyacetylene chain, by the use of the matching procedure, [27, 28] with nearest and next nearest neighbours elastic force constants. The scattering properties are calculated in full accordance with Landauer-Buttiker formalism for electron transport in mesoscopic systems. The defect is treated as the doped unit cell by substituting two hydrogen atoms by two potassium or two sodium atoms.

The paper is organized as follows. In Section 2, we present a brief description of the model used and the theoretical approach employed for studying the dynamic properties of undoped transpolyacetylene chain. In Section 3, we present the mathematical framework for determining the scattering properties due to given incorporated atoms or dopants in the unit cell of the undoped transpolyacetylene chain. Some numerical results and discussions are given in Section 4 for each case along with the conclusion of the paper.

2. Model and Its Dynamical Properties

The prototype of the conjugated polymers is polyacetylene . In its simplest form, a single-polymer chain consists of a zigzag backbone of carbon atoms. One hydrogen atom is bonded to each carbon atom, so that every carbon atom is threefold coordinated. In Figure 1(a), we show, within the box defined by the dotted lines, the unit cell for transpolyacetylene chain. This unit cell of an infinite one-dimensional chain contains two inequivalent carbon atoms A and B. The simplest system model used in this study is depicted in Figure 1(b).

(a)

(b)

Figure 1: (a) A transpolyacetylene chain structure, showing the unit cell bounded by a box defined by the dotted lines with two inequivalent Carbone atoms A and B in the unit cell. (b) A schematic representation of doped transpolyacetylene chain. The doped unit cell may contain different masses and force constants. The elastic constants in the boundary may also be different from the rest of the undoped transpolyacetylene chain. The doped unit cell, the incident waves (IW), the reflected waves (RW) and the transmitted waves (TW) are shown. The two perfect wave guides (LW) and (RW) and the two matching domains (LMD) and (RMD) are also indicated.

The atoms in the transpolyacetylene chain are held together by the chemical bonds between the atoms. These bonds are not rigid, but act like rings, which connect the atoms, creating a spring mass system with defined stiffness constants. Consider , , and and , , and to denote the elastic force constants between nearest and next nearest neighbor sites as shown in Figure 1(b). We define for convenience the ratio: ,  , and . The approximation that we make is to neglect the motion of side groups relative to the backbone by compressing the atoms in the groups into masses or . Each mass or corresponding to the mass of a (C–H) or (C–K or C–Na) group, respectively, is linked to its neighbors by harmonic springs.

When an atom is displaced in the chain, this displacement can travel as a wave through the medium. These waves can either be longitudinal, where the displacement of the atom is in the same direction as the propagation of the wave, or can be transverse, where the atomic displacement is perpendicular to the direction of propagation. The equation of motion for an atom located at site in the two semiinfinite transpolyacetylene chains is given by: where the indices denote Cartesian co-ordinates. Identical mass is considered for all the sites in the semiinfinite undoped polyacetylene chain. is the mass of the (CH) group and is the corresponding displacement vector of the atom along the Cartesian direction. The radius vector between and , has Cartesian components , and . Equation (1) is applied systematically for the transpolyacetylene chains, at the doped unit-cell boundary and in the domains to the left (LW) and the right (RW) of the doped unit cell. The same interatomic distance is taken between nearest neighbour sites throughout the system.

The solutions to (1) for the two sites in the unit cell of the semiinfinite undoped transpolyacetylene chain far from the doped unit cell of the structure illustrated in Figure 1(b) may hence be expressed as where is a dimensionless frequency given by . The frequency characteristic of the system is given by , for and nearest neighbors. is the dynamic matrix for the configuration described in Figure 1(a). is the corresponding unit matrix. The complex factors determine the phase changes of the solution of (2) along the undoped polyacetylene chain. We can distinguish between two different types of solutions: those corresponding to , which describes propagating modes, and to , which describes evanescent modes. The diagonalisation of the matrix for leads to the propagating eigenmodes and the corresponding eigenvectors for the undoped chain. The resulting phonon dispersion curves are usually given as function of the wave vector , where is the interatomic distance and the one-dimensional reciprocal lattice wave vectors along the direction in the first Brillouin zone.

The study of the phonon scattering in the presence of the doped unit cell requires a detailed analysis of both the propagating modes and the evanescent solutions of the system. In the undoped transpolyacetylene chain, these modes are described by the phase factor doublets , which are complex conjugate pairs that verify symmetrically the characteristic secular equations in each half space. The complex phase factors are the solution of (2). The determinant of the dynamical matrix vanishes leading to a characteristic secular equation, that is solved for a given system parameters at frequency . For the transpolyacetylene chain in Figure 1(b), only physical acceptable solutions are retained, namely, the propagating and the evanescent modes. The evanescent modes do not transport any energy but are necessary for a complete description of the scattering problem, and to match the solutions on either side of the defect. The evolution of the absolute values of these phase factors may be given as a function of the normalized frequency and the elastic force constant parameters of the system and . A more interesting alternative is to present the dispersion curves via 3D representation of as a function of the Real and Imaginary parts of the phase factors .

3. Scattering at the Doped Unit Cell

A number of methods exist to study the vibration dynamics and scattering properties of perturbed crystals. The matching procedure that we employ in this paper allows one to deal with both aspects of localized modes and of wave scattering, within the same mathematical framework. The reflection and transmission probabilities as well as the conductance of the doped unit cell in transpolyacetylene chain are calculated as a function of the phonon scattering frequency in full accordance with the Landauer-Buttiker formalism for electron transport.

In the absence of inelastic scattering, for an incoming propagating phonon eigenmodes , incident from the left to right of the doped unit cell along the direction, at a frequency , the resulting scattered fields are composed of the reflected and the transmitted parts. For sites inside the left semiinfinite transpolyacetylene chain (LW), the Cartesian components of the displacement field may be expressed as the sum of the incident phonon field and a superposition of the eigenmodes of the undoped polyacetylene chain, that are backward scattered in reflection at the same frequency denotes the reflection amplitudes for the backward scattering of incident phonon into the eigenmodes . The quantities denote the eigenvectors of the dynamic matrix at a frequency , which are associated with the atomic displacements for the mode . For sites inside the right semi-infinite transpolyacetylene chain (RW), the displacement field may be expressed by another appropriate superposition of the eigenmodes that are forward scattered at the same frequency denotes the transmission amplitudes for the forward scattering of the incident phonon into the eigenmodes . We denote by the basis vectors for the reflection and transmission in a constructed Hilbert space, and by the vector composed of the displacement for sites in the doped unit cell. The equations of motion for atoms of the doped unit-cell boundary, coupled to the rest of the polyacetylene chains via the matching domains (LMD) and (RMD), may be written as where , and is a square-matched scattering matrix, and the column vector , mapped appropriately onto the basis vectors, regroups the inhomogeneous terms describing the incident phonon. The solution of (5) yields the displacements of the irreducible set of atomic sites for the doped unit cell as well as the coherent reflection and transmission amplitudes and at the scattering frequency .

The scattering is usually described in terms of the scattering matrix. Its elements are given by the reflection and transmission amplitudes at the scattering frequency . When the scattering domain is connected to external reservoirs by phonon wave-guides with open channels, as is the case in the presence of doped transpolyacetylene chain, the reflection and transmission probabilities, and , are related to the elements of the scattering matrix, by

In order to satisfy the unitary requirements of the scattering matrix, the scattered waves have to be normalized with respect to their group velocity, . It is taken equal to zero for evanescent eigenmodes. Summing over all outgoing channels, we obtain the reflection and transmission probabilities for a given eigenmodes at the frequency

Further, in order to describe the overall phonon transmission or reflection of the doped transpolyacetylene chain, it is possible to sum over all input and output channels. This allows one to obtain the system phonon reflectance and the system phonon conductance at any scattering frequency , as

The transmission probabilities per eigenmodes and the phonon conductance are important quantities to calculate since each corresponds to an experimentally measurable observable. They play an essential role in the calculation, for example, of the net heat current across the defect between the two domains held at different temperatures.

4. Numerical Resultants and Discussion

Numerical calculations are carried out using the mathematical framework of the matching procedure in the harmonic approximation with a fixed ratio and for the undoped transpolyacetylene chain. Those model parameters used in this work are those generally chosen for transpolyacetylene [29].

By solving (2), one can determine the normalized frequency of the waves as a function of their wave number (also called wave vector) , where . The relationship between and may then be plotted in the form of a dispersion relation. One important parameter that can be found from the dispersion relation is the group velocity of the phonons, which is the slope of the branch on the dispersion relation.

The resulting phonon dispersion curves far from the doped unit cell in each half-space are given in Figure 2(a) over the first Brillouin zone, for . In this case, the eigenmodes labelled are propagating modes in the following respective frequency intervals: Note that the four modes are indexed from the bottom to the top in Figure 2(a). It should be emphasized that in the frequency intervals some of the modes are propagating phonons, the others are evanescent modes. Figure 2(a) shows two distinct phonon modes in the dispersion curves. The two lower modes are called the acoustic modes because the linear relationship for low frequency is similar to that for a sound wave. These modes correspond to two neighboring atoms moving in phase with one another. The two upper modes, called the optic modes, correspond to the case where two neighboring atoms are moving out of phase with one another; for low frequencies, these modes are characterized by a vanishing group velocity. Because the group velocity of the optic modes is small or equal to zero at low frequencies, these modes generally do not participate in the heat transport process.

(a)

(b)

Figure 2: (a) A typical Phonon dispersion curve for the considered undoped transpolyacetylene chain, with the ratio of and for the system. They are shown here over the interval of the first Brillouin zone. (b) A three-dimensional view of the vibration eigenmodes for the considered transpolyacetylene chain with nearest and next nearest neighbor interactions and for the specific case of and ; with the phase factors, is taken here with the values 1, 2, 3, and 4.

The behavior of the eigenmodes can be illustrated graphically by the curves on and inside the unit circle of the complex plane. The resulting dispersion curves in the complex plane are represented in Figure 2(b), with nearest and next nearest neighbor interactions and for the specific case of and for the system, showing typically the generalized functional behavior for the four eigenmodes. The propagating modes, namely, , are represented by the curves following the unit circle. Usually the latter is represented by the phonon dispersion curves as a function of the wave vector , rather than as a function of the parameter . The evanescent modes inside the unit circle are described for . The absolute values of the phase factors as a function of the normalized frequency are given specifically in Figure 3(a). For completeness, we also present in Figure 3(b) the group velocities of the phonon eigenmodes over their frequency intervals of propagation, under the same conditions.

(a)

(b)

Figure 3: (a) Evolution of the absolute values of the complex phase factors , as a function of the dimensionless frequency , for the above and values. (b) Group velocity for the different phonon eigenmodes as a function of the dimensionless frequency , for the specific case of and of the undoped transpolyacetylene chain.

The phonon scattering properties of the doped transpolyacetylene chain are calculated using (7) and (8) with reference to wave incident from the left to the right of the doped unit cell, in Figure 1(b). Numerical results are obtained for a number of choices concerning the system force constants in the scattering domain. We have investigated the evolution of these results and found that it is interesting to present those for scattering from a set of system parameters , , and in two cases of interest as regards the mass of atoms constituting the doped unit cell, namely, the potassium and sodium atoms.

In Figures 4(a)–4(d) and 5(a)–5(d), we present the phonon scattering frequency dependency of the calculated transmission , , , and and reflection , , , probabilities, owing to the scattering of phonon modes 1, 2, 3 and 4 respectively, in the cases of potassium- and sodium-doped unit cell in polyacetylene chain.

(a)

(b)

(c)

(d)

(e)

Figure 4: (a) The transmission and the reflection coefficients for phonon incident mode 1, (b) the transmission and the reflection coefficients for phonon incident mode 2, (c) the transmission , and the reflection coefficients for phonon incident mode 3, (d) the transmission and the reflection coefficients for phonon incident mode 4, and (e) the total transmission probability or phonon conductance of the doped transpolyacetylene chain and the total transmission probability or phonon conductance of the undoped transpolyacetylene chain as a function of the scattering frequency in the case of potassium-doped transpolyacetylene chain with .

(a)

(b)

(c)

(d)

(e)

Figure 5: Similar to Figure 4 for the case of sodium-doped transpolyacetylene chain with .

The results show that the transmission probabilities are a sensitive function of the scattering frequency . As expected, the influence of the dopants is smallest at low frequency. At the Brillouin zone center, where tends to zero, only acoustical modes 1 and 2 are excited. In this domain, the transmission probabilities approach the unity for these incident modes, which indicates that the doped unit cell does not influence the transmission for wavelengths much greater than the lattice constant. The presence of dopants in the transpolyacetylene chain leads, however, to a general decrease of the transmission probabilities with increasing the scattering energy and becomes zero at the BZ boundaries. For high frequencies, only the optical modes 3 and 4 are excited. We note that the frequency interval correspond to the propagating phonon modes 1, 2, 3, and 4, indicated in Figure 2(a). The upper and lower scattering frequency limits of the and spectra can be identified with reference to the dispersion curves of the undoped transpolyacetylene chain. The most striking consequence of these dopants is the appearance of resonance structures in the transmission spectra as a function of the scattering frequency, which can be understood as Fano resonances, usually invoked to describe the coherent coupling between the propagating modes of the system and the localized vibration modes induced by the potassium atoms. The transmission and the reflection cross-sections verify the unitarity condition for the scattering matrix, and this is used throughout as a check on the numerical calculation.

In Figures 4(e) and 5(e), the total transmission probability or phonon conductance of potassium- and sodium-doped transpolyacetylene chain, respectively, and the total transmission probability or phonon conductance of undoped transpolyacetylene chain, respectively, are presented as a function of the scattering frequency . The phonon conductance of the system is useful to calculate, as it corresponds to an experimentally measurable quantity, for example, for thermal transport across such doped polymer. The histogram step of in Figures 4(e) and 5(e) corresponds to the change in the number of propagating modes in the different frequency intervals.

In summary, we have studied theoretically the phonon dynamics and scattering properties in doped conjugated polymers such as potassium- and sodium-doped transpolyacetylene chain. The dispersion curves and group velocities as well as the phase factors of the propagating modes in undoped transpolyacetylene chain are calculated. Numerical results are obtained and graphically presented as function of the scattering frequency. We have also analyzed the dependence of the phonon transmission and reflection probabilities as well as the phonon conductance on the scattering frequency. This is done using the mathematical framework of the matching procedure in the harmonic approximation with nearest and next nearest neighbor elastic force constants. Our numerical results show that in spite of their special character, the scattering of phonons has some features in common with the scattering of electron waves, which can be described in terms of basically the same interference phenomena. The results show that the transmission probabilities and the phonon conductance, calculated in full accordance with the Landauer-Buttiker description of electron transport in mesoscopic systems are a strong function of the scattering frequency. The most striking consequence of this defect is the appearance of resonance structures in the transmission spectra as a function of the scattering frequency and system parameters. They are understood as Fano resonances, which are usually invoked to describe the coherent coupling between the propagating modes of the system and the localized vibration modes in the neighborhood of the dopants.

References

A. G. MacDiarmid, “Nobel lecture: “synthetic metals”: a novel role for organic polymers,” Reviews of Modern Physics, vol. 73, no. 3, pp. 701–712, 2001.

W. P. Su, J. R. Schrieffer, and A. J. Heeger, “Soliton excitations in polyacetylene,” Physical Review B, vol. 22, no. 4, pp. 2099–2111, 1980.

H. Takayama, Y. R. Lin-Liu, and K. Maki, “Continuum model for solitons in polyacetylene,” Physical Review B, vol. 21, no. 6, pp. 2388–2393, 1980.

M. J. Małachowski and J. Zmija, “Organic field-effect transistors,” Opto-Electronics Review, vol. 18, no. 2, pp. 121–136, 2010.

Y. Sun, Y. Liu, and D. Zhu, “Advances in organic field-effect transistors,” Journal of Materials Chemistry, vol. 15, no. 1, pp. 53–65, 2005.

J. Kovac, L. Peternai, and O. Lengyel, “Advanced light emitting diodes structures for optoelectronic applications,” Thin Solid Films, vol. 433, no. 1-2, pp. 22–26, 2003.

J. H. Burroughes, D. D. C. Bradley, A. R. Brown et al., “Light-emitting diodes based on conjugated polymers,” Nature, vol. 347, no. 6293, pp. 539–541, 1990.

M. H. Chen, J. Hou, Z. Hong et al., “Efficient polymer solar cells with thin active layers based on alternating polyfluorene copolymer/fullerene bulk heterojunctions,” Advanced Materials, vol. 21, no. 42, pp. 4238–4242, 2009.

P. W. M. Blom, V. D. Mihailetchi, L. J. A. Koster, and D. E. Markov, “Device physics of polymer: fullerene bulk heterojunction solar cells,” Advanced Materials, vol. 19, no. 12, pp. 1551–1566, 2007.

S. Günes, H. Neugebauer, and N. S. Sariciftci, “Conjugated polymer-based organic solar cells,” Chemical Reviews, vol. 107, no. 4, pp. 1324–1338, 2007.

S. Roth and W. Graupner, “Conductive polymers: evaluation of industrial applications,” Synthetic Metals, vol. 57, no. 1, pp. 3623–3631, 1993.

N. Basescu, Z. X. Liu, D. Moses, A. J. Heeger, H. Naarmann, and N. Theophilou, “High electrical conductivity in doped polyacetylene,” Nature, vol. 327, no. 6121, pp. 403–405, 1987.

A. J. Heeger, “Semiconducting and metallic polymers: the fourth generation of polymeric materials,” Journal of Physical Chemistry B, vol. 105, no. 36, pp. 8475–8491, 2001.

A. J. Heeger, S. Kivelson, J. R. Schrieffer, and W. P. Su, “Solitons in conducting polymers,” Reviews of Modern Physics, vol. 60, no. 3, pp. 781–850, 1988.

H. Zhao, Z. An, and C. Q. Wu, “Lattice vibrational analysis of polyacene,” The European Physical Journal B, vol. 43, no. 1, pp. 53–64, 2005.

H. Naaramann and N. Theophilou, “New process for the production of metal-like, stable polyacetylene,” Synthetic Metals, vol. 22, no. 1, pp. 1–8, 1987.

C. K. Chiang, C. R. Fincher, Y. W. Park et al., “Electrical conductivity in doped polyacetylene,” Physical Review Letters, vol. 39, no. 17, pp. 1098–1101, 1977.

D. S. Suh, T. J. Kim, Y. W. Park, K. Akagi, H. Shirakawa, and J. S. Brooks, “Low dimensional transport properties of iodine doped ultrathin polyacetylene of 100 nm thickness,” Synthetic Metals, vol. 119, no. 1–3, pp. 467–468, 2001.

D. B. Tanner, G. L. Doll, A. M. Rao, P. C. Eklund, G. A. Arbuckle, and A. G. MacDiarmid, “Optical properties of potassium-doped polyacetylene,” Synthetic Metals, vol. 141, no. 1-2, pp. 75–79, 2004.

B. Lesiak, A. Kosinski, M. Krawczyk et al., “Characterization of polyacetylene and polyacetylene doped with palladium,” Polish Journal of Chemistry, vol. 74, no. 6, pp. 847–865, 2000.

R. Landauer, “Conductance determined by transmission: probes and quantised constriction resistance,” Journal of Physics, vol. 1, no. 43, pp. 8099–8110, 1989.

M. Büttiker, “Four-terminal phase-coherent conductance,” Physical Review Letters, vol. 57, no. 14, pp. 1761–1764, 1986.

M. Belhadi, A. Khater, J. Hardy, and K. Maschke, “Phonon transmission via a three-terminal junction in low dimensional wave-guides,” The European Physical Journal—Applied Physics, vol. 35, no. 3, pp. 185–190, 2006.

M. Belhadi, A. Khater, O. Rafil, and J. Hardy, “Scattering of phonons in hexagonal 2D crystals on their domain boundaries,” Physica Status Solidi B, vol. 228, no. 3, pp. 685–693, 2001.

M. Belhadi, R. Chadli, A. Khater, and M. Abou Ghantous, “Spin dynamics across an inhomogeneous atomic boundary separating ultrathin Heisenberg ferromagnetic films,” The European Physical Journal—Applied Physics, vol. 37, no. 1, pp. 25–31, 2007.

M. Belhadi and R. Chadli, “Scattering and localization of spin waves in quasi-2D Heisenberg ferromagnets with extended flat interface boundaries,” Surface Review and Letters, vol. 11, no. 3, pp. 321–329, 2004.

F. Mila and J. Szeftel, “Calculation of phonon dispersion in superlattices using the matching procedure,” Physical Review B, vol. 38, no. 9, pp. 5931–5937, 1988.

J. Szeftel and A. Khater, “Calculation of surface phonons and resonances: the matching procedure revisited: I,” Journal of Physics C, vol. 20, no. 29, pp. 4725–4736, 1987.

W. P. Su, J. R. Schrieffer, and A. J. Heeger, “Solitons in polyacetylene,” Physical Review Letters, vol. 42, no. 25, pp. 1698–1701, 1979.