| 查看: 2150 | 回复: 5 | |||
我家小屁孩铜虫 (小有名气)
|
[求助]
量化中出错的log文件
|
|
虫虫们好： 有劳见多识广的虫虫帮俺分析一下此log文件出错的原因并怎样改正，单线态已顺利完成，在算五线态时出错了，多谢了！！！！ Gaussian 09: EM64L-G09RevA.01 8-May-2009 17-Aug-2011 ****************************************** %chk=/home/lyy/dqs/perpen.chk %mem=2000MB %nprocshared=4 Will use up to 4 processors via shared memory. --------------------------------------------------- # opt freq=noraman b3lyp/6-31g(d) geom=connectivity --------------------------------------------------- 1/14=-1,18=20,19=15,26=3,38=1,57=2/1,3; 2/9=110,12=2,17=6,18=5,40=1/2; 3/5=1,6=6,7=1,11=2,16=1,25=1,30=1,71=1,74=-5/1,2,3; 4//1; 5/5=2,38=5/2; 6/7=2,8=2,9=2,10=2,28=1/1; 7//1,2,3,16; 1/14=-1,18=20,19=15/3(2); 2/9=110/2; 99//99; 2/9=110/2; 3/5=1,6=6,7=1,11=2,16=1,25=1,30=1,71=1,74=-5/1,2,3; 4/5=5,16=3/1; 5/5=2,38=5/2; 7//1,2,3,16; 1/14=-1,18=20,19=15/3(-5); 2/9=110/2; 6/7=2,8=2,9=2,10=2,19=2,28=1/1; 99/9=1/99; ------ perpen ------ Symbolic Z-matrix: Charge = 0 Multiplicity = 5 （中间部分略） (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) The electronic state of the initial guess is 5-A. Initial guess 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. EnCoef did 100 forward-backward iterations EnCoef did 11 forward-backward iterations EnCoef did 3 forward-backward iterations EnCoef did 6 forward-backward iterations EnCoef did 5 forward-backward iterations EnCoef did 3 forward-backward iterations Rare condition: small coef for last iteration: 0.000D+00 Rare condition: small coef for last iteration: 0.000D+00 Restarting incremental Fock formation. Rare condition: small coef for last iteration: 0.000D+00 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(UB3LYP) = -4473.53934912 A.U. after 129 cycles Convg = 0.2980D-04 -V/T = 2.0026 Annihilation of the first spin contaminant: S**2 before annihilation 6.0847, after 6.0019 DiagDN has N= 676 LTot= 2060 but NE2= 3 cannot use DSYEVD. Convergence failure -- run terminated. Error termination via Lnk1e in /home/lyy/g09/l502.exe at Thu Aug 18 02:36:02 2011. Job cpu time: 0 days 18 hours 8 minutes 33.1 seconds. File lengths (MBytes): RWF= 427 Int= 0 D2E= 0 Chk= 11 Scr= 1 小女是实验为主，量化知识浅薄，希望能得到高手虫虫们的指点！ |
回复此楼» 收录本帖的淘帖专辑推荐
 高斯日常学习总结 |
