计算化学公社

 找回密码 Forget password
 注册 Register

使用Multiwfn通过ETS-NOCV方法深入分析片段间的轨道相互作用

查看数: 38237 | 评论数: 64 | 收藏 Add to favorites 22
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2021-7-31 13:11

正文摘要:

使用Multiwfn通过ETS-NOCV方法深入分析片段间的轨道相互作用In-depth analysis of orbital interactions between fragments using the ETS-NOCV method in Multiwfn 文/Sobereva@北京科音  2021-Jul-31 ...

回复 Reply

sobereva 发表于 Post on 2026-1-8 17:36:52
Dave666 发表于 2026-1-8 17:28
老师您好,请问使用Multiwfn+Gaussian进行ETS-NOCV分析的话,除了Multiwfn的两篇必引文献和Gaussian的引用 ...

Multiwfn手册里ETS-NOCV部分给了原文,引用之
Dave666 发表于 Post on 2026-1-8 17:28:02
老师您好,请问使用Multiwfn+Gaussian进行ETS-NOCV分析的话,除了Multiwfn的两篇必引文献和Gaussian的引用外,对于ETS-NOCV,应该引用什么文献呢?
smiling123 发表于 Post on 2025-11-28 08:39:46
谢谢sob老师的回复,多谢!
抱歉,忘记仔细研读Multiwfn手册2.9.2.2节中周期功能介绍的了。
期待以后可以。
sobereva 发表于 Post on 2025-11-28 06:27:53
smiling123 发表于 2025-11-27 20:52
sob老师您好,请问可以以ETS-NOCV方法计算周期性边界条件下两片段之间的轨道相互作用能吗?审稿人要求计算 ...

目前不支持周期性体系
Multiwfn目前什么功能支持周期性在手册2.9.2.2节写明了
smiling123 发表于 Post on 2025-11-27 20:52:39
本帖最后由 smiling123 于 2025-11-27 21:07 编辑

sob老师您好,请问可以以ETS-NOCV方法计算周期性边界条件下两片段之间的轨道相互作用能吗?审稿人要求计算,头大。我以gaussian试着优化了周期性边界条件下的晶体结构,再计算了两个碎片的单点能,再运行Multiwfn,按流程以ETS-NOCV方法计算pair energy,试试了好几次, 总是报错:
Input path of wavefunction file of fragment   
For example, D:\ll\A-Rise.mwfn
fragment1.fchk

Input path of wavefunction file of fragment   2
For example, D:\ll\A-Rise.mwfn
fragment2.fchk

Loading wavefunction information from fragment1.fchk
Loading wavefunction information from fragment2.fchk

Alpha and beta electrons of all fragments:    181    181

Reloading complex.fchk

Calculating overlap integral matrix between fragment orbitals ...

Intel MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

Lowdin orthogonalization between occupied fragment orbitals ...

Generating frozen state density matrix (P_frz) ...
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image              PC                Routine            Line        Source
Multiwfn           00000000020E85F3  Unknown               Unknown  Unknown
libpthread-2.17.s  00002B7F89426630  Unknown               Unknown  Unknown
Multiwfn           00000000005A0AE0  Unknown               Unknown  Unknown
Multiwfn           0000000000C3E146  Unknown               Unknown  Unknown
Multiwfn           000000000085AD7B  Unknown               Unknown  Unknown
Multiwfn           0000000000431822  Unknown               Unknown  Unknown
libc-2.17.so       00002B7F89655555  __libc_start_main     Unknown  Unknown
Multiwfn           0000000000431729  Unknown               Unknown  Unknown

已经用了最小的基组和最小K点取样了,还是报错,报错内容同上。
但如果用这台机器的Multiwfn程序以ETS-NOCV方法计算非周期性结构两片段的轨道相互作用能,或者测试 Multiwfn程序自带的ETS-NOCV,都成功运行.
不知道哪里出了问题?
f195 发表于 Post on 2025-11-24 10:08:07

谢谢老师!
sobereva 发表于 Post on 2025-11-24 02:39:35
f195 发表于 2025-11-23 19:53
sob老师您好,请问我用ETS-NOCV得出的Sum of pair energies是sobEDA得出的E_orb的两倍以上,是否可以将其 ...

可以
f195 发表于 Post on 2025-11-23 19:53:44
sobereva 发表于 2025-11-11 00:30
1 哪种片段定义方式ΔE_orb越小,暗示哪种片段状态越接近整个体系中其实际状态、在形成实际成键方式时越 ...

sob老师您好,请问我用ETS-NOCV得出的Sum of pair energies是sobEDA得出的E_orb的两倍以上,是否可以将其校正到E_orb来分析?
sobereva 发表于 Post on 2025-11-15 05:05:25
Dave666 发表于 2025-11-13 23:06
老师您好,在文献中,我看到有些做EDA-NOCV分析的,会把NOCV pair分解成两个片段上的轨道来表示,并且标上H ...

