计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: sobereva
打印 Print 上一主题 Last thread 下一主题 Next thread

[综合交流] 简单量子化学问题答疑专帖(旧版,已关帖)

   关闭 [复制链接 Copy URL]

91

帖子

0

威望

419

eV
积分
510

Level 4 (黑子)

4891#
发表于 Post on 2017-10-25 16:20:06 | 只看该作者 Only view this author
sob老师  HF用单电子积分描述波函数,其未考虑电子相关   请问这里的电子相关具体指些什么

259

帖子

0

威望

851

eV
积分
1110

Level 4 (黑子)

4892#
发表于 Post on 2017-10-25 21:31:03 | 只看该作者 Only view this author
老师,VMD做轨道图的时候一般是不是,先用multiwfn生成格点文件,然后就可以生成轨道图像,如http://sobereva.com/134使用Multiwfn绘制NBO及相关轨道所述(为什么要说NBO?)。VMD在量化计算方面还可以用来绘制哪些图像。官网的图例好多是生物高分子方面的,VMD的手册感觉实例很少,不知道有没有更好的手册,还是不需要学那么多啊。另外,量化计算还有哪些常用的绘图软件啊。谢谢老师!

100

帖子

0

威望

414

eV
积分
514

Level 4 (黑子)

4893#
发表于 Post on 2017-10-25 21:50:32 | 只看该作者 Only view this author
sob老师好,在解决SCF不收敛(L502错误)时,我在关键词里加了scf(novaracc,noincfock) int=acc2e=12, 单点能数据能直接拿来用吗?是否还需要读取波函数并去掉上述的关键词重新算一遍?

357

帖子

0

威望

2071

eV
积分
2428

Level 5 (御坂)

真 掘墓者

4894#
发表于 Post on 2017-10-26 05:48:16 | 只看该作者 Only view this author
SMD下边做的结构优化,感觉上能量是有下降的,开始的几步跑30几个cycle,然后就某一步超128步数停了,继续算,情况差不多,能量略微下降。
这种情况可以做回大菜鸟吗,改maxcyc?

微信截图_20171026084014.png (25.65 KB, 下载次数 Times of downloads: 26)

微信截图_20171026084014.png
圣诞刨坟忙

559

帖子

0

威望

6220

eV
积分
6779

Level 6 (一方通行)

4895#
发表于 Post on 2017-10-26 09:12:35 | 只看该作者 Only view this author
请问老师:
1.对于复合物计算的偶极矩值是否能比较复合物稳定性的大小?
2.文献中有很多计算A与B形成复合物的偶极矩值,解释相互作用,这样做有什么根据?
3.有些文献中计算偶极矩值采用的基组没加弥散,这样的结果可信吗?

113

帖子

0

威望

271

eV
积分
384

Level 3 能力者

4896#
发表于 Post on 2017-10-26 09:37:57 | 只看该作者 Only view this author
sobereva 发表于 2017-10-25 12:21
NBO较新版本里有个CMO关键词,可以给出分子轨道中各个NBO的贡献。NBO是可以区分sigma,sigma*,pi,n*等 ...

老师,分子轨道中NBO的贡献可以在multiwfn里面计算吗?还是必须要用新的NBO软件。

159

帖子

0

威望

1231

eV
积分
1390

Level 4 (黑子)

4897#
发表于 Post on 2017-10-26 10:25:50 | 只看该作者 Only view this author

如何用高斯计算乙醇每一个单键的键能?:

本帖最后由 XZDZ 于 2017-10-26 10:34 编辑

请问一下,已经对C2H5OH做了opt+ferq,要怎要我才可以得知C2H5OH当中具体的每一个单键的键能呢?是还需要做什么样的计算吗?谢谢
在追加一个问题
     因为做了opt+freq之后,得出了一个最稳定的结构,这个结构就是乙醇分子在空间转动的时候能量最低的一个位置结构,那么,在这一个结构的时候,每一个键角、键长我们可以得知并且标定出来,但是每一个键的能量为啥不能直接得知呢?  比如在计算乙醇与HO2的反应的时候,乙醇可能会有α,β等三个不同位置的脱氢,所以我们就需要知道这三处不同位置的键的能量啊,但是再对过渡态的位置进行opt+freq以及能量计算的时候,这个键能能否得知呢?谢谢