» 猜你喜欢
调剂招生
已经有5人回复
-
招收0703开头的学硕调剂生 线上复试
已经有5人回复
-
物理化学论文润色/翻译怎么收费?
已经有152人回复
-
厦门大学侯旭团队招聘青年教师、博士后及博硕士研究生
已经有2人回复
-
哈工大深圳 表界面物理化学方向（胶体与界面，光谱探测研究） 博士招生
已经有0人回复
-
元素硫的化合价和氧化态
已经有1人回复
-
江苏科技大学2026年招聘启事
已经有1人回复
-
International Journal of Nanoscience属于中科院4区吗？
已经有0人回复
-
澳大利亚南昆士兰大学量子点课题组招收博士生
已经有0人回复
-
香港中文大学（深圳）2026年招聘公告
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- ansys fluent并行计算许可文件试用【转载】
已经有11人回复
- 如何将Gaussian计算后的输出文件的类型.chk或.log改为*.mol2类型
已经有5人回复
- 关于化合物的log D PSA 等理化特征
已经有16人回复
- 【求助】 GAMESS 里log文件里描述原子轨道AG，B1G，B2G，B3G，B1U等等是什么意思啊？
已经有9人回复
- 【求助】量化计算投稿ACS JPCC Supporting Material的准备问题
已经有7人回复
- 【讨论】大家好，加水后跑完分子动力学，在log中发现了问题
已经有15人回复
- 【求助】vasp优化结构未结束，log文件出现下列问题停止计算，请高手指教
已经有6人回复
- 【求助】log和ln的困惑，是否一样？
已经有5人回复
- 【求助】log文件和fchk文件不一致
已经有4人回复
fatpig8832
铁杆木虫 (著名写手)
- 应助: 16 (小学生)
- 金币: 7420.2
- 散金: 31
- 红花: 28
- 帖子: 1438
- 在线: 610.9小时
- 虫号: 770403
- 注册: 2009-05-14
【答案】应助回帖
★ ★
zhou2009(金币+2): 2011-08-24 19:40:30
我家小屁孩(金币+10): 不好意思，我是小木虫新手，不怎么上，发金币也晚了点。很感谢你的！ 2011-09-02 17:39:48
zhou2009(金币+2): 2011-08-24 19:40:30
我家小屁孩(金币+10): 不好意思，我是小木虫新手，不怎么上，发金币也晚了点。很感谢你的！ 2011-09-02 17:39:48
|
无他，构型不好而已... 1. 改构型 2. 可用 scfcyc=500 增大scf计算的cycle，默认的是129 （SCF Done: E(UB3LYP) = -4473.53934912 A.U. after 129 cycles） 3. scf=qc，用了这个，速度很慢；所以一般用scf=qc算第一点即可，然后再读取第一点的chk文件继续往下做... |
2楼2011-08-24 11:06:23
bubuweiying
金虫 (著名写手)
- 应助: 4 (幼儿园)
- 金币: 2512.3
- 散金: 1447
- 红花: 5
- 沙发: 3
- 帖子: 1084
- 在线: 351.6小时
- 虫号: 968063
- 注册: 2010-03-11
- 性别: MM
- 专业: 有机合成
★ ★
zhou2009(金币+2): 2011-08-24 19:40:47
zhou2009(金币+2): 2011-08-24 19:40:47
|
如果还不行，只能拿出杀手锏了，就是使用qc，但不建议直接使用，而是使用xqc关键词，比如scf=(maxcycle=80,xqc)，意思是如果scf正常计算(dc)在80个循环之内不收敛才进行昂贵的qc计算，因为scf不收敛多数在几个优化的过程中出现，无法判断哪一步优化的时候会出现scf不收敛，所以用xqc比纯粹使用qc要省时的多。 使用scf=qc，这个计算会慢，而且需要用stable关键字来测试结果是否波函数稳定。如果这个还不收敛，会提示L508错误，请问如何测试波函数稳定啊，是加上关键词scf=xqc stable 就可以么，即# opt freq rb3lyp/6-31+g(d) pop=(nbo,savenbo,full) scf=xqc stable 谢谢啊 |
3楼2011-08-24 11:43:49
fatpig8832
铁杆木虫 (著名写手)
- 应助: 16 (小学生)
- 金币: 7420.2
- 散金: 31
- 红花: 28
- 帖子: 1438
- 在线: 610.9小时
- 虫号: 770403
- 注册: 2009-05-14
【答案】应助回帖
★ ★
zhou2009(金币+2): 2011-08-24 19:41:04
我家小屁孩(金币+5): 3ks 2011-09-12 12:50:20
zhou2009(金币+2): 2011-08-24 19:41:04
我家小屁孩(金币+5): 3ks 2011-09-12 12:50:20
|
测试波函数稳定否，用stable=opt，只需算一个单点即可，不用算freq之类 例如 # hf/6-31g* stable=opt 若用scf=qc，也可以写成 # hf/6-31g* stable=opt scf=qc 如要在波函数稳定的基础上继续优化或算freq，可用--link1-- # hf/6-31g* stable=opt scf=qc ... ... ... 空一行 --link1-- %chk=...跟上面的一致... # hf/6-31g* opt freq=noraman guess=read geom=allcheck |
4楼2011-08-24 12:10:41
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +459 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13336.5
- 散金: 27307
- 红花: 301
- 帖子: 4683
- 在线: 928.6小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
5楼2011-08-24 21:15:14
fwenling0506
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1686.4
- 红花: 3
- 帖子: 300
- 在线: 175.5小时
- 虫号: 866403
- 注册: 2009-10-09
- 专业: 理论和计算化学
6楼2011-09-04 10:34:01
20