计算化学公社

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

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

   关闭 [复制链接 Copy URL]

477

帖子

0

威望

887

eV
积分
1364

Level 4 (黑子)

16861#
发表于 Post on 2022-10-21 21:04:05 | 只看该作者 Only view this author
wzkchem5 发表于 2022-10-21 20:58
B4-10比文献长得太多了,估计是因为和文献理论级别不统一,比如忘加色散校正了

老师我用的是um062x 6-311**,所以没加。文献是MP2 cc-pVTZ;

1万

帖子

0

威望

8998

eV
积分
20752

Level 6 (一方通行)

16862#
发表于 Post on 2022-10-21 21:29:20 | 只看该作者 Only view this author
糖糖糖豆9988 发表于 2022-10-21 14:04
老师我用的是um062x 6-311**,所以没加。文献是MP2 cc-pVTZ;

那可以检查一下是不是格点太粗糙了,M06-2X对格点很敏感,如果格点不够大,势能面会坑坑洼洼的,可能你的分子会陷在一个小坑里
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?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

16863#
 楼主 Author| 发表于 Post on 2022-10-21 22:20:15 | 只看该作者 Only view this author
小翟w 发表于 2022-10-21 19:41
老师们,晚上好,请教个问题,这是输入文件和输出的截图(文件太大,传不上来),第一次没有opt=maxcycle=2 ...

都已经给你合并到这个答疑专贴里了,你又重新发个新帖子干嘛?
过于简单初级的问题不要单独发主贴,避免拉低论坛帖子平均水准。初级问题一律在这里以跟帖方式提问。别人早就回复过你了,好好看16853和16854楼
北京科音自然科学研究中心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

78

帖子

0

威望

244

eV
积分
322

Level 3 能力者

16864#
发表于 Post on 2022-10-22 08:15:54 | 只看该作者 Only view this author
好,第一次,现在知道了,我以为帖子没了又发了一遍,谢谢老师们的指导

78

帖子

0

威望

244

eV
积分
322

Level 3 能力者

16865#
发表于 Post on 2022-10-22 08:21:42 | 只看该作者 Only view this author
wzkchem5 发表于 2022-10-21 19:45
你贴的输出文件部分不包含有效的报错信息。仔细阅读http://bbs.keinsci.com/thread-4829-1-1.html,了解怎 ...

好的,谢谢老师的指导

78

帖子

0

威望

244

eV
积分
322

Level 3 能力者

