在5月22-26日期间,本论坛仅限等级≥level 4的成员发新的主题贴,带来的不便请谅解!(此期间请勿给管理员Sobereva发私信问学术问题)

计算化学公社

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

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

   关闭 [复制链接 Copy URL]

26

帖子

0

威望

159

eV
积分
185

Level 3 能力者

16186#
发表于 Post on 2022-8-7 19:05:53 | 只看该作者 Only view this author
关于QM/MM MD的疑问:
1.QM/MM MD 是如何迭代的?每次一QM/MM MD迭代是有哪些关键步骤?
2.QM/MM MD的能量是如何进行计算的?个人理解MD是通过对对构象内能通过概率转化为自由能,如果在一个case中QM区域没有MM力学参数,如何和MM部分的能量耦合起来,计算得到自由能?
以上问题尝试着自己去寻找答案但是无奈网上似乎相关话题实在是太少,也希望老师能给个小白系统学习理解QM/MM MD的途径。
谢谢老师!

1万

帖子

0

威望

1万

eV
积分
22586

Level 6 (一方通行)

16187#
发表于 Post on 2022-8-7 22:35:38 | 只看该作者 Only view this author
wanghaoyu 发表于 2022-8-7 12:05
关于QM/MM MD的疑问:
1.QM/MM MD 是如何迭代的?每次一QM/MM MD迭代是有哪些关键步骤?
2.QM/MM MD的能 ...

第一个问题:看http://sobereva.com/576,这个讲的是AIMD,但QM/MM MD可以认为是AIMD的一种,原理是一样的
第二个问题参见https://bdf-manual.readthedocs.i ... %20Guide.html#qm-mm
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
积分
127029

管理员

公社社长

16188#
 楼主 Author| 发表于 Post on 2022-8-8 04:01:34 | 只看该作者 Only view this author
wanghaoyu 发表于 2022-8-7 19:05
关于QM/MM MD的疑问:
1.QM/MM MD 是如何迭代的?每次一QM/MM MD迭代是有哪些关键步骤?
2.QM/MM MD的能 ...

QM/MM的理论介绍文章或书籍多了去了,比如Introduction to Computational Chemistry (3ed,Frank Jensen,2017)就有介绍
北京科音自然科学研究中心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

72

帖子

0

威望

2343

eV
积分
2415

Level 5 (御坂)

16189#
发表于 Post on 2022-8-8 09:04:29 | 只看该作者 Only view this author
老师好,我需要做一个离子液体稳定构型的搜索;
阳离子有30个原子左右(C、H、O、N组成),阴离子有7个原子(C、H、O组成);
我想用molclus+Gaussian做构型搜索;
请问老师,我这里用B3LYP-D3(BJ)/6-311G** 可以吗?

6

帖子

0

威望

17

eV
积分
23

Level 1 能力者

16190#
发表于 Post on 2022-8-8 09:48:26 | 只看该作者 Only view this author

如何运用高斯软件进行结合能的计算?

我想要通过高斯软件判断哪种单体与模板结合能力最好,现所遇到的问题是不太懂高斯软件如何进行结合能的计算。哪位老师可以给我提供一下计算教程和计算时所需注意事项吗?

6

帖子

0

威望

17

eV
积分
23

Level 1 能力者

16191#
发表于 Post on 2022-8-8 09:57:21 | 只看该作者 Only view this author
请问结合能的计算与单点能的计算有什么不同?

1016

帖子

0

威望

2128

eV
积分
3144

Level 5 (御坂)

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

16192#
发表于 Post on 2022-8-8 10:13:14 | 只看该作者 Only view this author
忘不掉的回忆 发表于 2022-8-8 09:04
老师好,我需要做一个离子液体稳定构型的搜索;
阳离子有30个原子左右(C、H、O、N组成),阴离子有7个原 ...

建议分步优化,先用低级别优化,再用高级别。
最后一步优化至少阴离子要加弥散。
第一步用B3LYP-D3(BJ)/6-311G**还是太高,建议第一步用半经验。
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

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

6万

帖子

99

威望

6万

eV
积分
127029

管理员

公社社长

16193#
 楼主 Author| 发表于 Post on 2022-8-8 10:14:43 | 只看该作者 Only view this author
XYWO 发表于 2022-8-8 09:48
我想要通过高斯软件判断哪种单体与模板结合能力最好,现所遇到的问题是不太懂高斯软件如何进行结合能的计算 ...

公社论坛首页google框搜 Gaussian 结合能 计算
花半个小时若还没搞明白,再在这里具体问、描述具体不会的地方

这俩至少要看
谈谈BSSE校正与Gaussian对它的处理
http://sobereva.com/46
在ORCA中做counterpoise校正并计算分子间结合能的例子
http://sobereva.com/542http://bbs.keinsci.com/thread-16585-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
积分
127029

管理员

公社社长

16194#
 楼主 Author| 发表于 Post on 2022-8-8 10:14:51 | 只看该作者 Only view this author
XYWO 发表于 2022-8-8 09:57
请问结合能的计算与单点能的计算有什么不同?

