计算化学公社

标题: 非绝热分子动力学拟合出两个纳米级寿命 [打印本页]

作者
Author: Shana 时间: 2024-10-28 16:52
标题: 非绝热分子动力学拟合出两个纳米级寿命
各位老师和同学. 请教大家一个关于非绝热分子动力学的问题. 我所采用的软件是vasp+hefei-namd. 模拟载流子的带间复合速率. 采样50个样本, 1000条轨迹,模拟10ns.

我所研究的体系是直接带隙,但是gap在(0.5 0.5 0.5)点. 我采用了两种方法进行AIMD的计算.

(1) 采用2*2*2扩胞方法 (40个原子, K-mesh 2*2*2), 将gap折叠到gamma点, 然后进行AIMD和SCF的计算, 得到非绝热耦合系数. 最后得到的电子布局数变化如下图所示. 可以看到能够很好的被单指数拟合,寿命约为20ns.
(, 下载次数 Times of downloads: 12)
(2) 采用3*3*3扩胞方法(135个原子, 单gamma点计算), 此情形下gap仍然在(0.5 0.5 0.5)点. 所以我先用单gamma点做AIMD,然后使用单(0.5 0.5 0.5)点进行SCF计算,得到非绝热耦合系数. 最后得到的电子布局数变化如下图所示. 此时出现了两个截然不同的变化趋势, 在4ns之前, 电子复合的比较慢. 4ns之后, 电子快速的复合.
(, 下载次数 Times of downloads: 8)

请教一下各位, 这种区别可能是哪些原因造成的?

作者
Author: huanghao 时间: 2024-11-19 10:25
您好楼主，我可以请教一下这个电子布局数变化是怎么处理出来的吗。我按照hefei-NAMD的教程最后使用poen_fssh.py处理出了平均电子能量随时间变化的图。不知道如何处理出电子布局数和时间变化曲线，之后怎么拟合出载流子寿命也没找到明确的教程。请问楼主有相关教程可以分享一下吗。

作者
Author: Shana 时间: 2024-11-26 11:16

huanghao 发表于 2024-11-19 10:25
您好楼主，我可以请教一下这个电子布局数变化是怎么处理出来的吗。我按照hefei-NAMD的教程最后使用poen_fss ...

褚老师的github里面有画布局数的例子。至于拟合寿命，用python的fit_curve就可以了

作者
Author: fany 时间: 2025-1-2 01:46

Shana 发表于 2024-11-26 11:16
褚老师的github里面有画布局数的例子。至于拟合寿命，用python的fit_curve就可以了

您好楼主，我正在计算材料中载流子的布局数。在褚老师的github里面我没有找到画布局数的python脚本，请问可以分享一下吗？非常感谢！

作者
Author: Shana 时间: 2025-1-3 09:51

fany 发表于 2025-1-2 01:46
您好楼主，我正在计算材料中载流子的布局数。在褚老师的github里面我没有找到画布局数的python脚本，请问 ...

这个脚本在fsshmp文件夹里面，poen_fssh.py

作者
Author: Lgd--NAMD 时间: 2025-4-11 15:15
你好，有几个问题想问一下可以加一下微信吗

欢迎光临计算化学公社 (http://ccc.keinsci.com/)