24小时热门版块排行榜    

查看: 1716  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 297求调剂 +5 戏精丹丹丹 2026-03-17 5/250 2026-03-18 12:46 by Linda Hu
[考研] 299求调剂 +5 △小透明* 2026-03-17 5/250 2026-03-18 11:49 by 尽舜尧1
[考研] 303求调剂 +4 睿08 2026-03-17 6/300 2026-03-18 11:01 by Iveryant
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +3 1孙悟空 2026-03-17 3/150 2026-03-18 10:28 by macy2011
[考研] 材料专硕306英一数二 +9 z1z2z3879 2026-03-16 12/600 2026-03-18 08:53 by zhukairuo
[考研] 085601专硕,总分342求调剂,地区不限 +4 share_joy 2026-03-16 4/200 2026-03-17 21:32 by hmn_wj
[考研] 考研化学学硕调剂,一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +7 Ncdx123456 2026-03-13 8/400 2026-03-16 12:15 by karry wen
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 材料与化工 323 英一+数二+物化,一志愿:哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
信息提示
请填处理意见