应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】为何高斯不报告错误，但是算到一定的步就停止了呢
51/1返回列表
查看: 810  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

Jldbaobei

木虫 (著名写手)

按时毕业找到好工作

[交流] 【求助】为何高斯不报告错误，但是算到一定的步就停止了呢

刚开始学Gaussian
自己写了一个C2与NO的反应做练习，寻找过渡态
使用OPT=QST2，UHF/6-31G(d)方法
设定maxcyc=300，但是在输出文件中发现计算总是显示
“Step number  **  out of a maximum of 100”
为什么最大不是300呢？
还有一个问题就是整个计算过程没有出现错误，但是Gaussian每次计算到第83步时便停止不动了，不报告有错误，一直是“Processing”的显示，但是计算就是不再进行了
算了好几次，总是这个样子
请问是为什么呢？
回复此楼

» 猜你喜欢
1楼 2008-10-22 14:47:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jldbaobei

木虫 (著名写手)

按时毕业找到好工作
我也是门外汉硬拿起这个工具
理论方面的比较差
现在纯是自己做练习摸索
大家有什么好的入门经验多教教我吧
赞一下 回复此楼
7楼2008-10-23 12:21:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

Jldbaobei

木虫 (著名写手)

按时毕业找到好工作
最终一直处于L103.exe is processing...

请大家帮忙解释一下吧
赞一下 回复此楼
2楼2008-10-22 15:37:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jldbaobei

木虫 (著名写手)

按时毕业找到好工作
输出文件的最后一部分如下：

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')
of initial guess= 1.0415
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

quantumor

金虫 (著名写手)

快乐兔子

“Step number  **  out of a maximum of 100”是Gaussian程序对计算体系结构优化完成所需结构变更数的内在设定，也就是说某结构优化大约能在100个结构变更内完成，所以，当实际计算的结构变更数超过该数时，计算回自行终止。但这不妨碍你从输出文件中取出最后的或你认为最好的结构，重写输入文件继续进行优化。
看了一下LZ的分子结构，似乎结构上不合理，计算的结构似乎离合理稳定结构相差很远，所以，可能是Gaussian被你的结构搞瞢了，面对这样奇特的结构，程序不知道该如何处理的好，在这种情况下，程序可能会莫名其妙地退出来——实际上可能是死机了。
要继续计算，还是修改一下结构吧。下图是你的结构。
赞一下 回复此楼
愿好运与快乐伴随你！
4楼2008-10-22 18:47:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