计算化学公社

 找回密码 Forget password
 注册 Register

使用Multiwfn预测晶体密度、蒸发焓、沸点、溶解自由能等性质

查看数: 40604 | 评论数: 42 | 收藏 Add to favorites 31
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2016-6-21 04:05

正文摘要:

使用Multiwfn预测晶体密度、蒸发焓、沸点、溶解自由能等性质Using Multiwfn to predict crystal density, heat of vaporization, boiling point and solvation free energy 文/Sobereva @北京科音First release: ...

回复 Reply

hydrogen2333 发表于 Post on 2025-6-27 15:24:23
sobereva 发表于 2025-6-25 20:03
对晶胞直接进行计算需要用第一性原理程序(Gaussian也支持周期性计算,只不过对于三维周期性体系效率很低 ...

明白了,谢谢老师
sobereva 发表于 Post on 2025-6-25 20:03:13
hydrogen2333 发表于 2025-6-24 10:23
老师,请问您文章中提到的叠氮化铵是用量子化学软件计算出来的吗?我有点困惑,比如叠氮化铵,维基百科上写 ...

对晶胞直接进行计算需要用第一性原理程序(Gaussian也支持周期性计算,只不过对于三维周期性体系效率很低,很少有人用罢了)。而当前是基于晶体中单个分子/离子的属性预测晶体的性质,Gaussian算的对象是单个分子或离子,所以不需要第一性原理程序或周期性计算
hydrogen2333 发表于 Post on 2025-6-24 10:23:17
老师,请问您文章中提到的叠氮化铵是用量子化学软件计算出来的吗?我有点困惑,比如叠氮化铵,维基百科上写得它是晶体,印象中对于晶体比如氯化钠用的是第一性原理的软件进行的计算,为什么对于叠氮化铵用的就是量子化学的方法呢?甚至更极端一点比如二氧化碳在低温下也可以形成晶体但它也用的是量子化学的计算方法呀。
njtech_wk 发表于 Post on 2024-6-17 00:58:52
sobereva 发表于 2024-6-16 05:16
不适合,原文的训练集里不包含这种情况
没有简单且准确的预测方法,顶多是凑合用这个预测公式,但不要指 ...

谢谢sob老师
sobereva 发表于 Post on 2024-6-16 05:16:46
njtech_wk 发表于 2024-6-15 22:53
请问sob老师,带重原子的比如I原子的有机分子,推荐用这个预测中性体系的晶体密度吗,不推荐的话,还有什么 ...

不适合,原文的训练集里不包含这种情况
没有简单且准确的预测方法,顶多是凑合用这个预测公式,但不要指望能有多准。严格的做法是做晶体结构预测,但复杂度和耗时远高于本文的经验预测方法。
njtech_wk 发表于 Post on 2024-6-15 22:53:08
请问sob老师,带重原子的比如I原子的有机分子,推荐用这个预测中性体系的晶体密度吗,不推荐的话,还有什么方法可以预测这种含有重原子中性体系晶体的密度呢?
黑色桃花 发表于 Post on 2023-9-24 13:32:17
sobereva 发表于 2023-9-24 13:30
2 逻辑极其混乱
跟chk文件毫无直接联系
热力学校正量和波函数相关的分析结果没任何直接联系

好的,谢谢老师,我再看看帖子好好学学
sobereva 发表于 Post on 2023-9-24 13:30:58
黑色桃花 发表于 2023-9-24 13:25
嗯嗯老师,我确实没说清楚,以下是我的补充:
2. 校正方式这个我不知道,说这句话的原因是因为我看了sher ...

2 逻辑极其混乱
跟chk文件毫无直接联系
热力学校正量和波函数相关的分析结果没任何直接联系

3 详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/60
按原子序号设置基组

该说的都说了
谈谈弥散函数和“月份”基组
http://sobereva.com/119

4 该说的都说了
谈谈量子化学中基组的选择
http://sobereva.com/336http://bbs.keinsci.com/thread-3545-1-1.html
黑色桃花 发表于 Post on 2023-9-24 13:25:59
sobereva 发表于 2023-9-24 12:58
1 是

2 说清楚什么校正

嗯嗯老师,我确实没说清楚,以下是我的补充:
2. 校正方式这个我不知道,说这句话的原因是因为我看了shermo和热力学校正因子帖子后,我看到了得出了chk文件会有没考虑非谐振频率的问题,我就是单纯地想问一下,这个问题会对上述计算结果有影响没,比如离子体积,电荷分布那些
3. 我不知道怎么单独给原子加弥散函数,目前我只是会给整体体系加,那么首先老师,我该看哪个帖子去学,其次,整体都加弥散可以吗,还是说对0,1的离子体系可以不用考虑弥散了
4. 这个问题是我不是按照上述的做法验证,而是假如现在我自己做,那么我用B3LYP/6-31+G*,是不是更合理,换句话说,我只是想根据我学帖子的知识,判断出上述所用基组参数的错误(①需要加弥散,②一个极化函数就行,不用两个,这样才是正确的)
5. 1号问题的回答解释了,这个没问题了
麻烦老师再看看我的问题,谢谢卢老师~~~
sobereva 发表于 Post on 2023-9-24 12:58:22
黑色桃花 发表于 2023-9-24 12:45
好的老师,我还有如下问题:
1.由于要进行频率分析,所以是不是在整个计算过程中(我自己算的体系),B3 ...

1 是

2 说清楚什么校正

3 带显著负电荷原子加弥散函数,其它无所谓,加了就是白浪费时间

4 用文中的经验预测方法,原文用什么级别你就必须用什么级别,什么可再解释的

5 看不懂你的问题的逻辑
黑色桃花 发表于 Post on 2023-9-24 12:45:52
sobereva 发表于 2023-9-24 12:20
作者就是在那个级别拟合的参数,显然自己算的时候也必须得用那个级别

好的老师,我还有如下问题:
1.由于要进行频率分析,所以是不是在整个计算过程中(我自己算的体系),B3LYP/6-311+G*就完全够用,我的意思是,在整哥过程中没有涉及到电子能量的计算,因此不需要采用高级别方法基组
2. 在fchk文件分析这些数据,是否用考虑校正情况,因为没有非谐振的频率(我不知道这个过程是从fchk文件分析出来的)
3.考虑弥散函数,阴离子体系必须加,但是如果是价态为0,自旋多重度为1的离子化合物体系,是否仍需把弥散函数当成必须考虑的因素,因为在整体看来其价态为0.属于中性,但是细分其实也有阴离子体系,所以我不明白这种情况怎么办
4.事实上,上述1.2章节如果正规的算(按照帖子)在对其阴离子体系计算时,是否必须加上弥散函数,并且我看上面加了两个极化函数,但是按照咱们的帖子,对于如上体系,一个极化函数便以足可(如上体系H原子不用考虑极化)
5. 在第三点默认0,1体系,不用特别考虑弥散函数的情况下,我在M062X/def2TZVPP下对一个离子化合物进行了opt/freq(仅仅为了得到好一点的能量),那么此时如果我想对阴阳离子进行如上的计算且我不想用额外的ma-TZVPP极化函数文件,此时我想在B3LYP/6-311+G* D3对阴阳离子体系(阳离子不用弥散,但是这样也行吧)计算,这样是否可行,整体逻辑是 好的方法基组对整体体系算了一遍opt/freq,然后对于单独阴阳离子体系用低级别的再算,这样的逻辑行吗
以上便是我想的问题,请老师批评指正,谢谢卢老师!
sobereva 发表于 Post on 2023-9-24 12:20:15
黑色桃花 发表于 2023-9-24 11:39
卢老师。在1.2章节中,对阴离子的计算中采用的是6-31G**,但是由于算了阴离子,为什么在这个时候不考虑弥散 ...

作者就是在那个级别拟合的参数,显然自己算的时候也必须得用那个级别
黑色桃花 发表于 Post on 2023-9-24 11:39:57
卢老师。在1.2章节中,对阴离子的计算中采用的是6-31G**,但是由于算了阴离子,为什么在这个时候不考虑弥散函数了呢,是因为这样算的数据比较吻合试验数据吗,请老师赐教
sobereva 发表于 Post on 2023-5-14 13:32:45
zjb-lufia 发表于 2023-5-14 10:47
我想找的粘度 这篇搜不到,需要DOI号

自己把期刊名页码那些输到Google学术里搜,或者到期刊官网上输入卷号、页码之类搜文章页面,基本学术检索常识

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

GMT+8, 2025-8-12 21:20 , Processed in 0.203431 second(s), 26 queries , Gzip On.

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