| 查看: 6597 | 回复: 12 | ||||||||||||
| 本帖产生 1 个 1ST强帖 ，点击这里进行查看 | ||||||||||||
xiao72379木虫 (正式写手)
宅男甲
|
[求助]
求助vasp做过电荷密度差分的文献以及答疑，谢谢！
|
|||||||||||
|
Sample Text 大家好 以前我分析成键都用DOS+ELF+Bader 最近的一个工作 我师姐建议做下电荷差分密度（charge density difference） 我没做过 请大家帮忙 谢了先！ 1）求用vasp做过电荷密度差分的文献 这个用ISI实在不好查 老板催工作我时间又不多- - 如果同学手上有现成的文献给我参考 感激不尽 2）查了下论坛原来的帖子 找到一个我认为比较好的参考贴[http://muchong.com/bbs/viewthread.php?tid=2129228]（我不知道如何给他金币，版主看到了教下我，呵呵），渔火江枫 的 先引用下： VASP与电荷密度差 一般分为：差分电荷密度图(deformation charge density)，二次差分电荷密度图(difference charge density) 两者的区别有很多种说法，其中一种定义认为： 差分电荷定义为成键后的电荷密度与对应的点的原子电荷密度之差。通过差分电荷密度的计算和分析，可以清楚地得到在成键和成键电子耦合过程中的电荷移动以及成键极化方向等性质。 “二次”是指同一个体系化学成分或者几何构型改变之后电荷的重新分布。 如何做对电荷密度做差，目前有两个版本的公式：  人认为公式(1)和差分电荷密度的定义比较吻合,但是这种方法会产生很多干扰信息，而文献中常用的是公式(2)。 如何用VASP做电荷密度差图？ ① 对优化后的结构做静态自洽计算 要注意的是，AB、A和B要分别放在同样大小的空间格子中并保证A和B与AB中相应坐标不变，计算时也要保证三次自洽计算所采用的FFT mesh 一致（NGXF,NGYF,NGZF）。 INCAR中几个注意的参数： IBRION = -1;NSW = 0; NGXF,NGYF,NGZF ② 做差 将第一步得到的三个CHGCAR文件依次命名为CHGCAR1、CHGCAR2、CHGCAR3，采用脚本程序按照公式②对三个文件做差，得到CHGdiff文件。 ③ 视图 3D显示可以将得到的CHGdiff文件添加vasp后缀(CHGdiff.vasp)，用VESTA读取，另外VESTA也可以做2D切面，缺点是没有标度。 2D的显示也可以用lev00去做切面，然后用origin作图，得到的2D面电荷密度差带有标度便于分析(具体做法可参阅lev00说明书)。 好了 我的问题就是 比如我的体系 基底是A材料 然后呢是多个B原子吸附或者是嵌入（这里你视为化学键合的作用就行），形成了A（B）n,对上面的公式2应用到我的体系如何处理？比如说步骤中的1,2步 B做静态自洽计算的时候 是n个B一起算（我理解是扣除A的所有坐标得n个B的坐标，不知对否）还是一个一个B分开算，如果是分开，那脚本做CHGdiff如何做呢？如果都不是，请教下我该如何做，谢谢！ [ Last edited by xiao72379 on 2011-5-25 at 10:32 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 材料计算模拟实用技巧 | 第一性原理 | 第一原理资料汇编 | VASP |
| vasp计算资料 | 第一性原理 | VASP相关 | 学习资料 |
| 第一性原理计算-tg |
» 猜你喜欢
环境工程专硕342求调剂
已经有9人回复
-
274求调剂
已经有9人回复
-
求调剂
已经有4人回复
-
085700资源与环境专硕290求调剂
已经有9人回复
-
一志愿合肥工业大学085701环境工程279求调剂
已经有10人回复
-
085700资源与环境，专硕290求调剂
已经有15人回复
-
材料与化工求调剂
已经有15人回复
-
求调剂
已经有16人回复
-
求070300化学调剂
已经有5人回复
-
化学工程 286分求调剂
已经有17人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 差分电荷密度图
已经有17人回复
- VASP结果分析求助
已经有6人回复
- vasp可以计算密立根电荷转移吗
已经有4人回复
- 求vasp输出文件的分析
已经有5人回复
- 关于用Multiwfn软件做电荷密度差分图遇到的一个问题
已经有3人回复
- 送红花N朵！MS怎么得出差分电荷密度图！
已经有5人回复
- 修改RWIGS 维格纳半径之后，差分电荷密度会不会变化
已经有4人回复
- VESTA做的差分电荷密度图分析
已经有20人回复
- 再次求问，差分电荷密度图的定性分析
已经有17人回复
- 【求教】电子密度，差分密度，原子电荷这三个之间有什么关系？
已经有7人回复
- vasp 那位高手做过四氧化三铁的计算，能否指点一二啊
已经有12人回复
- vasp如何计算差分电荷密度
已经有7人回复
- 差分电荷密度的计算
已经有11人回复
- 如何做电子密度差分图
已经有16人回复
- 在vaspview中画电荷密度等高线
已经有12人回复
- 如何让差分电荷密度图的等高线光滑呢？
已经有5人回复
- 【求助】做差分电荷密度
已经有5人回复
- 【求助】原子的电荷密度分布怎么求？
已经有3人回复
- 【求助】如何得到体系原子不同能级间的电荷转移？
已经有6人回复
- 【求助】电荷密度
已经有14人回复
2楼2011-05-25 19:46:56
xiao72379
木虫 (正式写手)
宅男甲
- 1ST强帖: 2
- 应助: 10 (幼儿园)
- 金币: 4314.3
- 散金: 303
- 红花: 5
- 帖子: 454
- 在线: 580.7小时
- 虫号: 638344
- 注册: 2008-10-27
- 性别: GG
- 专业: 半导体材料
3楼2011-05-25 19:58:56
4楼2011-05-26 05:43:58
漂泊四方
金虫 (小有名气)
- 1ST强帖: 2
- 应助: 4 (幼儿园)
- 金币: 1321.2
- 散金: 50
- 红花: 6
- 帖子: 198
- 在线: 115.2小时
- 虫号: 897145
- 注册: 2009-11-08
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
★ ★
xiao72379(金币+10): 感谢回帖 2011-05-26 09:11:34
youzhizhe(金币+2, 1ST强帖+1): 谢谢交流。 2011-05-26 19:27:15
xiao72379(金币+10): 感谢回帖 2011-05-26 09:11:34
youzhizhe(金币+2, 1ST强帖+1): 谢谢交流。 2011-05-26 19:27:15
|
好了 我的问题就是 比如我的体系 基底是A材料 然后呢是多个B原子吸附或者是嵌入（这里你视为化学键合的作用就行），形成了A（B）n,对上面的公式2应用到我的体系如何处理？比如说步骤中的1,2 步 B做静态自洽计算的时候 是n个B一起算（我理解是扣除A的所有坐标得n个B的坐标，不知对否）还是一个一个B分开算，如果是分开，那脚本做CHGdiff如何做呢？如果都不是，请教下我该如何做，谢谢！ 方案一里面：把B当做一个整体，主要考虑的是A整体和B的作用。 方案二里面：实际上考虑了A和B的作用，还考虑了B与B之间的作用。 所以重点看你关注什么，衬底与吸附原子之间的作用还是吸附原子之间的作用。两种作用之间的强弱关系？看你的意思。如果B形成了团簇，方案一可能是你关注的，如果B比较离散，作用力较弱，方案一也可以。如果BB之间作用力介于他们之间，和衬底的作用类似，分析起来就比较吃力了。 操作起来其实还不是太难，记得p4vasp好像可以做这个。 |
5楼2011-05-26 08:43:26
xiao72379
木虫 (正式写手)
宅男甲
- 1ST强帖: 2
- 应助: 10 (幼儿园)
- 金币: 4314.3
- 散金: 303
- 红花: 5
- 帖子: 454
- 在线: 580.7小时
- 虫号: 638344
- 注册: 2008-10-27
- 性别: GG
- 专业: 半导体材料
|
我的情况是：B在体系A表面不成团簇，B与A体系的作用从bader计算看是有一个|e|左右的电子从B转移到体系A（B是2价主族元素） DOS上看B的S电子与体系A的峰有杂化现象，呈现一定的共价成分，不过B的S电子在费米能级附近有明显裂分现象，一个峰出现在价带，一个峰出现在导带，还是很好的印证了bader电荷分析的结果。所以我初步判断了B与A体系的键合作用是主离子性辅共价性的混合键性。做电荷差分是进一步支持下我的论点。我在计算过程中除了考虑一个B与A的作用，也做了多个B与A的作用。当然一个B与A的作用做电荷差分现在看来我在步骤上已经是清楚了，但当多个的时候，情况应该是A与B是较强的化学作用，BB之间作用较弱，从ELF做的图看，没有化学键，顶多是静电排斥，我主要关心的是B与A的键合情况是不是我上面提到的情况，这个差分做出来也是为了多一个证据说明。写的有点乱，总结下就是：B不成团簇，B较分散，但不同B的化学环境是不同的，因为他们嵌入的位置不一样，不过成键性质是一样的（比如dos+ELF的结果显示定性是类似的，不过是定量上的不一样）。我确定B与A之间有化学键合，我希望做差分进一步研究成键情况，当然结果能匹配我上面提到的 我已经做过的一些成键分析结果是最好的了。如果是B一起做个差分（考虑到相同键性），我还能分析到每个B与A的电荷转移情况？还是该一个一个来？如你做过，能不能给个细致点的操作步骤，如何处理vasp文件进行计算得到差分结果？ |
6楼2011-05-26 09:51:12
xiao72379
木虫 (正式写手)
宅男甲
- 1ST强帖: 2
- 应助: 10 (幼儿园)
- 金币: 4314.3
- 散金: 303
- 红花: 5
- 帖子: 454
- 在线: 580.7小时
- 虫号: 638344
- 注册: 2008-10-27
- 性别: GG
- 专业: 半导体材料
7楼2011-05-26 09:56:52
8楼2011-05-26 17:34:51
漂泊四方
金虫 (小有名气)
- 1ST强帖: 2
- 应助: 4 (幼儿园)
- 金币: 1321.2
- 散金: 50
- 红花: 6
- 帖子: 198
- 在线: 115.2小时
- 虫号: 897145
- 注册: 2009-11-08
- 性别: GG
- 专业: 理论和计算化学
9楼2011-05-26 22:01:44
gloomy2004
金虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 1461.1
- 散金: 64
- 红花: 32
- 沙发: 1
- 帖子: 1177
- 在线: 167.5小时
- 虫号: 49848
- 注册: 2004-07-10
10楼2011-07-20 03:37:07
