求助：振动分析不收敛的原因和解决办法

yihanxu · 发表于 Post on 2020-1-24 15:30:23

大家新年好！遇到振动分析不收敛，向大家求助该怎么办呢？详见附件！
衷心祝愿大家2020新年万福！
感谢！

sobereva · 发表于 Post on 2020-1-24 17:43:03

输入文件也要上传，否则不知道gen对应的什么基组。当前体系也没有任何理由需要自定义基组。
scrf=read显得莫名其妙，意义完全不明
当前的结构怎么来的也需要说明。

naonao5205 · 发表于 Post on 2020-1-24 18:34:28

sobereva 发表于 2020-1-24 17:43
输入文件也要上传，否则不知道gen对应的什么基组。当前体系也没有任何理由需要自定义基组。
scrf=read显得 ...

sob老师，之前我也遇到过这个问题，推测可能的原因是弥散函数，在我砍掉弥散函数后6-31+g(d,p)到仅反应中心带弥散后面不带即可解决。或者int=acc2e=11（之前是10），也解决了。

sobereva · 发表于 Post on 2020-1-25 23:19:48

naonao5205 发表于 2020-1-24 18:34
sob老师，之前我也遇到过这个问题，推测可能的原因是弥散函数，在我砍掉弥散函数后6-31+g(d,p)到仅反应中 ...

用大量弥散函数的时候建议通过cphf=grid=x（x=fine、ultrafine等）增加DFT积分格点精度，有助于解决CPHF阶段不收敛。
当然，弥散函数能不用就尽量不用

yihanxu · 发表于 Post on 2020-1-30 07:46:13

sobereva 发表于 2020-1-24 03:43
输入文件也要上传，否则不知道gen对应的什么基组。当前体系也没有任何理由需要自定义基组。
scrf=read显得 ...

老师好，我是想算这个分子在溶剂环境中的RESP电荷，所以先优化了一下，然后想做freq看看有没有虚频。计算级别是b3lyp/def2-tzvp谢谢！

sobereva · 发表于 Post on 2020-1-30 22:30:42

yihanxu 发表于 2020-1-30 07:46
老师好，我是想算这个分子在溶剂环境中的RESP电荷，所以先优化了一下，然后想做freq看看有没有虚频。计算 ...

本来def2在你的G16 B01中就是内置的，何故用gen
用def2-TZVP做opt freq对这么大体系完全是浪费时间，def-TZVP或更次一点的6-311G*都已经足够了
另外你这体系里每个N上真的只连一个氢？难道不是氨基么？

yihanxu · 发表于 Post on 2020-1-31 07:36:04

sobereva 发表于 2020-1-30 08:30
本来def2在你的G16 B01中就是内置的，何故用gen
用def2-TZVP做opt freq对这么大体系完全是浪费时间，def ...

谢谢老师！
内置的话，是直接写def2TZVP关键词就可以吗
谢谢老师问的问题，我也是自己想像这个是=N-这种，双键那边连六元环，单键这边连氢，请问老师有什么不合理的地方吗？谢谢！

sobereva · 发表于 Post on 2020-2-1 02:13:56

yihanxu 发表于 2020-1-31 07:36
谢谢老师！
内置的话，是直接写def2TZVP关键词就可以吗
谢谢老师问的问题，我也是自己想像这个是=N-这 ...

y
我没见过你这种=N-H的基团。掉一个氢的氨基自由基才这模样

yihanxu · 发表于 Post on 2020-2-1 05:19:37

sobereva 发表于 2020-1-31 12:13
y
我没见过你这种=N-H的基团。掉一个氢的氨基自由基才这模样

谢谢老师的回复，我化学基础弱，当时确实也没细想，附件这两篇文献里有这样的结构

sobereva · 发表于 Post on 2020-2-1 05:59:51

yihanxu 发表于 2020-2-1 05:19
谢谢老师的回复，我化学基础弱，当时确实也没细想，附件这两篇文献里有这样的结构

倘若有证据表明（可以google搜更多相关信息，或者到诸如pubchem、chemspider里通过自己画的结构来搜索）那种结构就是对的，那就那么算，反正我是没见过。这和freq的CPHF阶段是否能收敛并没必然关系。

PS：pdf文件没必要压缩，压缩不压缩大小基本没变化

yihanxu · 发表于 Post on 2020-2-1 09:40:47

本帖最后由 yihanxu 于 2020-1-31 19:42 编辑

sobereva 发表于 2020-1-31 15:59
倘若有证据表明（可以google搜更多相关信息，或者到诸如pubchem、chemspider里通过自己画的结构来搜索） ...

谢谢老师，查到了这种分子，是只要在pubchem搜到就可以说这种含氮官能团是对的吗？https://pubchem.ncbi.nlm.nih.gov/compound/602819

sobereva · 发表于 Post on 2020-2-3 04:49:30

yihanxu 发表于 2020-2-1 09:40
谢谢老师，查到了这种分子，是只要在pubchem搜到就可以说这种含氮官能团是对的吗？https://pubchem.ncbi. ...

是

yihanxu · 发表于 Post on 2020-2-3 08:59:25

老师好，我的体系这种freq的CPHF阶段不收敛还可能是因为什么呢？我没用弥散，所以应该不是弥散的影响;虽然基组大但是这会导致不收敛吗？还是可能有其他的原因?

sobereva · 发表于 Post on 2020-2-4 10:10:54

yihanxu 发表于 2020-2-3 08:59
老师好，我的体系这种freq的CPHF阶段不收敛还可能是因为什么呢？我没用弥散，所以应该不是弥散的影响;虽然 ...

mutron · 发表于 Post on 2020-2-9 21:59:34

C=N-H是氨分子和酮羰基形成类似席夫碱结构的亚胺，一般芳香胺产物才比较稳定可以分离，氨产物比较不稳定

题主用的g16默认双电子积分精度acc2e=12，而且不是弥散函数，所以不是这个问题。如果把中间两个C=N键其中一个反转180度变成C2对称的话，在6-311G*基组、int=fine（CPHF grid对应CoarseGrid）下计算频率都没有这个问题

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 求助：振动分析不收敛的原因和解决办法

评分 Rate

评分 Rate

评分 Rate

评分 Rate

评分 Rate

评分 Rate

评分 Rate

评分 Rate

评分 Rate

浏览过的版块