| 查看: 3471 | 回复: 6 | ||||
| 【奖励】 本帖被评价3次，作者beefly增加金币 2.4 个 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[资源]
用数值方法计算偶极矩
|
||||
|
此文目的： 1. 介绍偶极矩计算公式 2. 对于没有密度矩阵的理论方法，通过数值差分计算偶极矩 通常情况下，偶极矩都是用解析方法得到的。以CO分子为例，键长取1.2 Angstrom，用Gaussian程序计算CCSD/3-21G级别偶极矩的输入如下： # ccsd/3-21g density=cc Dipole moment with CCSD density 0 1 C O 1 1.2 计算得到的偶极矩为-0.0743 Debye，方向为-Z（即，从O指向C）。很多初学者忘了加density=cc，这种情况下得到的是Hartree-Fock级别的偶极矩。 如果用数值方法，偶极矩该怎么算呢？偶极矩定义为能量对电场强度的一阶导数：(μx, μy, μz) = (dE/dFx, dE/dFy, dE/dFz)，因此可以把偶极矩计算分解为不同电场下同一理论级别的一系列单点能计算。 从很多数值分析书上可以查到3点数值差分公式， dE/dFx = [-E(-Dx) + E(+Dx)]/(2Dx) 其中Dx是X坐标轴方向上某个自定义的电场大小；E(-Dx)、E(+Dx)分别为给定电场下计算的分子总能量。还可以采用更精确的5点数值差分公式， dE/dFx = [E(-2Dx) -8E(-Dx)+8E(+Dx)-E(+2Dx)]/(12Dx) 类似地，可以得到dE/dFy和dE/dFz的数值差分公式。 从上面的公式可以看到，3点公式只需要算两个点，5点公式只需要4个点。这是因为在一阶导数的数值差分计算中，用不到电场为0的点。 有了公式，下面可以设计计算流程了。仍是CO分子，C-O键仍放在Z轴。因为在X和Y方向，分子是对称的，偶极矩为零，故只需要算μz = dE/dFz。Gaussian提供了几种定义电场的方法，其中一种的命令为Field=z-n或Field=z+n，z表示电场方向，-n/+n表示电场为-/+n * 0.0001 a.u.。以下是第一个点的输入，电场取-0.002 a.u.： # ccsd/3-21g Field=z-20 CCSD energy calculation with an electric dipole field in the Z direction of -0.002 au 0 1 C O 1 1.2 第二、三、四个点的电场分别为-0.001，+0.001，+0.002 a.u.，需要分别把Field=z-20替换为Field=z-10，Field=z+10，和Field=z+20。四个点的CCSD能量如下（单位Hartree）： -112.2970354 -112.2970453 -112.2971037 -112.2971523 代入5点公式（Dz=0.001），换算成Debey（换算因子2.542）后，为-0.0742 Debey，非常接近上面通过解析方法得到的-0.0743 Debey。如果用三点公式（取第2、3个点），仍然得到-0.0742 Debey。 应用 上面计算偶极矩的数值方法有什么用呢？我们知道，Gaussian程序目前还不支持CCSD(T)，BCCD(T)，QCISD(TQ)，等理论级别的密度矩阵。如果想得到这些理论级别的偶极矩，要么换程序（例如免费的CFour程序），要么就用上面的数值方法。 还有些理论方法是没有密度矩阵的，例如CBS，BSSE，Gn理论等。用上面的数值方法，原则上也可以计算这些理论级别的偶极矩。 [ Last edited by beefly on 2011-6-26 at 11:56 ] 补充： 全电子相对论方法（如ZORA，DKH2，等）的哈密顿需要做图景变换，因此一般情况下不能用上面的数值方法算偶极矩，除非在电场计算程序中加入图景变换。 [ Last edited by beefly on 2012-8-2 at 21:56 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 @个人收集@ | 我的收藏 |
» 猜你喜欢
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生
已经有11人回复
-
分子动力学模拟合作需求
已经有4人回复
-
物理化学论文润色/翻译怎么收费?
已经有192人回复
-
中国科学院成都山地灾害与环境研究所2026年特别研究助理（博士后）招聘启事
已经有0人回复
-
湖北文理学院材料与化工专硕还有调剂名额，欢迎大家推荐考生！
已经有1人回复
-
湖南大学刘巧玲课题组2026年博士研究生招生信息
已经有28人回复
-
南方科技大学表面催化反应动力学研究课题组招收2026届博士生
已经有1人回复
-
北京邮电大学物理学院老师招收博士生
已经有3人回复
-
求助：三乙酰丙酮铁（Fe(acac)₃）标准红外(IR)与XRD谱图/数据
已经有1人回复
-
Comsol 6.0 桶内肥料产热产湿过程模拟求助
已经有1人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 数值计算方法（浙江大学内部资料）文字版
已经有236人回复
- 有什么方法可以改变分子的偶极矩呢？
已经有10人回复
- 数值取值范围计算
已经有4人回复
- 数值计算的目的是什么？
已经有5人回复
- 如何计算水的介电常数，公式中的偶极矩怎么计算
已经有20人回复
- 振动分辨的电子光谱的计算
已经有30人回复
- 求助：如何用Gaussian计算激发态偶极矩？
已经有6人回复
- 一本关于c语言数值计算的书，比较全面
已经有7人回复
- gauss 计算电偶极矩
已经有22人回复
- 数值积分方法求助
已经有6人回复
- 关于G09计算激发态偶极矩，基态偶极矩，跃迁偶极矩的讨论
已经有16人回复
- 求助数值计算结果和仿真结果有什么不一样的
已经有6人回复
- 关于频率的计算
已经有6人回复
- 数值计算怎么学
已经有4人回复
- 大规模优化数值计算，梯度的计算方式讨论
已经有24人回复
- 已知偶氮苯的偶极矩及介电常数，怎样计算极性
已经有3人回复
- 关于偶极矩的计算公式
已经有3人回复
- 【求助】计算方法与数值方法的区别
已经有9人回复
- 【求助】数值计算和解析计算有什么不同
已经有3人回复
- 【求助】极化率、偶极矩、振动频率关系
已经有5人回复
送红花一朵 
6楼2017-08-03 11:21:16
2楼2011-06-26 13:08:43
3楼2011-06-26 13:09:09
4楼2017-08-02 16:47:22
