| 查看: 2521 | 回复: 12 | ||
flyingfish00铜虫 (正式写手)
|
[求助]
求支援--gaussian 柔性扫描总是莫名退出，找不到原因。
|
|
做一个小的计算，进行邻苯二酚的柔性扫描 （opt=modredundant）。分子结构是已经优化好的分子结构，导入进行柔性扫描。 开始使用小基组(uhf/3-21g)来做柔性扫描，还挺顺利的，顺利的找到了能量较高的点。 所使用的gjf文件内容为： —————————————————————————————— %chk=C:\Gaussian\phenol\04_phenol_0H_othro_redundant.chk %nprocshared=2 # opt=modredundant uhf/3-21g nosymm scf(maxcycle=300,qc) geom=connectivity phenol-OH ortho redundant 0 2 C -1.67101400 -1.02358000 0.23264100 C -1.96232700 0.33247600 -0.06700600 C -0.96971200 1.21545500 -0.39033700 C 0.48579000 0.84642900 -0.36921900 C 0.68967700 -0.62759100 -0.13072800 C -0.33582000 -1.48367900 0.18943500 O 1.99025800 -1.00483700 -0.24180800 H -2.46949000 -1.71202800 0.48972100 H -2.99664000 0.66750300 -0.04492700 H -1.19241400 2.25467200 -0.61602000 H -0.11817000 -2.53147100 0.39541100 H 2.06095800 -1.95751000 -0.06419300 H 0.98144000 1.12280900 -1.31073500 O 1.20543100 1.64142600 0.60539400 H 0.74924500 1.50625700 1.45334400 1 2 1.5 6 1.5 8 1.0 2 3 2.0 9 1.0 3 4 1.0 10 1.0 4 5 1.0 13 1.0 14 1.0 5 6 2.0 7 1.0 6 11 1.0 7 12 1.0 8 9 10 11 12 13 14 15 1.0 15 B 4 14 S 3 0.0500000 ———————————————————————————— 结果还挺顺利，顺利完成扫描，下面是结果： 分子结构：  能量变化曲线  小基组完了以后，就准备采用更高一点的基组(ub3lyp/6-31g(d))来重新进行柔性扫描。gjf内容还是上面的内容，就是基组换成了 ub3lyp/6-31g(d)。 结果问题来了，运行了一段时候后，gaussian就弹出个对话框： openprocess failed in startlink。也没说原因。 再查看out文件，out文件最后也没有说明是什么原因。仿佛就是自然中断了一样。 ———————————————————————————————— Harris functional with IExCor= 402 diagonalized for initial guess. ExpMin= 1.61D-01 ExpMax= 5.48D+03 ExpMxC= 8.25D+02 IAcc=1 IRadAn= 1 AccDes= 0.00D+00 HarFok: IExCor= 402 AccDes= 0.00D+00 IRadAn= 1 IDoV= 1 ScaDFX= 1.000000 1.000000 1.000000 1.000000 FoFCou: FMM=F IPFlag= 0 FMFlag= 100000 FMFlg1= 0 NFxFlg= 0 DoJE=T BraDBF=F KetDBF=T FulRan=T Omega= 0.000000 0.000000 1.000000 0.000000 0.000000 ICntrl= 500 IOpCl= 0 NMat0= 1 NMatS0= 1 NMatT0= 0 NMatD0= 1 NMtDS0= 0 NMtDT0= 0 I1Cent= 4 NGrid= 0. Symmetry not used in FoFCou. Requested convergence on RMS density matrix=1.00D-08 within 40 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. Initial convergence to 1.0D-05 achieved. Increase integral accuracy. SCF Done: E(UB3LYP) = -383.217669131 A.U. after 15 cycles Convg = 0.5646D-08 -V/T = 2.0094 Annihilation of the first spin contaminant: S**2 before annihilation 0.7826, after 0.7506 QCSCF skips out because SCF is already converged. Calling FoFJK, ICntrl= 2127 FMM=F ISym2X=0 I1Cent= 0 IOpClX= 1 NMat=1 NMatS=1 NMatT=0. ------------------------------------------------------------------- Center Atomic Forces (Hartrees/Bohr) Number Number X Y Z ------------------------------------------------------------------- 1 6 0.000552860 0.000261531 -0.000034388 2 6 -0.000137206 -0.000662756 0.000380411 3 6 0.003328125 0.001029937 0.001394776 4 6 0.010272746 0.012638895 0.019845318 5 6 0.001611720 0.003593828 0.001147491 6 6 -0.001022950 -0.000054718 0.000259065 7 8 -0.000318132 -0.000142463 -0.000093789 8 1 0.000010625 -0.000006488 -0.000002513 9 1 0.000045225 -0.000024698 -0.000178775 10 1 -0.000091426 -0.000127936 -0.000143429 11 1 0.000010046 0.000045639 -0.000135751 12 1 0.000034210 -0.000016400 -0.000039541 13 1 0.001273744 0.001644310 0.002561579 14 8 -0.014140392 -0.016929990 -0.024093930 15 1 -0.001429195 -0.001248691 -0.000866523 ------------------------------------------------------------------- Cartesian Forces: Max 0.024093930 RMS 0.006279918 ———————————————————————————————— 目前，我所采取的对策有： 1. 使用优化后的分子构型； 2. 采用scf=qc； 3. 减小 柔性扫描步长，原来步长是0.1，后来改成 0.05 4. 增加迭代次数 maxcyc=300; 5. 不考虑对称性，增加参数 nosymm； 但是都没有解决。而且用GaussianView打开chk文件，也打不开。说是非正常中断，不能打开。但是以前有柔性扫描中断什么的，用GaussianView也是可以打开的，比如扫了4-5步，这几步的数据就能打开。 6.现在我又尝试减少步数（steps）。开始是3步，成功；然后增加到5步，成功；再增加到10次，失败，莫名退出。估计这个步数有问题了。请问我如何解决？ 真是奇怪，还请各位大牛帮忙指点一下。还有什么更好的解决方案没。多谢！ [ Last edited by flyingfish00 on 2012-11-24 at 13:36 ] |
回复此楼» 猜你喜欢
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生
已经有11人回复
-
中国科学院成都山地灾害与环境研究所2026年特别研究助理（博士后）招聘启事
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有230人回复
-
招生材料与化工、环境科学调剂研究生
已经有5人回复
-
湖南大学刘巧玲课题组2026年博士研究生招生信息
已经有28人回复
-
南方科技大学表面催化反应动力学研究课题组招收2026届博士生
已经有1人回复
-
北京邮电大学物理学院老师招收博士生
已经有3人回复
-
Comsol 6.0 桶内肥料产热产湿过程模拟求助
已经有1人回复
-
微纳hub，微纳加工的宝库
已经有0人回复
-
微纳加工里的 lift-off 工艺
已经有0人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有37人回复
» 本主题相关商家推荐: (我也要在这里推广)
flyingfish00
铜虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 32.4
- 散金: 412
- 红花: 27
- 帖子: 942
- 在线: 402.6小时
- 虫号: 289541
- 注册: 2006-10-22
- 性别: GG
- 专业: 应用高分子化学与物理
2楼2012-11-24 20:01:44
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
flyingfish00: 金币+10, ★有帮助, 非常感谢！ 2012-11-25 15:39:51
感谢参与，应助指数 +1
flyingfish00: 金币+10, ★有帮助, 非常感谢！ 2012-11-25 15:39:51
|
本帖内容被屏蔽 |
3楼2012-11-25 13:38:10
flyingfish00
铜虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 32.4
- 散金: 412
- 红花: 27
- 帖子: 942
- 在线: 402.6小时
- 虫号: 289541
- 注册: 2006-10-22
- 性别: GG
- 专业: 应用高分子化学与物理
4楼2012-11-25 15:40:57
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
flyingfish00: 金币+40, ★★★很有帮助, 非常感谢！ 2012-11-25 20:57:45
flyingfish00: 金币+40, ★★★很有帮助, 非常感谢！ 2012-11-25 20:57:45
|
本帖内容被屏蔽 |
5楼2012-11-25 15:52:29
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
flyingfish00: 金币+60, ★★★很有帮助, 非常感谢！ 2012-11-25 21:02:02
flyingfish00: 金币+60, ★★★很有帮助, 非常感谢！ 2012-11-25 21:02:02
|
本帖内容被屏蔽 |
6楼2012-11-25 16:01:14
flyingfish00
铜虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 32.4
- 散金: 412
- 红花: 27
- 帖子: 942
- 在线: 402.6小时
- 虫号: 289541
- 注册: 2006-10-22
- 性别: GG
- 专业: 应用高分子化学与物理
7楼2012-11-25 21:01:28
flyingfish00
铜虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 32.4
- 散金: 412
- 红花: 27
- 帖子: 942
- 在线: 402.6小时
- 虫号: 289541
- 注册: 2006-10-22
- 性别: GG
- 专业: 应用高分子化学与物理
8楼2012-11-25 21:02:31
9楼2012-11-27 09:19:50
flyingfish00
铜虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 32.4
- 散金: 412
- 红花: 27
- 帖子: 942
- 在线: 402.6小时
- 虫号: 289541
- 注册: 2006-10-22
- 性别: GG
- 专业: 应用高分子化学与物理
10楼2012-11-27 16:11:24
20