| 查看: 7103 | 回复: 22 | |||||
[交流]
gauss 计算电偶极矩
|
|||||
| 优化了一个团簇，计算了电偶极矩，导师希望得到各个分子的电偶极矩，我就用matlab把各个分子的电荷与坐标相乘去计算各个分子的电偶极矩。现在又这样一个问题，我认为用matlab把所有的电荷与坐标相乘得到的应该是总的电偶极矩，但算完后发现与gauss输出文件中的dipole有误差，已经考虑了guass中dipole的单位是德拜。不知怎么回事。另外，我计算用的电荷是ATP charge。 |
回复此楼» 收录本帖的淘帖专辑推荐
 激发态计算 | QC Learning |
» 猜你喜欢
博士生招生, 材料模拟方向，University of South Carolina
已经有0人回复
-
Postdoc position, materials simulations/machine learning, U. of South Carolina
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有126人回复
-
祈福青基
已经有34人回复
-
牛来！米来！面来！
已经有8人回复
-
2027申博咨询
已经有11人回复
-
寻求CrC2Al的cif文件，感谢
已经有0人回复
-
求问：LATP煅烧应该使用什么容器？
已经有0人回复
-
ORCA 纯新手向入门指南0.7版
已经有0人回复
-
哈尔滨工业大学（深圳）绿色合成化学团队招收27年春季或秋季入学的博士生
已经有5人回复
-
吉林大学化学学院理论化学研究所刘靖尧教授课题组招收2027级申请考核制博士研究生
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 吸附计算可否人为添加化学键？
已经有10人回复
- 关于电偶极矩的一个简单问题
已经有5人回复
- bader电荷布局分析怎么运行啊？
已经有14人回复
- 关于频率的计算
已经有6人回复
- 请教压电原理的计算（有点类似电容）。
已经有7人回复
- 提交siesta任务计算失败
已经有11人回复
- 偶极矩问题，请大家帮帮忙，都糊涂了
已经有8人回复
- 滚动阻力矩如何计算？
已经有4人回复
- 请问Gauss中溶剂化的计算怎样指定溶剂
已经有7人回复
- 偶极矩分析
已经有20人回复
- 单掺或者共掺ZnO的电偶极矩怎么求啊
已经有3人回复
- MOLPRO优化计算为什么如此神奇？急啊！
已经有8人回复
- siesta计算dipole moment（电偶极矩）
已经有7人回复
- 激发态的跃迁偶极矩
已经有9人回复
- 关于SIESTA中如何计算功函数的一种做法探讨
已经有8人回复
- 【求助】vasp Berry phase计算自发极化
已经有4人回复
- 【求助】电荷转移积分的计算
已经有10人回复
» 抢金币啦！回帖就可以得到:
-
项目结题验收，希望一切顺利
+1/971
-
【推免】】中科院国家级人才团队招收27级推免硕士生
+1/72
-
福建农林大学王志刚教授团队诚聘青年教师、博士后、博士生 - 材料、化学、力学、高分�
+1/42
-
天津大学 刘浩业教授团队 招收27级博士生
+1/40
-
中山大学中法核工程与技术学院放射生物学课题组招收推免生
+1/38
-
山东征女友，坐标济南
+1/36
-
中科院大连化学物理研究所刘生忠研究员团队招收联合培养博士
+1/25
-
中国科大-合肥国家实验室冷原子量子网络团队招聘启事
+2/16
-
深圳大学2027级光电信息工程/物理学研究生推免
+1/14
-
生物功能应用的仿生超湿材料专题征稿（Scopus、DOAJ 收录）
+1/13
-
中国科学院院士成会明招收2027届高分子方向推免直博生（截止2026年9月19日）
+1/5
-
【生物材料与生物打印专题征稿】Scopus期刊，免APC
+1/5
-
中国科学院院士成会明招收2027届推免直博生（截止2026年9月19日）
+1/3
-
信息提供：中国科学院院士成会明招收2027届推免直博生（截止2026年9月19日）
+1/3
-
北理工集成电路杰青团队 | 诚招科助理
+1/3
-
口腔特刊客座招募客编Biofunctional Materials牙科生物功能材料专题，Scopus期刊
+1/3
-
水利工程博士点师资博士后招聘
+1/3
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/2
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/1
-
中国科学院院士成会明招收2027届推免直博生（截止2026年9月19日）
+1/1
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
“系统找不到指定路径”并不是出错，不用管它。 如果直接用Multiwfn自带的原子波函数文件就把atomwfn拷到上一级目录，Multiwfn就会自动用它而不会再去调用高斯去计算原子波函数文件。 如果没这样做，或者你的体系有超过第四周期的元素，Multiwfn就会先让你输入计算所用的理论方法和基组，然后输入高斯的路径（如果已经在settings.ini里的gaupath中设定了路径则不会再让你输入）。 你的输入不对。注意看程序提示。第一步让你输入理论方法和基组，例子已经给明了，如B3LYP/cc-pVDZ。第二步输入高斯可执行文件的路径，提示中也已经明确地给出了输入例子，如d:\study\g03w\g03.exe。 在Multiwfn中凡是需要用户输入的地方都会给出明确的提示，e.g.后是输入例子。一定要注意阅读和理解屏幕上的输出。如果仍不明白就看手册的相应章节的说明。 [ Last edited by sobereva on 2013-9-23 at 02:56 ] |
17楼2013-09-23 02:54:58
2楼2013-09-13 15:45:18
3楼2013-09-16 19:19:40
★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
gmy1990: 金币+3 2013-09-17 03:02:10
春华秋实2020: 金币+1 2013-09-17 21:38:49
春华秋实2020: 金币+2 2013-09-17 21:39:32
小木虫: 金币+0.5, 给个红包，谢谢回帖
gmy1990: 金币+3 2013-09-17 03:02:10
春华秋实2020: 金币+1 2013-09-17 21:38:49
春华秋实2020: 金币+2 2013-09-17 21:39:32
| 所谓的ATP charge有两种，一个是GAPT，Gaussian在做频率计算的时候会打印此电荷，但是GATP charge的计算方法没有理论依据，结果也不太好。另一个是意大利人提出的“IR charge”，不过这个名字误导人，其实并未用到红外强度信息。IR charge理论上比较严谨，能够重复电偶极矩，但是仅适用于平面分子。 |
4楼2013-09-16 20:18:08
5楼2013-09-17 02:57:07
★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
春华秋实2020: 金币+5, 太感谢了，受益匪浅 2013-09-17 21:42:37
小木虫: 金币+0.5, 给个红包，谢谢回帖
春华秋实2020: 金币+5, 太感谢了，受益匪浅 2013-09-17 21:42:37
|
用Mulliken电荷算偶极矩结果很差，就连整体分子的偶极矩都无法很好重现，对于团簇内某个分子的偶极矩也不可能很好表现。GAPT电荷也不推荐使用。ADCH电荷可以完全精确重现体系总偶极矩，也可以期望用ADCH电荷计算的分子偶极矩也是很合理的。ADCH电荷原文见J. Theor. Comput. Chem., 11, 163，各种原子电荷包括ADCH的全面对比见《原子电荷计算方法的对比》（物理化学学报, 28, 1)。Multiwfn支持ADCH电荷的计算。 LS的“电荷与坐标相乘是得不到总偶极矩的”这句话是不严格的。像ADCH这样能精确重现体系总偶极矩的方法这句话是错误的。通过电荷与坐标相乘再加和就能得到。 值得一提的是，团簇当中分子间会有电荷转移，因此其中的分子可以认为是有一定净电荷的。对于带有净电荷的体系，其偶极矩是受到原点选取的影响的，也就是说偶极矩并非是唯一确定的，要计算的话必须指明原点，比如取此分子质心作为原点。 |
6楼2013-09-17 08:32:48
7楼2013-09-17 21:45:13
8楼2013-09-17 23:06:45
9楼2013-09-18 03:43:11
10楼2013-09-18 06:27:57
11楼2013-09-21 08:48:16
12楼2013-09-21 08:57:08
★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
春华秋实2020: 金币+3, 多谢 2013-09-22 08:45:50
小木虫: 金币+0.5, 给个红包，谢谢回帖
春华秋实2020: 金币+3, 多谢 2013-09-22 08:45:50
|
手册(pdf文件)在Multiwfn主页http://multiwfn.codeplex.com上下载 Multiwfn是功能最全面的波函数分析程序，作图只是其用处的很小一部分。 .wfn/.wfx/.fch/.molden都可以用作计算ADCH电荷的输入文件。例子见手册4.7.2节。计算ADCH电荷需要原子波函数文件，最简单的近似处理办法就是把examples目录下的atomwfn目录挪到上一级目录（对于Windows版而言），Multiwfn就会直接用这些自带的原子波函数文件。 |
13楼2013-09-21 09:22:22
14楼2013-09-22 10:17:34
15楼2013-09-22 14:53:34
16楼2013-09-22 15:18:16
18楼2013-09-23 14:02:01
19楼2014-07-23 17:38:22
20楼2014-07-27 09:12:28
21楼2014-07-27 16:41:45
22楼2015-11-16 09:20:12
23楼2019-11-20 00:00:39