“准分子轨道”语义不明
如果是指片段的分子轨道,我有空时倒是可以考虑加入NOCV orbital向片段轨道投影的功能。目前版本无法给出

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
Dave666 + 1 谢谢sob老师,期待新功能

查看全部评分 View all ratings

Dave666 发表于 Post on 2025-11-13 23:06:58
老师您好,在文献中,我看到有些做EDA-NOCV分析的,会把NOCV pair分解成两个片段上的轨道来表示,并且标上HOMO,LUMO等。
比如这篇文章(Nat. Chem. 13, 89–93 (2021),doi:10.1038/s41557-020-00579-w)的图片
您的帖文中,我看到的功能是指出NOCV pair分解成壳层,另外还可以显示片段的准分子轨道。
请问ETS-NOCV是否有和图中类似的功能,能分辨某一NOCV pair是来源于哪个准分子轨道吗?
sobereva 发表于 Post on 2025-11-11 16:39:31
Dave666 发表于 2025-11-11 02:39
谢谢sob老师。
那如果我对同一分子的两个片段,按两种情况(Singlet+Singlet或Triplet+Triplet),只进行E ...

基本可以,但更建议直接用sobEDA给出的E_orb项说事

Dave666 发表于 Post on 2025-11-11 02:39:08
本帖最后由 Dave666 于 2025-11-11 03:33 编辑
sobereva 发表于 2025-11-11 00:30
1 哪种片段定义方式ΔE_orb越小,暗示哪种片段状态越接近整个体系中其实际状态、在形成实际成键方式时越 ...

谢谢sob老师。
那如果我对同一分子的两个片段,按两种情况(Singlet+Singlet或Triplet+Triplet),只进行ETS-NOCV分析后,比较两种情况下输出的Sum of pair energies的绝对值。
我是否可以在文章中通过两者的比较,说明绝对值更小的那种情况所对应的成键方式是更准确的?
另外,对于Sum of pair energies,在文章中我是否可以直接用ΔE_orb表示?(虽然您本文中有指出,Sum of pair energies是总轨道相互作用能的近似值
sobereva 发表于 Post on 2025-11-11 00:30:50
Dave666 发表于 2025-11-10 22:57
sob老师您好,在我所在的方向中,其他课题组为了分析某两原子成键是donor-acceptor或electron-share,会把 ...

1 哪种片段定义方式ΔE_orb越小,暗示哪种片段状态越接近整个体系中其实际状态、在形成实际成键方式时越不需要片段轨道相对于参考态发生变形

2 sobEDA结合ETS-NOCV分析 = EDA-NOCV
仔细看
谈谈分子间结合能的构成以及分解分析思想
http://sobereva.com/733http://bbs.keinsci.com/thread-50890-1-1.html
使用sobEDA和sobEDAw方法做非常准确、快速、方便、普适的能量分解分析
http://sobereva.com/685http://bbs.keinsci.com/thread-39446-1-1.html

有什么审稿人不接受的理由?Multiwfn是用户最多的、最流行的波函数分析程序,知名度极高,到现在都被引用了超过37000次,海量用户都使用了Multiwfn做ETS-NOCV分析
ADF那个功能诞生的早所以现存文献里才见得多



评分 Rate

参与人数
Participants 2
eV +4 收起 理由
Reason
Dave666 + 1 谢谢
student0618 + 3

查看全部评分 View all ratings

Dave666 发表于 Post on 2025-11-10 22:57:55
本帖最后由 Dave666 于 2025-11-10 23:01 编辑

sob老师您好,在我所在的方向中,其他课题组为了分析某两原子成键是donor-acceptorelectron-share,会把这个键切开得到两个片段,然后比较在singlet+singlet和triplet+triplet两种情况下的ΔE_orb,他们认为这个能量越小,越代表这种情况是能更好地描述成键情况的。
比如如果singlet+singlet的ΔE_orb较小,那么则偏向认为是donor-acceptor;反之,如果triplet+triplet的较小,那么认为是electron-share
我搞不太懂为什么ΔE_orb越小越认为是合适的,希望老师能点拨下。
另外,在我所在的方向中,比较多使用的是ADF软件进行的EDA-NOCV分析。想问问这和老师您文中介绍的Multiwfn作的ETS-NOCV有什么区别吗?Multiwfn的ETS-NOCV是否也能用ΔE_orb和同样的原理来判断成键类型,会被审稿人接受吗?

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

GMT+8, 2026-1-24 03:10 , Processed in 0.204603 second(s), 26 queries , Gzip On.

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