计算化学公社

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

[综合交流] 使用sobEDA和sobEDAw方法做非常准确、快速、方便、普适的能量分解分析

  [复制链接 Copy URL]

6万

帖子

99

威望

5万

eV
积分
120062

管理员

公社社长

61#
 楼主 Author| 发表于 Post on 2025-3-25 21:48:19 | 只看该作者 Only view this author
本文的FAQ里专门增加了一节“(9)SCF不收敛怎么办”,并且今日更新了官方sobEDA教程,还在教程包里新添加了额外脚本,使得sobEDA计算片段的过程中可以方便地读取其它级别已收敛的波函数当初猜。
北京科音自然科学研究中心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

139

帖子

0

威望

826

eV
积分
965

Level 4 (黑子)

62#
发表于 Post on 2025-3-26 04:31:36 | 只看该作者 Only view this author
用sobEDA又发了一篇ACS Catalysis,sobEDA的高效率使得我可以批量计算IRC曲线上的数十个结构中的相互作用能分解,观察各项在反应历程中的变化趋势,这在以前只有SAPT可用的时候是完全不可想象的。

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
sobereva + 2 赞!

查看全部评分 View all ratings

8

帖子

0

威望

210

eV
积分
218

Level 3 能力者

63#
发表于 Post on 2025-4-27 09:42:31 | 只看该作者 Only view this author
sob老师,请问此方法是否可以用于一个分子上两个片段间的相互作用?

6万

帖子

99

威望

5万

eV
积分
120062

管理员

公社社长

64#
 楼主 Author| 发表于 Post on 2025-4-27 16:43:22 | 只看该作者 Only view this author
ffcc 发表于 2025-4-27 09:42
sob老师,请问此方法是否可以用于一个分子上两个片段间的相互作用?

得恰当修改分子,饱和边缘、删除有严重位阻的原子,当成双分子对待
北京科音自然科学研究中心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

212

帖子

0

威望

887

eV
积分
1099

Level 4 (黑子)

65#
发表于 Post on 2025-5-6 12:39:41 | 只看该作者 Only view this author
老师 我看了bash脚本中的部分
是专门针对gaussian的输出写的
原理上是不是能对orca的输出作类似的处理 来做sobedaw能量分解?
orca原生的LED方法要求dlpno CCSD(T)方法 对大体系不友好
学生有g16版权 能用g16做计算,出于原理好奇问这个问题

6万

帖子

99

威望

5万

eV
积分
120062

管理员

公社社长

66#
 楼主 Author| 发表于 Post on 2025-5-6 23:51:59 | 只看该作者 Only view this author
18217265596 发表于 2025-5-6 12:39
老师 我看了bash脚本中的部分
是专门针对gaussian的输出写的
原理上是不是能对orca的输出作类似的处理 来 ...

ORCA目前缺乏sobEDA需要的Gaussian才具有的一些设计(至少是对普通用户不可见)

LED给出的能量项的物理意义在我来看不是很清楚
北京科音自然科学研究中心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

威望

5万

eV
积分
120062

管理员

公社社长

67#
 楼主 Author| 发表于 Post on 2025-6-1 20:06:32 | 只看该作者 Only view this author
给本文FAQ部分增加了一段
(10)为什么sobEDA.sh给出的相互作用能和我手动用整体能量减去各片段能量得到的结果不同?
原理上,只要计算得当、有可比性,sobEDA.sh脚本给出的相互作用能和手动计算的相互作用能应当严格相同。若发现不同,可能有以下原因
(1)计算级别、色散校正、几何结构、溶剂效应的考虑、积分格点精度等影响电子能量的因素不精确相同
(2)counterpoise的使用有差异
(3)数据读错了
(4)自己的计算或sobEDA.sh涉及的计算没有收敛到相同的波函数。自己手动计算片段和整体的能量的时候都可以通过Gaussian的stable关键词做波函数稳定性测试以及用stable=opt优化稳定的波函数,从而确保得到的片段和整体的能量对应稳定的波函数,因此手动计算相互作用能完全由可控性。若想让sobEDA.sh计算的总相互作用能能和自己手动计算的相互作用能严格对应上,考虑以下两点:
• 确保sobEDA.sh计算各个片段时收敛到和你手动计算片段时相同的波函数。将sobEDA.sh产生的fragment[序号].out里的能量和自己计算片段时得到的能量相对比,看是否一致可检验。若发现不一致,可以按上面提到的sobEDA官方教程里的第6节,让sobEDA.sh计算片段时读取你做片段计算产生的chk文件里的收敛的波函数当初猜。
• 确保sobEDA.sh做最终整体计算时收敛到了和你手动计算整体时相同的波函数。将sobEDA.sh产生的final.out里的能量和自己计算整体时得到的能量相对比,看是否一致,若发现不一致,说明sobEDA.sh基于片段组合波函数(promol.chk里的波函数)作为初猜没收敛到和你手动计算整体时相同的波函数。对于这种情况,可以把你自己算的整体能量与final.out中SCF过程第一次给出的能量(冻结态能量)相减作为ΔE_orb项,而其它sobEDA的能量项保持不变。
北京科音自然科学研究中心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

52

帖子

0

威望

271

eV
积分
323

Level 3 能力者

