一维FPU-β介电原子链中的热传导特性研究.pdfVIP

第29卷 湖北师范学院学报 (自然科学版) Vo1．29 第4期 JournalofHubeiNormalUniversity(NaturalScience) No．4，2009 一 维 FPU一 介 电原子链中的热传导特性研究 李金星，赵桂芳，石大庆，朱国良 (湖北师范学院 物理与电子科学学院，湖北 黄石 435002) 摘要：运用椎广的FKM方法对一维FPu一口介 电原子链的热传导特性进行 了研 究，在一定的条件下，分别 得到了热流、热导率和局域动能的近似表达式。通过数值模拟发现，不同的边界条件、非线性相互作用对 热传导特性有很大的影响。 关键词 ：热输运 ；FPU一口链 ；FKM方法 中图分类号：0145 文献标识码 ：A 文章编号：1009-2714(2009)04—0016—04 0 引言 自从Fermi，Pasta和Ulam(FPu)¨的原创性工作做 出以后 ，人们就开始 了对这个领域的广泛研 究 。一维FPU型晶格在非线性科学、分子动力学 以及热传导的研究中起 了很重要的作用 。通 过模拟计算，人们得到在系统趋近于无穷大时，FPU型晶格的热导率是发散的 ，这种模拟给出了热 导率随系统大小的变化关系k。cN 。这方面的工作一般都是基于耦合理论模型 ，而应用Peierls 方程来分析FPU型晶格的热传导只能局限于对周期性链的定性讨论。近年来大量的试验(如单壁纳 米管、绝缘纳米线)也表现出了不可预测的物理机制 J。因此，在量子的热传导中，对边界、接点，特 别是原子之间的非线性相互作用的理解变得极其重要。 在Bardeen的隧道理论中，有几种处理量子热传导的有效方法，像 Landauer公式 ，Keldysh技 术，量子 Langevin方程 和量子散射理论。后来 Ford，Kac和 Mazur(FKM)发展了一种新的方法 ， 这种方法最初用来处理粒子的布朗运动，随后扩展到处理粒子和热源的耦合问题。在处理开放性系 统的时候，FKM方法相对于其他方法更直接，它可以很容易得到热流和局域动能的一般表达式。 本文，我们运用推广的FKM平均场方法来描述与两热源相连接的FPU一 链 ，在导出热流、热导 率和局域动能的表达式后，利用不同的边界条件，对热流进行了数值分析。 1 理论模型 当FPU链组成的系统与左右热源通过导线相连时，总的哈密通量可 以表示为 ： Hr= ” +HL+HR一 ／xx1Xp— XPr (I) 中问N粒子系统的哈密通量可以写为： N 2 N-I 盟 +A 一 Ⅳ = + (2) 其中 和P，分别代表坐标和动量算符，m 为随机质量。第 1个粒子和第 JⅣ个粒子分别与左右热源相 连，这里给定边界条件 。= =0．对于每个热源，我们都用 个谐振子来描述，左边热源的哈密 [收稿 日期] 2oo9—03—24 [基金项 目] 湖北师范学院科研项 目(2007D22) [作者简介] 李金星(1980一 )，男，湖北大悟人 ，硕士，助教 · l6 · 通量可表示为 ： 村 一 一 = ∑ p~／2m+∑KXX／2=∑P／2+()￡霹／2 (3) ，．m 1 K 是弹性耦合的对称矩阵，{X，P}分别代表热源的坐标和动量算符，X，P分别对应相应的正则算 子，而且存在关系 =∑ ，这里的 满足本征