| 查看: 858 | 回复: 6 | |||
| 当前主题已经存档。 | |||
Jldbaobei木虫 (著名写手)
按时毕业找到好工作
|
[交流]
【求助】为何高斯不报告错误，但是算到一定的步就停止了呢
|
||
|
刚开始学Gaussian 自己写了一个C2与NO的反应做练习，寻找过渡态 使用OPT=QST2，UHF/6-31G(d)方法 设定maxcyc=300，但是在输出文件中发现计算总是显示 “Step number ** out of a maximum of 100” 为什么最大不是300呢？ 还有一个问题就是整个计算过程没有出现错误，但是Gaussian每次计算到第83步时便停止不动了，不报告有错误，一直是“Processing”的显示，但是计算就是不再进行了 算了好几次，总是这个样子 请问是为什么呢？ |
回复此楼» 猜你喜欢
调剂信息
已经有0人回复
-
河北农大生物与医药专业招收调剂，名额充足
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有81人回复
-
齐鲁工业大学 国家重点实验室-新能源、锂离子电池、材料、化学方向
已经有0人回复
-
杰青/长江团队招收调剂
已经有0人回复
-
中南林业科技大学罗勇锋教授团队2026年招收材料类调剂研究生
已经有0人回复
-
投稿期刊求助
已经有1人回复
-
电子科技大学材料学院王翀教授拟招收工程博士一名（化学、材料方向）
已经有1人回复
-
湖南农业大学能源材料创新团队招收2026化学专业研究生（学硕）调剂
已经有3人回复
-
江西师大化学工程学院调剂今晚0点开放！ 缺额100个，速冲！
已经有0人回复
-
奥尔堡Xinxin Xiao课题组招收一名岗位制博士（酶基生物电化学）
已经有0人回复
Jldbaobei
木虫 (著名写手)
按时毕业找到好工作
- 应助: 2 (幼儿园)
- 金币: 3294.4
- 散金: 310
- 帖子: 1019
- 在线: 458.6小时
- 虫号: 507695
- 注册: 2008-02-19
- 专业: 燃烧学
2楼2008-10-22 15:37:28
Jldbaobei
木虫 (著名写手)
按时毕业找到好工作
- 应助: 2 (幼儿园)
- 金币: 3294.4
- 散金: 310
- 帖子: 1019
- 在线: 458.6小时
- 虫号: 507695
- 注册: 2008-02-19
- 专业: 燃烧学
|
输出文件的最后一部分如下： Item Value Threshold Converged? Maximum Force 0.016586 0.000450 NO RMS Force 0.012317 0.000300 NO Maximum Displacement 0.001265 0.001800 YES RMS Displacement 0.001067 0.001200 YES Predicted change in Energy= 1.576758D-05 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Input orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 6 0 -0.080000 0.000000 -0.015798 2 6 0 0.278376 0.000000 1.313980 3 7 0 0.087982 0.000000 2.464664 4 8 0 -0.104614 0.000000 3.628657 --------------------------------------------------------------------- Distance matrix (angstroms): 1 2 3 4 1 C 0.000000 2 C 1.377222 0.000000 3 N 2.486144 1.166330 0.000000 4 O 3.644538 2.346149 1.179819 0.000000 Stoichiometry C2NO(2) Framework group CS[SG(C2NO)] Deg. of freedom 5 Full point group CS NOp 2 Largest Abelian subgroup CS NOp 2 Largest concise Abelian subgroup C1 NOp 1 Standard orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 6 0 0.443834 1.972433 0.000000 2 6 0 -0.126810 0.718995 0.000000 3 7 0 -0.126810 -0.447335 0.000000 4 8 0 -0.126810 -1.627153 0.000000 --------------------------------------------------------------------- Rotational constants (GHZ): 436.9977704 5.0576939 4.9998272 Standard basis: 6-31G(d) (6D, 7F) There are 44 symmetry adapted basis functions of A' symmetry. There are 16 symmetry adapted basis functions of A" symmetry. Integral buffers will be 262144 words long. Raffenetti 2 integral format. Two-electron integral symmetry is turned on. 60 basis functions, 112 primitive gaussians, 60 cartesian basis functions 14 alpha electrons 13 beta electrons nuclear repulsion energy 84.7413681490 Hartrees. NAtoms= 4 NActive= 4 NUniq= 4 SFac= 1.00D+00 NAtFMM= 60 Big=F One-electron integrals computed using PRISM. NBasis= 60 RedAO= T NBF= 44 16 NBsUse= 60 1.00D-06 NBFU= 44 16 Initial guess read from the read-write file: Initial guess orbital symmetries: Alpha Orbitals: Occupied (A') (A') (A') (A') (A') (A') (A') (A') (A") (A') (A") (A') (A') (A") Virtual (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') (A") (A') (A') (A") (A') (A') (A') (A") (A') (A') (A') (A') (A') (A') Beta Orbitals: Occupied (A') (A') (A') (A') (A') (A') (A') (A') (A') (A") (A') (A") (A') Virtual (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") (A') (A") (A') (A') (A") (A') (A') (A') (A") (A') (A') (A') (A') (A') (A') 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. Keep R1 and R2 integrals in memory in canonical form, NReq= 3921837. Rare condition: small coef for last iteration: 0.000D+00 Rare condition: small coef for last iteration: 0.000D+00 Rare condition: small coef for last iteration: 0.000D+00 Rare condition: small coef for last iteration: 0.000D+00 Restarting incremental Fock formation. SCF Done: E(UHF) = -204.785403018 A.U. after 29 cycles Convg = 0.6130D-08 -V/T = 2.0020 S**2 = 0.9608 Annihilation of the first spin contaminant: S**2 before annihilation 0.9608, after 0.7789 ***** Axes restored to original set ***** ------------------------------------------------------------------- Center Atomic Forces (Hartrees/Bohr) Number Number X Y Z ------------------------------------------------------------------- 1 6 0.008630697 0.000000000 0.012509389 2 6 -0.008316091 0.000000000 -0.029430588 3 7 -0.006059895 0.000000000 0.030178532 4 8 0.005745288 0.000000000 -0.013257333 ------------------------------------------------------------------- Cartesian Forces: Max 0.030178532 RMS 0.013911991 |
3楼2008-10-22 15:43:25
quantumor
金虫 (著名写手)
快乐兔子
- 应助: 13 (小学生)
- 金币: 8506.3
- 散金: 10
- 红花: 14
- 帖子: 1898
- 在线: 215.9小时
- 虫号: 630331
- 注册: 2008-10-19
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
 “Step number ** out of a maximum of 100”是Gaussian程序对计算体系结构优化完成所需结构变更数的内在设定，也就是说某结构优化大约能在100个结构变更内完成，所以，当实际计算的结构变更数超过该数时，计算回自行终止。但这不妨碍你从输出文件中取出最后的或你认为最好的结构，重写输入文件继续进行优化。 看了一下LZ的分子结构，似乎结构上不合理，计算的结构似乎离合理稳定结构相差很远，所以，可能是Gaussian被你的结构搞瞢了，面对这样奇特的结构，程序不知道该如何处理的好，在这种情况下，程序可能会莫名其妙地退出来——实际上可能是死机了。 要继续计算，还是修改一下结构吧。下图是你的结构。 |
4楼2008-10-22 18:47:05
Jldbaobei
木虫 (著名写手)
按时毕业找到好工作
- 应助: 2 (幼儿园)
- 金币: 3294.4
- 散金: 310
- 帖子: 1019
- 在线: 458.6小时
- 虫号: 507695
- 注册: 2008-02-19
- 专业: 燃烧学
5楼2008-10-23 09:22:09
6楼2008-10-23 10:19:39
Jldbaobei
木虫 (著名写手)
按时毕业找到好工作
- 应助: 2 (幼儿园)
- 金币: 3294.4
- 散金: 310
- 帖子: 1019
- 在线: 458.6小时
- 虫号: 507695
- 注册: 2008-02-19
- 专业: 燃烧学
7楼2008-10-23 12:21:37
100