计算化学公社

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

[Gaussian/gview] 如何定义胺分子umbrella motion的扫描坐标

[复制链接 Copy URL]

7

帖子

0

威望

41

eV
积分
48

Level 2 能力者

本帖最后由 Dwppppp 于 2018-5-16 23:45 编辑

大家好,我有一个关于柔性扫描胺分子umbrella motion如何定义坐标的问题。
我们知道三级胺分子离子态是平面结构,基态是三角锥,所以当被激发到激发态时会从三角锥转换为平面模式然后剩下的能量会导致氮原子和周围的三个碳沿着这个震动方向继续来回做震动,我想计算一下震动到不同位置的能量并且获得各个震动角度时最稳定的离子态结构。但我一直弄不明白怎么定义这个扫描的坐标。如图所示,我要研究该分子1,2,6,7四个原子的umbrella motion,我应该扫描什么呢?


Screen Shot 2018-05-16 at 11.33.33 AM.png (76.69 KB, 下载次数 Times of downloads: 50)

分子图片

分子图片

269

帖子

1

威望

3030

eV
积分
3319

Level 5 (御坂)

2#
发表于 Post on 2018-5-17 02:10:55 | 只看该作者 Only view this author
1-2-6-7的二面角?

6万

帖子

99

威望

6万

eV
积分
125153

管理员

公社社长

3#
发表于 Post on 2018-5-17 09:33:32 | 只看该作者 Only view this author
先尝试扫1-2-6-7二面角,如果发现达不到目的到时候再说。测试的时候可以先用半经验方法快速过一下
北京科音自然科学研究中心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

7

帖子

0

威望

41

eV
积分
48

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2018-5-19 01:19:38 | 只看该作者 Only view this author
谢谢楼上两位的回答! 确实是这个二面角没错,我之前好像扫错了顺序导致扫出奇怪的结构了。

我现在还有一个问题。这个二面角的初始角度是 -7.68 度, 我往 -2.5 度每步扫 没有问题,但是往+2.5度扫 2个结构之后就报错停了, 改成1.5度还是两个结构后就报错。是因为接近0了所以报错了吗?怎么能让我完成对于往另一个方向(+)的扫描呢?
报错信息是: >>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UB3LYP) =  -287.393093532     A.U. after  129 cycles
             Convg  =    0.3952D-02             -V/T =  1.9967

我大概知道这是因为结构不稳定所以SCF不收敛,但是0度附近应该是稳定的结构啊, 还是说二面角从负扫到正中间需要做什么操作才能跳过去。谢谢大家帮助!

6万

帖子

99

威望

6万

eV
积分
125153

管理员

公社社长

5#
发表于 Post on 2018-5-19 22:35:04 | 只看该作者 Only view this author
直接提取出现SCF不收敛对应的结构,看到底是怎么结构不合理
北京科音自然科学研究中心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

7

帖子

0

威望

41

eV
积分
48

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2018-5-21 05:39:25 | 只看该作者 Only view this author
我看了结构并没有什么不合理, 能量也还好, 我尝试用了scf=xqc和scf=qc,成功的算出了往后的10个结构,但还是报错了。现在从正常离子态的-7.8°扫到了+10° (二面角1267),停止的结构如图,实在不知道看着挺正常的结构为什么总报错。

Screen Shot 2018-05-20 at 5.35.01 PM.png (60.76 KB, 下载次数 Times of downloads: 62)

Screen Shot 2018-05-20 at 5.35.01 PM.png

6万

帖子

99

威望

6万

eV
积分
125153

管理员

公社社长

7#
发表于 Post on 2018-5-23 18:38:36 | 只看该作者 Only view this author
Dwppppp 发表于 2018-5-21 05:39
我看了结构并没有什么不合理, 能量也还好, 我尝试用了scf=xqc和scf=qc,成功的算出了往后的10个结构,但 ...

别太依赖scf=qc,又慢能解决问题几率又不高
尝试其它做法
解决SCF不收敛问题的方法
http://sobereva.com/61

这种结构SCF不会有多难收敛
北京科音自然科学研究中心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

7

帖子

0

威望

41

eV
积分
48

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2018-5-23 23:23:20 | 只看该作者 Only view this author
谢谢sob君,已经解决了收敛的问题。但是这次报错的是
Iteration 85 RMS(Cart)=  0.00000409 RMS(Int)=  0.02033888
Iteration 86 RMS(Cart)=  0.00000381 RMS(Int)=  0.02034269
Iteration 87 RMS(Cart)=  0.00000355 RMS(Int)=  0.02034624
Iteration 88 RMS(Cart)=  0.00000331 RMS(Int)=  0.02034955
Iteration 89 RMS(Cart)=  0.00000309 RMS(Int)=  0.02035263
Iteration 90 RMS(Cart)=  0.00000288 RMS(Int)=  0.02035551
Iteration 91 RMS(Cart)=  0.00000268 RMS(Int)=  0.02035819
Iteration 92 RMS(Cart)=  0.00000250 RMS(Int)=  0.02036068
Iteration 93 RMS(Cart)=  0.00000233 RMS(Int)=  0.02036301
Iteration 94 RMS(Cart)=  0.00000217 RMS(Int)=  0.02036518
Iteration 95 RMS(Cart)=  0.00000203 RMS(Int)=  0.02036721
Iteration 96 RMS(Cart)=  0.00000189 RMS(Int)=  0.02036909
Iteration 97 RMS(Cart)=  0.00000176 RMS(Int)=  0.02037085
Iteration 98 RMS(Cart)=  0.00000164 RMS(Int)=  0.02037249
Iteration 99 RMS(Cart)=  0.00000153 RMS(Int)=  0.02037402
Iteration100 RMS(Cart)=  0.00000143 RMS(Int)=  0.02037544
New curvilinear step not converged.
这个计算力场的迭代好像差一点点就可以了,有没有方法能增加最大步数呢。或者我还有什么问题导致这一步出现迭代不收敛的情况。

抱歉我太笨了 这么多蠢问题 QAQ

6万

帖子

99

威望

6万

eV
积分
125153

管理员

公社社长

9#
发表于 Post on 2018-5-23 23:29:54 | 只看该作者 Only view this author
Dwppppp 发表于 2018-5-23 23:23
谢谢sob君,已经解决了收敛的问题。但是这次报错的是
Iteration 85 RMS(Cart)=  0.00000409 RMS(Int)=  0. ...


据我所知没有可以控制这个的关键词,或许某不知名IOp能实现
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

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

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