计算化学公社

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

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

   关闭 [复制链接 Copy URL]

87

帖子

0

威望

539

eV
积分
626

Level 4 (黑子)

15046#
发表于 Post on 2022-4-2 11:25:50 | 只看该作者 Only view this author
slope 发表于 2022-4-2 10:10
建议把输入输出文件附上

这是输入文件和输出文件
链接:https://pan.baidu.com/s/17JFANlcSoQRI7x3v_DH8HQ?pwd=za4l
提取码:za4l

24

帖子

0

威望

175

eV
积分
199

Level 3 能力者

15047#
发表于 Post on 2022-4-2 15:37:39 | 只看该作者 Only view this author
本帖最后由 zzxxcc 于 2022-4-2 15:39 编辑
sobereva 发表于 2022-4-1 01:17
AM1不支持过渡金属

谢谢sob老师!那我想先用低精度猜测过渡态 那我的oniom(pm6:ram1)中应该把ram1换成什么呢?关键词写成  # opt=(calcfc,ts,noeigen) freq oniom(pm6/genecp:???) temperature=793
赝式写成下边这样可以吗? 谢谢老师的解答!



202204021539109793..png (9.48 KB, 下载次数 Times of downloads: 13)

202204021539109793..png

6万

帖子

99

威望

6万

eV
积分
125181

管理员

公社社长

15048#
 楼主 Author| 发表于 Post on 2022-4-3 20:31:49 | 只看该作者 Only view this author
zzxxcc 发表于 2022-4-2 15:37
谢谢sob老师!那我想先用低精度猜测过渡态 那我的oniom(pm6:ram1)中应该把ram1换成什么呢?关键词写成  # ...

没有基本常识,甚至都不知道半经验没法指定基组
甭用ONIOM了,否则之后还得犯不知多少错误,白瞎折腾

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
zzxxcc + 5 谢谢sob老师!

查看全部评分 View all ratings

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

173

帖子

0

威望

1036

eV
积分
1209

Level 4 (黑子)

15049#
发表于 Post on 2022-4-3 21:30:32 | 只看该作者 Only view this author
各位老师,我想请问一下,在进行TDDFT计算后,能否根据S1的振子强度比较不同体系的PL谱峰强度?

1万

帖子

0

威望

9892

eV
积分
22140

Level 6 (一方通行)

15050#
发表于 Post on 2022-4-3 22:37:26 | 只看该作者 Only view this author
din5g 发表于 2022-4-3 14:30
各位老师,我想请问一下,在进行TDDFT计算后,能否根据S1的振子强度比较不同体系的PL谱峰强度?

PL谱峰强度正比于体系吸收的光子数乘以荧光量子产率,体系吸收的光子数正比于垂直吸收的振子强度,荧光量子产率等于 辐射跃迁速率常数/(辐射跃迁速率常数+非辐射跃迁速率常数) ,辐射跃迁速率常数正比于垂直发射的振子强度,非辐射跃迁速率常数和普通TDDFT算出来的任何量都没有简单关系,需要用MOMAP等软件单独计算。
由以上分析可知,当且仅当你知道这些体系垂直吸收的振子强度相似,且非辐射跃迁速率常数也相似时,才能根据垂直发射的振子强度比较PL谱峰强度。但是这个条件几乎不可能成立,因为很少能遇到辐射跃迁速率常数差很多但是非辐射跃迁速率常数一样的体系,如果有,也是你算出来非辐射跃迁速率常数确实一样你才知道,这个是没办法从普通TDDFT计算结果或者分子结构看出来的。而如果你已经算出来了非辐射跃迁速率常数,你就可以直接算荧光量子产率,用荧光量子产率比较PL谱峰强度了,这个比用振子强度靠谱得多。
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

68

帖子

0

威望

1646

eV
积分
1714

Level 5 (御坂)

15051#
发表于 Post on 2022-4-4 15:25:13 | 只看该作者 Only view this author
请问老师,我利用ORCA5.0.1的STEOM-DLPNO-CCSD计算多环芳烃的三重态激发能,发现较实验明显的低估
是不是跟我下面的提示有关呢

Percentage Active Character     95.91

Warning:: the state may have not converged with respect to active space
-------------------- Handle with Care --------------------

这是我的关键词
%mdci
DoTriplet true
maxiter 500
nroots 3
NDav 50
OThresh 0.00005
VThresh 0.00005
end

1万

帖子

0

威望

9892

eV
积分
22140

Level 6 (一方通行)

15052#
发表于 Post on 2022-4-4 16:27:25 | 只看该作者 Only view this author
心向暖阳 发表于 2022-4-4 08:25
请问老师,我利用ORCA5.0.1的STEOM-DLPNO-CCSD计算多环芳烃的三重态激发能,发现较实验明显的低估
是不是 ...

