计算化学公社

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

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

   关闭 [复制链接 Copy URL]

44

帖子

0

威望

193

eV
积分
237

Level 3 能力者

6871#
发表于 Post on 2018-6-4 11:25:38 | 只看该作者 Only view this author
sobereva 发表于 2018-6-4 10:35
没法变
非要更大的,就导出数据到sigmaplot里重新绘图,届时随意调整。或者自行修改源代码。或者自行p ...

好的,谢谢老师

15

帖子

0

威望

89

eV
积分
104

Level 2 能力者

6872#
发表于 Post on 2018-6-4 11:37:27 | 只看该作者 Only view this author
Sob老师,我有个问题想请教您。模拟计算一个结构的光电子能谱时出现了负值,这种情况应该怎么解决啊?我算的是TD(50-50,nstates30))第一个激发态算出来是负数
Excited State   1:      Triplet-B1     -0.1698 eV    -7302.90 nm  f=-0.0000  <S**2>=2.000
Excited State   2:      Triplet-A2     0.8880 eV 1396.19 nm  f=0.0000  <S**2>=2.000
应该怎样正确的处理这个问题呢

146

帖子

0

威望

1329

eV
积分
1475

Level 4 (黑子)

6873#
发表于 Post on 2018-6-4 12:39:51 | 只看该作者 Only view this author
老师,绘制体系减去片段EDD时,能用“opt时的molden” 减去 “片段的自洽的molden”吗?这个总体系还需要算sp吗?

6万

帖子

99

威望

6万

eV
积分
125190

管理员

公社社长

6874#
 楼主 Author| 发表于 Post on 2018-6-4 13:57:42 | 只看该作者 Only view this author
刘小英雄steven 发表于 2018-6-4 12:39
老师,绘制体系减去片段EDD时,能用“opt时的molden” 减去 “片段的自洽的molden”吗?这个总体系还需要算 ...


不需要
北京科音自然科学研究中心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
积分
125190

管理员

公社社长

6875#
 楼主 Author| 发表于 Post on 2018-6-4 13:58:47 | 只看该作者 Only view this author
yu947937636 发表于 2018-6-4 11:11
sob老师,看了您的博文写到6-311G*基组有时会比6-31+G*得到的结果更好,在网上也看到一篇博文计算的磷氧键 ...


什么时候该加弥散这里说了
谈谈弥散函数和“月份”基组
http://sobereva.com/119

凡是这里没说非得加弥散的情况,6-31+G*的结果绝大多数情况都不会比6-311G*更好,而计算还更贵

别盲目效仿文献,50%以上的文献用的基组都不是最理想的基组
北京科音自然科学研究中心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
积分
125190

管理员

公社社长

6876#
 楼主 Author| 发表于 Post on 2018-6-4 13:59:15 | 只看该作者 Only view this author
小尼姑 发表于 2018-6-4 11:37
Sob老师,我有个问题想请教您。模拟计算一个结构的光电子能谱时出现了负值,这种情况应该怎么解决啊?我算 ...


北京科音自然科学研究中心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
积分
125190

管理员

公社社长

6877#
 楼主 Author| 发表于 Post on 2018-6-4 14:03:46 | 只看该作者 Only view this author
haha 发表于 2018-6-4 11:17
sob老师,我的分子是带有ONO基团的烷氧自由基,根据tddft计算激发态得到第四态符合实验波长,对应C-O成键轨 ...

如置顶的新人必读贴所示,上传较大文本型文件前一定要先压缩

本帖只用于回复简单问题,涉及较多输入输出文件的情况请单独发帖
北京科音自然科学研究中心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

11

帖子

0

威望

336

eV
积分
347

Level 3 能力者

6878#
发表于 Post on 2018-6-4 14:16:51 | 只看该作者 Only view this author
sobereva 发表于 2018-6-4 13:58
什么时候该加弥散这里说了
谈谈弥散函数和“月份”基组
http://sobereva.com/119

好的,谢谢老师的解答

15

帖子

0

威望

89

eV
积分
104

Level 2 能力者

