漫谈分子动力学计算热导率的五种方法

2023-08-08

漫谈分子动力学计算热导率的五种方法

本仓库提供了一个资源文件,详细介绍了使用分子动力学(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功能提出。我们非常欢迎您的反馈和贡献,以帮助我们不断改进和完善本资源文件。

下载链接

漫谈分子动力学计算热导率的五种方法分享