24小时热门版块排行榜    

查看: 1497  |  回复: 15
当前主题已经存档。

chrinide

木虫 (正式写手)

把输入文件贴出来瞧瞧
未知必须借由更深的未知去获悉,隐晦必须借由更深的隐晦来明了
11楼2009-10-11 00:57:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lihb734

铁杆木虫 (职业作家)

站在计算化学入门的门槛上

★ ★
lei0736(金币+1,VIP+0):谢谢 10-11 08:55
gongxd325(金币+1,VIP+0): 10-11 17:14
建议楼主在两个测试中的路径行加 5D 7F进行比照一下。
看看是否一样?
测试结果可以和大家说一声。
前途光明,出路难觅!
12楼2009-10-11 08:34:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gongxd325

至尊木虫 (文坛精英)

这是输入
% chk=1
# B3LYP/GENECP OPT GFinput optcyc=200

Test1

0  1
Ag    0.000000    0.973585    0.000000
O    -0.035321   -1.043566    0.000000
N     0.003644   -2.322076    0.000000
C     0.042843   -3.525907    0.000000

C O N 0
CEP-121G
****
Ag 0
CEP-121G
****

Ag 0
CEP-121G

--Link1--
% chk=2
# B3LYP/CEP-121G OPT GFinput optcyc=200

Test2

0  1
Ag    0.000000    0.973585    0.000000
O    -0.035321   -1.043566    0.000000
N     0.003644   -2.322076    0.000000
C     0.042843   -3.525907    0.000000

测试了5D,对结果没有影响。
用LanL2DZ,两种输入方式的计算结果完全一样,能量都是-313.697668au。
但CEP-121G的结果就不一样,basis functions分别是67和70,能量是-281.819101和-177.65283au。
似乎是基组的问题?
一个人做点好事并不难难的是一辈子做好事!!!
13楼2009-10-11 17:20:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lihb734

铁杆木虫 (职业作家)

站在计算化学入门的门槛上


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by gongxd325 at 2009-10-11 17:20:
这是输入
% chk=1
# B3LYP/GENECP OPT GFinput optcyc=200

Test1

0  1
Ag    0.000000    0.973585    0.000000
O    -0.035321   -1.043566    0.000000
N     0.003644   -2.322076    0.000000
C  ...

从这看不出来你是否加入了5D,7F啊。
前途光明,出路难觅!
14楼2009-10-11 19:34:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天空空

木虫 (正式写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
gongxd325(金币+2,VIP+0):Thanks! 10-12 09:48
哈哈,不好意思,给弄错了!

[ Last edited by 天空空 on 2009-10-12 at 09:53 ]
好好学习,天天向上!嘿嘿
15楼2009-10-11 21:38:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gongxd325

至尊木虫 (文坛精英)

★ ★ ★
yjcmwgk(金币+3,VIP+0):问题解决 10-12 13:48
现在搞清楚了,如下输入,两种方式才等价:
# B3LYP/GENECP OPT 6D 10F

Test1

0  1
Ag    0.000000    0.973585    0.000000
O    -0.035321   -1.043566    0.000000
N     0.003644   -2.322076    0.000000
C     0.042843   -3.525907    0.000000

C O N 0
CEP-121G
****
Ag 0
CEP-121G
****

C O N 0
CEP-121G
Ag 0
CEP-121G

--Link1--
# B3LYP/CEP-121G OPT

Test2

0  1
Ag    0.000000    0.973585    0.000000
O    -0.035321   -1.043566    0.000000
N     0.003644   -2.322076    0.000000
C     0.042843   -3.525907    0.000000

用GENECP时CON也要读ECP,还要加关键词6D和10F,但LanL2DZ基组不需要这样。
一个人做点好事并不难难的是一辈子做好事!!!
16楼2009-10-12 09:43:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gongxd325 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见