计算化学公社

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

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

   关闭 [复制链接 Copy URL]

1万

帖子

0

威望

9886

eV
积分
22130

Level 6 (一方通行)

12031#
发表于 Post on 2020-11-14 11:19:38 | 只看该作者 Only view this author
CayJ 发表于 2020-11-14 11:05
ORCA4.2.0, MPI10.0.12498.5,根据http://sobereva.com/542计算分子间的相互作用能,
任务1:输入文件.inp ...

第一个问题按http://sobereva.com/61的方法解决,所有方法都尝试过仍然不奏效再来问
第二个问题,错误信息提示SOSCF迭代发散,加NoSOSCF关键字把SOSCF关掉就行了,如果仍然不收敛的话,和第一个问题同样处理
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
本团队长期招收研究生,有意者可私信联系

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

12032#
 楼主 Author| 发表于 Post on 2020-11-14 23:47:19 | 只看该作者 Only view this author
CayJ 发表于 2020-11-14 11:05
ORCA4.2.0, MPI10.0.12498.5,根据http://sobereva.com/542计算分子间的相互作用能,
任务1:输入文件.inp ...

提示得很清楚SCF NOT CONVERGED AFTER 125 CYCLES         ,即SCF没收敛
论坛首页google框搜 ORCA 收敛,就能找到很多解决办法
北京科音自然科学研究中心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

175

帖子

0

威望

1057

eV
积分
1232

Level 4 (黑子)

12033#
发表于 Post on 2020-11-16 09:52:44 | 只看该作者 Only view this author
sob老师您好,隐式模型下算溶质在溶剂环境下的自由能,是不是任何情况下的模型都要加上1.89kcal/mol?如果都加这个,再算能量差的时候,相当于所有的都没加吧?

1万

帖子

0

威望

9886

eV
积分
22130

Level 6 (一方通行)

12034#
发表于 Post on 2020-11-16 14:52:53 | 只看该作者 Only view this author
329678001 发表于 2020-11-16 09:52
sob老师您好,隐式模型下算溶质在溶剂环境下的自由能,是不是任何情况下的模型都要加上1.89kcal/mol?如果 ...

每个分子加1.89kcal/mol,所以当反应前后分子数不一致的时候,就会影响能量差
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
本团队长期招收研究生,有意者可私信联系

12

帖子

0

威望

532

eV
积分
544

Level 4 (黑子)

12035#
发表于 Post on 2020-11-16 16:46:20 | 只看该作者 Only view this author
sob老师您好:能不能通过计算来比较下图两个分子的亲电或亲核加成反应的难易程度?如果可以,需要计算哪些量?

1万

帖子

0

威望

9886

eV
积分
22130

Level 6 (一方通行)

12036#
发表于 Post on 2020-11-16 22:29:53 | 只看该作者 Only view this author
xxu若白 发表于 2020-11-16 16:46
sob老师您好:能不能通过计算来比较下图两个分子的亲电或亲核加成反应的难易程度?如果可以,需要计算哪些 ...

可以参见http://sobereva.com/484
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
本团队长期招收研究生,有意者可私信联系

75

帖子

0

威望

632

eV
积分
707

Level 4 (黑子)

12037#
发表于 Post on 2020-11-16 23:17:42 | 只看该作者 Only view this author
老师好,采用molclus对下图中各个构型做团簇构型搜索时,有一些问题:1. 批量产生初始结构:xtb程序/genmer程序/...,用哪一种程序是最佳的,且耗时较短。
2. 对初始结构预优化:xtb程序/MOPAC/...,用哪一种程序是最佳的,且耗时较短。




ACS Sustainable Chem. Eng. 2020, 8, 16241−16251.png (308.01 KB, 下载次数 Times of downloads: 40)

ACS Sustainable Chem. Eng. 2020, 8, 16241−16251中的图

ACS Sustainable Chem. Eng. 2020, 8, 16241−16251中的图

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

12038#
 楼主 Author| 发表于 Post on 2020-11-17 04:47:17 | 只看该作者 Only view this author
xxu若白 发表于 2020-11-16 16:46
sob老师您好:能不能通过计算来比较下图两个分子的亲电或亲核加成反应的难易程度?如果可以,需要计算哪些 ...

