Modern Physics
Vol.07 No.05(2017), Article ID:21768,6 pages
10.12677/MP.2017.75019

Theoretical Study on the Heat Conductivities of Cyclic Molecules in Molecular Machines

Jianping Huang1*, Jing Tang2

1College of Information Science and Technology, Hunan Normal University, Changsha Hunan

2College of Physics and Information Science, Hunan Normal University, Changsha Hunan

Received: Aug. 3rd, 2017; accepted: Aug. 17th, 2017; published: Aug. 22nd, 2017

ABSTRACT

The formulas for energy flux of atomic vibrations and linewidths of phonons of cyclic molecules in molecular machines were derived based on the lattice dynamics, and then based on those formulas and Green-Kubo formula, the formula for heat conductivities of cyclic molecules was derived. Finally the numerical calculations were carried out. The numerical results show that the heat conductivity of a short cyclic molecule increases when its length increases and will tend to infinity when its length tends to infinity, so there is size effect of the heat conductivity in a cyclic molecule with any length. The numerical results also show that the main contributions to the heat conductivity of a long cyclic molecule are made by phonons with short wave vectors.

Keywords:Molecular Machine, Cyclic Molecule, Heat Conductivity, Linewidth of Phonon, Lattice Dynamics

分子机器中环状分子热传导性质的理论研究

黄建平1*,唐婧2

1湖南师范大学信息科学与工程学院,湖南 长沙

2湖南师范大学物理与信息科学学院,湖南 长沙

收稿日期:2017年8月3日;录用日期:2017年8月17日;发布日期:2017年8月22日

摘 要

运用晶格动力学理论推导了分子机器中环状分子振动的能量通量公式,在此基础上再应用格林-久保公式推导了环状分子的热传导系数公式。由于环状分子的热传导系数与声子的谱线宽度有关,因此还推导了其声子谱线宽度公式。数值计算结果表明,环状分子的热传导系数随着环状分子长度的增加而增加,并且在环状分子无限增长时趋于无限,因此任何长度的环状分子其热传导性质都存在尺寸效应。数值计算结果还表明,长环状分子的热传导系数主要来自于短波矢声子的贡献。

关键词 :分子机器,环状分子,热传导系数,声子谱线宽度,晶格动力学

Copyright © 2017 by authors and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

1. 引言

诺贝尔奖得主费曼于1959年美国物理学会年会上所做的演讲中,提出了纳米科技的概念。分子机器属于纳米科技范畴,是最小尺度的微纳机电系统,具有小尺寸、多样性、自适应,以及仅依靠热能、光照或化学能驱动等性能,在生物医药、新材料、传感器等领域具有广阔的应用前景。由于在分子机器研究方面做出了杰出贡献 [1] [2] [3] ,索维奇、斯图达特和费林加共同获得2016年诺贝尔奖。

分子机器常用热能驱动。只有通过高效的热传导,才能使热量到达指定位置及部件而进行高效的热驱动。也只有通过高效的热传导,才能避免分子机器工作期间聚集过多的热量,而使部件乃至整个分子机器因过热而失效。因此研究环状分子这种分子机器中最常用构件的热传导性质,显得十分重要。由于分子机器尺寸过于微小,因此目前通过实验手段对其构件的热传导性质进行研究还存在较多困难。分子动力学方法被广泛运用于材料的热传导性质模拟 [4] ,由于它是一种建立在牛顿力学基础上的纯粹的数值计算方法,因而无法通过解析公式分析物理机制,也无法考虑量子效应而得到准确的计算结果。我们曾运用晶格动力学研究了纳米晶体薄膜热传导性质 [5] 。由于晶格动力学考虑了晶格振动的量子效应,并且基于该理论能推导得到材料的热传导系数的解析公式,因此晶格动力学克服了分子动力学的缺点,本文仍将使用该方法对环状分子热传导性质进行计算。

2. 环状分子的晶格动力学与热传导理论

