24小时热门版块排行榜    

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

明-明-

金虫 (小有名气)

[求助] cif文件中的原子电荷分配 已有2人参与

我在数据库中找到了伊利石的cif文件,导入MS后,先点击calculate bonds,每个原子的电荷为零,然后没有进行手动分配电荷,构建超胞后,根据文献,采用UFF力场,Qeq算法来对其进行优化,然后就出现报错了。  
This periodic system bears a net charge (-64.000000e) and the
QEq calculation cannot be performed.
This periodic system bears a net charge (-64.000000e) and the
QEq calculation cannot be performed.
Unable to calculate energy and gradients.
Geometry optimization finished with an error condition.
查看报错后的晶体文件,每个原子的力场类型都分配了,但是charge还是为零,用compass力场就不会出现这种现象,优化后原子的电荷、力场均分配了。但查文献,一般是用UFF+Qeq来进行优化以及后续的吸附计算,请问这是什么原因?新人求助。。。是找来的cif文件构建晶体模型时都必须先手动分配原子电荷信息吗?
回复此楼
互相交流,互相学习
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

橙子学术

银虫 (小有名气)

你好,这个问题你解决了吗?我也是构建超晶胞,然后采用UFF力场,Qeq算法来对其进行优化,主要是模型电荷也不平衡,所以QE电荷平衡法会出错!
6楼2016-01-10 11:00:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 32 个回答

xie158795

金虫 (小有名气)

【答案】应助回帖

★
感谢参与,应助指数 +1
明-明-: 金币+1, ★有帮助 2015-11-19 16:22:22
你好!应该先初始化既先赋予电荷,再去优化!你有没有学过雅可比迭代,也要先赋予初始值的!同是想问问你的cif文件是在哪找到的???

发自小木虫Android客户端
愿同志同道和的人前行!
2楼2015-11-19 15:00:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

明-明-

金虫 (小有名气)

引用回帖:
2楼: Originally posted by xie158795 at 2015-11-19 15:00:15
你好!应该先初始化既先赋予电荷,再去优化!你有没有学过雅可比迭代,也要先赋予初始值的!同是想问问你的cif文件是在哪找到的???

我再问一下,赋予电荷是需要我查文献,找到各个原子的电荷信息后,直接在properties里面的charge中直接改吗?还是用某一task直接赋予电荷,看文章里(关于蒙脱石的)是说引用的电荷参数。 但优化时选择compass,以及forcefield assigned,正常运行,优化后电荷,力场都有了。  晶体文件是在这个网站上找的http://www.crystallography.net/index.php
另外,Findit软件中也可以找到,直接生成cif文件。我装的08版
互相交流,互相学习
3楼2015-11-19 16:21:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xie158795

金虫 (小有名气)

直接改听我师兄说,好像是不行的,应该用某种规折赋予,不然后期会模拟不下去

发自小木虫Android客户端
愿同志同道和的人前行!
4楼2015-11-20 08:18:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-15 11:14 by xSrQNSSjcgoX
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-15 09:35 by xSrQNSSjcgoX
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 SX5mJmPWDfOO 2026-11-13 3/150 2026-11-15 07:56 by xSrQNSSjcgoX
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 uKOSCocCb6Fm 2026-11-12 7/350 2026-11-15 02:48 by xSrQNSSjcgoX
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 uKOSCocCb6Fm 2026-11-12 8/400 2026-11-15 02:26 by xSrQNSSjcgoX
[文学芳草园] 厉娜 +4 张丹枫1991 2026-11-09 6/300 2026-11-15 00:48 by Equinoxhua
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +8 zSbvVX9JFl9J 2026-11-09 9/450 2026-11-13 11:11 by Equinoxhua
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 C6CRn66RYYXS 2026-11-10 6/300 2026-11-13 11:04 by Equinoxhua
[考博] 27博士申请 +5 白白白ly 2026-11-09 5/250 2026-11-13 05:01 by wlf1982
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 16:55 by zzy933268
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 cdqGD6ftHmbl 2026-11-11 4/200 2026-11-12 12:26 by ECLTXnJnRSXf
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 czAiolzR0IJC 2026-11-10 7/350 2026-11-12 11:09 by HajghKOMynmz
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 8/400 2026-11-12 11:00 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-12 10:38 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 YHiw1gou1Eda 2026-11-10 6/300 2026-11-12 10:08 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 FChKX8meu1ds 2026-11-09 11/550 2026-11-12 09:57 by HajghKOMynmz
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 FChKX8meu1ds 2026-11-09 10/500 2026-11-12 09:52 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 08:15 by HajghKOMynmz
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cdqGD6ftHmbl 2026-11-11 3/150 2026-11-12 07:53 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 gRND4NVAZ8VG 2026-11-11 3/150 2026-11-12 07:42 by HajghKOMynmz
信息提示
请填处理意见