计算化学公社

标题: 请问要如何优化出局域极小值点 [打印本页]

作者
Author: 真步 时间: 2022-11-20 09:54
标题: 请问要如何优化出局域极小值点
(, 下载次数 Times of downloads: 2) 刚开始学优化，尝试着用高斯view构造的初始结构进行优化，但是每次优化最终都收敛在Anti构象。#p b3lyp/6-31g gfinput opt=(maxcyc=500) scf=(maxcyc=500,conver=8)
这是优化出来的结构，请问如果想优化到Syn构象，需要学习或者注意些什么呢，求大佬指导。
(, 下载次数 Times of downloads: 2)

作者
Author: luchen 时间: 2022-11-20 10:11
建模的时候最大限度让模型接近于想得到的结构，如果想找能量最小点的结构可以对于柔性分子可以进行构象搜索

作者
Author: exity 时间: 2022-11-20 10:53
一个E式一个Z式
建模的时候通过二面角调整初始结构，一般情况下如果有极小值点的话，opt过程会落到这个点的。

作者
Author: wjc404 时间: 2022-11-20 12:39
本帖最后由 wjc404 于 2022-11-20 12:55 编辑

https://www.sciencedirect.com/sc ... i/S0022286001005932对二苯乙二酮的OCCO二面角进行了扫描，只有一个最小点，是介于顺式和反式之间的扭曲式。

用B3LYP-D3/def2-SVP对二苯乙二酮的OCCO二面角做柔性扫描，得到下图，只有一个最低点。
(, 下载次数 Times of downloads: 3)
这也比较好理解，二苯乙二酮的两个PhCO基团之间的共轭本来就不强（平面化倾向不强），顺式时两个苯环空间排斥太强，反式时一个PhCO的2号位氢和另一个PhCO的氧有空间排斥。

作者
Author: sobereva 时间: 2022-11-20 13:52
甭用这种垃圾基组，仔细看下文之前不要做任何计算，白浪费时间
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）

甭写scf=(maxcyc=500,conver=8) gfinput 这种垃圾关键词，好好看
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331（http://bbs.keinsci.com/thread-3460-1-1.html）

Syn型分子内色散作用不可忽视，用B3LYP必须带上D3，仔细看下文
谈谈量子化学研究中什么时候用B3LYP泛函优化几何结构是适当的
http://sobereva.com/557（http://bbs.keinsci.com/thread-17899-1-1.html）
谈谈“计算时是否需要加DFT-D3色散校正？”
http://sobereva.com/413（http://bbs.keinsci.com/thread-9772-1-1.html）

建议把下文也看了，免得老犯低级错误
谈谈学量子化学如何正确地入门
http://sobereva.com/355（http://bbs.keinsci.com/thread-4447-1-1.html）

作者
Author: 真步 时间: 2022-11-20 19:56

sobereva 发表于 2022-11-20 13:52
甭用这种垃圾基组，仔细看下文之前不要做任何计算，白浪费时间
谈谈量子化学中基组的选择
http://soberev ...

谢谢sob老师!

作者
Author: 真步 时间: 2022-11-20 20:17

wjc404 发表于 2022-11-20 12:39
https://www.sciencedirect.com/sc ... i/S0022286001005932对二苯乙二酮的OCCO二面角进行了扫描，只有一个 ...

非常感谢，但是还有一个问题，我在G16上使用B3LYPD3/def2-SVP关键词进行尝试，结果在LINK1就出现报错，用b3lyp/6-31g却能正常计算。
我特别想重复一下大佬的柔性扫描，想请教一下是哪个环节出现了问题，像这样的计算会使用多少计算资源。

%chk=pi2.chk
%mem=4GB
%nprocs=4
Will use up to 4 processors via shared memory.
-----------------------
#p B3LYPD3/def2-SVP opt
-----------------------
QPErr --- A syntax error was detected in the input line.
#p B3LYPD3/def2-SVP opt

Last state= "GCL"
TCursr=    3918 LCursr=       3
Error termination via Lnk1e in /home/software/g16A03/g16/l1.exe at Sun Nov 20 19:55:50 2022.
Job cpu time:    0 days  0 hours  0 minutes  0.1 seconds.
Elapsed time:    0 days  0 hours  0 minutes  0.0 seconds.
File lengths (MBytes):  RWF=    6 Int=    0 D2E=    0 Chk=    1 Scr=    1

作者
Author: wjc404 时间: 2022-11-20 20:48
(1)我不是大佬，我只是个学生

。
(2)在http://sobereva.com/210的后半段，社长非常详细地写了怎么在高斯中规范地使用DFT-D。
(3)在https://gaussian.com/basissets/上面有高斯能识别的各个基组关键词的写法。

作者
Author: rugals 时间: 2022-11-20 20:51
本帖最后由 rugals 于 2022-11-20 20:52 编辑

真步发表于 2022-11-20 20:17
非常感谢，但是还有一个问题，我在G16上使用B3LYPD3/def2-SVP关键词进行尝试，结果在LINK1就出现报错，用 ...

你这么写关键词高斯不认识
D3是写成em=gd3 def2-svp的横杠去掉

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