设有个相同原子构成的环状分子,原子间距为。原子振动位移和动量可表示为 [6] :

(1)

(2)

其中,为原子质量,是波矢,为环状分子波动的频率,它们之间有以下色散关系:

(3)

其中,为原子间线性相互作用力常数,代表代表声子产生与消灭算符,波矢取值,其中,为整数,且

环状分子中的能量通量可表示为 [7]

(4)

将(1)式和(2)式代入上式,得环状分子的能量通量公式

(5)

其中表示,为声子数算符,为环状分子中的波动群速度。

将(5)式代入以下的Green-Kubo公式 [8]

(6)

并按以下公式解耦 [9]

(7)

可得

(8)

考虑声子衰减,可得以下公式 [10] :

(9)

(10)

其中,是波矢为的声子谱线宽度。将(9)式和(10)式代入(8)式,得

(11)

上式中的每个求和项就代表波矢为的声子对环状分子热传导系数的贡献

3. 环状分子的声子格林函数与声子寿命

环状分子中原子间有非和谐相互作用势能项,忽略四次以上非和谐势能项,哈密顿可表示为 [11]

(12)

定义声子Green函数 [10]

(13)

其中,为单位阶梯函数,该式右边也可记为

将(13)式对时间求1阶导数,得

(14)

其中

(15)

将(15)式代入(14)式,并进行傅里叶变换,得

(16)

用同样的方法,将对时间求导,再进行傅里叶变换,得

(17)

由(16)式可知,上式中的可近似为。将上式代入(16)式,得

(18)

其中

(19)

利用(9)式和(10)式,并考虑到,由上式得

(20)

其中,

(21)

为三阶非和谐相互作用力常数。当等于0或倒格矢的整数倍时,等于1,其它情形为0。将中的取值记为分解为实部和虚部之和。由于声子格林函数的极点代表声子频率,因此为声子谱线宽度,可表示为

(22)

取值为时,(22)式转化为

(23)

上式与Maradudine等人 [11] 得到的公式是一致的,这也说明推导(22)式的过程是正确的。Maradudine等人 [11] 在使用(23)式时,用替代进行数值计算,人为取值的对计算结果有很大影响,影响了结果的准确性。而基于(22)式进行迭代计算时,能快速收敛,声子谱线宽度初值对计算结果影响很小,能保证计算结果准确。

4. 数值计算结果与分析

虽然氩原子构成的环形构件在分子机器中没有具体应用,但由于氩原子间的相互作用能用L-J模型准确描述 [12] ,因此在本文中通过对它进行计算来了解环状分子的热膨胀性质。设温度T = 60 K,此时远大于所有声子能量,计算可采用高温近似。

不同长度的环状分子的声子谱线宽度与波矢之间的关系,如图1所示。环状分子声子中,最短波矢为,随长度增加,最短波矢减小,由图1可知,其谱线宽度也明显减小,当长度趋于无穷大时,最短波矢趋于无穷小,其对应的声子谱线宽度也趋于无穷小。

环状分子的热传导系数与长度之间的关系如图2所示。由图2可知,环状分子的热传导系数随长度的增加而增加,并随长度无限增加而趋于无穷。

为了分析环状分子的热传导系数随长度变化的物理机制,我们计算了不同长度环状分子中的各种波矢的声子对热传导的贡献,如图3所示。由图3可知,随环状分子的长度增加,关系曲线的重心向短波矢方向移动,短波矢的声子对热传导系数的贡献变大,长环状分子热传导系数主要来自于短波矢声子的贡献。当N趋于无穷时,最短波矢也趋于无穷,由图1知其声子谱线宽度趋于零,最后由(11)式可知热传导系数也趋于无穷,因此任何长度的环状分子其热传导系数都存在尺寸效应。

Figure 1. Phonon linewidth vs. wave vector in cyclic molecule

图1. 环状分子声子谱线宽度与波矢之间的关系

Figure 2. Heat conductivity vs. length of cyclic molecule

