24小时热门版块排行榜    

查看: 547  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 326求调剂 +3 上岸的小葡 2026-03-15 4/200 2026-03-15 18:50 by 无际的草原
[考研] 化学工程321分求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-15 18:49 by a不易
[考研] 085600材料与化工 326 求调剂 +5 热爱生活ing 2026-03-09 5/250 2026-03-14 02:39 by JourneyLucky
[考研] 化学工程321分求调剂(南京工业,浙江工业) +3 大米饭! 2026-03-09 4/200 2026-03-14 02:34 by JourneyLucky
[考研] 云南财经大学信息学院计算机学硕专硕学位点 +3 zjptai 2026-03-10 5/250 2026-03-14 01:23 by 飞行琦
[考研] 307求调剂 +7 超级伊昂大王 2026-03-10 7/350 2026-03-14 00:49 by JourneyLucky
[考研] 一志愿华中农业大学071010,总分三百二,求调剂 +3 困困困困坤坤 2026-03-10 3/150 2026-03-14 00:35 by JourneyLucky
[考研] 311求调剂 +5 牛乳糖的卡卡 2026-03-10 5/250 2026-03-14 00:05 by JourneyLucky
[考研] 材料专硕288分求调剂 一志愿211 +4 在家想你 2026-03-11 4/200 2026-03-13 22:49 by JourneyLucky
[考研] 材料与化工304求B区调剂 +5 邱gl 2026-03-11 6/300 2026-03-13 22:37 by JourneyLucky
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 5/250 2026-03-13 22:06 by 星空星月
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
[考研] 282分材料专业求调剂院校 +18 枫桥ZL 2026-03-09 25/1250 2026-03-13 10:47 by 白夜悠长
[考研] 一志愿华中师范071000,325求调剂 +5 RuitingC 2026-03-12 5/250 2026-03-13 10:43 by hyswxzs
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考研] 一志愿河海大学085900土木水利专硕279求调剂不挑专业 +4 SunWwWwWw 2026-03-10 8/400 2026-03-13 02:23 by SunWwWwWw
[考研] 一志愿山东大学,总分327,英语二79,有论文,有竞赛,已过四六级 +3 木木目目1 2026-03-09 3/150 2026-03-09 19:52 by yuningshan
信息提示
请填处理意见