计算用的基组是什么?低估多少eV?计算考虑溶剂了吗?
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

68

帖子

0

威望

1646

eV
积分
1714

Level 5 (御坂)

15053#
发表于 Post on 2022-4-4 16:30:54 | 只看该作者 Only view this author
wzkchem5 发表于 2022-4-4 16:27
计算用的基组是什么?低估多少eV?计算考虑溶剂了吗?

谢谢您的回复,我是想算detEST的,三重态激发能低估了0.5eV以上,所以效果很不好

用的def2-SVP基组,因为再大些内存就不够了,但是以我TDDFT的经验,基组对激发能影响并不是很大

我测试的结果是溶剂考虑与否,影响不大的

1万

帖子

0

威望

9892

eV
积分
22140

Level 6 (一方通行)

15054#
发表于 Post on 2022-4-4 18:29:04 | 只看该作者 Only view this author
心向暖阳 发表于 2022-4-4 09:30
谢谢您的回复,我是想算detEST的,三重态激发能低估了0.5eV以上,所以效果很不好

用的def2-SVP基组, ...

STEOM-DLPNO-CCSD对基组还是比较敏感的,毕竟是post-HF方法,比DFT要敏感得多。建议可以试试把基组加大到def2-TZVP或def2-TZVP(-f)试试
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

68

帖子

0

威望

1646

eV
积分
1714

Level 5 (御坂)

15055#
发表于 Post on 2022-4-4 21:16:22 | 只看该作者 Only view this author
wzkchem5 发表于 2022-4-4 18:29
STEOM-DLPNO-CCSD对基组还是比较敏感的,毕竟是post-HF方法,比DFT要敏感得多。建议可以试试把基组加大到 ...

我再减小核数试试
谢谢~

200

帖子

1

威望

1325

eV
积分
1545

Level 5 (御坂)

15056#
发表于 Post on 2022-4-4 23:28:04 | 只看该作者 Only view this author
请问下各位老师,我笔记本电脑上同时安装了G16 B01,GaussView 6和G09 D01,GaussView 5两个版本的高斯和GaussView。请问我用GaussView 5打开GaussView 6产生的gif文件,或者用GaussView 6打开GaussView 5产生的gjf文件,会出现问题(比如文本内容变化等)吗?现在有点虚......
计算小白 水平很菜
随时可喷 万望指点

1562

帖子

0

威望

5005

eV
积分
6567

Level 6 (一方通行)

15057#
发表于 Post on 2022-4-5 08:46:14 | 只看该作者 Only view this author
WaterMirror 发表于 2022-4-4 23:28
请问下各位老师,我笔记本电脑上同时安装了G16 B01,GaussView 6和G09 D01,GaussView 5两个版本的高斯和Gau ...

http://sobereva.com/439

  1. 应当注意Gaussian和GaussView的版本兼容性。根据笔者经验,Gaussian 09建议搭配GaussView 5.0.9使用,Gaussian 16建议搭配GaussView 6.0.16(目前的最新版本)使用,这样兼容性是最好的,否则会有很多问题。比如GaussView 6.0.16打开Gaussian 09的IRC任务的帧号顺序是错乱的,而GaussView 5.0.9无法正确载入Gaussian 16的振动分析输出文件等。
复制代码
又菜又爱玩

410

帖子

0

威望

2519

eV
积分
2929

Level 5 (御坂)

15058#
发表于 Post on 2022-4-5 16:02:40 | 只看该作者 Only view this author
请问一下,有哪些方法可以计算电负性

6万

帖子

99

威望

6万

eV
积分
125181

管理员

公社社长

15059#
 楼主 Author| 发表于 Post on 2022-4-5 21:46:27 | 只看该作者 Only view this author
mfdsrax2 发表于 2022-4-5 16:02
请问一下,有哪些方法可以计算电负性

看具体什么体系
如果是孤立体系,直接按标准Mulliken电负性基于IP和EA算。
北京科音自然科学研究中心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

68

帖子

0

威望

1646

eV
积分
1714

Level 5 (御坂)

15060#
发表于 Post on 2022-4-6 20:44:25 | 只看该作者 Only view this author
wzkchem5 发表于 2022-4-4 18:29
STEOM-DLPNO-CCSD对基组还是比较敏感的,毕竟是post-HF方法,比DFT要敏感得多。建议可以试试把基组加大到 ...

老师,我将基组提升到def2-TZVP(-f),三重态的激发能提升了0.05eV
我还要加大基组么,但是现在提示Percentage Active Character     80.30

本版积分规则 Credits rule

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

GMT+8, 2026-2-26 05:14 , Processed in 0.183146 second(s), 25 queries , Gzip On.

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