计算化学公社

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

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

   关闭 [复制链接 Copy URL]

45

帖子

0

威望

95

eV
积分
140

Level 2 能力者

2281#
发表于 Post on 2016-10-12 20:43:46 | 只看该作者 Only view this author
sobereva 发表于 2016-10-12 14:22
实现不了。只记录最后一步的。
你只能自己提取第n步的结构算个单点

好的,谢谢

45

帖子

0

威望

95

eV
积分
140

Level 2 能力者

2282#
发表于 Post on 2016-10-12 20:47:16 | 只看该作者 Only view this author
在用cas方法的时候,最后可视化轨道,并不显示能量,我记得好像是cas迭代过程中并不产生fock矩阵。
那样的话,在scf迭代不收敛的时候,关于调试fock矩阵的关键词,比如能级提升法vshift等等,是不是用了也不会对cas方法里的scf迭代产生有利的影响了。

720

帖子

0

威望

1252

eV
积分
1972

Level 5 (御坂)

2283#
发表于 Post on 2016-10-12 21:12:17 | 只看该作者 Only view this author
sobereva 发表于 2016-10-12 14:20

原子电荷计算方法的对比
http://www.whxb.pku.edu.cn/CN/abstract/abstract27818.shtml

Sob老师,您给我推荐的那篇文章我看了一下,但是我还是有疑问。
弥散函数会使重叠积分失准,某些原本只属于单个原子的轨道被误划分至其它原子的轨道中。Mulliken方法把重叠的部分平均分配给两个原子,从而造成原子电荷的计算错误。
其它利用重叠积分的划分来计算原子电荷的方法是不是也不能用弥散函数啊?
谢谢老师

407

帖子

0

威望

856

eV
积分
1263

Level 4 (黑子)

2284#
发表于 Post on 2016-10-12 21:22:24 | 只看该作者 Only view this author
老师,混合基组的频率矫正因子,应该怎么给啊?

6万

帖子

99

威望

6万

eV
积分
125179

管理员

公社社长

2285#
 楼主 Author| 发表于 Post on 2016-10-12 23:00:49 | 只看该作者 Only view this author
小书童 发表于 2016-10-12 21:22
老师,混合基组的频率矫正因子,应该怎么给啊?


用大多数原子的基组的频率校正因子
北京科音自然科学研究中心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
积分
125179

管理员

公社社长

2286#
 楼主 Author| 发表于 Post on 2016-10-12 23:01:53 | 只看该作者 Only view this author
shangfang 发表于 2016-10-12 20:47
在用cas方法的时候,最后可视化轨道,并不显示能量,我记得好像是cas迭代过程中并不产生fock矩阵。
那样的 ...

SCF不收敛那些处理方法大多无法用在CASSCF迭代上
北京科音自然科学研究中心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
积分
125179

管理员

公社社长

2287#
 楼主 Author| 发表于 Post on 2016-10-12 23:04:02 | 只看该作者 Only view this author
lastzealot 发表于 2016-10-12 21:12
Sob老师,您给我推荐的那篇文章我看了一下,但是我还是有疑问。
弥散函数会使重叠积分失准,某些原本只 ...

“弥散函数会使重叠积分失准”这句话是错的,对基础理论没有透彻理解
不管有没有弥散函数,基函数间重叠积分都是精确计算的,只不过导致Mulliken布居方法在划分交叉项时候不合理,缺乏化学意义。
文中仔细比较了不同原子电荷计算方法对基组的依赖性,明显看得出诸如NPA电荷,虽然也是直接基于波函数,但不怕弥散函数。
北京科音自然科学研究中心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

45

帖子

0

威望

95

eV
积分
140

Level 2 能力者

2288#
发表于 Post on 2016-10-13 08:20:19 | 只看该作者 Only view this author
sobereva 发表于 2016-10-12 23:01
SCF不收敛那些处理方法大多无法用在CASSCF迭代上

老师,那请问有什么资料推荐,可以解决cas迭代的不收敛问题呢。

407

帖子

0

威望

856

eV
积分
1263

Level 4 (黑子)

2289#
发表于 Post on 2016-10-13 08:26:01 | 只看该作者 Only view this author
sobereva 发表于 2016-10-12 23:00
用大多数原子的基组的频率校正因子

恩恩,多谢老师

720

帖子

0

威望

1252

eV
积分
1972

Level 5 (御坂)

2290#
发表于 Post on 2016-10-13 08:47:17 | 只看该作者 Only view this author
sobereva 发表于 2016-10-12 23:04
“弥散函数会使重叠积分失准”这句话是错的,对基础理论没有透彻理解
不管有没有弥散函数,基函数间重叠 ...

非常感谢sob老师的解答

97

帖子

0

威望

337

eV
积分
434

Level 3 能力者

2291#
发表于 Post on 2016-10-13 13:35:22 | 只看该作者 Only view this author
请问老师,

我优化了激发态, 然后想再计算激发态优化结构下的频率, 但是不知道写的对不对, 想请教一下

优化激发态时如下,

%mem=20
%nprocshared=8
%chk=***.chk
# opt td=(singlets, nstates=10) b3ylp/6-31g(d,p)

我在gview里打开优化后的log文件, 计算频率设置如下,

%mem=20
%nprocshared=8
%chk=***.chk
# freq td=(singlets, nstates=10) b3ylp/6-31g(d,p)

不知道这样写对不对? 想问要不要写上geom=check guess=read? 如果写上的话和不写计算结果有没有什么区别?

ps, 没有写root=1是不是没有影响,因为默认就是这个?

6万

帖子

99

威望

6万

eV
积分
125179

管理员

公社社长

2292#
 楼主 Author| 发表于 Post on 2016-10-13 18:22:49 | 只看该作者 Only view this author
shangfang 发表于 2016-10-13 08:20
老师,那请问有什么资料推荐,可以解决cas迭代的不收敛问题呢。

调整初猜轨道、活性空间设置,或者更好的CASSCF程序如molcas/molpro
北京科音自然科学研究中心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
积分
125179

管理员

公社社长

2293#
 楼主 Author| 发表于 Post on 2016-10-13 18:34:01 | 只看该作者 Only view this author
ccccn345 发表于 2016-10-13 13:35
请问老师,

我优化了激发态, 然后想再计算激发态优化结构下的频率, 但是不知道写的对不对, 想请教一下


b3ylp是b3lyp
root=1不用写
nstates=10也没必要,默认算3个态够了
写geom=check的话就不用写坐标了,如果直接写了坐标就别写这个
guess=read写了之后从chk中读取初猜可以把第一次的SCF部分的时间省下来,不过相对于频率计算总耗时来说微乎其微。
北京科音自然科学研究中心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

97

帖子

0

威望

337

eV
积分
434

Level 3 能力者

2294#
发表于 Post on 2016-10-14 00:52:59 | 只看该作者 Only view this author
sobereva 发表于 2016-10-13 18:34
b3ylp是b3lyp
root=1不用写
nstates=10也没必要,默认算3个态够了

老师要求计算10个, 就加了nstates=10
b3lyp笔误, 非常感谢

2

帖子

0

威望

21

eV
积分
23

Level 1 能力者

2295#
发表于 Post on 2016-10-14 11:28:16 | 只看该作者 Only view this author

读取chk

请问,上一步的计算有分碎片,但接下来的一步计算没分碎片,那下一步的计算能读取上一步是chk文件吗

本版积分规则 Credits rule

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

GMT+8, 2026-2-25 16:36 , Processed in 0.594951 second(s), 26 queries , Gzip On.

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