| 查看: 6222 | 回复: 22 | |||||
[交流]
gauss 计算电偶极矩
|
|||||
| 优化了一个团簇，计算了电偶极矩，导师希望得到各个分子的电偶极矩，我就用matlab把各个分子的电荷与坐标相乘去计算各个分子的电偶极矩。现在又这样一个问题，我认为用matlab把所有的电荷与坐标相乘得到的应该是总的电偶极矩，但算完后发现与gauss输出文件中的dipole有误差，已经考虑了guass中dipole的单位是德拜。不知怎么回事。另外，我计算用的电荷是ATP charge。 |
回复此楼» 收录本帖的淘帖专辑推荐
 激发态计算 | QC Learning |
» 猜你喜欢
国内青年破解费曼难题，从分子层面阐释液体蒸发（气液相变）微观过程
已经有5人回复
-
研究发现一个新的临界系数，独立于临界压缩因子
已经有7人回复
-
物理化学论文润色/翻译怎么收费?
已经有288人回复
-
200年来人们一直无法破解的难题： 气液临界点状态方程
已经有7人回复
-
细胞培养，这22个细节一定要注意！（一）
已经有0人回复
-
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
分子动力学模拟合作需求
已经有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人回复
» 抢金币啦！回帖就可以得到:
-
【2026/2027 哈工大计算机类博士招生】
+1/78
-
[求助] 半导体物理PPT课件 东华理工机电学院 彭新村老师
+1/76
-
【2026/2027 哈工大计算机类博士/硕士招生】
+1/75
-
深圳大学张雷教授课题组诚聘博士后
+1/74
-
广东工业大学自动化学院鲁仁全教授团队刘勇华老师招收2026年博士研究生（申请制）
+1/74
-
关于本子打包
+1/70
-
西交利物浦大学黄彪院士招收26年全奖博士生1名（工业智能方向）
+1/29
-
中国科学院上海光学精密机械研究所 特种强激光薄膜课题组
+1/29
-
2026湖南科技大学化学化工学院招新能源电池方向博士生
+1/28
-
【青岛大学】2026年生物与医药博士生招生（含少数民族骨干人才）
+1/28
-
法国斯特拉斯堡大学有机光伏全奖博士招聘
+1/20
-
美国密苏里大学堪萨斯城分校（UMKC）生物材料诚聘全奖博士
+1/12
-
博士后, 博士招生 美国大纽约地区 NJIT
+1/11
-
顺磁计算交流
+1/6
-
澳科大诚招2026年秋季生物材料全奖博士研究生（今日16:30线上宣讲会）
+1/6
-
华南理工大学宋波教授招收2026年博士生（二氧化碳转化方向优先）
+1/4
-
中国科学技术大学 招特任副研究员 有机-钙钛矿方向
+1/4
-
顾敏院士课题组招收2026级光学工程专业博士研究生-上海理工大学智能科技学院
+1/3
-
换工作
+1/2
-
中南民族大学-国家级领军人才团队超支化聚合物方向2026年博士研究生招生
+1/2
★
小木虫: 金币+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