第11届北京科音分子动力学与GROMACS培训班将于8月13~16日于北京举办,是系统性学习分子动力学与GROMACS使用的最佳机会,请点击此链接查看培训详情,欢迎参加!

计算化学公社

 找回密码 Forget password
 注册 Register
Views: 1016|回复 Reply: 13

[新手求助] 关于有机分子介电常数的计算问题

[复制链接 Copy URL]

33

帖子

0

威望

129

eV
积分
162

Level 3 能力者

发表于 Post on 2022-5-25 11:16:25 | 显示全部楼层 Show all |阅读模式 Reading model
请教各位老师,我结合Clausius-Mossotti方程在进行几个分子的介电常数计算时,下面这几个结构类似分子的极化率和分子体积总是出现较大误差,想问下为什么会这样?下面我附上了gjf文件和其中一个分子结构。
202205251113477310..png
202205251115402603..png

5508

帖子

0

威望

3231

eV
积分
8739

Level 6 (一方通行)

发表于 Post on 2022-5-25 16:32:56 | 显示全部楼层 Show all
误差是和什么比的?和实验的极化率/分子体积比的?还是算出介电常数以后和实验的介电常数比的?

3万

帖子

99

威望

4万

eV
积分
84403

管理员

公社社长+计算化学玩家

发表于 Post on 2022-5-26 03:47:19 | 显示全部楼层 Show all
说清楚体积怎么算的
而且分子体积也根本不是可实验观测的,至少没有唯一实验上的定义

当前基组不带弥散函数,算极化率完全离谱,仔细看
谈谈量子化学中基组的选择
http://sobereva.com/336http://bbs.keinsci.com/thread-3545-1-1.html

注意看此文
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚
http://sobereva.com/620http://bbs.keinsci.com/thread-25787-1-1.html
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班。这些培训是计算化学快速入门以及全面系统性提升研究水平的最佳途径,培训各种相关信息见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人,讨论范畴相同
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

33

帖子

0

威望

129

eV
积分
162

Level 3 能力者

 楼主 Author| 发表于 Post on 2022-5-27 10:47:01 | 显示全部楼层 Show all
sobereva 发表于 2022-5-26 03:47
说清楚体积怎么算的
而且分子体积也根本不是可实验观测的,至少没有唯一实验上的定义

老师您好,分子体积的部分我是参照您的博文《谈谈分子体积的计算》中MT算法由multiwfn直接查看的。极化率我也是利用multiwfn直接查看的。
之前我用的基组是6-311g*,我算的有机分子体系大概在20-80个原子范围内不等,在gaussian16里采用b3lyp进行优化和频率分析,要加弥散函数的话,是不是改为6-31+g*会好些?但是感觉会蛮耗时的,如果只是想通过这个部分计算介电常数,之后利用此步的数值在隐式溶剂模型下计算电离能,是不是只优化结构就好了,频率分析感觉不是必要的,因为后面我还会用其他方法优化结构和计算能量。请老师指教。

33

帖子

0

威望

129

eV
积分
162

Level 3 能力者

 楼主 Author| 发表于 Post on 2022-5-27 10:49:53 | 显示全部楼层 Show all
wzkchem5 发表于 2022-5-25 16:32
误差是和什么比的?和实验的极化率/分子体积比的?还是算出介电常数以后和实验的介电常数比的?

我是和之前师兄师姐的计算结果对比了一下,发现这种类似体系下我的和他们的计算结果误差很大。

33

帖子

0

威望

129

eV
积分
162

Level 3 能力者

 楼主 Author| 发表于 Post on 2022-5-27 11:02:13 | 显示全部楼层 Show all
其实,对于这部分的计算我还有一个小小的问题,我刚才看了老师《谈谈弥散函数和“月份”基组》的博文,在计算极化率时需要加弥散函数,但是为什么必须要加啊?我有些不太能理解。如果说对于阴离子体系,对外层电子的束缚能力较弱,这个时候必须加弥散函数,而阳离子体系,并不会出现这种问题,所以不加,这个我是可以理解的。

5508

帖子

0

威望

3231

eV
积分
8739

Level 6 (一方通行)

发表于 Post on 2022-5-27 15:21:47 | 显示全部楼层 Show all
王若曦 发表于 2022-5-27 03:49
我是和之前师兄师姐的计算结果对比了一下,发现这种类似体系下我的和他们的计算结果误差很大。

你是体系和他们不一样,还是理论级别和他们不一样?具体是哪里不一样?
你的计算和前人计算不一样,那么显然要既告诉我们你的计算细节,也告诉我们前人的计算细节,我们才能判断原因。这是问问题的常识

5508

