| 查看: 6408 | 回复: 22 | |||||
[交流]
gauss 计算电偶极矩
|
|||||
| 优化了一个团簇，计算了电偶极矩，导师希望得到各个分子的电偶极矩，我就用matlab把各个分子的电荷与坐标相乘去计算各个分子的电偶极矩。现在又这样一个问题，我认为用matlab把所有的电荷与坐标相乘得到的应该是总的电偶极矩，但算完后发现与gauss输出文件中的dipole有误差，已经考虑了guass中dipole的单位是德拜。不知怎么回事。另外，我计算用的电荷是ATP charge。 |
回复此楼» 收录本帖的淘帖专辑推荐
 激发态计算 | QC Learning |
» 猜你喜欢
085404（计算机技术）293分 希望调剂到材料相关专业
已经有0人回复
-
ELISA试验中不容忽视的细节盘点(二)
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有64人回复
-
细胞培养，这22个细节一定要注意！（一）
已经有0人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
七嗪类物质合成求助
已经有0人回复
-
分子动力学模拟合作需求
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 吸附计算可否人为添加化学键？
已经有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/183
-
淮北师范大学资源与环境硕士研究生调剂名额充足（26人），欢迎调剂。
+2/148
-
华北理工大学 矿业工程 研究生调剂
+3/127
-
聊城大学省重点实验室收化学调剂/网上面试无笔试/名额充足/可调专业广泛
+1/86
-
2026年赣南师范大学接收 电子信息工程+化学+材料与化工专业 调剂学生
+1/83
-
武汉纺织大学全国重点实验室招收材料、化学、纺织相关研究生
+1/64
-
青岛科技大学调剂招生
+1/55
-
中国石油大学（北京）理学院2026年硕士研究生接收调剂生：物理、化学、材料等
+1/43
-
昆明理工大学余富朝团队2026年诚招有机化学/生物化学等方向“申请-审核”博士研究生
+1/39
-
江西理工大学稀土学院招材料类硕士研究生，稀土光功能材料方向，非诚勿扰
+1/37
-
河南中部某本科高校 资源与环境（专硕）研究生 招生信息
+1/37
-
诚聘博士后 | 26年华东师大海优课题组——多相催化方向：CO₂资源化与废塑料转化
+1/31
-
可调声波驱动仪
+1/20
-
哈尔滨理工大学材料与化工学院刘刚老师团队招募“申请-考核”制博士
+1/8
-
中国农业大学资源与环境学院 - 招聘杰出人才（教授-专技四级）、青年研究员（副教授-�
+1/8
-
杭州师范大学-浙江省湿地智慧监测与生态修复重点实验室团队硕士调剂公告
+1/8
-
2026年西华大学材料学院-材料近净成形与表面技术研究团队-招收调剂研究生
+1/6
-
双一流天津工业大学电信学院李鸿强教授招收2026年申请审核制博士（4月3号截止）
+1/4
-
西南石油大学招收2026年“申请-考核“博士研究生 方向：不对称催化氢化
+1/4
-
河南城建学院生命科学与工程学院（085700资源与环境方向）有调剂名额
+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
