| 查看: 1910 | 回复: 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 小女是实验为主，量化知识浅薄，希望能得到高手虫虫们的指点！ |
回复此楼» 收录本帖的淘帖专辑推荐
 高斯日常学习总结 |
» 猜你喜欢
汕头大学化学招调剂，相关内容如下。
已经有4人回复
-
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有11人回复
-
物理化学论文润色/翻译怎么收费?
已经有197人回复
-
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有10人回复
-
求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340
已经有2人回复
-
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料
已经有0人回复
-
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建
已经有0人回复
-
半透性容器进行高湿或者低湿研究
已经有0人回复
-
求ENDNOTE破解软件安装包
已经有4人回复
-
河海大学 （211.双一流高校）钙钛矿光电实验室丁勇团队欢迎有志青年！
已经有3人回复
-
黑龙江省原子能研究院化学课题组招收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
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +458 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 12988.7
- 散金: 27307
- 红花: 300
- 帖子: 4672
- 在线: 926.9小时
- 虫号: 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
5