有别人回复之前若需要对帖子进行修改、补充,应直接编辑原帖,不要通过回帖进行补充,这点在置顶的新社员必读贴里明确说了。
北京科音自然科学研究中心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
积分
127029

管理员

公社社长

16195#
 楼主 Author| 发表于 Post on 2022-8-8 10:17:14 | 只看该作者 Only view this author
忘不掉的回忆 发表于 2022-8-8 09:04
老师好,我需要做一个离子液体稳定构型的搜索;
阳离子有30个原子左右(C、H、O、N组成),阴离子有7个原 ...

我的所有molclus的教程/例子里都不是一步得到结果,都是分步进行,以达到最佳性价比。把例子都仔细看了、思想充分领会了,提问时说具体哪个环节用哪个级别
B3LYP-D3(BJ)/6-311G**对优化浪费,对单点不够
北京科音自然科学研究中心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

72

帖子

0

威望

2343

eV
积分
2415

Level 5 (御坂)

16196#
发表于 Post on 2022-8-8 10:27:33 | 只看该作者 Only view this author
chands 发表于 2022-8-8 10:13
建议分步优化,先用低级别优化,再用高级别。
最后一步优化至少阴离子要加弥散。
第一步用B3LYP-D3(BJ) ...

好的,谢谢老师,我这里没有表达清楚;
我这里想问的是做完半经验优化后再用Gaussian做最终优化

72

帖子

0

威望

2343

eV
积分
2415

Level 5 (御坂)

16197#
发表于 Post on 2022-8-8 10:38:55 | 只看该作者 Only view this author
sobereva 发表于 2022-8-8 10:17
我的所有molclus的教程/例子里都不是一步得到结果,都是分步进行,以达到最佳性价比。把例子都仔细看了、 ...

谢谢老师的回复,下次提问问题时我一定注意;
我这里想问的是根据您的这篇文章(http://bbs.keinsci.com/thread-577-1-1.html)通过使用molclus+MOPAC做完初步优化后,
再用molclus+Gaussian做优化时泛函和基组的选择。
那么请问老师既然B3LYP-D3(BJ)/6-311G** 对优化浪费,那我可不可以先用一下小基组做一下优化,比如6-31G** ,
优化排序后筛选出几个构型再用大一点的基组优化一下。
谢谢老师

1万

帖子

0

威望

1万

eV
积分
22586

Level 6 (一方通行)

16198#
发表于 Post on 2022-8-8 14:57:11 | 只看该作者 Only view this author
忘不掉的回忆 发表于 2022-8-8 03:38
谢谢老师的回复,下次提问问题时我一定注意;
我这里想问的是根据您的这篇文章(http://bbs.keinsci.com ...

6-31G**优化以后直接换更高的理论级别算单点就行了
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
本团队长期招收研究生,有意者可私信联系

72

帖子

0

威望

2343

eV
积分
2415

Level 5 (御坂)

16199#
发表于 Post on 2022-8-8 16:09:23 | 只看该作者 Only view this author
wzkchem5 发表于 2022-8-8 14:57
6-31G**优化以后直接换更高的理论级别算单点就行了

好的,谢谢老师

47

帖子

0

威望

403

eV
积分
450

Level 3 能力者

16200#
发表于 Post on 2022-8-8 17:05:23 | 只看该作者 Only view this author
您好老师,我做计算时出现了以下内容导致link dead,请问这是什么原因?
%mem=4GB
%nprocshared=8
Will use up to    8 processors via shared memory.
%chk=E:\test.chk
----------------------------------------------------------------------
#p opt=calcfc casscf(10,8)/3-21g* scrf=(solvent=chloroform) guess=read
  int=ultrafine scf=maxcyc=200 iop(1/8=3)
----------------------------------------------------------------------
1/8=3,10=4,14=-1,18=20,19=15,38=1/1,3;
2/9=110,12=2,17=6,18=5,40=1/2;
3/5=5,7=1,16=1,25=1,32=1,70=2201,71=2,72=7,75=-5,116=101,140=1/1,2,3;
4/5=1,17=10,18=8/1,5;
5/5=2,7=200,17=200,38=6,53=7/10;
8/6=4,10=90,11=11/1;
11/42=1,45=1/1;
10/6=1,13=1/3;
6/7=2,8=2,9=2,10=2,28=1/1;
7/10=1,18=20,25=1/1,2,3,16;
1/8=3,10=4,14=-1,18=20,19=15/3(2);
2/9=110/2;
99//99;
2/9=110/2;
3/5=5,7=1,16=1,25=1,32=1,70=2205,71=1,72=7,75=-5,116=101/1,2,3;
4/5=5,16=3,17=10,18=8,69=1/1,5;
5/5=2,7=200,23=1,38=5,53=7/10;
7//1,2,3,16;
1/8=3,14=-1,18=20,19=15/3(-5);
2/9=110/2;
6/7=2,8=2,9=2,10=2,19=2,28=1/1;
99/9=1/99;
Leave Link    1 at Mon Aug 08 17:00:36 2022, MaxMem=   268435455 cpu:         0.0

本版积分规则 Credits rule

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

GMT+8, 2026-5-24 06:04 , Processed in 1.015941 second(s), 21 queries , Gzip On.

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