用户工具

站点工具

本页面的其他翻译:
  • zh

adf:dftbmdexample

这是本文档旧的修订版!


使用DFTB-MD模拟银表面甲醇分子吸附的动力学过程

DFTB-MD适用于三维、二维、一维周期性体系,本例是一个二维体系。本例包含32个银原子,一个甲醇分子,使用精确度较高的SCC-DFTB方法,模拟步数为3000步,在双核Mac Book Pro计算耗费时间约为19分钟。

建模

点击ADFinput右上角的🔍符号,输入silver,搜索软件自带的Ag晶体。

转为常用单胞:

点击AMSinput窗口底部的四边形按钮(该按钮用于切割表面),作为演示,我们设置切面米勒指数为111,厚度为2层。切出这个表面之后,体系的周期性变成了二维,原子沿着xy平面无限延伸,z方向上下为半无限大的真空,因此Main面板中Periodicity自动变为Slab(Bulk为三维周期性,Chain为一维周期性)。

我们希望做一个大一些大超胞,例如2*2:点击ADFinput窗口底部大❄️状按钮 → Generate Super Cell → 对角线上两个数字修改为2(BAND也支持非整数超胞,点击Preset即可选择)。这样就得到一个单胞包含32个银原子的二位无限大银薄膜。在表面画一个甲醇分子。

参数设置

如下图:

说明:

  • 精度方面,DFTB3 > SCCDFTB > DFTB
  • 对分子晶体而言,K-space可以选择Normal即可,与选择Good区别不大,对半导体等晶体结构选择Good、Very Good更好。
  • 如果是分子晶体优化,出现ERROR: PeriodicDFTBModule: error with occful的报错,可以考虑降低k点个数,改为basic
  • Parameter Director选择DFTB参数。会自动根据体系元素、Method,列出所有可选参数,用户需要自行测试哪个参数对自己的体系而言,精度更高(点击下拉窗口后面的?按钮,可以显示该参数的详细介绍)
  • 分子晶体可以使用色散修正,例如D3-BJ,但这会让结构收敛变得更困难。不存在分子间相互作用的体系,不需要使用色散修正。有的参数、方法不支持色散修正

点击Task:Molecular Dynamics后面的 > 图标,可以对分子动力学进行详细设置,例如:设置模拟30000步,步长为0.25fs,每个10步保存一次原子轨迹:

点击Thermostat后面的 > 按钮,设置NVT系综(点击Barostat,则设置系综关于压强的部分,适用于NPT系综): 其中Thermostat是实现NVT系综的算法,可以选择NHC。

如果要实现温度变化:

温度如果只设置了一个,则不需要设置Duration(s),整个模拟过程都是该温度;如果需要改变温度,则如上图所示设置多个温度即可,配合Duration(s)的设置,可以实现控温。例如图中所示表示:298K运行5000步,然后348K运行25000步,然后回到298K直到结束。所以Duration(s)的数字个数比温度的数字个数少1个。

结果查看

SCM LOGO → Movie,点击窗口底部的播放按钮,可以看到原子的运动轨迹:

如果不是表面体系吸附,而是纯粹的分子混合物体系,则Movie中的MD Properties选项会多很多,包括Molecules可以显示分子数量变化曲线:

Movie - MD properties - Reaction Event Detection,可以分析单步反应(基元反应):

  • Start at Timestep,从多少步开始分析
  • Process steps,分析接下来的多少步
  • Recrossing Filter,该数值设置为保存轨迹频率的10倍,例如这里设置为100

之后点击Process,处理完毕之后,点击Browse即可在网页中看到分析信息:

  • 首页是排名前五的“反应物”、产物,以及发生次数最多的前五个反应
  • 分子结构都可以点击而得到该分子的“来源”(参与过的基元反应的反应物有哪些)与“去向”(参与过的基元反应的产物有哪些)
  • 右上角有其他汇总信息,例如所有反应、所以物种、时间线(每种物质首次出现的时间线)等等
adf/dftbmdexample.1606220887.txt.gz · 最后更改: 2020/11/24 20:28 由 liu.jun

© 2014-2022 费米科技(京ICP备14023855号