请选择 进入手机版 | 继续访问电脑版
第13届北京科音初级量子化学培训班将于10月5~8日于北京举办,请点击此链接查看详情。这是新人一次性正确、完整学习量子化学计算的最好、最快机会,能少走无数弯路,欢迎参加并相互转告!(已报满)

计算化学公社

 找回密码
 现在注册!
查看: 2164|回复: 11

[波函数分析交流] 最理想的用于波函数分析的波函数文件格式mwfn已发表!

[复制链接]

2万

帖子

25

威望

3万

eV
积分
58518

管理员

公社社长

发表于 2020-2-22 00:02:12 | 显示全部楼层 |阅读模式
之前已经有很多记录波函数的文件格式,比如wfn、wfx、fch、molden等等,然而这些文件都有这样或那样的缺点,给波函数分析程序的开发带来极为严重的不便,甚至有时导致结果错误,作为波函数分析程序Multiwfn的开发者我对此有极为深刻的感受。这个问题长久以来一直没有解决,程序开发者们总是将就着用现有的不理想的格式,带来很多局限性、混乱、不兼容问题。

为了彻底解决缺乏理想的波函数记录格式给波函数信息的储存与交换,以及给波函数分析程序的开发带来的严重不便,我提出了一种新的格式叫mwfn,全称Multiwfn wavefunction file,此格式已经被http://sobereva.com/multiwfn上的最新版Multiwfn所支持。此格式解决了现有各种格式存在的各种缺点,并且格式紧凑、易于读写、定义严格、高度可扩展。此格式的介绍文章已经发表在ChemRxiv上,欢迎查看:https://doi.org/10.26434/chemrxiv.11872524.v1。文中包含格式提出的思想、详细定义、目前常见格式的简介和对比,还专门说明如何测试文件内容的合理性。

非常欢迎广大量子化学程序开发者支持此格式,只要支持了这个格式,并且能通过文中说的检测方法,波函数就可以被Multiwfn进行分析。我相信此格式终究会彻底代替wfn、wfx、molden等格式,成为波函数分析与量子化学领域储存与交换波函数信息的新的公认的标准。

评分

参与人数 22eV +109 收起 理由
Aminus + 5 我很赞同,大统一
void + 5 牛!
yjmaxpayne + 5 牛!
wuy069 + 5 赞!牛🐂
alonewolfyang + 5 精品内容
hewei8912 + 5 赞!
小红豆 + 5 好物!
bnulk + 5 制定标准!
ggdh + 5 金标准!
dreamyeye + 5 GJ!
fallleave + 5 不明觉厉
朙天儿 + 5 mwfn终究会成为波函数分析与量子化学领域储.
asdf + 5 牛!
小小范 + 5 GJ!
阿锋001 + 5 牛!
dalaomao2008 + 5 赞!
终_焉 + 5 牛!
zsu007 + 5 牛!
猴子好坏 + 5 GJ!
pika02 + 4 撒花

查看全部评分

