拉曼光谱加弥散时，怎么计算极化率导数

wangchunchun · 发表于 Post on 2020-9-7 12:17:35

在老师量子化学中关于基组的选择方法这篇文章中，我看到了关于拉曼光谱的一段话，“J. Chem. Theory Comput., 7, 3323 (2011)提出计算Raman、ROA的两步方法，在优化和振动分析时不用弥散函数（低至6-31G*高至def-TZVP就够），但是计算极化率导数的时候再考虑弥散函数（建议>= aug-cc-pVDZ）”

其中我不是很懂怎么计算极化率导数的时候考虑弥散函数，我所做的需要计算分子的拉曼光谱
我是先进行了结构优化
#p b3lyp/6-311g** opt(gdiis,maxstep=3,notrust,recalc=5) EmpiricalDispersion=GD3BJ out=wfn
然后需要进行拉曼计算，需要加弥散函数，但是看文章说振动时方法必须和优化一致，所以拉曼也只能用b3lyp/6-311g**，即#p b3lyp/6-311g** empiricaldispersion=gd3bj freq=raman
请问老师我想在算拉曼时加弥散函数应该怎么做呢

sobereva · 发表于 Post on 2020-9-7 13:07:07

至少截止到G16 C.01为止，两步法算Raman的功能都有bug，没法直接用。但两步法算ROA的功能正常，同时也会输出Raman

wangchunchun · 发表于 Post on 2020-9-7 19:42:15

谢谢老师，十分感谢sob老师，祝老师工作顺利

indec · 发表于 Post on 2020-9-15 01:04:34

sobereva 发表于 2020-9-7 13:07
至少截止到G16 C.01为止，两步法算Raman的功能都有bug，没法直接用。但两步法算ROA的功能正常，同时也会输 ...

社长，ORCA 能用类似的方法两步法算拉曼吗？最近有个体系需要算拉曼，体系偏大，高斯要算很久很久。不知道ORCA 能否降低一些耗时。

sobereva · 发表于 Post on 2020-9-16 09:24:52

indec 发表于 2020-9-15 01:04
社长，ORCA 能用类似的方法两步法算拉曼吗？最近有个体系需要算拉曼，体系偏大，高斯要算很久很久。不知 ...

ORCA应该不行。而且ORCA的振动分析速度远不如它算能量和受力，到了算Hessian的场合，相对于Gaussian的速度优势会失去大半，甚至有时还可能更慢。另外，ORCA用DFT算Raman应该是基于有限差分算的Raman活性（我没仔细证实过），因为这涉及到三阶导数，从这点上也不如Gaussian，Gaussian可以解析地计算。

wangchunchun · 发表于 Post on 2020-10-20 09:16:13

sobereva 发表于 2020-9-7 13:07
至少截止到G16 C.01为止，两步法算Raman的功能都有bug，没法直接用。但两步法算ROA的功能正常，同时也会输 ...

老师请问在第一步中先优化完然后再计算振动freq可以吗然后再进行第二步

sobereva · 发表于 Post on 2020-10-21 11:38:03

wangchunchun 发表于 2020-10-20 09:16
老师请问在第一步中先优化完然后再计算振动freq可以吗然后再进行第二步

可

wangchunchun · 发表于 Post on 2020-11-2 11:05:26

sobereva 发表于 2020-9-7 13:07
至少截止到G16 C.01为止，两步法算Raman的功能都有bug，没法直接用。但两步法算ROA的功能正常，同时也会输 ...

老师，我用您的ROA例子试着算了一下，加了色散和溶剂，但是在加溶剂的情况下就算画不出拉曼的图像，请问老师这个我应该怎么解决呢
optfreq：#p B3LYP/6-31G* opt freq EmpiricalDispersion=GD3BJ scrf=(smd,solvent=n-hexane)
polar：#p B3LYP/aug-cc-pVDZ polar=ROA guess=read geom=check EmpiricalDispersion=GD3BJ scrf=(smd,solvent=n-hexane)

sobereva · 发表于 Post on 2020-11-2 13:19:51

wangchunchun 发表于 2020-11-2 11:05
老师，我用您的ROA例子试着算了一下，加了色散和溶剂，但是在加溶剂的情况下就算画不出拉曼的图像，请问 ...

这和色散校正没关系。
看看输出文件中拉曼活性有没有输出

wangchunchun · 发表于 Post on 2020-11-2 20:48:45

sobereva 发表于 2020-11-2 13:19
这和色散校正没关系。
看看输出文件中拉曼活性有没有输出

老师我用gview打开out文件是有拉曼活性的 raman activity ，但是用Multiwfn画图画不出来，不加溶剂计算的结果就能画出来

sobereva · 发表于 Post on 2020-11-4 14:13:48

wangchunchun 发表于 2020-11-2 20:48
老师我用gview打开out文件是有拉曼活性的 raman activity ，但是用Multiwfn画图画不出来，不加 ...

仔细看到底运行正常结束了没有。
说不清楚就把输出文件压缩后上传，我看看

大头攒毛 · 发表于 Post on 2023-9-8 16:34:26

请教一下，拉曼活性和极化率以及极化率导数之间有什么关系吗？有没有一组公式可以说明？谢谢

sobereva · 发表于 Post on 2023-9-8 22:16:49

大头攒毛发表于 2023-9-8 16:34
请教一下，拉曼活性和极化率以及极化率导数之间有什么关系吗？有没有一组公式可以说明？谢谢

参考北京科音初级量子化学培训班（http://www.keinsci.com/workshop/KEQC_content.html）ppt

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 拉曼光谱加弥散时，怎么计算极化率导数