Abstract
We present molecular dynamics simulation of phonon thermal conductivity of semiconductor nanoribbons with an account for phonon quantum statistics. In our semiquantum molecular dynamics simulation, dynamics of the system is described with the use of classical Newtonian equations of motion where the effect of phonon quantum statistics is introduced through random Langevinlike forces with a specific power spectral density (color noise). The color noise describes interaction of the molecular system with the thermostat. The thermal transport of silicon and germanium nanoribbons with atomically smooth (perfect) and rough (porous) edges are studied. We show that the existence of rough (porous) edges and the quantum statistics of phonon change drastically the lowtemperature thermal conductivity of the nanoribbon in comparison with that of the perfect nanoribbon with atomically smooth edges and classical phonon dynamics and statistics. The roughedge phonon scattering and weak anharmonicity of the considered lattice produce a weakly pronounced maximum of thermal conductivity of the nanoribbon at low temperature.
Keywords:
Thermal conductivity; Molecular dynamics simulation; Nanoribbon; Silicon; Germanium; Isotopic effectBackground
It has been recently shown [1] that silicon and germanium nanowires can give a figure of merit of over 2 at 800 K due to strong reduction of phonon thermal conductivity in nanowires as compared with their equivalent bulk material, i.e., the reduction is caused not only by the alloy disorder, but also by the decrease of the phonon mean free path by reduceddimensional effects. The effect of temperature on the thermal conductivity of silicon and germanium may be quite different since the Debye temperature of silicon almost doubles that of germanium. The purpose of the present work is to analyze quantum statistic effects on thermal phonon conductivity in silicon and germanium nanoribbons with the use of the novel semiquantum molecular dynamics simulation [2].
Molecular dynamics is a method of numerical modeling of molecular systems based on classical Newtonian mechanics. It does not allow for the description of pure quantum effects such as the freezing out of highfrequency oscillations at low temperatures and the related decrease to zero of heat capacity for T→0. On the other hand, because of its complexity, a pure quantummechanical description does not allow, in general, the detailed modeling of the dynamics of manybody systems. To overcome these obstacles, different semiclassical methods, which allow to take into account quantum effects, have been proposed [39].
The most convenient method for numerical modeling is to use the Langevin equations with colornoise random forces [5,7]. In this approximation, the dynamics of the system is described with the use of classical Newtonian equations of motion while the quantum effects are introduced through random Langevinlike forces with a specific power spectral density (the color noise), which describes the interaction of the molecular system with the thermostat. Here, we apply such semiquantum approach to the simulation of heat transport in lowdimensional nanostructures such as semiconductor nanoribbons with atomically smooth (perfect) and porous (rough) edges. Our previous analytical studies [10] and molecular dynamics simulations [11] have revealed the dramatic decrease of phonon thermal conductivity in quasionedimensional nanostructures with rough (porous) surface and edge layers.
Methods
In the semiquantum molecular dynamics approach, the dynamics of the system is described
with the use of the classical Newtonian equations of motion while the effects of phonon
quantum statistics are introduced through random Langevinlike forces with a specific
power spectral density (the color noise). If the random forces are deltacorrelated
in a time domain, this corresponds to the white noise with a flat power spectral density.
This situation corresponds to highenough temperatures, when k_{B}T is larger than the quantum of the highest phonon frequency mode in the system,
Results and discussion
We consider a system which consists of K parallel atomic chains in one plane [11]. To model the diamondlike lattice, we assume that each atom has four nearest neighbors. In this connection, we would like to mention that the considered model cannot be applied directly to the predicted [1619] and recently grown [20,21] twodimensional lattice with graphenelike structure, made from Si or Ge atoms, the silicene. Our main goal is to provide semiquantum modeling of the heat transport and effective ‘isotopic effect’ on phonon heat transport in lowdimensional structures made from Si or Ge atoms, arranged in lattices, which reflect the symmetry of corresponding bulk materials. Since the lattice structure (the number of nearest neighbors) of the considered quasitwodimensional nanoribbons reflects the bulk one, our model can also be applied to the quasithreedimensional nanowires with bulklike structure. The isotopic effect on phonon heat transport can be used for the understanding and prediction of the trends in the changes of thermal conductivity in lowdimensional nanostructures caused by the essential change in ion masses accompanied by less strong change in interion force constants.
The Hamiltonian of the system describes the kinetic energy and harmonic interparticle interaction potentials. The characteristic energy of the nearestneighbor interaction energy E_{0} can be related with the energy of the LO phonon mode in the semiconductor, which is approximately 15 THz in Si and approximately 9 THz in Ge. The ratio of these maximal frequencies is close to the ratio of the Debye temperatures, T_{D }= 645 K in Si and T_{D }= 374 K in Ge, and to the ratio of the inverse square root of Si and Ge atomic masses, which reflect the approximate isotopic effect in phonon properties of Si and Ge lattices when the materials can be described approximately with the same force constants and different atomic masses (see [22]). The particle mass (M) and lattice constant (a) are determined by the mass and characteristic period of the corresponding bulk semiconductor material, a = 5.43 Å and a = 5.658 Å for Si and Ge, respectively.
We consider a ribbon which consists of K = 18 atomic chains. To model the roughness of the ribbon edges, we delete with probability (porosity) p = 1− d some atoms from K_{1} chains adjacent to each ribbon edge. Here, K_{1} is a width of the rough edges, and d, 0 ≤ d ≤ 1, is a fraction of the deleted atoms in the edge atomic chains. In our simulations, we take K_{1} = 4 and d = 0.80. In Figure 1, we show an example of the nanoribbon with porous edges, cut from the twodimensional diamondlike lattice in which each atom has four nearest neighbors.
Figure 1. Nanoribbon with porous edges cut from twodimensional diamondlike lattice where each atom has four nearest neighbors.
We computed the thermal conductivity κ(NT) for the nanoribbons with the length of N = 500 unit cells. In Figure
2, we plot the results of the semiquantum molecular dynamics simulation of thermal
conductivity of Si and Ge nanoribbons. As one can see in this figure, the thermal
conductivities of both Si and Ge nanoribbons have a weakly pronounced maximum at low
temperatures, T_{max }= 85 K for Si and T_{max }= 91 K for Ge. This property of thermal conductivity temperature dependence is a consequence
of roughedge scattering as the main phonon scattering mechanism at elevated temperatures
and the absence of (or weak) anharmonicity of the lattice potential and correspondingly
the absence of (or weak) anharmonic (Umklapp) scattering. The latter causes a clear
peak in the thermal conductivity versus temperature both in finite bulk crystals of
pure silicon
[23] and in lowdimensional nanoribbons
[2]. The values of thermal conductivities of the Si and Ge nanoribbons for T > T_{max} approximately reproduce an isotopic effect because
Figure 2. Thermal conductivity κ of roughedge nanoribbon versus temperature for ribbon length of N = 500 unit cells. Thermal conductivity κ of roughedge nanoribbon (ribbon width K = 18 atomic chains, rough edges widths K_{1} = 4 atomic chains, porosity of rough edges p = 0.20) versus temperature T for ribbon length of N = 500 unit cells of the twodimensional diamondlike lattice of Ge (blue circles, line 1) or Si (red diamonds, line 2) atoms.
Conclusions
Semiquantum molecular dynamics simulations with random Langevinlike forces with a
specific power spectral density show that quantum statistics of phonons and porosity
of edge layers dramatically change the thermal conductivity of Si and Ge nanoribbons
at low and room temperatures in comparison with that of the nanoribbons with perfect
edges and classical phonon dynamics and statistics. Phonon scattering by the rough
edges and weak anharmonicity of the considered lattice produce weakly pronounced maximum
of the thermal conductivity of the nanoribbon at low temperature. The approximate
isotopic effect is manifested in the
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
This work was finished through the collaboration of all authors. YAK proposed the model for the lattice and isotopic effect. AVS has been working on the MD simulation. YAK and AC have participated in the interpretation of results and in revising the manuscript. All authors read and approved the final manuscript.
Acknowledgements
The authors thank the financial support given by the project CSD20100044, which belongs to the ‘Consolider Ingenio’ Programme of the Spanish Ministry of Finances and Competitiveness.
References