如果用全局指数衡量,亲电加成看亲核指数,亲核加成看亲电指数,用Multiwfn能直接算出来,看http://sobereva.com/484
如果具体到反应位点特征,亲电反应可以用局部亲核指数、平均局部离子化能(ALIE),亲核反应可用Hirshfeld电荷、局部亲电指数

这些指数的具体定义看
概念密度泛函综述和重要文献合集 http://bbs.keinsci.com/thread-384-1-1.html
ALIE相关:
静电势与平均局部离子化能综述合集 http://bbs.keinsci.com/thread-219-1-1.html
北京科音自然科学研究中心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
积分
125155

管理员

公社社长

12039#
 楼主 Author| 发表于 Post on 2020-11-17 04:51:09 | 只看该作者 Only view this author
CayJ 发表于 2020-11-16 23:17
老师好,采用molclus对下图中各个构型做团簇构型搜索时,有一些问题:1. 批量产生初始结构:xtb程序/genmer ...

你已经会用xtb的话,用xtb跑MD产生初始构型没问题。不会用的话,genmer产生初始结构也可以,只不过由于产生时没有考虑相互作用,所以需要纳入优化的构型数需要比较多来确保优化后最终能得到对应最稳定结构。

预优化用GFN2-xTB或MOPAC跑PM7都可以,对此类体系效果半斤八两,但前者耗时更低
北京科音自然科学研究中心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

10

帖子

0

威望

465

eV
积分
475

Level 3 能力者

12040#
发表于 Post on 2020-11-17 08:07:48 | 只看该作者 Only view this author

这种结果算是振荡吗,优化是否还有必要继续

大家好,我用b3lyp/6-31+g*计算级别对一个较大分子进行优化+振动分析。已经优化了60多步,快一个星期了还没收敛(),能量变化图是这样的,这算是发生了震荡吗,是否需要停止工作改变优化方法。

搜狗截图20201117080741.png (123.77 KB, 下载次数 Times of downloads: 24)

搜狗截图20201117080741.png

1245

帖子

1

威望

3860

eV
积分
5125

Level 6 (一方通行)

12041#
发表于 Post on 2020-11-17 08:57:37 | 只看该作者 Only view this author
显然是震荡了

47

帖子

0

威望

1071

eV
积分
1118

Level 4 (黑子)

12042#
发表于 Post on 2020-11-17 09:21:47 | 只看该作者 Only view this author
建议大分子优化的时候不要加弥散。

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

12043#
发表于 Post on 2020-11-17 09:25:41 | 只看该作者 Only view this author
中间确实震荡了,但是应该已经跳出去了

10

帖子

0

威望

465

eV
积分
475

Level 3 能力者

12044#
发表于 Post on 2020-11-17 09:30:22 | 只看该作者 Only view this author
喵星大佬 发表于 2020-11-17 09:25
中间确实震荡了,但是应该已经跳出去了

我也是想再等等看,毕竟已经算了这么久了

12

帖子

0

威望

532

eV
积分
544

Level 4 (黑子)

12045#
发表于 Post on 2020-11-17 19:54:18 | 只看该作者 Only view this author
本帖最后由 xxu若白 于 2020-11-17 19:57 编辑
sobereva 发表于 2020-11-17 04:47
如果用全局指数衡量,亲电加成看亲核指数,亲核加成看亲电指数,用Multiwfn能直接算出来,看http://sober ...

谢谢sob老师,还要向您请教一下,比如说研究这两个分子(a和b)与亲核试剂(Nu-)进行加成反应。1. 对于全局亲电指数a<b,是不是说明b分子比a分子与Nu-更容易发生亲核加成反应?
2. 比较a分子的C13和b分子的C12位点对于亲核试剂Nu-的反应难易,Hirshfeld电荷和简缩福井函数如下表,q(N)这一列的值b>a,是不是就足够说明b分子的C12位点比a分子的C13位点更容易与Nu-发生反应?如果不够还需要计算哪些量?另外这种比较是合理的吗?
3. 局部亲电指数中的f+(r)是如何计算的?






本版积分规则 Credits rule

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

GMT+8, 2026-2-24 04:56 , Processed in 0.195430 second(s), 22 queries , Gzip On.

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