请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 1801|回复: 32

[VASP] 关于vasp计算两态之间的跃迁偶极矩

[复制链接]

121

帖子

0

威望

303

eV
积分
424

Level 3 能力者

发表于 2018-6-6 10:39:50 | 显示全部楼层 |阅读模式
请问,vasp如何计算两态之间的跃迁偶极矩以及非绝热耦合的大小呢?多谢了

197

帖子

4

威望

1640

eV
积分
1917

Level 5 (御坂)

发表于 2018-6-6 19:49:13 | 显示全部楼层
参见科大赵瑾老师组郑奇靖博士的程序:
https://github.com/QijingZheng/vaspwfc

评分

参与人数 1eV +5 收起 理由
obaica + 5

查看全部评分

121

帖子

0

威望

303

eV
积分
424

Level 3 能力者

 楼主| 发表于 2018-6-14 13:35:11 | 显示全部楼层
万里云 发表于 2018-6-6 19:49
参见科大赵瑾老师组郑奇靖博士的程序:
https://github.com/QijingZheng/vaspwfc

十分感谢您。我试过了,通过该程序可以获得wfc_i.vasp和wfc_r.vasp,那请问,如何得到Hamiltonian的real和imaginary部分呢?

197

帖子

4

威望

1640

eV
积分
1917

Level 5 (御坂)

发表于 2018-6-15 10:28:18 | 显示全部楼层
Theor_Comp 发表于 2018-6-14 13:35
十分感谢您。我试过了,通过该程序可以获得wfc_i.vasp和wfc_r.vasp,那请问,如何得到Hamiltonian的real ...

哈密顿量的对角元是能级,非对角元是波函数内积(<t,i|t+1,j>-<t+1,i|t,j>)/2dt。

翻一翻他的程序就能找到。

121

帖子

0

威望

303

eV
积分
424

Level 3 能力者

 楼主| 发表于 2018-6-15 14:38:09 | 显示全部楼层
万里云 发表于 2018-6-15 10:28
哈密顿量的对角元是能级,非对角元是波函数内积(-)/2dt。

翻一翻他的程序就能找到。

您好!能否稍微详细介绍下呢,我想分别输出elelctronic hamiltonian的real part和imaginary part,麻烦了。

197

帖子

4

威望

1640

eV
积分
1917

Level 5 (御坂)

发表于 2018-6-16 10:02:18 | 显示全部楼层
Theor_Comp 发表于 2018-6-15 14:38
您好!能否稍微详细介绍下呢,我想分别输出elelctronic hamiltonian的real part和imaginary part,麻烦了 ...

https://github.com/QijingZheng/VaspBandUnfolding
看那个nac.py。

121

帖子

0

威望

303

eV
积分
424

Level 3 能力者

 楼主| 发表于 2018-6-16 10:36:57 | 显示全部楼层
万里云 发表于 2018-6-16 10:02
https://github.com/QijingZheng/VaspBandUnfolding
看那个nac.py。

谢谢。这个程序我试了,运行之后得到的是eig.txt和nac.txt,不知如何输出hamiltonian矩阵呢?

197

帖子

4

威望

1640

eV
积分
1917

Level 5 (御坂)

发表于 2018-6-16 15:25:59 | 显示全部楼层
Theor_Comp 发表于 2018-6-16 10:36
谢谢。这个程序我试了,运行之后得到的是eig.txt和nac.txt,不知如何输出hamiltonian矩阵呢?

nac.txt中的数据乘以-ihbar,对角元加上能量,就是哈密顿量了。

121

帖子

0

威望

303

eV
积分
424

Level 3 能力者

 楼主| 发表于 2018-6-16 15:58:49 | 显示全部楼层
万里云 发表于 2018-6-16 15:25
nac.txt中的数据乘以-ihbar,对角元加上能量,就是哈密顿量了。

不太清楚具体如何操作,能否麻烦您再介绍下吗?实在是麻烦您了
另外,怎么单独获得real和imaginary部分的输出呢?

197

帖子

4

威望

1640

eV
积分
1917

Level 5 (御坂)

发表于 2018-6-16 20:17:59 | 显示全部楼层
本帖最后由 万里云 于 2018-6-16 20:38 编辑
Theor_Comp 发表于 2018-6-16 15:58
不太清楚具体如何操作,能否麻烦您再介绍下吗?实在是麻烦您了
另外,怎么单独获得real和imaginary部分 ...
nac_vasp.tar.bz2 (8.28 KB, 下载次数: 11)

评分

参与人数 1eV +3 收起 理由
sobereva + 3

查看全部评分

121

帖子

0

威望

303

eV
积分
424

Level 3 能力者

 楼主| 发表于 2018-6-16 23:54:31 | 显示全部楼层

十分感谢您!

121

帖子

0

威望

303

eV
积分
424

Level 3 能力者

 楼主| 发表于 2018-6-17 10:15:29 | 显示全部楼层

十分谢谢。通过使用以上代码后,确实可以输出0与1之间。不知是否存在这样的代码,通过vasp的MD之后得到的xdatcar,里面包含了多个坐标,分别对每个坐标的构型进行计算,可以得到多个wavecar,之后对它们均进行以上处理,同时输出多个_re.dat和多个_im.dat?

197

帖子

4

威望

1640

eV
积分
1917

Level 5 (御坂)

发表于 2018-6-17 15:06:21 | 显示全部楼层
Theor_Comp 发表于 2018-6-17 10:15
十分谢谢。通过使用以上代码后,确实可以输出0与1之间。不知是否存在这样的代码,通过vasp的MD之后得到的 ...

用shell脚本就能完成,不算很复杂。

121

帖子

0

威望

303

eV
积分
424

Level 3 能力者

 楼主| 发表于 2018-6-17 16:43:56 | 显示全部楼层
万里云 发表于 2018-6-17 15:06
用shell脚本就能完成,不算很复杂。

自己不是太熟悉脚本这块。所有有两个问题:

1. 不知是否有可以参考的脚本,可实现以上目的?

2. 另外,通过您的代码得到的_re.dat和_im.dat 与 通过郑博士nac.py的脚本得到的eig.txt和nac.txt有什么区别呢?

197

帖子

4

威望

1640

eV
积分
1917

Level 5 (御坂)

发表于 2018-6-17 22:43:53 | 显示全部楼层
Theor_Comp 发表于 2018-6-17 16:43
自己不是太熟悉脚本这块。所有有两个问题:

1. 不知是否有可以参考的脚本,可实现以上目的?

1. PYXAID这个程序中有类似的脚本,可以拿来修改;
2. _re.dat和_im.dat是PYXAID格式,程序本身也是拿郑博士的程序改的。联系就是上面提到的,nac.txt中数据乘以-ihbar,对角元再加上eig.txt,就是_re.dat和_im.dat。
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2018-9-25 21:58 , Processed in 0.244318 second(s), 27 queries .

快速回复 返回顶部 返回列表