计算化学公社

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

[综合交流] 关于SOS计算超极化率的问题

[复制链接 Copy URL]

40

帖子

1

威望

882

eV
积分
942

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
最近做第一超极化率的计算,碰到几个小问题,有点不明白,特来请教。
对于六十个原子左右(含CNHF)的共轭配合物体系,如果利用SOS方法计算其第一超极化率,选取CAM-B3LYP 泛函,按照培训讲义推荐的基组选用aug-cc-pVDZ。但是计算的体系中如若包含一个金属原子 In,似乎aug-cc-pVDZ并没有定义它,该换成什么基组?另外,对于不同的分子,结合Multiwfn计算出的β在偶极矩方向上的分量β||,有的是负值有的是正值比如-1000a.u.和300a.u.,这个正负号表示的是什么含义?如果是要比较其大小,是比较数值还是绝对值大小?
兵火有馀烬,贫村才数家。
无人争晓渡,残月下寒沙。

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

2#
发表于 Post on 2018-6-26 22:55:04 | 只看该作者 Only view this author
用aug-cc-pVDZ-PP
谈谈赝势基组的选用
http://sobereva.com/373

根据超极化率导数形式的定义理解。比绝对值

PS:Multiwfn 3.5的SOS功能有bug,如果用的3.5,改用最新的3.6(dev)或3.5之前的版本
北京科音自然科学研究中心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 2018-6-27 21:21:33 | 只看该作者 Only view this author
sobereva 发表于 2018-6-26 22:55
用aug-cc-pVDZ-PP
谈谈赝势基组的选用
http://sobereva.com/373

哦哦,好的谢谢老师!
我正是用的Multiwfn3.5版本计算的β值,那我重新计算一遍。
另,发现计算过程中生成S0S.txt文件这一步骤极度耗费时间,对于正在计算的体系,笔记本要跑三四十小时甚至更久才能完成,除了增大settings.ini中的nthreads,不知有没有其他加快生成文件速度的方法?
兵火有馀烬,贫村才数家。
无人争晓渡,残月下寒沙。

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

4#
发表于 Post on 2018-6-28 02:07:09 | 只看该作者 Only view this author
cornerj 发表于 2018-6-27 21:21
哦哦,好的谢谢老师!
我正是用的Multiwfn3.5版本计算的β值,那我重新计算一遍。
另,发现计算过程中 ...


我刚刚更新了主页上的3.6(dev)版,这个版本在生成SOS.txt的过程中把非常耗时的一步(产生分子轨道间的跃迁偶极矩)的耗时降低了甚巨,可以试试。如果还慢,那就只能用更小点基组,或者把IOp(9/40=x)的x设更小点(但不能小于3),代码已经没法进一步优化了。
北京科音自然科学研究中心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 2018-6-28 23:30:33 | 只看该作者 Only view this author
sobereva 发表于 2018-6-28 02:07
我刚刚更新了主页上的3.6(dev)版,这个版本在生成SOS.txt的过程中把非常耗时的一步(产生分子轨道间的 ...

十分感谢老师的指点!今天尝试了下用3.6(dev)版重新生成S0S.txt文件,速度快的飞起!太爽了!
大致耗时仅仅半小时至一小时而已,极度省时了!
麻烦再问点问题,对比了一下3.5版和3.6(dev)版得到的数据,对于同样的.fch和.log文件,体系的极化率α数值是可重现的,一模一样,
但β|| 的数值却有较大差别,如有略微增大几百a.u.的,也有暴减的,从之前+300a.u.变为-20a.u,对这个变化很是疑惑,不知是何原因?
兵火有馀烬,贫村才数家。
无人争晓渡,残月下寒沙。

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

6#
发表于 Post on 2018-6-29 00:12:35 | 只看该作者 Only view this author
cornerj 发表于 2018-6-28 23:30
十分感谢老师的指点!今天尝试了下用3.6(dev)版重新生成S0S.txt文件,速度快的飞起!太爽了!
大致耗 ...

3.5版的alpha没问题,beta和delta是错的
北京科音自然科学研究中心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 2018-6-29 19:40:37 | 只看该作者 Only view this author
sobereva 发表于 2018-6-29 00:12
3.5版的alpha没问题,beta和delta是错的

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

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 06:26 , Processed in 0.156301 second(s), 20 queries , Gzip On.

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