24小时热门版块排行榜    

查看: 546  |  回复: 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的回帖

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的回帖
查看全部 5 个回答

hbnuwangshuling

铜虫 (小有名气)

有人用半经验算法计算的吗?能优化带电体系吗?
2楼2009-12-12 18:11:59
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 289求调剂 +5 步川酷紫123 2026-03-11 5/250 2026-03-15 00:45 by kruisytel
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[教师之家] 焦虑 +5 水冰月月野兔 2026-03-13 7/350 2026-03-14 15:14 by 农药害害
[考研] 290求调剂 +4 @将就将就看 2026-03-10 8/400 2026-03-14 14:23 by 千千运气
[考研] 331求调剂(0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 一志愿华中农业大学071010,总分三百二,求调剂 +3 困困困困坤坤 2026-03-10 3/150 2026-03-14 00:35 by JourneyLucky
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 341求调剂 +4 番茄头--- 2026-03-10 4/200 2026-03-13 23:12 by JourneyLucky
[考研] 材料专硕288分求调剂 一志愿211 +4 在家想你 2026-03-11 4/200 2026-03-13 22:49 by JourneyLucky
[考研] 332求调剂 +3 zjy101327 2026-03-11 6/300 2026-03-13 22:48 by JourneyLucky
[考研] 一志愿西南交大,材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 【0856】化学工程(085602)313 分,本科学科评估A类院校化学工程与工艺,诚求调剂 +7 小刘快快上岸 2026-03-11 7/350 2026-03-13 16:06 by ruiyingmiao
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 一志愿河海大学085900土木水利专硕279求调剂不挑专业 +4 SunWwWwWw 2026-03-10 8/400 2026-03-13 02:23 by SunWwWwWw
[考博] 2026年博士申请 +3 QwQwQW10 2026-03-11 3/150 2026-03-12 17:58 by gxch43
[考研] 085600 材料与化工 295 求调剂 +10 dream…… 2026-03-10 12/600 2026-03-12 13:46 by dream……
[考博] 26申博求助 +3 跳跃饼干 2026-03-10 4/200 2026-03-10 21:15 by Tntcnn
信息提示
请填处理意见