24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1745  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 279求调剂 +8 j的立方 2026-03-29 8/400 2026-03-30 17:21 by 无际的草原
[考研] 26考研-291分-厦门大学(085601)-柔性电子学院材料工程专业求调剂 +4 min3 2026-03-24 5/250 2026-03-30 16:41 by 梨子~
[考研] 材料工程专硕求调剂 +6 hyl3153942 2026-03-29 6/300 2026-03-30 12:04 by 1018329917
[考研] 南京大学化学调剂 +10 景随风 2026-03-29 15/750 2026-03-30 11:21 by limeifeng
[考研] 296求调剂 +10 彼岸t 2026-03-29 10/500 2026-03-30 10:50 by 探123
[考研] 285求调剂 +5 AZMK 2026-03-29 8/400 2026-03-30 10:41 by AZMK
[考研] 085404求调剂,总分309,本科经历较为丰富 +6 来财aa 2026-03-25 6/300 2026-03-30 09:48 by 青海小西牛
[考研] 求调剂 +4 QiMing7 2026-03-25 5/250 2026-03-29 21:10 by 唐沐儿
[考研] 289求调剂 +5 BrightLL 2026-03-29 5/250 2026-03-29 17:24 by zhyzzh
[考研] 340求调剂 +6 Amber00 2026-03-26 6/300 2026-03-29 12:06 by 无际的草原
[考研] 0856求调剂 +13 zhn03 2026-03-25 14/700 2026-03-29 08:13 by fmesaito
[考研] 081200-11408-276学硕求调剂 +6 崔wj 2026-03-26 6/300 2026-03-29 01:11 by hanserlol
[考研] 数一英一271专硕(085401)求调剂,可跨 +7 前行必有光 2026-03-28 8/400 2026-03-28 23:22 by 小木虫tim
[考研] 一志愿太原理工安全工程300分,求调剂 +5 0857求调剂. 2026-03-24 6/300 2026-03-28 22:04 by zhq0425
[考研] 11408软件工程求调剂 +3 Qiu学ing 2026-03-28 3/150 2026-03-28 21:50 by zhq0425
[考研] 083000学硕274求调剂 +8 Li李鱼 2026-03-26 8/400 2026-03-28 20:33 by 加油向未来啊
[考研] 070300求调剂306分 +4 26要上岸 2026-03-27 4/200 2026-03-28 13:06 by 唐沐儿
[考研] 材料277求调剂 +7 min3 2026-03-24 7/350 2026-03-28 11:39 by xuxiang
[考研] 309求调剂 +4 gajsj 2026-03-25 5/250 2026-03-26 00:27 by Dyhoer
[考研] 【2026考研调剂】制药工程 284分 求相关专业调剂名额 +4 袁奂奂 2026-03-25 8/400 2026-03-25 14:32 by lbsjt
信息提示
请填处理意见