24小时热门版块排行榜    

查看: 3597  |  回复: 28
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

starrywang

新虫 (正式写手)

[求助] Electron density difference 的sets of atoms 什么意思 已有1人参与

如题,我知道如何设置sets of atoms ,但是不知道这个是什么意思,算出来的结果我也看不懂,好像有两个单独的castep文件。

另外求教一下下面的图怎么画。该文计算了表面有无吸附物前后的电荷差,然后画出图形。请问,这个电荷是通过什么得到的,光是population的数据是画不出来的吧。
Electron density difference 的sets of atoms 什么意思
Electron density difference 的sets of atoms 什么意思-1
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★ ★ ★
franch: 金币+4, 谢谢回帖交流。。。 2014-10-06 23:25:29
差分电荷密度(二次差分,两片段间相互作用前后的电荷密度差)需要指定两个片段具体是那些原子。

差分电荷的定义是△ρ(AB)=ρ(AB)-ρ(A)-ρ(B),其中ρ(A)、ρ(B)分别都是维持体系AB中A和B的构型,但是体系中,A具体由哪些原子,B具体由哪些原子构成,这个就要指定,set就是指定这个。

如果我们知道了上述定义,那么其实手工计算差分电荷也是可以实现的。只是本人对CASTEP输出的check文件实在是破解不能。。。有个工具叫check2xsf似乎可以转化成cub或是其他格式,这样也许可以输出文本文件用origin作图。如果不限定作图方式,可以通过在MS下建立Slice调整大小后直接截图即可。
不一定挂在论坛,计算问题问题欢迎留言。
2楼2014-09-30 21:05:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★ ★
franch: 金币+3, 谢谢回帖交流。。。 2014-10-06 23:25:35
引用回帖:
3楼: Originally posted by starrywang at 2014-09-30 22:56:13
我知道了,sets of atoms定义的叫分子片段,剩余的部分肯能认为是另一个部分。计算electron density difference时,可以计算出这两个或多个部分组合和分离前后的电荷密度差,也就是我前面给的图。是这个意思吧。不 ...

对。

差分电荷主要描述的是两个片段间相互作用中电子态的变化,不包括核结构变化导致电子的弛豫,体现的是两片段相互作用前后电荷的聚集与发散,以及电荷的偏向。当然也可以和布居分析联系起来。
不一定挂在论坛,计算问题问题欢迎留言。
5楼2014-10-01 00:01:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★ ★ ★
franch: 金币+4, 谢谢回帖交流。。。 2014-10-06 23:25:39
引用回帖:
4楼: Originally posted by starrywang at 2014-09-30 23:04:30
另外还有两个小问题,1~构型比较立体,无法单看某个面,这样怎么可以很好的表现出小范围的体的电荷情况。2~是不是用population计算直接比较反应物和产物的原子电荷,得到的是反应前后的电荷移动;用electron den ...

1、你可尝试做差分电荷的等值面,类似于电荷密度或前线轨道的等值面那样,然后选择合适的角度截图。
2、electron density difference就可以直接大致判断共价、离子性的特征。布居分析是分析原子电荷的。估计电荷偏向除了布居分析也可以通过电负性均衡原理。

布居分析多种多样,没啥统一性。常规的布居分析操作从简单的例子来看,两原子成键,一部分电荷几乎维持原来的状态,一部分游离于两原子间,这些被共用的电子按照布局分析按一定人为形式划分归属到某个原子,有时候未必完全具有合理性,但容易定量化,对应着电荷密度分布。
不一定挂在论坛,计算问题问题欢迎留言。
6楼2014-10-01 00:22:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
8楼: Originally posted by dxcharlary at 2014-10-02 19:19:38
在化学里面有个词叫basis set翻译过来叫基组。这里面也是这个意思,指的是单原子的轨道函数。
楼主说的就是电荷密度的一些加加减减,这个现成的MS是做不了的。推荐vasp。处理一下CHGCAR就行了。

前辈,CASTEP能直接做,两个片段需指定一个片段叫"DensityDifference"即可,然后一次差分和二次差分都能做。
Dmol3能做一次差分,二次差分通过手工对grd文件也可完成。
不一定挂在论坛,计算问题问题欢迎留言。
9楼2014-10-02 20:42:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★
franch: 金币+2, 谢谢回帖交流。。。 2014-10-06 23:26:10
引用回帖:
10楼: Originally posted by dxcharlary at 2014-10-02 21:30:18
我不清楚CASTEP能对电荷密度进行操作,有可能是新版本功能。不过原来只能做单个原子的差分电荷密度。...

现在可以了。以前的话应该原则上应该可以通过check2xsf插件转化成其他格式后对格点文件做差,不过超麻烦。。
不一定挂在论坛,计算问题问题欢迎留言。
11楼2014-10-02 23:35:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★
franch: 金币+2, 谢谢回帖交流。。。 2014-10-06 23:26:22
引用回帖:
12楼: Originally posted by starrywang at 2014-10-03 01:10:40
弱问一次差分和二次差分分别指什么。
...

形变电荷和差分电荷
不一定挂在论坛,计算问题问题欢迎留言。
13楼2014-10-03 08:28:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★
franch: 金币+2, 谢谢回帖交流。。。 2014-10-06 23:26:32
引用回帖:
14楼: Originally posted by starrywang at 2014-10-03 10:07:38
这两个在castep中分别对应的是什么选项?
...

一次对应atomic densities
二次对应sets of atoms
不一定挂在论坛,计算问题问题欢迎留言。
15楼2014-10-03 10:43:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
16楼: Originally posted by starrywang at 2014-10-03 11:11:44
原来中文可以这么叫。感谢指点。最近做了几个电荷密度和差分的分布图,固体表面小分子吸附的。希望可以继续得到您指点。
...

提不上指点,大家多交流
不一定挂在论坛,计算问题问题欢迎留言。
17楼2014-10-03 11:44:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
18楼: Originally posted by 一川烟草 at 2014-10-03 16:57:48
请问怎么设置sets of atoms呢?因为选这项后好像还要set一下,但找不到地方…
...

楼上正解,指定片段必须叫DensityDifference
不一定挂在论坛,计算问题问题欢迎留言。
21楼2014-10-03 19:47:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★ ★
franch: 金币+3, 谢谢回帖交流。。。 2014-10-06 23:27:41
引用回帖:
23楼: Originally posted by 一川烟草 at 2014-10-03 20:53:01
灰常感谢!

还有个问题求教大侠,就是怎么把*.check转换成cube格式以备后续作图呢?...

木虫论坛有个叫check2xsf的插件,http://muchong.com/bbs/viewthread.php?tid=5406515&fpage=1&target=blank
你可以下载后(源码的话要编译),之后按照使用说明转换格式即可。对于*.check格式本身是加密的,我就不了解了。另外check文件本身是隐藏的(没记错的话)。
不一定挂在论坛,计算问题问题欢迎留言。
24楼2014-10-03 21:02:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 starrywang 的主题更新
信息提示
请填处理意见