Lee EK, Yin L, Lee Y, Lee JW, Lee SJ, Lee J, Cha SN, Whang D, Hwang GS, Hippalgaonkar K, Majumdar A, Yu C, Choi BL, Kim JM, Kim K: Large thermoelectric figureofmerits from SiGe nanowires by simultaneously measuring electrical and thermal transport properties.

Savin AV, Kosevich Yu A, Cantarero A: Semiquantum molecular dynamics simulation of thermal properties and heat transport in lowdimensional nanostructures.

Wang JS: Quantum thermal transport from classical molecular dynamics.
Phys Rev Lett 2007, 99:160601. PubMed Abstract  Publisher Full Text

Donadio D, Galli G: Thermal conductivity of isolated and interacting carbon nanotubes: comparing results from molecular dynamics and the Boltzmann transport equation.
Phys Rev Lett 2007, 99:255502. PubMed Abstract  Publisher Full Text

Heatwole EM, Prezhdo OV: Secondorder Langevin equation in quantized Hamilton dynamics.
J Physical Soc Jpn 2008, 77:044001. Publisher Full Text

Buyukdagli S, Savin AV, Hu B: Computation of the temperature dependence of the heat capacity of complex molecular systems using random color noise.

Ceriotti M, Bussi G, Parrinello M: Nuclear quantum effects in solids using a colorednoise thermostat.
Phys Rev Lett 2009, 103:030603. PubMed Abstract  Publisher Full Text

