24小时热门版块排行榜    

查看: 534  |  回复: 4
当前主题已经存档。

hbnuwangshuling

铜虫 (小有名气)

[交流] 【求助】am1能就算带电体系吗?

我计算聚乙炔体系时 用AM1方法计算,优化出的分子体系结构。但总是出错,不知道为什么?am1能就算带电体系吗?
输入文件是:AM1  opt
Charge =  1 Multiplicity = 2
C                    -0.05022  -11.71397   0.
C                    -0.65811  -10.51136   0.
C                     0.05022  -9.24884   0.
C                    -0.56499  -8.03702   0.
C                     0.13286  -6.77996   0.
C                    -0.48743  -5.56625   0.
C                     0.20653  -4.31113   0.
C                    -0.4156   -3.0964    0.
C                     0.27682  -1.84216   0.
C                    -0.34599  -0.62699   0.
C                     0.34599   0.62699   0.
C                    -0.27682   1.84216   0.
C                     0.4156    3.0964    0.
C                    -0.20653   4.31113   0.
C                     0.48743   5.56625   0.
C                    -0.13286   6.77996   0.
C                     0.56499   8.03702   0.
C                    -0.05022   9.24884   0.
C                     0.65811   10.51136   0.
C                     0.05022   11.71397   0.
H                     1.03295  -11.80111   0.
H                    -0.6169   -12.63887   0.
H                    -1.7475   -10.465    0.
H                     1.14113  -9.2892    0.
H                    -1.65644  -8.00246   0.
H                     1.22433  -6.81211   0.
H                    -1.57909  -5.5357    0.
H                     1.2982   -4.34075   0.
H                    -1.50734  -3.06743   0.
H                     1.36857  -1.87076   0.
H                    -1.43776  -0.59856   0.
H                     1.43776   0.59856   0.
H                    -1.36857   1.87076   0.
H                     1.50734   3.06743   0.
H                    -1.2982    4.34075   0.
H                     1.57909   5.5357    0.
H                    -1.22433   6.81211   0.
H                     1.65644   8.00246   0.
H                    -1.14113   9.2892    0.
H                     1.7475    10.465    0.
H                    -1.03295   11.80111   0.
H                     0.6169    12.63887   0.
报错信息是  """""""""""""UNABLE TO ACHIEVE SELF-CONSISTENCE


          DELTAE=  -0.4907E+05     DELTAP=   0.3673E+02



Error termination via Lnk1e in /export/home/gauss/g03/l402.exe at Wed Dec  9 14:27:05 2009.
Job cpu time:  0 days  0 hours  5 minutes  3.1 seconds.
File lengths (MBytes):  RWF=     13 Int=      0 D2E=      0 Chk=      7 Scr=      1

[ Last edited by yjcmwgk on 2009-12-12 at 13:37 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hbnuwangshuling

铜虫 (小有名气)

有人用半经验算法计算的吗?能优化带电体系吗?
2楼2009-12-12 18:11:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

★ ★ ★
erylingjet(金币+3,VIP+0):感谢交流 12-24 15:18
引用回帖:
Originally posted by hbnuwangshuling at 2009-12-12 18:11:
有人用半经验算法计算的吗?能优化带电体系吗?

我给你运行了一下
23秒后正常结束
输出文件看附件吧
搞不懂你为什么会遇到scf不收敛的问题
3楼2009-12-13 13:32:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hbnuwangshuling

铜虫 (小有名气)

谢谢!但是仍然有问题!

我对比了一下,我的输出和您的输出结果,对最初轨道对称性分析时,我的结果是A',而您的输出结果中是 There are    41 symmetry adapted basis functions of AG  symmetry.
There are    10 symmetry adapted basis functions of BG  symmetry.
There are    10 symmetry adapted basis functions of AU  symmetry.
There are    41 symmetry adapted basis functions of BU  symmetry.
还有就是点群对称性我的结果是:Stoichiometry    C20H22(1+,2)
Framework group  CS[SG(C20H22)]
Deg. of freedom    81
Full point group                 CS
Largest Abelian subgroup         CS      NOp   2
Largest concise Abelian subgroup C1      NOp   1而您的结果是
Deg. of freedom    41
Full point group                 C2H     NOp   4
Largest Abelian subgroup         C2H     NOp   4
Largest concise Abelian subgroup C2      NOp   2
咱们两个的输入文件是一样的,而输出就不同了,这是什么原因呢?
4楼2009-12-14 09:23:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
erylingjet(金币+2,VIP+0):感谢交流 12-24 15:18
引用回帖:
Originally posted by hbnuwangshuling at 2009-12-14 09:23:
我对比了一下,我的输出和您的输出结果,对最初轨道对称性分析时,我的结果是A',而您的输出结果中是 There are    41 symmetry adapted basis functions of AG  symmetry.
There are    10 symmetry adapted ba ...

呵呵,我们用的当然不是同一个文件。你的输入文件被我调整成了新对称性了。
在out文件的开头的分子说明部分,你是很容易看出这一点的,所以……我猜你肯定没好好看我给你的out文件,呵呵

[ Last edited by yjcmwgk on 2009-12-14 at 10:06 ]
5楼2009-12-14 10:01:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hbnuwangshuling 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见