计算化学公社

 找回密码 Forget password
 注册 Register
Views: 6326|回复 Reply: 3

[Multiwfn功能更新] Multiwfn的SOS功能现已可以对动态超极化率beta(w1,w2)的频率做二维扫描

[复制链接 Copy URL]

4万

帖子

99

威望

4万

eV
积分
89888

管理员

公社社长+计算化学玩家

发表于 Post on 2019-11-21 01:55:26 | 显示全部楼层 Show all |阅读模式 Reading model
之前在《使用Multiwfn基于完全态求和(SOS)方法计算极化率和超极化率》(http://sobereva.com/232)里笔者详细介绍了Multiwfn的SOS模块的功能和用法。今天在思想家公社群里有人问怎么绘制ACS Appl. Nano Mater. 2019, 2, 1648−1654这篇文章里下面这种图,即对beta(-(w1+w2);w1,w2)的外场频率w1和w2进行二维扫描。利用这种图的峰的位置,可以考察体系的第一超极化率特征。

beta二维频率扫描例子(ACS Appl Nano Mater,2019)-Spiral Graphene Nanoribbons with .jpg
图中X=0和Y=0的位置对应单光子共振,w1与w2符号相同的区域对应合频生成,符号相反区域对应差频生成。
如果w1(w2)为0时在某个w2(w1)频率下出现了峰,体现出相应频率下的electro-optical Pockels效应;如果w1=w2为某个频率时出现了峰,体现出相应频率下的SHG效应;如果w1=-w2为某个频率时出现了峰,体现出相应频率下的optical rectification效应。

为了能绘制这种图,笔者把Multiwfn的SOS功能更新了一下,现已可以在官网http://sobereva.com/multiwfn下载。

用Multiwfn照常载入SOS模块计算超极化率要用的输入文件,之后在SOS计算模块中,可以看到选项19 Scanning w1 and w2 of beta(-(w1+w2);w1,w2)。进入之后,输入w1的初值、终值、扫描的步数,再输入输入w2的初值、终值、扫描的步数,经过计算后(w1和w2都取100步其实也花不了太多时间),就可以在当前目录下得到beta_w.txt和beta_w_comp.txt,如屏幕上的提示所示,内容是各种形式的beta以及beta的各个分量随w1和w2的变化。之后把其中之一载入比如Sigmaplot、Origin之类的程序,将w1和w2作为x和y轴,再将某个beta作为z轴绘制地形图即可,例如:


Clipboard01.png


一个具体计算例子看最新版手册4.200.8.1节。

评分 Rate

参与人数
Participants 5
eV +24 收起 理由
Reason
leebo + 5 好物!谢谢社长这么快速的更新这个功能
Novice + 4 とてもいい!
朙天儿 + 5 牛!
zsu007 + 5 赞!
liyuanhe211 + 5 好物!

查看全部评分 View all ratings

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班,内容介绍以及往届资料购买请点击链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的最佳途径。培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人,讨论范畴相同
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

23

帖子

0

威望

2713

eV
积分
2736

Level 5 (御坂)

发表于 Post on 2019-12-12 12:16:30 | 显示全部楼层 Show all
本帖最后由 leebo 于 2019-12-12 12:20 编辑

社长,要是能直接输出electro-optical Pockels效应,SHG效应,optical rectification效应发生的具体频率就好了。

4万

帖子

99

威望

4万

eV
积分
89888

管理员

公社社长+计算化学玩家

 楼主 Author| 发表于 Post on 2019-12-15 10:31:38 | 显示全部楼层 Show all
leebo 发表于 2019-12-12 12:16
社长,要是能直接输出electro-optical Pockels效应,SHG效应,optical rectification效应发生的具体频率就好 ...

这些效应无处不在,只是不同频率下的大小的问题,而且并非只有一个极值点
通过excel产生对应的扫描的w1,w2频率列表,在Multiwfn里做个频率的扫描即可
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班,内容介绍以及往届资料购买请点击链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的最佳途径。培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人,讨论范畴相同
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

100

帖子

0

威望

1348

eV
积分
1448

Level 4 (黑子)

PIG

发表于 Post on 2020-1-19 18:16:43 | 显示全部楼层 Show all
sobereva 发表于 2019-12-15 10:31
这些效应无处不在,只是不同频率下的大小的问题,而且并非只有一个极值点
通过excel产生对应的扫描的w1, ...

同意,只要有这功能,搞Hyper-Raman, CARS, SRS都是有可能的

本版积分规则 Credits rule

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

GMT+8, 2023-2-2 21:50 , Processed in 0.215877 second(s), 25 queries .

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