计算化学公社

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

[综合交流] 简单量子化学问题答疑专帖

  [复制链接 Copy URL]

38

帖子

0

威望

267

eV
积分
305

Level 3 能力者

1711#
发表于 Post on 2023-6-13 19:49:43 | 只看该作者 Only view this author
sobereva 发表于 2023-6-13 16:50
0 1 0 1 0 2 0 2设得不对
第二个和第三个片段的alpha都比beta电子数多一个,显然整体的自旋多重度不可能 ...

sob老师,我想问一下片段初猜首先是构建片段组合波函数,然后再生成自旋极化单重态,我是不是可以直接将第二步替换成柔性扫描或直接优化过渡态?

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

1712#
 楼主 Author| 发表于 Post on 2023-6-13 22:31:05 | 只看该作者 Only view this author
hehe321123 发表于 2023-6-13 19:49
sob老师,我想问一下片段初猜首先是构建片段组合波函数,然后再生成自旋极化单重态,我是不是可以直接将 ...

柔性扫描、优化过渡态等各种任务都可以用guess=read从构造片段组合波函数任务产生的chk文件里读取波函数当初猜
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

109

帖子

0

威望

1730

eV
积分
1839

Level 5 (御坂)

1713#
发表于 Post on 2023-6-13 23:30:23 | 只看该作者 Only view this author
老师好,我在寻找一个反应的过渡态时,先通过柔性扫描确定了初猜结构,但在构建过渡态的计算输入文件时,我忘了将上一步中扫描的变量删掉(即输入文件最后的B 10 111 S 4 -0.100000),这样会对过渡态的计算有影响吗(比如说优化过程中键长会被固定)?我从输出的结果来看,似乎是没有影响的(扫描的键长在优化过程中是变化的)。
这个是我的过渡态计算输入文件
  1. %chk=aaa.chk
  2. #p opt=(calcfc,ts,noeigen) freq ub3lyp def2svp em=gd3bj

  3. Title Card Required

  4. 0 2
  5. 坐标

  6. B 10 111 S 4 -0.100000

复制代码

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

1714#
 楼主 Author| 发表于 Post on 2023-6-14 09:55:34 | 只看该作者 Only view this author
MercuryLamp 发表于 2023-6-13 23:30
老师好,我在寻找一个反应的过渡态时,先通过柔性扫描确定了初猜结构,但在构建过渡态的计算输入文件时,我 ...

没影响
当前没有任何关键词去读那部分信息
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

38

帖子

0

威望

267

eV
积分
305

Level 3 能力者

1715#
发表于 Post on 2023-6-14 14:30:06 | 只看该作者 Only view this author
各位老师,我想问一下怎么用ef方法追踪特定的振动模式?以及ef方法的写法以及格式?

38

帖子

0

威望

267

eV
积分
305

Level 3 能力者

1716#
发表于 Post on 2023-6-14 15:30:12 | 只看该作者 Only view this author
请问各位老师,这个得到的是自旋极化单重态吗?

202306141527547128..png (56.09 KB, 下载次数 Times of downloads: 13)

202306141527547128..png

202306141529296398..png (143.04 KB, 下载次数 Times of downloads: 11)

202306141529296398..png

202306141530097252..png (29.76 KB, 下载次数 Times of downloads: 12)

202306141530097252..png

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

1717#
发表于 Post on 2023-6-14 15:54:29 | 只看该作者 Only view this author
hehe321123 发表于 2023-6-14 08:30
请问各位老师,这个得到的是自旋极化单重态吗?

是。由自旋密度非零,以及alpha、beta电子数相同这两点就足以判断了
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

1718#
发表于 Post on 2023-6-14 15:54:59 | 只看该作者 Only view this author
hehe321123 发表于 2023-6-14 07:30
各位老师,我想问一下怎么用ef方法追踪特定的振动模式?以及ef方法的写法以及格式?

必须告诉我们你用的是什么计算软件,我们才可能回答。支持EF的程序不唯一
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

38

帖子

0

威望

267

eV
积分
305

Level 3 能力者

