求助高斯怎么计算拉曼极化率导数和红外跃迁振动

飒飒萨达

老师们好，本人最近正在自己摸索高斯计算的问题，但找了相关的文章也没搞明白，特来求助，提前感谢！
1. 我要计算的是Raman polarizability Derivatives, 但始终没搞懂这个东西到底是什么意思，是一个vector还是一个tensor，还有请问这个东西高斯可以算吗，具体怎么计算呢？
2. 另外一个问题是想计算红外的跃迁振动（vibrational transition dipole moment)，想问这个东西怎么用高斯怎么计算呢？ 拜读了老师的教程，仔细研究electronic transition dipole moment
但不是很明白”振动“跃迁偶极矩，这该怎么操作？

我自己是研究SFG（sum frequency generation) 所以想通过计算红外的振动和拉曼的极化率乘积来对比实验上的二级实验数据，但量化新手实在搞不懂怎么计算，希望可以得到老师的说明和帮助

sobereva

建议你仔细把这个文档看了，写得很清楚
(Amos,书,1987)-MOLECULAR PROPERTY DERIVATIVES.pdf (2.73 MB, 下载次数 Times of downloads: 29)

2023-10-13 12:36 上传 Uploaded

点击下载
Click to download

模拟拉曼谱需要拉曼强度，拉曼强度通过拉曼活性转化而来（参考http://sobereva.com/578介绍的我这篇文章Chem Asian J. 2021, 16, 56–63的computational details部分），计算拉曼活性需要算极化率对正则坐标的一阶导数。拉曼活性可以直接由Gaussian做振动分析产生，freq=raman。

Gaussian做非谐振计算（freq=anharm）会给出振动态之间的跃迁偶极矩，比如：

     ==================================================
               Anharmonic Transition Moments
     ==================================================

Electric dipole : Fundamental Bands
------------------------------------------------------------------------
   Mode(n)                            X               Y               Z
      1(1)                       0.929015D-16   -0.633311D-14    0.279965D-02
      2(1)                      -0.364575D-16    0.771845D-15   -0.524004D-01
      3(1)                      -0.855190D-16   -0.156384D-01    0.918100D-15

Electric dipole : Overtones
------------------------------------------------------------------------
   Mode(n)                            X               Y               Z
      1(2)                      -0.174609D-15    0.870374D-14    0.318513D-02
      2(2)                       0.658918D-16    0.151073D-15    0.344299D-02
      3(2)                       0.567367D-16    0.766023D-14    0.104821D-02

Electric dipole : Combination Bands
------------------------------------------------------------------------
   Mode(n)     Mode(n)                X               Y               Z
      2(1)        1(1)           0.116730D-15   -0.196251D-14   -0.291637D-02
      3(1)        1(1)           0.137761D-15   -0.640374D-02   -0.121039D-14
      3(1)        2(1)          -0.135943D-15    0.619726D-02   -0.584640D-14

飒飒萨达

sobereva 发表于 2023-10-13 12:37
建议你仔细把这个文档看了，写得很清楚

谢谢老师的指导，太感激了，我好好看看