北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训,是提升计算化学研究水平的最佳选择。欢迎加入“北京科音”公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》(http://bbs.keinsci.com/thread-5098-1-1.html)。
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群:思想家公社QQ群1号:18616395,2号:466017436。合计6000人。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一定会被拒绝加入。
思想家公社的门口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!

236

帖子

1

威望

4008

eV
积分
4264

Level 6 (一方通行)

发表于 2020-2-22 07:43:32 | 显示全部楼层
社长威武!

168

帖子

0

威望

1083

eV
积分
1251

Level 4 (黑子)

发表于 2020-2-22 10:44:14 | 显示全部楼层
社长威武
上海交通大学计算化学与分子生物信息实验室
Shanghai JiaoTong University
Computational Chemistry and Molecular Bioinformatics Laboratory

168

帖子

0

威望

1083

eV
积分
1251

Level 4 (黑子)

发表于 2020-2-22 11:07:47 | 显示全部楼层
已经拜读,mwfn格式确实最全面

Overview and comparison of formats for recording wavefunction

Overview and comparison of formats for recording wavefunction
上海交通大学计算化学与分子生物信息实验室
Shanghai JiaoTong University
Computational Chemistry and Molecular Bioinformatics Laboratory

59

帖子

0

威望

1408

eV
积分
1467

Level 4 (黑子)

发表于 2020-2-22 12:21:29 | 显示全部楼层
.存与交换波函数信息的新的公认的标准。
—— Sob神


723

帖子

0

威望

1731

eV
积分
2454

Level 5 (御坂)

发表于 2020-2-22 13:24:54 来自手机 | 显示全部楼层
lijiayisjtu 发表于 2020-2-22 11:07
已经拜读,mwfn格式确实最全面

单看这个表,我个人感觉未必是包含信息越全越好,因为这些信息之间并非完全独立,例如,(对于单行列式波函数)知道了轨道系数自然就能得到密度矩阵,Fock矩阵和轨道占据数。所以fchk包含的都是必要信息,可以避免存储冗余信息从而节省文件体积。

2万

帖子

25

威望

3万

eV
积分
58518

管理员

公社社长

 楼主| 发表于 2020-2-22 17:58:56 | 显示全部楼层
granvia 发表于 2020-2-22 13:24
单看这个表,我个人感觉未必是包含信息越全越好,因为这些信息之间并非完全独立,例如,(对于单行列式波 ...

密度矩阵虽然可以通过轨道占据数和轨道系数简单得到,但不是所有程序自己都可以做这件事,诸如NBO程序就是直接去读密度矩阵,而自己没有密度矩阵产生功能。而且,mwfn格式明确说明密度矩阵是optional的,想存就存,不想存就不存,给出程序开发者足够的余地自行决定,而不是像诸如molden格式根本不给出这个可能性。

另外,光知道系数矩阵更没法简单得到Fock矩阵,有几个波函数分析程序同时自带单/双电子积分、交换相关势计算代码?

"fchk包含的都是必要信息"的观点不对。Fock矩阵对于涉及轨道能量的分析非常重要,诸如用于获得LMO、AdNDP轨道的能量,而fchk却没有自带。对于后HF波函数,自然轨道占据数很重要,fchk甚至都没有明确给出能够记录这样信息的专门的字段(除非像Multiwfn这样提供特殊的功能基于密度矩阵去得到自然轨道,但显然给其它波函数程序开发者带来不便)。反倒是fch里有很多字段,诸如Coordinates of each shell之类是真正多余的。而且fch还包含一大堆零七八碎的Gaussian私有的对波函数分析毫无意义的信息,第三方开发者在载入、导出的时候都不知道哪些是必须的,没个统一标准,因此都称不上是一个通用的记录、交换波函数信息的格式。

mwfn格式在定义的时候是结合我10年来开发Multiwfn大量功能的过程中积攒的经验,从各个角度多方面非常全面地考虑程序开发的便利性来定义的。没有多余的信息,而只有最必须、最重要的信息,并且留出可扩展的余地。如果你也长期开发功能各异的波函数分析程序,想必会对mwfn格式的重要性和价值深有体会。
北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训,是提升计算化学研究水平的最佳选择。欢迎加入“北京科音”公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》(http://bbs.keinsci.com/thread-5098-1-1.html)。
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群:思想家公社QQ群1号:18616395,2号:466017436。合计6000人。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一定会被拒绝加入。
思想家公社的门口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!

723

帖子

0

威望

1731

eV
积分
2454

Level 5 (御坂)

发表于 2020-2-22 18:23:16 来自手机 | 显示全部楼层
本帖最后由 granvia 于 2020-2-23 09:58 编辑
sobereva 发表于 2020-2-22 17:58
密度矩阵虽然可以通过轨道占据数和轨道系数简单得到,但不是所有程序自己都可以做这件事,诸如NBO程序就 ...

Fock矩阵可由分子轨道能级构成的对角阵通过系数矩阵逆变换得到,具体就是:F = inv(C')*E*inv(C),其中E是对角阵,对角元为各CMO的能量,C是系数矩阵,C’是其转置共轭

2万

帖子

25

威望

3万

eV
积分
58518

管理员

公社社长

 楼主| 发表于 2020-2-22 22:38:14 | 显示全部楼层
granvia 发表于 2020-2-22 18:23
Fock矩阵可由分子轨道能级构成的对角阵通过系数矩阵逆变换得到,具体就是:F = C*E*C’,其中E是对角阵, ...

以这种方式得到Fock矩阵还需牵扯到正交化变换,对于基函数存在线性依赖的时候还得有特殊考虑,让波函数分析程序做这个太麻烦还容易因考虑不周导致出错。mwfn储存Fock/KS矩阵对应的矩阵名是1-e Hamiltonian,是个广义化的名字,实质上不限于Fock或者KS矩阵,因此此字段还可以有额外用处,诸如后HF计算下有的程序也可以定义特殊形式的有效单电子算符记录到这个字段里。

北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训,是提升计算化学研究水平的最佳选择。欢迎加入“北京科音”公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》(http://bbs.keinsci.com/thread-5098-1-1.html)。
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群:思想家公社QQ群1号:18616395,2号:466017436。合计6000人。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一定会被拒绝加入。
思想家公社的门口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!

723

帖子

0

威望

1731

eV
积分
2454

Level 5 (御坂)

发表于 2020-2-23 10:02:13 | 显示全部楼层
sobereva 发表于 2020-2-22 22:38
以这种方式得到Fock矩阵还需牵扯到正交化变换,对于基函数存在线性依赖的时候还得有特殊考虑,让波函数分 ...

8L式子打错了,已更正。 不用正交化,但需取逆。当基函数存在线性依赖时,确实有麻烦,可以取伪逆或求解线性方程,但可能存在问题。

668

帖子

19

威望

2797

eV
积分
3845

Level 5 (御坂)

发表于 2020-3-6 19:04:24 | 显示全部楼层
能不能弄一个二进制的mwfn,高度压缩,使得其体积更小,方便储存和传输。就好像xtc相对于trr那种?
然后可以用Multwfn解压

2万

帖子

25

威望

3万

eV
积分
58518

管理员

公社社长

 楼主| 发表于 2020-3-7 10:33:15 | 显示全部楼层
ggdh 发表于 2020-3-6 19:04
能不能弄一个二进制的mwfn,高度压缩,使得其体积更小,方便储存和传输。就好像xtc相对于trr那种?
然后可 ...

我觉得必要性不高。一般波函数文件也不大,以mwfn格式来存的话,撑死了也就两三百兆,再大的话现有量化程序也算不动。而且文本文件油水很大,rar、7z一压缩也没多大了(比如molden文件压成rar后就几分之一了)。而且弄成二进制之后也不好人工修改、读取。

xtc比trr小是因为用了有损方式储存数据,但波函数储存要求精度比较高,以牺牲精度减小尺寸的话得不偿失,弄不好还导致分析结果有误导性之类就麻烦了。
北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训,是提升计算化学研究水平的最佳选择。欢迎加入“北京科音”公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》(http://bbs.keinsci.com/thread-5098-1-1.html)。
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群:思想家公社QQ群1号:18616395,2号:466017436。合计6000人。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一定会被拒绝加入。
思想家公社的门口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!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2020-9-20 02:45 , Processed in 0.181503 second(s), 31 queries .

快速回复 返回顶部 返回列表