24小时热门版块排行榜    

查看: 2052  |  回复: 7

Nibelungen

新虫 (初入文坛)

[求助] 如何在GAMESS输入文件中指定某个原子带电?

计算一个铜离子和一个蛋白质结合后的能量,其中铜离子为+2价,请问要如何在GAMESS的输入文件中指定这种带电的特性?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

★ ★ ★
gmy1990: 金币+3 2012-07-03 13:10:43
1,正确方法:由计算结果决定。在保证物理正确的前提下,可以做一些设置来获得所需的价态。

2,强行设置Cu++。在$DATA部分使用下面的输入

...(其它原子的坐标和基组)

Cu  -27.00    (x, y, z坐标)
(Cu的基组。不要用基组名称,因为程序不知道27个负电荷代表什么元素)

pcharge  29.00    (同样的x, y, z坐标)

这样就定义了x,y,z位置有29个正电荷和27个电子,其中电子用Cu的某种基组来描述。不过这个计算没有什么物理意义。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
2楼2012-07-03 11:28:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Nibelungen

新虫 (初入文坛)

谢谢!

不过在楼上的回答中我还有两处不太明白:

1. “由计算结果决定”的意思是这样的么:即使不在GAMESS的输入文件中指定某个离子的价态,程序也能根据计算的过程中也能算出此离子应有的价态?

2. 我在GAMESS的输入参数中看到在contrl group里有ICHARG这个参数,是否可以通过设置这一参数达到“在保证物理正确的前提下,可以做一些设置来获得所需的价态”这一目标?

我是CS出身的,对于量化所需的化学和物理知识虽然也在学,但是感觉距离初步理解量化还有好长的路要走,问的问题也许很肤浅,见笑了~
3楼2012-07-03 22:34:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ldzh0531

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
Nibelungen: 金币+7, ★★★很有帮助 2012-07-06 15:00:37
分碎片优化,分别指定电荷和多重度。(本人没这样做过,仅供参考)
4楼2012-07-04 21:37:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

引用回帖:
3楼: Originally posted by Nibelungen at 2012-07-03 22:34:37
谢谢!

不过在楼上的回答中我还有两处不太明白:

1. “由计算结果决定”的意思是这样的么:即使不在GAMESS的输入文件中指定某个离子的价态,程序也能根据计算的过程中也能算出此离子应有的价态?

2. 我在G ...

1,是

2,ICHARG是蛋白质+Cu的总电荷。Cu的电荷未必等于ICHARG
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
5楼2012-07-05 01:18:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

引用回帖:
4楼: Originally posted by ldzh0531 at 2012-07-04 21:37:26
分碎片优化,分别指定电荷和多重度。(本人没这样做过,仅供参考)

FMO方法都是这么做的,最典型的就是ADF程序。但是这只是提供一个参照系,整个体系中Cu的电荷几乎不受分子片电荷的影响。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
6楼2012-07-05 01:20:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuhuisfp

金虫 (小有名气)

★ ★
gmy1990: 金币+2 2012-07-15 20:22:14
引用回帖:
4楼: Originally posted by ldzh0531 at 2012-07-04 21:37:26
分碎片优化,分别指定电荷和多重度。(本人没这样做过,仅供参考)

进行LMOEDA能量分解的时候可以设定
7楼2012-07-15 20:02:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

童生

银虫 (小有名气)

可以用LMOEDA算,语句如下(其中n1是蛋白质的原子数):
$LMOEDA matom(1)=1 n1 mcharg(1)=2 0 mmult(1)=2 1 supbas=.t. $end
ANYONEWHOISNOTSHOCKEDBYQUANTUMTHEORYHASNOTUNDERSTOODIT.--NIELSBOHR(1885-1962)
8楼2013-07-26 16:21:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Nibelungen 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:34 by ZvyPK8n6Nfki
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:18 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-28 23:19 by b8fCuTHqckEt
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:45 by ez6fGg9abYaj
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +3 tianxiaoxian 2026-09-28 6/300 2026-09-28 22:44 by tianxiaoxian
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:28 by ez6fGg9abYaj
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-28 22:21 by ez6fGg9abYaj
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 17:30 by ZYXdhzDAy9ZX
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +7 zju2000 2026-09-26 8/400 2026-09-28 16:42 by beefly
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 15:32 by ZYXdhzDAy9ZX
信息提示
请填处理意见