请问大家能不能导出差分电荷的密度数据呢？哪位大牛帮帮忙 - 第一性原理 - 热点前沿问题讨论

努力学习必有所获

11楼2014-05-26 20:03:40

Joannaouc

木虫 (著名写手)

应助: 46 (小学生)
贵宾: 2.2
金币: 1686.5
散金: 894
红花: 10
帖子: 1262
在线: 266.1小时
虫号: 305787
注册: 2006-12-10
专业: 理论和计算化学

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

9楼: Originally posted by clown9111014 at 2014-05-26 16:21:07
是这样的，castep不能直接输出差分后的电荷密度，也不能输出单个原子的密度文件，所以要用lev000做差只能用vasp输出的结果，这样差分电荷密度图就是vasp得出的，而电荷密度分布是由castep得出的，会不会不准确呀！ ...

我还是没有理解。电荷密度分布已经由castep得出，做个简单的数学计算就能得到差分，为什么还要用lev000做差分？

12楼2014-05-27 10:28:55

clown9111014

银虫 (小有名气)

应助: 15 (小学生)
金币: 140.8
散金: 3
红花: 3
帖子: 255
在线: 155小时
虫号: 2098016
注册: 2012-10-31
性别: GG
专业: 金属材料的合金相、相变及

引用回帖:

12楼: Originally posted by Joannaouc at 2014-05-27 10:28:55
我还是没有理解。电荷密度分布已经由castep得出，做个简单的数学计算就能得到差分，为什么还要用lev000做差分？...

你好，不知你能否得出差分电荷密度数据，不是图，是数据，如果你会请教教我，先谢啦！！

努力学习必有所获

13楼2014-05-27 20:15:11

clown9111014

银虫 (小有名气)

应助: 15 (小学生)
金币: 140.8
散金: 3
红花: 3
帖子: 255
在线: 155小时
虫号: 2098016
注册: 2012-10-31
性别: GG
专业: 金属材料的合金相、相变及

引用回帖:

是这样的，我有电荷密度的三维数据，xyz坐标和电荷密度，做差分要减去对应原子未成键状态电荷，可是要怎么获得这一数据呢，castep好像不能直接给出，而vasp可以给出的。

努力学习必有所获

14楼2014-05-27 20:22:44

Joannaouc

木虫 (著名写手)

应助: 46 (小学生)
贵宾: 2.2
金币: 1686.5
散金: 894
红花: 10
帖子: 1262
在线: 266.1小时
虫号: 305787
注册: 2006-12-10
专业: 理论和计算化学

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

14楼: Originally posted by clown9111014 at 2014-05-27 21:22:44
是这样的，我有电荷密度的三维数据，xyz坐标和电荷密度，做差分要减去对应原子未成键状态电荷，可是要怎么获得这一数据呢，castep好像不能直接给出，而vasp可以给出的。

...

我明白了，你要做的是deformation density而不是大家所认为的density difference。 CASTEP是不能做这个任务，dmol可以，我不知道VASP也可以做这个。

这样的话可以用VASP来计算deformation density，结果不会差太多。

15楼2014-05-28 10:03:18

clown9111014

银虫 (小有名气)

应助: 15 (小学生)
金币: 140.8
散金: 3
红花: 3
帖子: 255
在线: 155小时
虫号: 2098016
注册: 2012-10-31
性别: GG
专业: 金属材料的合金相、相变及

引用回帖:

15楼: Originally posted by Joannaouc at 2014-05-28 10:03:18
我明白了，你要做的是deformation density而不是大家所认为的density difference。 CASTEP是不能做这个任务，dmol可以，我不知道VASP也可以做这个。

这样的话可以用VASP来计算deformation density，结果不会差太 ...

对，我感觉也是，密度图看不出来

努力学习必有所获

16楼2014-05-28 14:58:37

May.Be

新虫 (初入文坛)

应助: 1 (幼儿园)
金币: 39
帖子: 42
在线: 5.7小时
虫号: 3080324
注册: 2014-03-22
专业: 凝聚态物性 II ：电子结构

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

我就是用VASP算的，然后在VESTV中画图的

17楼2014-05-28 16:15:10

clown9111014

银虫 (小有名气)

应助: 15 (小学生)
金币: 140.8
散金: 3
红花: 3
帖子: 255
在线: 155小时
虫号: 2098016
注册: 2012-10-31
性别: GG
专业: 金属材料的合金相、相变及

引用回帖:

17楼: Originally posted by May.Be at 2014-05-28 16:15:10
我就是用VASP算的，然后在VESTV中画图的

真的啊，请问你是否得到了电荷密度数据，就是xyz和对应电荷密度值

努力学习必有所获

18楼2014-05-28 17:36:31

无所谓K

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1777.5
帖子: 9
在线: 18小时
虫号: 3327922
注册: 2014-07-17
专业: 无机非金属类光电信息与功

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

楼主，最终答案是什么呢，如何导出差分电荷密度的数据？求解

19楼2014-09-23 21:12:56

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

专家经验: +224
1ST强帖: 2
应助: 47 (小学生)
金币: 3480.8
散金: 14
红花: 223
帖子: 2477
在线: 1253.7小时
虫号: 1369423
注册: 2011-08-16
性别: GG
专业: 金属材料的磨损与磨蚀
管辖: 第一性原理

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

16楼: Originally posted by clown9111014 at 2014-05-28 14:58:37
对，我感觉也是，密度图看不出来

...

castep可以做形变电荷密度（一次差分），atomic density直接得到的就是，差分电荷密度（二次差分）需要指定sets of atoms。

你看看check2xsf将castep的check文件导出是否能够导出为vasp格式之后再进行lev00作图？http://muchong.com/html/201301/5406515.html

另外还有就是把MS下做到的图贴到Origin下面的Matrix绘制成矢量图，但是我觉得效果不好。

不一定挂在论坛，计算问题问题欢迎留言。

20楼2014-09-24 11:05:18