计算化学公社

 找回密码 Forget password
 注册 Register
Views: 15835|回复 Reply: 15
打印 Print 上一主题 Last thread 下一主题 Next thread

[新手求助] 从只有一个虚频的初猜开始寻找过渡态,为什么收敛不了而且终点没有了虚频?

[复制链接 Copy URL]

235

帖子

1

威望

1176

eV
积分
1431

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
我为了寻找双芳基镍络合物转金属的过渡态,首先做了一个限制性优化,文件freez1。优化之后频率分析只有一个虚频,应该是对应着我需要的过渡态,用它作为初猜进行过渡态优化。但是一直无法收敛,而别结构逐渐偏离过渡态的结构,我取了最后一帧做了一个频率分析,发现虚频消失了。过渡态优化不是在保留虚频的情况下优化吗?为什么会出现这种情况?本人新手,第一次求助,这个结构已经能够优化两个星期了,希望大家有时间看一下。

freez1.zip

2.93 MB, 下载次数 Times of downloads: 53

235

帖子

1

威望

1176

eV
积分
1431

Level 4 (黑子)

16#
 楼主 Author| 发表于 Post on 2020-5-21 08:55:23 | 只看该作者 Only view this author
sobereva 发表于 2020-5-21 07:43
显然得自定义基组,仔细看
详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/ ...

好的,社长,我试一试

235

帖子

1

威望

1176

eV
积分
1431

Level 4 (黑子)

15#
 楼主 Author| 发表于 Post on 2020-5-21 08:54:59 | 只看该作者 Only view this author
sobereva 发表于 2020-5-21 07:43
你若没写genecp显然还是用的6-311G*算的

谢谢社长,我查看了一下输入文件,确实是这样,按照6-311*计算的

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

14#
发表于 Post on 2020-5-21 07:43:56 | 只看该作者 Only view this author
星斗如盘 发表于 2020-5-20 21:44
我使用的G09貌似不支持社长推荐的Lanl2TZ,显示无法识别Lanl2TZ。。。

显然得自定义基组,仔细看
详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/60
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

13#
发表于 Post on 2020-5-21 07:43:27 | 只看该作者 Only view this author
星斗如盘 发表于 2020-5-20 21:49
想再问您一个小问题,我定义了赝势基组,但是在输入计算方法和基组的时候错误的输入了b3lyp/6-311g(d), ...

你若没写genecp显然还是用的6-311G*算的
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

235

帖子

1

威望

1176

eV
积分
1431

Level 4 (黑子)

12#
 楼主 Author| 发表于 Post on 2020-5-20 21:49:58 | 只看该作者 Only view this author
shenchaoren 发表于 2020-5-20 09:19
**** 本内容被作者隐藏 ****

想再问您一个小问题,我定义了赝势基组,但是在输入计算方法和基组的时候错误的输入了b3lyp/6-311g(d),这样后面计算的结果是在赝势基组下计算的,还是全电子基组6-311g(d)下计算的呀?

235

帖子

1

威望

1176

eV
积分
1431

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2020-5-20 21:44:08 | 只看该作者 Only view this author
shenchaoren 发表于 2020-5-20 09:19
**** 本内容被作者隐藏 ****

我使用的G09貌似不支持社长推荐的Lanl2TZ,显示无法识别Lanl2TZ。。。

19

帖子

2

威望

1672

eV
积分
1731

Level 5 (御坂)

10#
发表于 Post on 2020-5-20 09:19:04 | 只看该作者 Only view this author
星斗如盘 发表于 2020-5-20 08:58
Br用什么赝势比较合理呢?和镍用同一种可以吗?

http://sobereva.com/373 看社长写的《谈谈赝势基组的选用》里面的 “2 什么时候该用赝势基组?”和“3 常用的几种赝势基组”,社长也特意标红了,“Lanl2DZ曾被用得非常广泛,即便是现在还用得很多。但笔者强烈不推荐用Lanl2DZ,因为在如今看来其精度明显偏低。”。 里面也写了关于计算配合物时赝势基组与全电子基组的匹配问题

235

帖子

1

威望

1176

eV
积分
1431

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2020-5-20 08:58:38 | 只看该作者 Only view this author
shenchaoren 发表于 2020-5-18 18:46
楼主的基组搭配不合理.既然Ni已经用了赝势,那么原子序数更大的Br为什么不用赝势而要用全电子基组6-31G(d)? ...

Br用什么赝势比较合理呢?和镍用同一种可以吗?

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

8#
发表于 Post on 2020-5-18 23:38:40 | 只看该作者 Only view this author
星斗如盘 发表于 2020-5-18 09:16
好的,谢谢社长,我尝试一下。我想在补充问一下,如果使用calcall,因为计算的特别慢,还有必要使用maxst ...

先尝试maxstep=10结合recalc=3,实在不行再calcall结合maxstep=10
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

19

帖子

2

威望

1672

eV
积分
1731

Level 5 (御坂)

7#
发表于 Post on 2020-5-18 18:52:47 | 只看该作者 Only view this author
星斗如盘 发表于 2020-5-18 09:16
好的,谢谢社长,我尝试一下。我想在补充问一下,如果使用calcall,因为计算的特别慢,还有必要使用maxst ...

http://sobereva.com/164 楼主可以看看社长写的

19

帖子

2

威望

1672

eV
积分
1731

Level 5 (御坂)

6#
发表于 Post on 2020-5-18 18:46:05 | 只看该作者 Only view this author
楼主的基组搭配不合理.既然Ni已经用了赝势,那么原子序数更大的Br为什么不用赝势而要用全电子基组6-31G(d)?建议两个Br上也用赝势

235

帖子

1

威望

1176

eV
积分
1431

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2020-5-18 09:16:25 | 只看该作者 Only view this author
本帖最后由 星斗如盘 于 2020-5-18 09:34 编辑
sobereva 发表于 2020-5-17 22:50
没什么新鲜的,初始结构满足这个条件又不代表优化过程中必定能保持这个条件,更何况你之前还对结构做了冻 ...

好的,谢谢社长,我尝试一下。我想在补充问一下,如果使用calcall,因为计算的特别慢,还有必要使用maxstep之类的关键字吗?

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

4#
发表于 Post on 2020-5-17 22:50:38 | 只看该作者 Only view this author
星斗如盘 发表于 2020-5-17 18:06
初猜做振动分析有一个虚频,但是优化之后就没有了

没什么新鲜的,初始结构满足这个条件又不代表优化过程中必定能保持这个条件,更何况你之前还对结构做了冻结,引入了人为因素

结合calcall或者recalc=x再试

写5d关键词毫无意义,何况genecp的时候默认就是用的5d
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

235

帖子

1

威望

1176

eV
积分
1431

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2020-5-17 18:06:45 | 只看该作者 Only view this author
liyuanhe211 发表于 2020-5-17 17:58
某结构是过渡态的条件有两个:有且只有一个虚频并且(无冻结的前提下)【梯度为0】。另外不要过于迷信限制 ...

初猜做振动分析有一个虚频,但是优化之后就没有了

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-21 05:44 , Processed in 0.178919 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list