588

帖子

0

威望

3310

eV
积分
3898

Level 5 (御坂)

4898#
发表于 Post on 2017-10-26 10:36:52 | 只看该作者 Only view this author
个人觉得不是太好计算,通过柔性力常数考察键的强度http://sobereva.com/364

159

帖子

0

威望

1231

eV
积分
1390

Level 4 (黑子)

4899#
发表于 Post on 2017-10-26 11:07:06 | 只看该作者 Only view this author
xujian 发表于 2017-10-26 10:36
个人觉得不是太好计算,通过柔性力常数考察键的强度http://sobereva.com/364

这个好像只能比较不同的键,它们的键强度差异呀?不能够得到具体的键能的数值么?

6万

帖子

99

威望

6万

eV
积分
125245

管理员

公社社长

4900#
 楼主 Author| 发表于 Post on 2017-10-26 11:22:24 | 只看该作者 Only view this author
erya123 发表于 2017-10-26 09:37
老师,分子轨道中NBO的贡献可以在multiwfn里面计算吗?还是必须要用新的NBO软件。

在NBO里用NBOMO关键词输出MO中NBO的系数,把系数求平方就是NBO对MO的贡献
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

6万

eV
积分
125245

管理员

公社社长

4901#
 楼主 Author| 发表于 Post on 2017-10-26 11:23:35 | 只看该作者 Only view this author
dingniu2 发表于 2017-10-26 09:12
请问老师:
1.对于复合物计算的偶极矩值是否能比较复合物稳定性的大小?
2.文献中有很多计算A与B形成复合 ...

1 不能
2 复合物的偶极矩对于讨论复合物内的相互作用没有直接用处。如果是两个单体都是有偶极的分子,讨论两个单体间的偶极-偶极相互作用对于讨论复合物的稳定性才有用
3 不可信
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

6万

eV
积分
125245

管理员

公社社长

4902#
 楼主 Author| 发表于 Post on 2017-10-26 11:24:57 | 只看该作者 Only view this author
kyuu 发表于 2017-10-26 05:48
SMD下边做的结构优化,感觉上能量是有下降的,开始的几步跑30几个cycle,然后就某一步超128步数停了,继续 ...

加大SCF=maxcyc几乎总是没用的
当前情况根据这些信息不好判断,从第9个点开始,用calcfc再继续优化看看
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

6万

eV
积分
125245

管理员

公社社长

4903#
 楼主 Author| 发表于 Post on 2017-10-26 11:25:31 | 只看该作者 Only view this author
不懂计算 发表于 2017-10-25 21:50
sob老师好,在解决SCF不收敛(L502错误)时,我在关键词里加了scf(novaracc,noincfock) int=acc2e=12, 单点能 ...

可以用

但如果和其它任务的单点能横向对比,要把int都设为相同才有可比性
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

6万

eV
积分
125245

管理员

公社社长

4904#
 楼主 Author| 发表于 Post on 2017-10-26 11:25:55 | 只看该作者 Only view this author
yangjianfeng 发表于 2017-10-25 16:20
sob老师  HF用单电子积分描述波函数,其未考虑电子相关   请问这里的电子相关具体指些什么

http://bbs.keinsci.com/forum.php?mod=viewthread&tid=4969
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

6万

eV
积分
125245

管理员

公社社长

4905#
 楼主 Author| 发表于 Post on 2017-10-26 11:27:42 | 只看该作者 Only view this author
lufenghxu 发表于 2017-10-25 21:31
老师,VMD做轨道图的时候一般是不是,先用multiwfn生成格点文件,然后就可以生成轨道图像,如http://sobere ...

VMD可以载入格点数据绘制等值面图
你如果要用VMD显示轨道图形,显然得有轨道波函数的格点数据。轨道波函数的格点数据可以用Multiwfn计算出来。
你看到的那个帖子是专门讲怎么用Multiwfn绘制NBO程序产生的轨道,因为很多人不懂怎么绘制NBO产生的各种轨道。当然量化界还有其它各种类型的轨道
更多介绍看
使用Multiwfn观看分子轨道
http://sobereva.com/269
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

本版积分规则 Credits rule

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

GMT+8, 2026-2-28 18:16 , Processed in 0.198811 second(s), 22 queries , Gzip On.

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