计算化学公社

标题: 求助CFOUR中如何从偶极矩对核坐标一阶导数以及偶极矩计算红外强度 [打印本页]

作者
Author: liurui 时间: 2023-10-24 09:08
标题: 求助CFOUR中如何从偶极矩对核坐标一阶导数以及偶极矩计算红外强度
本帖最后由 liurui 于 2023-10-24 19:37 编辑

各位老师！学生在研究红外强度，得知红外强度与偶极矩与核坐标的一阶偏微分的平方成正比，以最简单的水二聚体为例，总偶极矩是沿着OO方向的(如下图)，如果想研究形成氢键的OH振动拉伸的红外强度如何得到的，是不是就是将总偶极矩投影到OH方向的分量，也就是总偶极矩乘以cos夹角，再与OH键长之间进行拟合多项式，然后进行求导？
亦或者根据CFOUR输出文件中Total dipole moment derivatives，但是这个貌似是一个矩阵，如何得到红外强度的呢

作者
Author: ahxb 时间: 2023-10-24 13:04
你的描述中有如下问题：
（1）”红外强度与偶极矩与核坐标的一阶偏微分的平方成正比“不准确，实际上是与偶极矩对相应本征振动模式的导数的平方成正比，或者写成 I 正比于 |dμ/dq|^2，其中 q 是 3N 维矢量。
（2）从对称性可知总偶极矩并非沿着O-O方向，但确实很接近。
（3）对于小分子而言，并不一定存在对应某个键伸缩的振动模式，通常而言是分子整体的振动，因此H2O二聚体中并不一定存在对应于O-H...O氢键中OH的振动拉伸，即使找到了，通常也并非只对应于O…H键的伸缩。所以你提到的将偶极矩投影再拟合曲线求导是不可行的做法。

作者
Author: liurui 时间: 2023-10-24 19:45

ahxb 发表于 2023-10-24 13:04
你的描述中有如下问题：
（1）”红外强度与偶极矩与核坐标的一阶偏微分的平方成正比“不准确，实际上是与 ...

老师，主要是之前看了最后一张照片中一篇jpcl的文章，看到OH键长和偶极矩之间存在线性关系，说红外强度偶极矩与OH键长的一阶偏微分的平方成正比，所以还是存在疑问，谢谢老师！！！非常感谢！！！！

作者
Author: liurui 时间: 2023-10-24 19:46

ahxb 发表于 2023-10-24 13:04
你的描述中有如下问题：
（1）”红外强度与偶极矩与核坐标的一阶偏微分的平方成正比“不准确，实际上是与 ...

文章是 J. Phys. Chem. Lett. 2019, 10, 918−924

作者
Author: ahxb 时间: 2023-10-25 10:01

liurui 发表于 2023-10-24 19:46
文章是 J. Phys. Chem. Lett. 2019, 10, 918−924

不是程序中的通常做法的话，我没法提供更多帮助。

作者
Author: liurui 时间: 2023-11-1 22:06

ahxb 发表于 2023-10-25 10:01
不是程序中的通常做法的话，我没法提供更多帮助。

收到，谢谢老师！！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)