24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1757  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +6 熊二想上岸 2026-04-06 6/300 2026-04-07 08:14 by houyaoxu
[考研] 071000生物学,一志愿深圳大学296分,求调剂 +5 TIckLw 2026-04-06 5/250 2026-04-07 07:56 by dongzh2009
[考研] 材料科学与工程320求调剂,080500 +7 黄瓜味薯片 2026-04-06 7/350 2026-04-06 22:38 by qlm5820
[考研] 304求调剂 +4 luoye0105 2026-04-05 4/200 2026-04-06 21:05 by 木子君1218
[考研] 338求调剂 +4 我想上岸ii 2026-04-05 4/200 2026-04-06 21:04 by 木子君1218
[考研] 生物学求调剂 一志愿沪9,326分 +6 刘墨墨 2026-04-06 6/300 2026-04-06 19:36 by lijunpoly
[考研] 本科211,293分请求调剂 +8 莲菜就是藕吧 2026-04-03 9/450 2026-04-05 19:12 by 蓝云思雨
[考研] 材料调剂 +13 一样YWY 2026-04-03 14/700 2026-04-05 18:20 by 蓝云思雨
[考研] 085600调剂 +9 东照照照 2026-04-04 9/450 2026-04-05 13:44 by ujn_zhuj
[考研] 工科08专硕机械275求调剂 +3 AaAa7420 2026-04-02 3/150 2026-04-05 13:26 by jp9609
[考研] 一志愿江南大学085501机械工程专硕326分,本科佳木斯大学 +5 顾若浮生 2026-04-03 9/450 2026-04-05 09:57 by 1753564080
[考研] 341求调剂 +3 学无止境,冲 2026-04-05 3/150 2026-04-05 09:40 by lbsjt
[考研] 考研调剂 +11 小sun要好运 2026-04-04 11/550 2026-04-05 08:02 by qlm5820
[考研] 化工求调剂 +11 荔香芝士椰奶 2026-04-03 11/550 2026-04-03 22:06 by 啵啵啵0119
[考研] 工科 267求调剂 +5 wanwan00 2026-04-02 7/350 2026-04-03 14:14 by zhangdingwa
[考研] 283求调剂 +3 jiouuu 2026-04-03 4/200 2026-04-03 13:28 by jiouuu
[考研] 重庆大学材料与化工085600,初试370+,求求调剂建议 +8 shzhou_ 2026-04-01 9/450 2026-04-03 09:31 by 蓝云思雨
[考研] 材料调剂 +11 一样YWY 2026-03-31 11/550 2026-04-01 22:25 by zhouyuwinner
[考研] 296求调剂 +4 汪!?! 2026-03-31 7/350 2026-04-01 22:04 by 客尔美德
[考研] 0855机械初试280求调剂 +3 kazenotori 2026-03-31 3/150 2026-04-01 10:08 by JourneyLucky
信息提示
请填处理意见