计算化学公社
标题: 求助在高斯计算中用极化率算介电常数的问题 [打印本页]
作者Author: hanmingxi 时间: 2020-4-27 00:45
标题: 求助在高斯计算中用极化率算介电常数的问题
请问下面图片公式中求介电常数时用到的极化率，是在高斯频率计算时得到的 Exact polarizability: 369.565 33.035 194.449 44.511 -77.532 286.211
通过所列出的值是对应标准坐标的下三角型格式xx，xy，yy，xz，yz，zz
计算平均极化率的公式如下：（a)=1/3(axx+ayy+azz)
这么求得的吗？如果是还想问下为什么我求的值比图中大的多，是还需要单位换算吗？
作者Author: sobereva 时间: 2020-4-27 00:52
应当用Multiwfn处理这类数据，方便又省心，还能给你许多额外的重要信息，输出单位还可以直接指定。仔细看
使用Multiwfn分析Gaussian的极化率、超极化率的输出
http://sobereva.com/231
为了避免你对单位理解有误，这里给你我培训班里的一张幻灯片，单位都注明了
(, 下载次数 Times of downloads: 68)
作者Author: 王寓于 时间: 2020-7-28 10:11
请问老师，这里的分子极化率就是根据《使用Multiwfn分析Gaussian的极化率、超极化率的输出》算出来的Isotropic average polarizability 这一项吗？因为我看文献α有用下面公式计算的，但我查了一些文献，感觉他说的不正确，但是又不确定对于α是不是有新的定义。谢谢。
where α′, μ, ε0, kB, and T are the isotropic polarizability (at zero frequency), dipole moment, electrical permittivity of vacuum, Boltzmann constant, and temperature, respectively.
作者Author: zhficcas 时间: 2020-7-28 20:33
文献中的极化率 就是把用gaussian计算的输出文件用Multiwfn计算出的第二项 ，我记得带有单位的那个数值（A^3).就是。
作者Author: sobereva 时间: 2020-7-29 03:17
别信他的
直接把数据代入我2L贴的公式就得到了
作者Author: 王寓于 时间: 2020-7-29 09:34
谢谢老师。
作者Author: xxzj 时间: 2021-7-12 09:07
老师，我想计算一个分子的介电常数，可以通过Multiwfn获得分子极化率，想再请教一下老师，分子数密度应该如何获得？
作者Author: sobereva 时间: 2021-7-12 12:42
做分子动力学模拟
作者Author: Marc0 时间: 2021-11-11 09:36
请问老师，通过分子动力学模拟计算分子数密度，是否有相关帖子？用什么软件最合适？
作者Author: sobereva 时间: 2021-11-12 00:22
普通有机分子用GROMACS跑动力学是首选
弄个盒子，填充分子，跑NPT的动力学，最后盒子尺寸会自发调整合适。里面有多少分子是你已知的，除以已平衡区间的盒子平均体积就是分子的数密度
动力学模拟的具体操作不是几句话能讲清楚的，北京科音分子动力学与GROMACS培训班（http://www.keinsci.com/workshop/KGMX_content.html）里有详细讲解
作者Author: ayanamiyao 时间: 2021-11-22 16:51
sob老师，请教一下，您这页ppt中极化率α，应该用Isotropic average polarizability volume值来算吧。感觉一个阿佛加德罗常数下去，就全成2了，请您指导一下。
作者Author: sobereva 时间: 2021-11-23 06:46
看具体例子
(, 下载次数 Times of downloads: 37)
作者Author: ayanamiyao 时间: 2021-11-23 10:39
谢谢sob老师指点~OrZ
作者Author: 哎呀妈呀 时间: 2021-12-9 11:34
您好，看到您的帖子讨论了Clausius−Mossotti。很感兴趣，您能把您的文章和引用的Clausius−Mossotti的文章发给我吗？万分感谢
作者Author: zzffzz33 时间: 2022-3-4 20:53
10.1021/acs.jpclett.7b00688
这篇的supporting里有
作者Author: yangjunfang 时间: 2022-6-7 14:10
sob老师，您好，我想通过极化率计算晶体的介电常数，请问极化率通过Gaussian对单分子进行DFT计算就可以吗，然后再按照您说的计算公式计算得到介电常数？
作者Author: sobereva 时间: 2022-6-9 07:24
没法用于晶体
作者Author: yangjunfang 时间: 2022-6-9 16:00
老师，那两种不同极性的分子组成的无定形聚集体的介电常数有办法预测吗
作者Author: sobereva 时间: 2022-6-10 06:56
至少不是那个预测公式能预测的
作者Author: yangjunfang 时间: 2022-6-27 21:59
sob老师，您好，请问您说的凝聚相体系的介电常数通过分子动力学计算是为了得到体系的密度，从而得到数密度，然后再通过这个公式计算吗
作者Author: sobereva 时间: 2022-6-29 01:47
分子动力学直接就能得到介电常数
(, 下载次数 Times of downloads: 20)
作者Author: yangjunfang 时间: 2022-6-29 09:40
好的，谢谢sob老师
作者Author: yangjunfang 时间: 2022-6-29 19:43
老师，您好我想问一下这个公式计算的是同种分子组成的凝聚相才可以使用吗？不同分子组成的凝聚相是不是不可以用
作者Author: wzkchem5 时间: 2022-6-29 20:22
不是的，这个公式推导的时候没有用到所有分子都是同一种这个条件。但是不同分子组成的凝聚相通常需要更长时间的MD才能平衡
作者Author: yangjunfang 时间: 2022-6-30 09:59
好的老师，我明白了
作者Author: yangjunfang 时间: 2023-7-4 18:44
本帖最后由 yangjunfang 于 2023-7-4 18:46 编辑
老师，您好，我用分子动力学模拟计算了不同掺杂比例的凝聚相体系的介电常数，但是算出来的值基本一致，都是一点多。请问这个原因有可能是什么啊？
作者Author: sobereva 时间: 2023-7-5 02:30
问题要问完整，都不说什么和什么掺杂
诸如你用环己烷和水掺杂，1:99和99:1，算出来的介电常数不可能没有显著差异
作者Author: yangjunfang 时间: 2023-7-5 08:31
本帖最后由 yangjunfang 于 2023-7-5 08:34 编辑
老师 是1wt%和20wt%的客体-主体掺杂 形成凝聚相
作者Author: cheng1234 时间: 2024-5-8 16:13
本帖最后由 cheng1234 于 2024-5-8 16:16 编辑
Sob老师您好！我按照您这页PPT里面给出的例子计算了CS2的极化率，但是不知道是不是我关键字写的有问题，我用Multiwfn处理得到如下的结果，烦请老师指点一二！谢谢您！
1.优化输入：
(, 下载次数 Times of downloads: 5)
D:\Structure\Polar\OPT
2.计算极化率：
(, 下载次数 Times of downloads: 6)
D:\Structure\Polar\Polar
3.Multiwfn处理极化率输出结果：
(, 下载次数 Times of downloads: 4)
D:\Structure\Polar\Multiwfn处理结果
作者Author: sobereva 时间: 2024-5-8 20:05
你当前用的是MP2算的超极化率，显然应该选2，选项上的提示都写明了
作者Author: cheng1234 时间: 2024-5-9 09:22
老师不好意思，我没认真仔细看，感谢您的指点！
| 欢迎光临 计算化学公社 (http://bbs.keinsci.com/) |
Powered by Discuz! X3.3 |