| 查看: 6405 | 回复: 22 | |||||
[交流]
gauss 计算电偶极矩
|
|||||
| 优化了一个团簇，计算了电偶极矩，导师希望得到各个分子的电偶极矩，我就用matlab把各个分子的电荷与坐标相乘去计算各个分子的电偶极矩。现在又这样一个问题，我认为用matlab把所有的电荷与坐标相乘得到的应该是总的电偶极矩，但算完后发现与gauss输出文件中的dipole有误差，已经考虑了guass中dipole的单位是德拜。不知怎么回事。另外，我计算用的电荷是ATP charge。 |
回复此楼» 收录本帖的淘帖专辑推荐
 激发态计算 | QC Learning |
» 猜你喜欢
200年来人们一直无法破解的难题： 气液临界点状态方程
已经有7人回复
-
细胞培养，这22个细节一定要注意！（一）
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有299人回复
-
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
已知面上挂了，小学校实验方向无平台，感觉太难了
已经有26人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
求助各位大佬，球球了，在这科萨头
已经有0人回复
-
分子动力学模拟合作需求
已经有3人回复
-
招生材料与化工、环境科学调剂研究生
已经有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人回复
» 抢金币啦！回帖就可以得到:
-
河北大学招收生物与医药专业调剂考生
+3/107
-
安徽宿州学院085600材料与化工(专业学位)2026年研究生招生调剂信息
+1/88
-
欢迎调剂北京石油化工学院聚合物功能材料组
+1/84
-
浙师大化与材学院2026年有大量硕士研究生调剂名额 （线上 调剂）
+1/43
-
2026年河北大学无机化学专业招收调剂生
+1/40
-
西南民族大学广招硕士调剂生！快到碗里来！
+1/39
-
山东理工大学资源与环境工程学院招收调剂考生
+1/23
-
武汉工程大学化学与环境工程学院2026年硕士研究生调剂公告
+1/15
-
★上海★上海工程技术大学·环境与资源创新中心（ICER）接收化工材料环境等学科调剂生
+1/12
-
西北工业大学机电学院 复合材料加工制造方向招聘博士后
+1/8
-
[杭州电子科技大学]材料学院电池课题组-2026年硕士研究生调剂信息
+1/8
-
【博士招生】天津理工大学国家杰青王铁课题组招收2026年博士研究生
+1/8
-
齐齐哈尔大学李莉课题组诚招2026级考研调剂生（学硕和专硕）
+1/7
-
苏州大学医学部纳米生物医学方向招收调剂生1名
+1/6
-
天津大学光电子能源团队招收博士生、硕士生
+1/5
-
26
+1/5
-
国家大人才课题组招收化学、材料化工类2026年调剂生
+1/5
-
海水鱼类遗传育种课题组招收硕士调剂生
+1/5
-
河南城建学院资源与环境招调剂生
+1/3
-
广东省五邑大学MOFs催化功能材料课题组招收化学/材料背景调剂生
+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