图2. 环状分子热传导系数与长度关系

Figure 3. The contribution made by phonon to heat conductivity vs. wave vector in cyclic molecule

图3. 环状分子中声子对热传导系数贡献与波矢之间关系

5. 结论

本文运用晶格动力学理论推导了环状分子的热传导系数公式和声子谱线宽度公式,并进行了数值计算。我们发现环状分子的热传导系数的贡献主要来自于短波矢声子;环状分子越长,短波矢声子的波矢和谱线宽度越小,这些短波矢声子对热传导系数的贡献也越大;当分子长度趋于无穷时,这些短波矢声子的波矢和谱线宽度也趋于零,这些短波矢声子对热传导系数的贡献趋于无穷,这是环状分子的热传导系数随其长度的增加而增加并最终发散的根本原因,也是环状分子热传导系数的尺寸效应存在于任何长度的环状分子中的根本原因。

文章引用

黄建平,唐 婧. 分子机器中环状分子热传导性质的理论研究
Theoretical Study on the Heat Conductivities of Cyclic Molecules in Molecular Machines[J]. 现代物理, 2017, 07(05): 169-174. http://dx.doi.org/10.12677/MP.2017.75019

参考文献 (References)

  1. 1. Carina, R.F., Dietrich-Buchecker, C. and Sauvage J.P. (2016) Molecular Composite Knots. Journal of the American Chemical Society, 118, 9110-9116. https://doi.org/10.1021/ja961459p

  2. 2. Badjic, J.D., Ronconi, C.M., Stoddart, J.F., et al. (2015) Operating Molecular Elevators. Journal of the American Chemical Society, 128, 1489-1499. https://doi.org/10.1021/ja0543954

  3. 3. Ruangsupapichat, N., Pollard, M.M., Harutyunyan, S.R., et al. (2011) Reversing the Direction in a Light-Driven Rotary Molecular Motor. Nature Chemistry, 3, 53-60. https://doi.org/10.1038/nchem.872

  4. 4. Wu, G. and Dong, J. (2005) Anomalous Heat Conduction in a Carbon Nanowire: Molecular Dynamics Calculations. Physical Review B, 71, Article ID: 115410. https://doi.org/10.1103/PhysRevB.71.115410

  5. 5. 李圣怡, 黄建平. 基于晶格动力学的纳米薄膜热特性研究[J]. 微纳电子技术, 2008, 45(5): 249-254.

  6. 6. Bottger, H. (1983) Principles of the Theory of Lattice Dynamics. Physik-Verlag, Weinheim, 85-90.

  7. 7. Frizzera, W., Viliani, G., Montagna, M., et al. (1997) One-Dimensional Dynamics of Lattice Heat Transport. Journal of Physics: Condensed Matter, 9, 10867-10879. https://doi.org/10.1088/0953-8984/9/49/007

  8. 8. Wang, J. and Li, B. (2004) Mode-Coupling Theory and Molecular Dynamics Simulation for Heat Conduction in a Chain with Transverse Motions. Physical Review E, 70, Article ID: 021204. https://doi.org/10.1103/PhysRevE.70.021204

  9. 9. Wilson, R.S. and Kim, S.K. (1962) Theory of Heat Conductivity of Anharmonic Crystals. Physical Review, 128, 4674-4677.

  10. 10. Semwal, B.S. and Sharma, P.K. (1972) Heat Conductivity of an Anharmonic Crystal. Physical Review B, 5, 3909-3913.

  11. 11. Maradudin, A.A. and Fein, A.E. (1973) Scattering of Neutrons by an Anharmonic Crystal. Physical Review B, 7, 2589-2608.

  12. 12. Hyzorek, K. and Tretiakov, K.V. (2016) Thermal Conductivity of Liquid Argon in Nanochannels from Molecular Dynamics Simulations. The Journal of Chemical Physics, 144, Article ID: 194507. https://doi.org/10.1063/1.4949270

  13. NOTES

    *通讯作者。

期刊菜单