应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 热点前沿问题讨论 » 请问大家能不能导出差分电荷的密度数据呢？哪位大牛帮帮忙
202/2返回列表上一页12
查看: 3153  |  回复: 19
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

clown9111014

银虫 (小有名气)
引用回帖:
10楼: Originally posted by jpchou at 2014-05-26 15:29:17
同意

差分只是後處理
vasp or castep 都能得到 charge density

得到 charge density 就可以做差分
vasp 可以用 lev00 做差分
castep 要用... ? 這我就不知了

若是能自己寫 code，就沒有這種困擾嘍...

好吧，谢谢你们了!!!
赞一下 回复此楼
努力学习 必有所获
11楼2014-05-26 20:03:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Joannaouc

木虫 (著名写手)
★
小木虫: 金币+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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

clown9111014

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

你好，不知你能否得出差分电荷密度数据，不是图，是数据，如果你会请教教我，先谢啦！！
赞一下 回复此楼
努力学习 必有所获
13楼2014-05-27 20:15:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

clown9111014

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

是这样的，我有电荷密度的三维数据，xyz坐标和电荷密度，做差分要减去对应原子未成键状态电荷，可是要怎么获得这一数据呢，castep好像不能直接给出，而vasp可以给出的。
赞一下 回复此楼
努力学习 必有所获
14楼2014-05-27 20:22:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Joannaouc

木虫 (著名写手)
★
小木虫: 金币+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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

clown9111014

银虫 (小有名气)
引用回帖:
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

May.Be

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
我就是用VASP算的，然后在VESTV中画图的
赞一下 回复此楼
17楼2014-05-28 16:15:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

clown9111014

银虫 (小有名气)
引用回帖:
17楼: Originally posted by May.Be at 2014-05-28 16:15:10
我就是用VASP算的，然后在VESTV中画图的

真的啊，请问你是否得到了电荷密度数据，就是xyz和对应电荷密度值
赞一下 回复此楼
努力学习 必有所获
18楼2014-05-28 17:36:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

无所谓K

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
楼主，最终答案是什么呢，如何导出差分电荷密度的数据？求解
赞一下 回复此楼
19楼2014-09-23 21:12:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
★
小木虫: 金币+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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 clown9111014 的主题更新
202/2返回列表上一页12
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 s3fFTmArrBt6 2026-09-14 6/300 2026-09-16 16:48 by i3KBRD19ERkG
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 5/250 2026-09-16 16:12 by i3KBRD19ERkG
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +6 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-16 15:24 by i3KBRD19ERkG
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +8 两三岁ss 2026-09-14 8/400 2026-09-16 09:46 by 水士口
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-16 08:51 by SjQobnC04j06
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-16 08:27 by SjQobnC04j06
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-16 07:12 by LH5NkK5Ud4bh
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-16 07:03 by LH5NkK5Ud4bh
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-16 06:11 by LH5NkK5Ud4bh
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-16 05:03 by LH5NkK5Ud4bh
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-16 03:35 by DgNGHc3h5tPl
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 19:59 by CgyNCDVNhGVg
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-15 19:23 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 23jxep3nCNZb 2026-09-14 4/200 2026-09-15 18:11 by CgyNCDVNhGVg
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:23 by CgyNCDVNhGVg
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-15 15:14 by CgyNCDVNhGVg
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +7 QUjhNVAcOSff 2026-09-13 9/450 2026-09-15 14:10 by VMDqgnlKimtZ
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[基金申请] 求教各位大神：2026教育部人文社科青年基金项目何时公示呀？ +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
信息提示
关闭
请填处理意见
关闭