1719#
发表于 Post on 2023-6-14 15:55:53 | 只看该作者 Only view this author
wzkchem5 发表于 2023-6-14 15:54
必须告诉我们你用的是什么计算软件,我们才可能回答。支持EF的程序不唯一

高斯

2

帖子

0

威望

45

eV
积分
47

Level 2 能力者

1720#
发表于 Post on 2023-6-14 16:13:16 | 只看该作者 Only view this author
【求助】请问各位老师尿素过氧化氢配合物(UHP)中O-O间的键能怎么算?

DFT RB3LYP/6-31++G(d,p) opt+freq 找到UHP的能量,然后断开O-O键生成的还是一个片段,计算能量,发现最后键能为0。这种断开键之后还是一个片段的,应该如何计算键能?

202306141612439618..png (50.8 KB, 下载次数 Times of downloads: 13)

202306141612439618..png

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

1721#
发表于 Post on 2023-6-14 16:27:55 | 只看该作者 Only view this author
LWS 发表于 2023-6-14 09:13
【求助】请问各位老师尿素过氧化氢配合物(UHP)中O-O间的键能怎么算?

DFT RB3LYP/6-31++G(d,p) opt+f ...

参见http://sobereva.com/522。这篇讲的是氢键,但是举一反三不难知道类似情况下的共价键键能怎么算。
另一种做法是拉长O-O键以后算三重态而非自旋破缺单重态,这样可以避免两个OH自由基复合
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

2

帖子

0

威望

45

eV
积分
47

Level 2 能力者

1722#
发表于 Post on 2023-6-14 16:30:44 | 只看该作者 Only view this author
wzkchem5 发表于 2023-6-14 16:27
参见http://sobereva.com/522。这篇讲的是氢键,但是举一反三不难知道类似情况下的共价键键能怎么算。
...

谢谢老师,我马上去试试

147

帖子

2

威望

1226

eV
积分
1413

Level 4 (黑子)

计算化学公社外门弟子

1723#
发表于 Post on 2023-6-15 00:23:07 | 只看该作者 Only view this author
卢老师和各位老师,我想问一下,在初级班计算激发态的那一部分,有些有机体系是先用 B3LYP、wB97XD 以及 M06-2X 其中之一的泛函来优化基态后,再用 PBE0 泛函计算 TDDFT。但是我还看见有些有机体系是基态和激发态的计算都使用了 PBE0,学生不太理解为什么优化基态和计算激发态都使用 PBE0,或者说这和先用 B3LYP 等泛函优化基态后,再用 PBE0 计算激发态有什么区别呢。
KimariYB, Postgraduate, XiaMen University
School of Electronic Science and Engineering (National Model Microelectronics College)
Research interests: theoretical and computational chemistry and machine learning

109

帖子

0

威望

1730

eV
积分
1839

Level 5 (御坂)

1724#
发表于 Post on 2023-6-15 01:46:57 | 只看该作者 Only view this author
sobereva 发表于 2023-6-14 09:55
没影响
当前没有任何关键词去读那部分信息

谢谢老师,还有一个问题想请教一下您。我在结构优化和振动分析时使用的是B3LYP/def2-SVP,现在想用Shermo计算分子的热力学数据,但找不到该计算级别的频率校正因子,我看到您在博文《谈谈谐振频率校正因子》(http://sobereva.com/221)中提到“各类校正因子对于所用基组都不敏感,可以直接用此方法在相近的基组下的校正因子”,那我用B3LYP/6-31G**的频率校正因子来代替进行计算是否可行?

4

帖子

0

威望

29

eV
积分
33

Level 2 能力者

1725#
发表于 Post on 2023-6-15 10:45:46 | 只看该作者 Only view this author
老师好,我有一个关于VASP的问题。在ENCUT的使用方面,无论是大师兄的教程,还是VASP官方使用指南都是写了用ENMAX*1.3作为ENCUT的值。请问为什么大师兄使用的ENCUT值更高呢
大师兄教程中:Fe的ENCUT 默认值是:ENMAX = 267.882 eV, 前面单点计算设置的 450 eV。结构优化设置成600 eV

本版积分规则 Credits rule

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

GMT+8, 2024-11-24 09:08 , Processed in 0.306750 second(s), 30 queries , Gzip On.

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