16866#
发表于 Post on 2022-10-22 08:22:51 | 只看该作者 Only view this author
wxsxjj 发表于 2022-10-21 19:48
向上检索  Optimization stop 得到opt stop原因
不如先去看李大博士整理的帖子(http://bbs.k ...

好的,谢谢老师的指导

477

帖子

0

威望

887

eV
积分
1364

Level 4 (黑子)

16867#
发表于 Post on 2022-10-22 08:57:40 | 只看该作者 Only view this author
wzkchem5 发表于 2022-10-21 15:53
按http://sobereva.com/278的方法消虚频,如果消不掉,说明是数值误差问题,把优化收敛限设严、格点加大 ...

#p opt=(tight,recalc=5) freq um062x 6-311g(d,p) nosymm收敛,谢谢老师

477

帖子

0

威望

887

eV
积分
1364

Level 4 (黑子)

16868#
发表于 Post on 2022-10-22 11:12:37 | 只看该作者 Only view this author
请问老师,对于图片里的①环丙烷氢提取反应 cC3H6 + O = cC3H5 + OH 的高压极限速率常数HPL如何计算呢?②因为第一步涉及无能垒的反应前复合物RC-abs的形成,所以高压极限速率常数HPL的计算需要用变分过渡态理论吗?③文献里该反应的HPL计算 “ The high-pressure limiting (HPL) rate constants were calculated by conventional transition state theory (TST) over a temperature range of 500–2000 K.  ” 以及“All rate constant calculations were carried out by the KiSThelP2019 program ”,作者的意思是直接采用kisthelp里的双分子反应,输入过渡态的.kinp/.out文件、环丙烷和H原子的.out文件计算吗?麻烦老师们看看,谢谢

202210221048106351..png (56.1 KB, 下载次数 Times of downloads: 7)

202210221048106351..png

449

帖子

0

威望

1559

eV
积分
2008

Level 5 (御坂)

16869#
发表于 Post on 2022-10-22 12:09:37 | 只看该作者 Only view this author
本帖最后由 cokie 于 2022-11-7 11:15 编辑

各位大佬老师好,关于Gview6默认的自旋多重度问题,要请教各位。

我的分子式为:C46 H46 N2 O4 Br2,虽然是季铵盐,但很明显体系的电子数是偶数,但是导入Gview里后,默认的自旋多重度只有2,4等等可以选择,都没有1这个选项。所以有两个问题想请教各位:

1. 会让Gview默认我输入的分子自旋多重度为2,是因为我在绘制结构亦或是其他地方出了问题吗?

2. 如果只是Gview自身默认的情况不准确,我是否可以或者说必须在gjf文件里将“0 2”改成“0 1”呢?

谢谢各位!

1016

帖子

0

威望

2096

eV
积分
3112

Level 5 (御坂)

天上地下三山亖海五湖六合八荒九州之泛粤大典編輯委員會編外委員

16870#
发表于 Post on 2022-10-22 12:41:19 | 只看该作者 Only view this author
cokie 发表于 2022-10-22 12:09
各位大佬老师好,关于Gview6默认的自旋多重度问题,要请教各位。

我的分子式为:C46 H46 N2 O4 Br2,虽 ...

这种情况多半是什么地方多H或缺H,如果你改成 0 1,Gaussian马上说电荷数和自旋多重度不配,马上报错。
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

山河自落蒼梧月,風雨猶驅草木兵。
san huɔɞ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ

449

帖子

0

威望

1559

eV
积分
2008

Level 5 (御坂)

16871#
发表于 Post on 2022-10-22 14:06:45 | 只看该作者 Only view this author
chands 发表于 2022-10-22 12:41
这种情况多半是什么地方多H或缺H,如果你改成 0 1,Gaussian马上说电荷数和自旋多重度不配,马上报错。

明白了,谢谢答疑!我重新构建一个图仔细检查一下

85

帖子

0

威望

1463

eV
积分
1548

Level 5 (御坂)

16872#
发表于 Post on 2022-10-22 15:08:16 | 只看该作者 Only view this author
老师请问,对于基态不同的极小点结构来算激发态,得到的振子强度变化大么?

1万

帖子

0

威望

8998

eV
积分
20752

Level 6 (一方通行)

16873#
发表于 Post on 2022-10-22 15:59:00 | 只看该作者 Only view this author
liang666 发表于 2022-10-22 08:08
老师请问,对于基态不同的极小点结构来算激发态,得到的振子强度变化大么?

可能大也可能不大。举两个简单的例子:
(1)如果这个分子只是一个烷基取代基转了一下,那肯定几乎不影响振子强度;
(2)如果这个分子有两个相同的共轭基团,当两者平行时跃迁偶极矩刚好符号相反抵消,当两者垂直时跃迁偶极矩没办法抵消,那么平行和垂直时的振子强度一个为0,一个很大
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?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

30

帖子

0

威望

228

eV
积分
258

Level 3 能力者

16874#
发表于 Post on 2022-10-22 17:47:04 | 只看该作者 Only view this author
老师们好,我在用高斯优化一个平面分子2,3,6,7,10,11-六氨基三苯的时候发现优化着优化着就变成一个碗装的分子了。我用的关键词是# opt freq scrf=smd b3lyp 6-31g* em=gd3bj。我想请问老师是哪里出了问题呢。还有就是在遇到未知的结构的时候如果像这样明显偏离了,但优化出的结构收敛了又没有虚频的话我怎么判断他就一定是正确的呢?谢谢老师

202210221743411439..png (88.67 KB, 下载次数 Times of downloads: 9)

202210221743411439..png

202210221743502135..png (86.47 KB, 下载次数 Times of downloads: 7)

202210221743502135..png

202210221743552875..png (68.95 KB, 下载次数 Times of downloads: 8)

202210221743552875..png

1万

帖子

0

威望

8998

eV
积分
20752

Level 6 (一方通行)

16875#
发表于 Post on 2022-10-22 18:14:47 | 只看该作者 Only view this author
white-album 发表于 2022-10-22 10:47
老师们好,我在用高斯优化一个平面分子2,3,6,7,10,11-六氨基三苯的时候发现优化着优化着就变成一个碗装的分 ...

苯环上的氢漏画了。画上氢原子,重新算
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?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

本版积分规则 Credits rule

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

GMT+8, 2025-8-15 16:59 , Processed in 0.160075 second(s), 22 queries , Gzip On.

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