6879#
发表于 Post on 2018-6-4 15:33:24 | 只看该作者 Only view this author
谢谢Sob老师的解答。老师,我还想问下一起算单三激发态和分两步算单重态和三重态结果一样吗?就是写成下边这样计算
TD(50-50,NSTATES=30)
TD(NSTATES=30) 先算单重态再算三重态这样的结果是不是该保持一致啊?或者一般都是一起算的?

15

帖子

0

威望

89

eV
积分
104

Level 2 能力者

6880#
发表于 Post on 2018-6-4 15:36:46 | 只看该作者 Only view this author
小尼姑 发表于 2018-6-4 15:33
谢谢Sob老师的解答。老师,我还想问下一起算单三激发态和分两步算单重态和三重态结果一样吗?就是写成下边 ...

Excited State   1:  2.917-?Sym   0.3163 eV 3919.85 nm  f=0.0000  <S**2>=1.877
Excited State   2:  2.992-?Sym    0.9368 eV 1323.42 nm  f=0.0000  <S**2>=1.988
用波函数稳定性去掉负值一起算单三激发态后都没有算出哪个是单重态哪个是三重态,已经变成一般的数据了,所以怀疑是不是算的有问题?

78

帖子

0

威望

814

eV
积分
892

Level 4 (黑子)

6881#
发表于 Post on 2018-6-4 17:46:44 | 只看该作者 Only view this author
请问CCSD(T)/jul-cc-PVTZ计算SCF迭代不收敛,可以尝试用cc-pVTZ收敛的波函数当初猜,再用jul-cc-PVTZ基组去读取波函数初猜继续计算单点能吗?具体是不是这样:
%chk=1.chk
#p CCSD(T)/cc-PVTZ 算完单点能以后用得到的.chk文件续算。

%chk=1.chk
#p CCSD(T)/jul-cc-PVTZ guess=read
此时坐标还是用第一步的坐标,还是将第一步的cc-PVTZ算的另存为gjf文件的坐标呢。算单点能应该是不改变结构的坐标吧?十分谢谢Sob老师
自己决定自己的DNA

6万

帖子

99

威望

6万

eV
积分
125190

管理员

公社社长

6882#
 楼主 Author| 发表于 Post on 2018-6-5 06:27:43 | 只看该作者 Only view this author
小尼姑 发表于 2018-6-4 15:33
谢谢Sob老师的解答。老师,我还想问下一起算单三激发态和分两步算单重态和三重态结果一样吗?就是写成下边 ...

结果一样(早期G09版本不一样,是程序bug所致)
如果单三重态都算,建议用50-50一起算
北京科音自然科学研究中心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
积分
125190

管理员

公社社长

6883#
 楼主 Author| 发表于 Post on 2018-6-5 06:28:33 | 只看该作者 Only view this author
小尼姑 发表于 2018-6-4 15:36
Excited State   1:  2.917-?Sym   0.3163 eV 3919.85 nm  f=0.0000  =1.877
Excited State   2:  2.992 ...

参考态是开壳层时候用singlet,triplet,50-50没有意义,没法得到自旋纯态。这里说了
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314
北京科音自然科学研究中心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
积分
125190

管理员

公社社长

6884#
 楼主 Author| 发表于 Post on 2018-6-5 06:29:26 | 只看该作者 Only view this author
QuantumHan 发表于 2018-6-4 17:46
请问CCSD(T)/jul-cc-PVTZ计算SCF迭代不收敛,可以尝试用cc-pVTZ收敛的波函数当初猜,再用jul-cc-PVTZ基组去 ...

关键词都没错
坐标还是第一步的坐标。单点不改变结构
北京科音自然科学研究中心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

95

帖子

0

威望

237

eV
积分
332

Level 3 能力者

6885#
发表于 Post on 2018-6-5 08:53:14 | 只看该作者 Only view this author
sobereva 发表于 2018-6-4 05:02
图中体系很简单,直接对R(N-O)进行扫描,就能得到两个态的势能曲线
而对于复杂情况,需要对优化/搜索M ...

老师,您指的某几个几何变量是单独扫描键长变化吗?最后的图怎么做啊。好多个键长变化和能量的对应关系。

本版积分规则 Credits rule

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

GMT+8, 2026-2-26 15:48 , Processed in 0.624499 second(s), 21 queries , Gzip On.

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