“第10届量子化学波函数分析与Multiwfn程序培训班将于5月4-8日于北京举办,这是一次性完整、系统学习波函数分析的各种理论知识和全面掌握强大的Multiwfn波函数分析程序使用的最不可错过的机会!请点击此链接查看详情和报名方式,欢迎参加!

“第18届北京科音分子动力学与GROMACS培训班” 将于5月23-26日于北京举办。这是一次性全面、系统学习分子动力学模拟知识和最流行的分子动力学程序GROMACS的关键机会!报名正在进行中,请点击此链接查看详情,欢迎参加!

计算化学公社

 找回密码 Forget password
 注册 Register
Views: 183|回复 Reply: 11
打印 Print 上一主题 Last thread 下一主题 Next thread

[ORCA] 求助:优化从分子HOMO所在区提取出的模型分子得到的HOMO与原分子中的一样吗?

[复制链接 Copy URL]

15

帖子

0

威望

45

eV
积分
60

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
以下是对问题的更详细描述:
1.原分子优化结构的结果不合理,于是想提取原分子HOMO所在区域的结构作为模型分子,只考察这一个模型分子的HOMO;
2.假设原分子中,其它部分对HOMO的影响并不大,那优化模型分子所得的HOMO能量及形态与原分子(真实的原分子,不是优化结构不合理的"原分子")差异如何?模型分子的HOMO的信息作为原分子的HOMO信息的猜测是否严谨?

15

帖子

0

威望

45

eV
积分
60

Level 2 能力者

12#
 楼主 Author| 发表于 Post on 2026-3-22 13:52:45 | 只看该作者 Only view this author
wzkchem5 发表于 2026-3-22 09:03
你算一下不就知道了
再者从有机常识也能知道,离去性能多强的离去基团的E1cb消除是无垒的,这都得有概念

好的,感谢您的解答!

1万

帖子

0

威望

1万

eV
积分
22374

Level 6 (一方通行)

11#
发表于 Post on 2026-3-22 09:03:43 | 只看该作者 Only view this author
ekidona 发表于 2026-3-21 09:55
老师,关于您最后提到的那个两性离子的稳定性,我查了一些资料,都没提到您所说的情况。请问这个两性离子 ...

你算一下不就知道了
再者从有机常识也能知道,离去性能多强的离去基团的E1cb消除是无垒的,这都得有概念
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
本团队长期招收研究生,有意者可私信联系

15

帖子

0

威望

45

eV
积分
60

Level 2 能力者

10#
 楼主 Author| 发表于 Post on 2026-3-21 09:55:00 | 只看该作者 Only view this author
wzkchem5 发表于 2026-3-19 12:54
为何确保能成键,有实验证据吗?
有的时候,能画出C-P成键的结构式,不代表这样的结构真的存在平衡结构 ...

老师,关于您最后提到的那个两性离子的稳定性,我查了一些资料,都没提到您所说的情况。请问这个两性离子的稳定性有什么实验数据支撑吗?

1万

帖子

0

威望

1万

eV
积分
22374

Level 6 (一方通行)

9#
发表于 Post on 2026-3-19 12:54:30 | 只看该作者 Only view this author
ekidona 发表于 2026-3-18 10:21
按照我的设想,原分子(有100多个原子)的其中一个碳原子和其中一个磷原子是成键的,但两次优化(第一次 ...

为何确保能成键,有实验证据吗?
有的时候,能画出C-P成键的结构式,不代表这样的结构真的存在平衡结构。
举例:两性离子Ph3P(+)-CH2-CH2(-)会自发分解成Ph3P和C2H4,不能指望优化出Ph3P(+)-CH2-CH2(-)的结构,尽管Ph3P(+)-CH2-CH2(-)的结构式是满足有机化学的规则的
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
本团队长期招收研究生,有意者可私信联系

15

帖子

0

威望

45

eV
积分
60

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2026-3-18 16:08:49 | 只看该作者 Only view this author
Stardust0831 发表于 2026-3-18 15:00
DFT-D4几乎可以在所有场景中替代DFT-D3(BJ),其中对过渡金属体系的精度提升最为显著。
不会矫枉过正,M0 ...

感谢您的解答!

477

帖子

3

威望

2916

eV
积分
3453

Level 5 (御坂)

7#
发表于 Post on 2026-3-18 15:00:02 | 只看该作者 Only view this author
ekidona 发表于 2026-3-18 14:22
确实有可能我最初建的模型就有问题,我会从这个研究一下原模型的合理性的。但我还是想问一句:D4会不会“ ...

DFT-D4几乎可以在所有场景中替代DFT-D3(BJ),其中对过渡金属体系的精度提升最为显著。
不会矫枉过正,M06-2X是因为已经表现出了中程相关所以再用BJ阻尼会double-counting,B3LYP没有这种问题。

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
wal + 5

查看全部评分 View all ratings

15

帖子

0

威望

45

eV
积分
60

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2026-3-18 14:22:26 | 只看该作者 Only view this author
wal 发表于 2026-3-18 13:16
算算键级,RDG比较侧面
还有一个问题,B3LYP-D4/def2TZVP这个级别已经不算小了,它的结果我觉得天然有一 ...

确实有可能我最初建的模型就有问题,我会从这个研究一下原模型的合理性的。但我还是想问一句:D4会不会“矫枉过正”?(我对量化计算学得还不是很深入)

793

帖子

15

威望

3382

eV
积分
4475

Level 6 (一方通行)

鸩羽

5#
发表于 Post on 2026-3-18 13:16:57 | 只看该作者 Only view this author
ekidona 发表于 2026-3-18 10:21
按照我的设想,原分子(有100多个原子)的其中一个碳原子和其中一个磷原子是成键的,但两次优化(第一次 ...

算算键级,RDG比较侧面
还有一个问题,B3LYP-D4/def2TZVP这个级别已经不算小了,它的结果我觉得天然有一定合理性,这种情况下如何证明你的设想才是对应真实原分子的?
某不知名实验组从苞米地里长出来的计算选手

211

帖子

0

威望

2363

eV
积分
2574

Level 5 (御坂)

4#
发表于 Post on 2026-3-18 10:40:49 | 只看该作者 Only view this author
你要是能确保你HOMO一定在你关注的区域倒也没什么大问题。
HOMO的核心标准是能量最低,而不是在哪个区域,受多大的影响。
如果你的体系中有前沿轨道能量比较接近的情况,难保不会稍微优化就会导致能量顺序发生变化,HOMO位置也发生变化。

15

帖子

0

威望

45

eV
积分
60

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2026-3-18 10:21:14 | 只看该作者 Only view this author
本帖最后由 ekidona 于 2026-3-18 11:21 编辑
wzkchem5 发表于 2026-3-18 09:52
原分子优化的结构具体哪里不合理?

按照我的设想,原分子(有100多个原子)的其中一个碳原子和其中一个磷原子是成键的,但两次优化(第一次没约束,第二次加上了键约束,都是方法B3LYP,D4,基组用def2-TZVP)的结果都是C与P之间的是断开的。后续用Multiwfn做了RDG分析,发现C与P间的作用大约为范德华力(按照手册操作,选用中等质量的模型)。

1万

帖子

0

威望

1万

eV
积分
22374

Level 6 (一方通行)

2#
发表于 Post on 2026-3-18 09:52:30 | 只看该作者 Only view this author
原分子优化的结构具体哪里不合理?
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
本团队长期招收研究生,有意者可私信联系

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

GMT+8, 2026-4-13 22:26 , Processed in 0.322045 second(s), 22 queries , Gzip On.

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