| 查看: 6459 | 回复: 22 | |||||
[交流]
gauss 计算电偶极矩
|
|||||
| 优化了一个团簇，计算了电偶极矩，导师希望得到各个分子的电偶极矩，我就用matlab把各个分子的电荷与坐标相乘去计算各个分子的电偶极矩。现在又这样一个问题，我认为用matlab把所有的电荷与坐标相乘得到的应该是总的电偶极矩，但算完后发现与gauss输出文件中的dipole有误差，已经考虑了guass中dipole的单位是德拜。不知怎么回事。另外，我计算用的电荷是ATP charge。 |
回复此楼» 收录本帖的淘帖专辑推荐
 激发态计算 | QC Learning |
» 猜你喜欢
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有257人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有1人回复
-
实际气体内压强的微观机制或将改写--------从一个基本的力学佯谬说起
已经有1人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有5人回复
-
招收0703开头的学硕调剂生 线上复试
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 吸附计算可否人为添加化学键？
已经有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人回复
» 抢金币啦！回帖就可以得到:
-
2026年江西师范大学人文地理与低碳经济方向拟招收1名“申请入学制”学术博士研究生
+2/132
-
太行山徒步之凌水河（二）
+1/92
-
北京-89175-事业单位-诚征女友
+1/62
-
双一流高校-南京林业大学-化学工程学院-国家海外优青团队招2026级博士（5月15号截止）
+1/47
-
Google ai pro 会员分享
+1/29
-
海外优青课题组26年招收学术博士2名
+1/28
-
诚征女友，坐标西安
+1/28
-
海南大学药学院黄玲教授2026招收药理（神经药理）方向博士生一名
+1/27
-
智合健物课题组2026年博士生招生(5月6日—2026年5月13日)
+1/25
-
南京邮电大学李巍教授招收2026博士生（5月12日前有效）
+1/21
-
北京大学新材料学院招聘柔性电子方向博士后
+2/12
-
南京林业大学-国家级青年人才团队 招2026级博士 （合成化学）
+1/11
-
深圳理工大学-湖南大学项目招收钙钛矿方向博士生
+1/10
-
招收2026年秋季入学新能源材料方向博士生
+1/7
-
2026年北京化工大学化学工程学院程元徽教授团队博士研究生招生-考核
+1/6
-
招收2026年秋季入学博士生1名（河北工业大学/北京科技大学联合 增材制造/生物材料）
+1/6
-
港科大广州 - 量子计算全额奖学金博士招聘（申请制，2026和2027均有名额）
+1/6
-
科研服务：病毒包装、细胞系构建、蛋白表达、基因编辑
+1/4
-
上海交通大学复合材料研究所诚聘科研助理
+1/4
-
美国博士后招聘，代谢的转录调控方向，干+湿结合实验室，欢迎Email咨询。
+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