Dammak H, Chalopin Y, Laroche M, Hayoun M, Greffet JJ: Quantum thermal bath for molecular dynamics simulation.
Phys Rev Lett 2009, 103:190601. PubMed Abstract  Publisher Full Text

Wang JS, Ni X, Jiang JW: Molecular dynamics with quantum heat baths: application to nanoribbons and nanotubes.

Kosevich YuA: Multichannel propagation and scattering of phonons and photons in lowdimension nanostructures.

Kosevich Yu A, Savin AV: Reduction of phonon thermal conductivity in nanowires and nanoribbons with dynamically rough surfaces and edges.
Europhys Lett 2009, 88:14002. Publisher Full Text

Turney JE, McGaughey AJH, Amon CH: Assessing the applicability of quantum corrections to classical thermal conductivity predictions.

Mingo N: Calculation of Si nanowire thermal conductivity using complete phonon dispersion relations.

Martin P, Aksamija Z, Pop E, Ravaioli U: Impact of phononsurface roughness scattering on thermal conductivity of thin Si nanowires.
Phys Rev Lett 2009, 102:125503. PubMed Abstract  Publisher Full Text

Zhang W, Mingo N, Fisher TS: Simulation of phonon transport across a nonpolar nanowire junction using an atomistic Green’s function method.

Roethlisberger U, Andreoni W, Parrinello M: Structure of nanoscale silicon clusters.
Phys Rev Lett 1994, 72:665. PubMed Abstract  Publisher Full Text

Takeda K, Shiraishi K: Theoretical possibility of stage corrugation in Si and Ge analogs of graphite.
Phys Rev B 1994, 50:14916. Publisher Full Text

Cahangirov S, Topsakal M, Aktuerk E, Seahin H, Ciraci S: Two and onedimensional honeycomb structures of silicon and germanium.
Phys Rev Lett 2009, 102:236804. PubMed Abstract  Publisher Full Text

Houssa M, Pourtois G, Afanasiev VV, Stesmans A: Can silicon behave like graphene? A firstprinciples study.
Appl Phys Lett 2010, 97:112106. Publisher Full Text

Vogt P, De Padova P, Quaresima C, Avila J, Frantzeskakis E, Asensio MC, Resta A, Ealet B, Le Lay G: Silicene: compelling experimental evidence for graphenelike twodimensional silicon.
Phys Rev Lett 2012, 108:155501. PubMed Abstract  Publisher Full Text

Fleurence A, Friedlein R, Ozaki T, Kawai H, Wang Y, YamadaTakamura Yu: Experimental evidence for epitaxial silicene on diboride thin films.
Phys Rev Lett 2012, 108:245501. PubMed Abstract  Publisher Full Text

Ziman JM: Electrons and Phonons. Oxford: Oxford University Press; 1960.

Klitsner T, VanCleve JE, Fisher HE, Pohl RO: Phonon radiative heat transfer and surface scattering.
Phys Rev B 1988, 38:7576. Publisher Full Text

Lim J, Hippalgaonkar K, Andrews SC, Majumdar A, Yang P: Quantifying surface roughness effects on phonon transport in silicon nanowires.