请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 468|回复: 2

[Multiwfn资源与经验] 使用Multiwfn计算与超瑞利散射(HRS)实验相关的量

[复制链接]

1万

帖子

25

威望

2万

eV
积分
44706

管理员

公社社长

发表于 2019-8-7 12:27:23 | 显示全部楼层 |阅读模式
使用Multiwfn计算与超瑞利散射(HRS)实验相关的量

文/Sobereva@北京科音  2019-Aug-7


摘要:超瑞利散射(HRS)是考察分子第一超极化率的一种实验手段,本文介绍如何使用Multiwfn基于Gaussian的polar任务的输出文件计算出与HRS实验相关的量。首先介绍基础知识,然后给出具体例子,并重复文献里的数据。


1 基本知识

之前笔者在《使用Multiwfn分析Gaussian的极化率、第一超极化率的输出》(http://sobereva.com/231)中已经对(超)极化率的相关知识和计算有过基本的介绍。实验上测定beta(第一超极化率)常用的方法是EFISHG(电场诱导二次谐波产生)实验,这个测的是平行于分子偶极矩的beta(-2w;w,w)的分量,这里w是外场频率。还有一种实验上考察beta(-2w;w,w)的方法是超瑞利散射(Hyper-Rayleigh Scattering, HRS),相关介绍见Acc. Chem. Res., 31, 675 (1998)(此文有个别错误)。HRS实验示意图如下。
1.png

根据入射光的强度和测定的垂直散射出的倍频光的强度,可以得到beta_HRS值,表达式如下,它是对分子的beta(-2w;w,w)的一种展现,是各向平均属性,不受分子朝向的影响。其中<beta^2 ZZZ>和<beta^2 XZZ>的具体定义看Multiwfn手册3.200.7节或PCCP,10,6223(2008),是基于beta(-2w;w,w)张量的各个分量计算的。还有个量叫退偏比(depolarization ratio, DR),定义也在下面给出了,理论最小值为1.5(对应Td点群的分子)。
2.png

<beta^2 ZZZ>和<beta^2 XZZ>可分解为偶极(J=1)和八极(J=3)成份的加和,如下所示。并进而可以定义各向异性参数ρ,以及偶极和八极对beta的相对贡献φ(beta J=1)和φ(beta J=3)。
3.png

对于小分子,其偶极矩大小与上述很多量有密切联系,本文就不多说了,读者可参看J. Chem. Phys., 136, 024506 (2012)里的例子。

入射光的极化态通过(Ψ,δ)两个参数确定。当入射光以X方向射入、相位延迟量δ=π/2时,顺着Y方向散射出的沿Z方向极化的光的强度与入射光极化角Ψ的关系为:
4.png

很多文章里都考察了散射光强度随Ψ的变化关系,比如JCP,136,024506(2012)、JPCA,116,10249(2012)。


2 用法&例子

下面就结合例子说一下用Multiwfn具体怎么计算前述的与HRS相关的量。注意从2019年7月30日更新的Multiwfn 3.7(dev)开始才有相关功能,绝对不要用在此日期之前发布的版本!!!Multiwfn可在其主页http://sobereva.com/multiwfn免费下载。

计算HRS相关的量需要提供Gaussian的polar=DCSHG任务的输出文件。此处我们重复一下JCP,136,024506(2012)中对CH2Cl2计算的结果。文中给了很多数据,我们这里将要对照文中在HF/d-aug-cc-pVTZ级别下结合IEFPCM描述CH2Cl2溶剂环境时算的数据,入射波长取的是1064nm。Gaussian的输入文件如下。

#p HF/daug-cc-pVTZ polar=DCSHG CPHF=rdfreq scrf=solvent=CH2Cl2
[空行]
Title Card Required
[空行]
0 1
C                  0.00000000    0.00000000    0.76822500
H                 -0.89768500    0.00000000    1.37593000
H                  0.89768500    0.00000000    1.37593000
Cl                 0.00000000    1.49574000   -0.21650600
Cl                 0.00000000   -1.49574000   -0.21650600
[空行]
1064nm

这个例子里,#P必须有,否则Multiwfn无法解析输出文件。d-aug-cc-pVTZ在Gaussian里写作daug-cc-pVTZ。由于HRS展现的是分子的beta(-2w;w,w),显然必须用polar=DCSHG关键词来计算这种形式的beta。CPHF=rdfreq结合末尾空一行处的1064nm代表要算的不仅有静态极限w=0的情况,还有w=1064nm的情况(由于polar=DCSHG关键词默认了CPHF=rdfreq,因此CPHF=rdfreq也可以不写)。scrf=solvent=CH2Cl2代表用默认的IEFPCM隐式模型描述CH2Cl2溶剂环境。上面输入文件里的坐标是笔者事先使用与JCP文章里相同的级别优化好的。

上面的输入文件和G09 E.01计算产生的输出文件可以直接在此下载:http://sobereva.com/attach/499/file.rar

启动Multiwfn,载入上述任务的输出文件CH2Cl2_solv.out,然后依次输入
200  //主功能200
7  //解析Gaussian的polar任务的输出文件里的(超)极化率信息并计算相关数据
-1   //切换为载入含频(超)极化率的模式
1   //当前用的计算级别是有解析三阶导数的级别,包括HF、普通DFT泛函、半经验
2   //载入对应w=1064nm的数据
2   //考察的形式是beta(-2w;w,w)

此时在屏幕上看到下面的输出

[beta各个分量略...]
Beta_X=        0.00000  Beta_Y=        0.00000  Beta_Z=       41.77407
Magnitude of first hyperpolarizability:       41.774067
Projection of beta on dipole moment:       41.774067
Beta ||     :       25.064440
Beta ||(z)  :       25.064440
Beta _|_(z) :        8.334480

Note: Kleinman's symmetry condition is not employed for below quantities:
<beta_ZZZ^2>:  3.32226140E+02
<beta_XZZ^2>:  1.05032972E+02
Hyper-Rayleigh scattering (beta_HRS):        20.911
Depolarization ratio (DR):   3.163
|<beta J=1>|:        32.374
|<beta J=3>|:        46.322
Nonlinear anisotropy parameter (rho):   1.431
Dipolar contribution to beta, phi_beta(J=1):     0.411
Octupolar contribution to beta, phi_beta(J=3):   0.589
< (beta_ZXZ+beta_ZZX)^2 - 2*betaZZZ*betaZXX >:  4.03004664E+02

前文提到的beta_HRS值、<beta_ZZZ^2>、DR等各个量都已经明确给出来了。由于计算前没要求程序切换成其它单位输出,所以单位用的都是默认的a.u.。

在前述的JCP文章的表5中的最后一行我们看到在当前计算条件下beta_HRS值为20.90,DR值为3.12,我们自己算的结果分别是20.911和3.16,可见与文献高度相符。

如果你想考察散射光强度随Ψ的变化,此时输入y,然后输入一个起始角度。如果输入的值为r,那么就会给出Ψ从r到r+360度范围内每隔1度对应的散射光强度值。这里我们输入-180,然后屏幕上提示HRS_angle.txt已经导出到了当前目录下。内容为:
-180.0  105.032972
-179.0  105.091742
-178.0  105.268017
...略
177.0  105.561698
178.0  105.268017
179.0  105.091742

之后可以把这个文件拖入到Origin里并绘制成曲线图。由于当前的数据是360度周期性的,很多研究者通过图像展现这类数据的时候喜欢绘制成极坐标的图。具体做法是选择Plot - Polar theta(X)r(Y),然后对图像做一番调整,就可以得到下图的效果。相应的Origin的.opj文件在上面提到的文件包里直接提供了。
5.png

此图很容易理解,角度对应Ψ值,每个特定角度下红线距离原点的距离就是这个Ψ值对应的散射强度值。在前述JCP文章的图9里,给了与上图相对应的图,曲线图也一并给出,如下所示。其中右图中绿色长虚线对应的是我们上面用的计算级别,可见我们的图和文献里的图完全一致。

6.png


顺带一提,之前密度泛函小卒写过一个基于Python的程序NLO Calculator也可以计算上述这些量,Multiwfn计算的结果和NLO Calculator是一样的,但输出格式明显不同。

评分

参与人数 4eV +20 收起 理由
aqhuangry + 5 赞!
asdf + 5 赞!
Novice + 5 社长真是无所不能啊!
zsu007 + 5 牛!

查看全部评分

北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

81

帖子

1

威望

420

eV
积分
521

Level 4 (黑子)

发表于 前天 17:35 | 显示全部楼层
老师,能根据 https://doi.org/10.1002/anie.198406901这篇文献中的取向气体模型近似用超极化率β计算宏观二阶系数,在multiwfn中实现计算吗?

1万

帖子

25

威望

2万

eV
积分
44706

管理员

公社社长

 楼主| 发表于 昨天 22:21 | 显示全部楼层
tjuptz 发表于 2019-8-19 17:35
老师,能根据 https://doi.org/10.1002/anie.198406901这篇文献中的取向气体模型近似用超极化率β计算宏观 ...

目前无此功能
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2019-8-21 20:40 , Processed in 0.247801 second(s), 28 queries .

快速回复 返回顶部 返回列表