【亲测免费】 漫谈分子动力学计算热导率的五种方法

本仓库提供了一个资源文件，详细介绍了使用分子动力学（Molecular Dynamics, MD）计算材料热导率的五种方法。这些方法包括平衡态的GK法和非平衡态的NEMD/RNEMD（速度重标法、局部热浴法、动量交换法）。资源文件中包含了动态group设置、脚本获取、后处理分析和结果对比等内容，适用于不同工况下的热导率研究。

内容概述

1. 平衡态分子动力学（EMD）

   • GK法（Green-Kubo方法）
   • 方法说明
   • 报道结果

2. 非平衡态分子动力学（NEMD）

   • 速度重标法
   • 局部热浴法
   • 动量交换法
   • 方法说明
   • 后处理分析
   • 计算结果

3. 非平衡态分子动力学（RNEMD）

   • 速度重标法
   • 局部热浴法
   • 动量交换法
   • 方法说明
   • 后处理分析
   • 计算结果

使用说明

1. 动态group设置

   • 针对流体体系，热源和热汇需要设置为动态group。

2. 代码说明

   • 以固态60K下的Argon为例进行测试，脚本分为两个：Run in以及几种不同方法。

3. 代码获取

   • 全部in文件已经开源在博客，可以直接复制。
   • Python绘图脚本以及In文件需要打包文件请下载。

后处理分析

• 温度梯度拟合
• 能量输出
• 温度拟合
• 报道结果

结果对比

• 五种方法结果对比
• 注意动态group设置

适用范围

本资源文件适用于需要进行热导率计算的研究人员，特别是使用LAMMPS进行分子动力学模拟的科研工作者。通过本资源文件，用户可以了解并应用不同的热导率计算方法，以满足不同工况下的研究需求。

贡献与反馈

如果您在使用过程中有任何问题或建议，欢迎通过GitHub的Issue功能提出。我们非常欢迎您的反馈和贡献，以帮助我们不断改进和完善本资源文件。