计算化学公社

标题: 高斯优化含有Eu的结构时不收敛的问题请教 [打印本页]

作者
Author: learn-GAUSSIAN 时间: 2023-4-17 19:18
标题: 高斯优化含有Eu的结构时不收敛的问题请教
我在优化含有Eu的体系时，程序总是报错。输出文件的末尾提示：
The polynomial fit failed. Using point 1.
A contracting polynomial of degree 16 produced 0.0000
Search did not lower the energy significantly.
No lower point found -- run aborted.

现在我将输入输出文件提交上，请求各位指点！

作者
Author: Boer 时间: 2023-4-17 20:49
有些元素你都没有定义基组；
建议给色散关键词加上BJ阻尼函数。

作者
Author: sobereva 时间: 2023-4-18 00:43
另外，不要轻易用xqc
SCF不收敛怎么解决看此文
解决SCF不收敛问题的方法
http://sobereva.com/61

专门针对镧系配合物计算的一些关键细节看此文
使用Gaussian做镧系金属配合物的量子化学计算
http://sobereva.com/581（http://bbs.keinsci.com/thread-21274-1-1.html）

作者
Author: wzkchem5 时间: 2023-4-18 03:51
分子结构极度不合理。建议查一下DMF的键长键角，把初始结构的键长键角调成和实验值差不多（键长差别小于0.05埃，键角差别小于10度）。同时系统学习如何不经查阅文献、靠化学常识判断常见分子的键长键角范围。

作者
Author: learn-GAUSSIAN 时间: 2023-4-19 22:12

sobereva 发表于 2023-4-18 00:43
另外，不要轻易用xqc
SCF不收敛怎么解决看此文
解决SCF不收敛问题的方法

谢谢！尝试修改后，这次.log末尾报错行提示如下，实在不知道如何解决，烦请不吝赐教
Symmetry turned off by external request.
Stoichiometry C14H23Cl2EuN2O8
Framework group  C1[X(C14H23Cl2EuN2O8)]
Deg. of freedom 144
Full point group                C1    NOp 1
Rotational constants (GHZ):    0.2142136    0.0571928    0.0508672
Leave Link  202 at Wed Apr 19 21:19:16 2023, MaxMem=  6291456000 cpu:       0.2
(Enter /opt/software/gaussian/g09.E01/g09/l301.exe)
General basis read from cards:  (5D, 7F)
Centers:    8    9    11    15    19    23    27    31    34    35
Centers:    36    37    38    39    44    48    2    5    21    33
Centers:    45    46    49    50    10    22    3    4    6    7
Centers:    12    13    14    16    17    18    20    24    25    26
Centers:    28    29    30    32    40    41    42    43    47
6-31G*
****
Centers:    1
mwb52
****
======================================================================================================
                                    Pseudopotential Parameters
======================================================================================================
  Center    Atomic    Valence    Angular    Power
  Number    Number    Electrons    Momentum    of R    Exponent       Coefficient SO-Coeffient
======================================================================================================
EOF while reading ECP pointer card.
Error termination via Lnk1e in /opt/software/gaussian/g09.E01/g09/l301.exe at Wed Apr 19 21:19:16 2023.
Job cpu time:    0 days  0 hours  0 minutes 10.4 seconds.
File lengths (MBytes):  RWF=    27 Int=    0 D2E=    0 Chk=    2 Scr=    2

作者
Author: AndrewX 时间: 2023-4-20 00:08

learn-GAUSSIAN 发表于 2023-4-19 22:12
谢谢！尝试修改后，这次.log末尾报错行提示如下，实在不知道如何解决，烦请不吝赐教
Symmetry turned of ...

末尾多打几个回车

作者
Author: sobereva 时间: 2023-4-20 04:31

learn-GAUSSIAN 发表于 2023-4-19 22:12
谢谢！尝试修改后，这次.log末尾报错行提示如下，实在不知道如何解决，烦请不吝赐教
Symmetry turned of ...

上传完整输入文件，否则别人没法准确判断你怎么写的

作者
Author: learn-GAUSSIAN 时间: 2023-4-24 16:44

Boer 发表于 2023-4-17 20:49
有些元素你都没有定义基组；
建议给色散关键词加上BJ阻尼函数。

谢谢提醒，刚注意到

作者
Author: learn-GAUSSIAN 时间: 2023-4-24 16:45

AndrewX 发表于 2023-4-20 00:08
末尾多打几个回车

按照您的办法已经把问题解决了，谢谢！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)