【摘要】 作者自己做了一次氩气的计算，推荐同学们在学习的时候自己先把这个Ar示例也跑一次，没问题了再去跑自己实际体系。

Forcite计算热导率的脚本原始地址是：

https://community.3dsbiovia.com/Communities_Topics?

id=90650000000L6FDAA0&name=Materials+Studio&tid=09a500000004QfAAAU#/feedtype=SINGLE_QUESTION_DETAIL

 

*Forcite计算热导率的脚本资源文件的获取方式见文末

 

现在该页面已经从网站移除或者转移到其他路径了。

 

脚本的作者是Reinier Akkermans博士，一位“普普通通”的达索系统 BIOVIA 的首席软件开发科学家，主要研究与开发方向是数学建模、化学、生物信息学和计算化学方面。

 

作者在发布脚本时候给出了具体的研究说明手册，具体见下方图片或者附件中的PDF文件。

 

作者自己做了一次氩气的计算，推荐同学们在学习的时候自己先把这个Ar示例也跑一次，没问题了再去跑自己实际体系。

 

以下把作者的几个叮嘱个再唠叨一下：

 

1 推荐使用MS的Amorphous Cell模块搭建结构（其实第三方工具也好使）

2 建议建模时候以C方向为延展方向（方便脚本改写和分析数据）

3 建模时候长边边长至少是短边边长的3倍

4 计算热导率的结构应该是在计算热导率要求的平均温度下充分弛豫的结构。

5 全原子力场会高估热导率（原理上导致的，无法修正，具体解释看PDF文档中Known issues部分）

 

该脚本基于以下两篇文章

 

1.定流量法

Jund, P.; Jullien, R., "Molecular dynamics calculation of the thermal conductivity of vitreous silica.", Phys. Rev. B, 59, 13707-13711 (1999).

 

2.可变流量法

Müller-Plathe, F., "A simple nonequilibrium molecular dynamics method for calculating the thermal conductivity.", J. Chem. Phys., 106, 6082-6085 (1997).

 

根据以下两个公式便可计算得到 z 方向上的热导率，J为 z 方向上的热流，dT/dz是z方向上的温度梯度。这个Z方向就是我们刚刚说的C方向。

 

由于通量的方向与梯度相反，热导率总是正的。

 

A是通量的横截面面积。

 

 

 

如果还想了解关于Perl脚本语言在MS中的应用可点击链接查看详情~

 

文件获取

科学指南针一模拟计算联盟

后台回复：0910

 

本文所有内容文字、图片和音视频资料，版权均属科学指南针网站所有，任何媒体、网站或个人未经本网协议授权不得以链接、转贴、截图等任何方式转载。