帖子

0

威望

3231

eV
积分
8739

Level 6 (一方通行)

发表于 Post on 2022-5-27 15:27:22 | 显示全部楼层 Show all
本帖最后由 wzkchem5 于 2022-5-27 08:29 编辑
王若曦 发表于 2022-5-27 04:02
其实,对于这部分的计算我还有一个小小的问题,我刚才看了老师《谈谈弥散函数和“月份”基组》的博文,在计 ...

极化率顾名思义算的是分子的电子云被电场极化的能力,极化的过程中电子会往电势更正的方向整体偏移,所以原本不是那么靠外的电子就会变得比平时更加靠外,需要用弥散函数描述。如果不加弥散函数,那么因为离原子核较远的地方缺乏基函数,电子没有办法极化得那么厉害,所以会导致低估极化率。
当然如果较真的话,对于带特别多正电的离子,即使算极化率,电子也不会跑多远,所以也不需要加弥散函数,比如假如你研究Al3+在气相里的极化率,那就不需要加弥散函数。但是绝大多数实际研究会涉及到的体系,都是需要加弥散函数才能把极化率算准的。

33

帖子

0

威望

129

eV
积分
162

Level 3 能力者

 楼主 Author| 发表于 Post on 2022-5-27 19:41:14 | 显示全部楼层 Show all
wzkchem5 发表于 2022-5-27 15:21
你是体系和他们不一样,还是理论级别和他们不一样?具体是哪里不一样?
你的计算和前人计算不一样,那么 ...

我是完全按照他们的实验步骤计算的,方法基组、计算内容完全相同,均是在b3lyp/6-311g*下对分子进行几何优化和频率分析,计算极化率,然后应用multiwfn查看log文件和fchk文件的。(当然这是我们都没加弥散函数之前的做法),然后对照了我和他们的极化率和介电常数计算数值,发现有很大差异。并且很巧合的是,出现这个问题的分子体系基本都是我在第一个帖子里附上的这一类型的分子,所以我感到很困惑。想问下是不是因为未考虑什么作用或者矫正项导致的计算错误呢?
202205251115402603..png

5508

帖子

0

威望

3231

eV
积分
8739

Level 6 (一方通行)

发表于 Post on 2022-5-27 20:21:42 | 显示全部楼层 Show all
王若曦 发表于 2022-5-27 12:41
我是完全按照他们的实验步骤计算的,方法基组、计算内容完全相同,均是在b3lyp/6-311g*下对分子进行几何 ...

是不是所选的构象和他们不一样?比如环己烷选的初始构象不合理,卡在了一个高能量的构象上

11

帖子

0

威望

172

eV
积分
183

Level 3 能力者

发表于 Post on 2022-5-27 20:30:04 | 显示全部楼层 Show all
王若曦 发表于 2022-5-27 19:41
我是完全按照他们的实验步骤计算的,方法基组、计算内容完全相同,均是在b3lyp/6-311g*下对分子进行几何 ...

什么材料还带环己烷的?

33

帖子

0

威望

129

eV
积分
162

Level 3 能力者

 楼主 Author| 发表于 Post on 2022-5-27 22:10:46 | 显示全部楼层 Show all
风起~ 发表于 2022-5-27 20:30
什么材料还带环己烷的?

就是选取的一些有机共轭分子算电离能

33

帖子

0

威望

129

eV
积分
162

Level 3 能力者

 楼主 Author| 发表于 Post on 2022-5-27 22:18:27 | 显示全部楼层 Show all
wzkchem5 发表于 2022-5-27 20:21
是不是所选的构象和他们不一样?比如环己烷选的初始构象不合理,卡在了一个高能量的构象上

那我再去看一下初始结构,再计算一下吧,谢谢老师

3万

帖子

99

威望

4万

eV
积分
84403

管理员

公社社长+计算化学玩家

发表于 Post on 2022-5-28 01:57:43 | 显示全部楼层 Show all
王若曦 发表于 2022-5-27 10:47
老师您好,分子体积的部分我是参照您的博文《谈谈分子体积的计算》中MT算法由multiwfn直接查看的。极化率 ...

6-31+g*如果你是指计算极化率,明显是很烂的选择

谈谈量子化学中基组的选择
http://sobereva.com/336http://bbs.keinsci.com/thread-3545-1-1.html
明确说了该用什么算极化率
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班。这些培训是计算化学快速入门以及全面系统性提升研究水平的最佳途径,培训各种相关信息见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人,讨论范畴相同
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )

GMT+8, 2022-8-12 02:26 , Processed in 0.769764 second(s), 27 queries .

快速回复 返回顶部 返回列表 Return to list