68#
发表于 Post on 2025-6-26 16:18:47 | 只看该作者 Only view this author
想请教一下老师,我目前在进行氯仿溶液中150个原子的有机分子笼(不带电)和阴离子体系的弱相互作用(C-H...X氢键)分析。体系的示意图如图1显示。
我目前使用您的sobEDA对我的弱相互作用体系进行能量分解分析,刚开始几次因为我的内存及其他设置的问题始终没有正常算完。这一次正常算完了,但是算出来的结果却显示相互作用能全部等于0.我也找不出问题,我请师兄帮我看了一下,他因为没用过也不太清楚,但是他看了输出文件以后说高斯计算应该没问题,应该是后续数据抓取的时候出问题了,因此我想请您帮我看一下
相关的输入输出文件我已经上传到论坛了,输出文件里面,sobEDA.out,sobEDA.err,屏幕上的输出(1.out),final.out我都上传了,final.chk文件因为很大,有500多MB,因此我就先不上传

图片2025.6.20.png (205.1 KB, 下载次数 Times of downloads: 37)

图片2025.6.20.png

sobEDA.806517.out

2.02 KB, 下载次数 Times of downloads: 2

sobEDA.806517.err

2.28 KB, 下载次数 Times of downloads: 1

1.out

2.71 KB, 下载次数 Times of downloads: 1

final.out

119.76 KB, 下载次数 Times of downloads: 0

fragment.txt

36 Bytes, 下载次数 Times of downloads: 0

system.xyz

9.81 KB, 下载次数 Times of downloads: 0

template.gjf

127 Bytes, 下载次数 Times of downloads: 0

6万

帖子

99

威望

5万

eV
积分
120062

管理员

公社社长

69#
 楼主 Author| 发表于 Post on 2025-6-27 07:40:32 | 只看该作者 Only view this author
peter123 发表于 2025-6-26 16:18
想请教一下老师,我目前在进行氯仿溶液中150个原子的有机分子笼(不带电)和阴离子体系的弱相互作用(C-H.. ...

.err文件里可见
/share/home/lics/chenyiz/sobEDA-cageallylchloCl-2/sobEDA.sh: line 32: dos2unix: command not found
/share/home/lics/chenyiz/sobEDA-cageallylchloCl-2/sobEDA.sh: line 105: bc: command not found
/share/home/lics/chenyiz/sobEDA-cageallylchloCl-2/sobEDA.sh: line 112: bc: command not found
...

明显dos2unix没装、bc命令用不了
认真阅读下文,该说的都说了
使用sobEDA和sobEDAw方法做非常准确、快速、方便、普适的能量分解分析
http://sobereva.com/685http://bbs.keinsci.com/thread-39446-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

8

帖子

0

威望

81

eV
积分
89

Level 2 能力者

70#
发表于 Post on 2025-7-1 16:57:16 | 只看该作者 Only view this author
sob老师你好,请问sobEDA和sobEDAW可以进行铀酰配合物的能量分解计算吗?我看的相关文献都是使用ADF在PBE/TZP/ZORA理论水配上进行能量分解的,但全网关于ADF的使用及讨论很少,并且学校未购买ADF,请问若sobEDA和sobEDAW不支持,可以用ORCA进行替代计算吗?

6万

帖子

99

威望

5万

eV
积分
120062

管理员

公社社长

71#
 楼主 Author| 发表于 Post on 2025-7-3 03:53:35 | 只看该作者 Only view this author
慢醇似梦 发表于 2025-7-1 16:57
sob老师你好,请问sobEDA和sobEDAW可以进行铀酰配合物的能量分解计算吗?我看的相关文献都是使用ADF在PBE/T ...

sobEDA和sobEDAw可以用。如果片段划分截断化学键,只能用前者
北京科音自然科学研究中心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

50

帖子

0

威望

353

eV
积分
403

Level 3 能力者

72#
发表于 Post on 2025-7-14 14:23:01 | 只看该作者 Only view this author
ALMO-EDA中distortion energy+Pauli repulsion认为是位阻效应。请问在sobEDA和sobEDAw中,交换-互斥项ΔE_xrep+变形能可以认为是位阻效应吗?

6万

帖子

99

威望

5万

eV
积分
120062

管理员

公社社长

73#
 楼主 Author| 发表于 Post on 2025-7-15 04:21:46 | 只看该作者 Only view this author
typ_gauss 发表于 2025-7-14 14:23
ALMO-EDA中distortion energy+Pauli repulsion认为是位阻效应。请问在sobEDA和sobEDAw中,交换-互斥项ΔE_x ...

北京科音自然科学研究中心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

278

帖子

0

威望

4434

eV
积分
4712

Level 6 (一方通行)

74#
发表于 Post on 2025-7-28 15:01:13 | 只看该作者 Only view this author
sobereva 发表于 2023-10-10 21:31
等于
这一般叫结合能

sob老师,这个形变能应该叫Distortion Energy吗?或者叫Deformation Energy?

6万

帖子

99

威望

5万

eV
积分
120062

管理员

公社社长

75#
 楼主 Author| 发表于 Post on 2025-7-29 01:13:56 | 只看该作者 Only view this author
dreamyeye 发表于 2025-7-28 15:01
sob老师,这个形变能应该叫Distortion Energy吗?或者叫Deformation Energy?

叫哪个都可以
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2025-8-12 21:30 , Processed in 0.197949 second(s), 23 queries , Gzip On.

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