Gaussian几何优化后的键长不合理 - 量子化学 - Gaussian - 第4页小木虫论坛-学术科研互动平台

一个人做点好事并不难难的是一辈子做好事!!!

31楼2011-07-23 14:35:55

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

jianqin

金虫 (正式写手)

应助: 7 (幼儿园)
金币: 3729.9
红花: 2
帖子: 627
在线: 394小时
虫号: 863555
注册: 2009-10-05
专业: 质谱分析

引用回帖:

Originally posted by gongxd325 at 2011-07-23 14:35:55:
设置内存时别超过总内存的80%，否则会很慢

这是win7下算的，机子3GB内存，intel Core2 duo CPU

32楼2011-07-23 14:39:13

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

jianqin

金虫 (正式写手)

应助: 7 (幼儿园)
金币: 3729.9
红花: 2
帖子: 627
在线: 394小时
虫号: 863555
注册: 2009-10-05
专业: 质谱分析

引用回帖:

Originally posted by fatpig8832 at 2011-07-18 17:20:13:
还有一种可能性，这东西说不定是三重态，不知有没有算过，看看和单重态哪个才是ground state

几何优化未完成，计算就结束了。我把最后几行出错的粘贴如下：
Petite list used in FoFCou.
Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on          energy=1.00D-06.
No special actions if energy rises.
Integral accuracy reduced to 1.0D-05 until final iterations.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UHF) =  -1289.08035487    A.U. after  129 cycles
         Convg  = 0.2150D-04          -V/T =  2.0010
= 0.0000 = 0.0000 = 1.0000 = 2.1029 S= 1.0339
= 0.000000000000E+00
Annihilation of the first spin contaminant:
S**2 before annihilation    2.1029, after    2.0035
Convergence failure -- run terminated.
Error termination via Lnk1e in /home/kss/g09/l502.exe at Sat Jul 23 17:29:37 2011.

33楼2011-07-23 18:10:32

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT

QC强帖: 5
应助: 99 (初中生)
金币: 6961.8
散金: 10406
红花: 49
帖子: 1761
在线: 763.2小时
虫号: 880392
注册: 2009-10-22
性别: GG
专业: 理论和计算化学

★ ★
zhou2009(金币+2): 2011-08-06 15:21:17

>>>>>>>>>> Convergence criterion not met.
SCF Done: E(UHF) = -1289.08035487 A.U. after 129 cycles
以上表明没能收敛。。参照置顶的本版FAQ，给出了很多解决方案，根据你的具体情况，有选择地修正。

一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭，其实这是修养不足学识浅薄的一种体现，也是可恶的偏执和愚蠢的自以

34楼2011-07-23 20:03:56

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)

应助: 16 (小学生)
金币: 7420.2
散金: 31
红花: 28
帖子: 1438
在线: 610.9小时
虫号: 770403
注册: 2009-05-14

引用回帖:

Originally posted by jianqin at 2011-07-23 18:10:32:
几何优化未完成，计算就结束了。我把最后几行出错的粘贴如下：
Petite list used in FoFCou.
Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested convergence on MAX d ...

没关系，你就对比下 -1289.08035487 和单重态的能量哪个高就行了...

35楼2011-07-24 12:10:11

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

jianqin

金虫 (正式写手)

应助: 7 (幼儿园)
金币: 3729.9
红花: 2
帖子: 627
在线: 394小时
虫号: 863555
注册: 2009-10-05
专业: 质谱分析

引用回帖:

Originally posted by fatpig8832 at 2011-07-24 12:10:11:
没关系，你就对比下 -1289.08035487 和单重态的能量哪个高就行了...

谢谢！虽然我不知道这个E(uhf)=-1289.08035487代表的物理意义，但是这个值应该比单重态的小。
单重态，rhf 6-31g能量结果：
Zero-point correction=                         0.358864 (Hartree/Particle)
Thermal correction to Energy=                   0.379415
Thermal correction to Enthalpy=                0.380359
Thermal correction to Gibbs Free Energy=       0.305788
Sum of electronic and zero-point Energies=       -1288.405788
Sum of electronic and thermal Energies=          -1288.385237
Sum of electronic and thermal Enthalpies=          -1288.384293
Sum of electronic and thermal Free Energies=       -1288.458864

36楼2011-07-25 08:42:08

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

B4-449

铁杆木虫 (著名写手)

应助: 62 (初中生)
贵宾: 0.218
金币: 7842.6
散金: 330
红花: 37
沙发: 2
帖子: 1675
在线: 787小时
虫号: 1140943
注册: 2010-11-07
性别: GG
专业: 生物力学与生物流变学

好复杂，正常情况计算机都是在缺省条件设置的，和事实不符

37楼2011-07-25 10:45:11

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)

应助: 16 (小学生)
金币: 7420.2
散金: 31
红花: 28
帖子: 1438
在线: 610.9小时
虫号: 770403
注册: 2009-05-14

★ ★
zhou2009(金币+2): 2011-08-06 15:19:21

引用回帖:

Originally posted by jianqin at 2011-07-25 08:42:08:
谢谢！虽然我不知道这个E(uhf)=-1289.08035487代表的物理意义，但是这个值应该比单重态的小。
单重态，rhf 6-31g能量结果：
Zero-point correction= 0.358864 (Hartree/Particle) ...

这个E(UHF)是指electronic energy，每优化一步，就会给出一个E(UHF)。对大部分寻找minimum的优化来说，每一步的E(UHF)总体趋势上是一个降低的过程，

在你的单重态输出中，electronic energy应该是
=-1288.405788-0.358864
=-1288.385237-0.379415
=...
=...
=-1288.764652

你可以再搜一下三重态的输出文件的其他E(UHF)，如果都在-1289.08035487的附近，都比单重态能量低的话，那么基本可以断定，它的ground state是三重态(也可以试试五重，我想机会不大)。然后，你现在做的单重态的这些东西就没有太大意义了...

» 本帖已获得的红花（最新10朵）

jianqin

38楼2011-07-25 13:15:03

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

jianqin

金虫 (正式写手)

应助: 7 (幼儿园)
金币: 3729.9
红花: 2
帖子: 627
在线: 394小时
虫号: 863555
注册: 2009-10-05
专业: 质谱分析

送鲜花一朵

引用回帖:

38楼: Originally posted by fatpig8832 at 2011-07-25 13:15:03:
这个E(UHF)是指electronic energy，每优化一步，就会给出一个E(UHF)。对大部分寻找minimum的优化来说，每一步的E(UHF)总体趋势上是一个降低的过程，

在你的单重态输出中，electronic energy应该是
=-1288.40 ...

谢谢！

39楼2011-08-06 10:45:30

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

tsb1980

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1.5
帖子: 11
在线: 16小时
虫号: 590665
注册: 2008-09-01
性别: GG
专业: 理论和计算化学

★ ★
zhou2009(金币+2): 2011-08-06 15:18:16

主要原因：
1 计算方法，PM3太粗糙
2 RHF计算时应调好开始的结构，
3 基组应选择恰当。