计算化学公社

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

[综合交流] 态求和计算150states还是没有呈现收敛趋势,需要继续增加态数嘛?

[复制链接 Copy URL]

40

帖子

1

威望

882

eV
积分
942

Level 4 (黑子)


请问态求和算β0值,计算150states,用Multiwfn导出来来的β0与态数的关系图,还是没有明显的收敛趋势,还有必要继续增加态数吗?

使用计算水平CAM-B3LYP/aug-cc-pVDZ,体系原子数60多个,再增加态数计算耗时太大估计算不动了。请问该如何处理这种情况?

兵火有馀烬,贫村才数家。
无人争晓渡,残月下寒沙。

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

2#
发表于 Post on 2019-4-15 09:43:02 | 只看该作者 Only view this author
只能再增加,没办法。这是SOS方法的缺点。估计到了两百多个态总应该能收敛了
北京科音自然科学研究中心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

40

帖子

1

威望

882

eV
积分
942

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2019-4-15 21:52:15 | 只看该作者 Only view this author
sobereva 发表于 2019-4-15 09:43
只能再增加,没办法。这是SOS方法的缺点。估计到了两百多个态总应该能收敛了

哦哦那我老老实实算个200试试看。谢谢老师!
还有个小问题,从图上看出来第3,4个激发态有很明显的正、负贡献,
怎么用Multiwfn导出这几个激发态对应的轨道跃迁组成?我在TD的输出文件里找不到轨道跃迁信息、
兵火有馀烬,贫村才数家。
无人争晓渡,残月下寒沙。

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

4#
发表于 Post on 2019-4-16 00:08:08 | 只看该作者 Only view this author
cornerj 发表于 2019-4-15 21:52
哦哦那我老老实实算个200试试看。谢谢老师!
还有个小问题,从图上看出来第3,4个激发态有很明显的正、负 ...

轨道跃迁信息肯定有,不然Multiwfn都没法产生用于做SOS的文件
如果你嫌输出量太大不好找主要的,可以用这个程序处理
简化用了IOp(9/40=x)的Gaussian的CIS/TDDFT任务的输出文件的程序:simpIOp940
http://sobereva.com/467http://bbs.keinsci.com/thread-12188-1-1.html

在Multiwfn里也可以直接输出某个电子激发中主要的轨道跃迁,见
电子激发任务中轨道跃迁贡献的计算
http://sobereva.com/230
北京科音自然科学研究中心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

40

帖子

1

威望

882

eV
积分
942

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2019-4-16 12:52:20 | 只看该作者 Only view this author
sobereva 发表于 2019-4-16 00:08
轨道跃迁信息肯定有,不然Multiwfn都没法产生用于做SOS的文件
如果你嫌输出量太大不好找主要的,可以用 ...

尝试了下,输出特别迅速!谢谢老师!

Some MO transitions sorted by absolute contributions:
#  43327    197B ->  208B   Coeff.:   0.53203   Contri.:   28.3056%
#  23226    210A ->  213A   Coeff.:  -0.37279   Contri.:   13.8972%

比如开壳层体系,对应的轨道跃迁197B ->  208B  210A ->  213A
   我想知道这两个跃迁涉及到的分子轨道 197B 208B 210A 213A 各自的分子轨道图是哪几个?
   因为用Multiwfn导出来的分子轨道的序号 :206 is HOMO of alpha spin, orbital  1045 is HOMO of beta spin
   请问这两种序号是存在对应关系的嘛?谢谢!


兵火有馀烬,贫村才数家。
无人争晓渡,残月下寒沙。

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

6#
发表于 Post on 2019-4-17 00:45:47 | 只看该作者 Only view this author
cornerj 发表于 2019-4-16 12:52
尝试了下,输出特别迅速!谢谢老师!

Some MO transitions sorted by absolute contributions:

使用Multiwfn观看分子轨道
http://sobereva.com/269http://bbs.keinsci.com/thread-462-1-1.html

这里说了,主功能0里面在文本框里输入负号就代表选择beta轨道,正值就代表alpha
北京科音自然科学研究中心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

40

帖子

1

威望

882

eV
积分
942

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2019-4-17 07:43:27 | 只看该作者 Only view this author
sobereva 发表于 2019-4-17 00:45
使用Multiwfn观看分子轨道
http://sobereva.com/269(http://bbs.keinsci.com/thread-462-1-1.html)

...

明白这个操作了,感谢老师深夜解答!
有个莫名其妙的问题,优化和SOS一低一高俩计算级别 都产生了fchk文件,
想看SOS中激发态对应的轨道跃迁,用哪一个fchk?这俩fchk包含的轨道信息有没有区别?
兵火有馀烬,贫村才数家。
无人争晓渡,残月下寒沙。

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

8#
发表于 Post on 2019-4-17 13:10:33 | 只看该作者 Only view this author
cornerj 发表于 2019-4-17 07:43
明白这个操作了,感谢老师深夜解答!
有个莫名其妙的问题,优化和SOS一低一高俩计算级别 都产生了fchk文 ...

显然用你做电子激发计算任务产生的fch文件
北京科音自然科学研究中心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

40

帖子

1

威望

882

eV
积分
942

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2019-4-17 13:17:45 | 只看该作者 Only view this author
sobereva 发表于 2019-4-17 13:10
显然用你做电子激发计算任务产生的fch文件

了解了,感谢老师的解答!
兵火有馀烬,贫村才数家。
无人争晓渡,残月下寒沙。

3

帖子

0

威望

214

eV
积分
217

Level 3 能力者

10#
发表于 Post on 2021-1-13 17:40:39 | 只看该作者 Only view this author
本帖最后由 Yuanyaya 于 2021-1-13 17:47 编辑

sob老师您好,我看了您使用Multiwfn基于完全态求和(SOS)方法计算极化率和超极化率(http://sobereva.com/232)的帖子,用到自己体系分析γ的时候,有时在100个态左右就收敛的很好,有时会出现算200个态仍不收敛的情况(如图),请问对于不收敛的可以怎么修改呢?
计算td用的cam-b3lyp/genecp em=gd3 td(nstates=200) IOp(9/40=5) nosymm,C H 用6-31+g(d,p),Sn 用 Lanl2DZ,大概90个原子。C:\Users\Yuanyaya\Desktop

Graph1.jpg (110.76 KB, 下载次数 Times of downloads: 32)

Graph1.jpg

1245

帖子

1

威望

3860

eV
积分
5125

Level 6 (一方通行)

11#
发表于 Post on 2021-1-13 19:39:51 | 只看该作者 Only view this author
Yuanyaya 发表于 2021-1-13 17:40
sob老师您好,我看了您使用Multiwfn基于完全态求和(SOS)方法计算极化率和超极化率(http://sobereva.com/23 ...

不收敛就加多态数目继续算

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

12#
发表于 Post on 2021-1-16 08:02:37 | 只看该作者 Only view this author
Yuanyaya 发表于 2021-1-13 17:40
sob老师您好,我看了您使用Multiwfn基于完全态求和(SOS)方法计算极化率和超极化率(http://sobereva.com/23 ...

6-31+g(d,p)算(超)极化率烂得一塌糊涂。用aug-cc-pVDZ还勉勉强强
多少态能收敛和体系关系很大,算到几百个态才收敛都有可能。
北京科音自然科学研究中心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-24 07:13 , Processed in 0.172411 second(s), 23 queries , Gzip On.

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