请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 449|回复: 13

[其它] 用于原子芯电子密度的EDF数据库

[复制链接]

188

帖子

10

威望

1441

eV
积分
1829

Level 5 (御坂)

发表于 2017-3-18 20:35:48 | 显示全部楼层 |阅读模式
本帖最后由 beefly 于 2017-3-18 23:34 编辑

Gaussian09产生的WFX波函文件加入了EDF字段,用于描述被ECP取代的芯电子的电子密度。但是Gaussian是有版权的,EDF数据无法移植到其他程序中。


最近给molden2aim做了一个EDF数据库。请sobereva考虑是否也加入到multiwfn里
https://github.com/zorkzou/Molden2AIM/blob/master/src/edflib.f90

两个EDF数据库的差异:
1. Gaussian的原子芯密度数据是在标量相对论(DKH2)的DFT/UGBS级别产生的,新的原子芯密度数据是在二分量相对论(X2C+AMFI)HF/Dyall-CV4Z级别产生的。其中HF采用分数平均占据以保持波函的球对称性。
2. 用于拟合的高斯指数,采用同一个指数生成公式。但是Gaussian稍微改了公式里的参数。
3. Gaussian对每个壳层分别拟合,导致源代码冗长,光是拟合数据就有4MB,好处是比较灵活,原则上可以处理任意芯电子数的ECP。而这里只对芯电子的总密度进行拟合,可以看到数据代码只有300多KB,不足之处是只支持常见的芯电子数的ECP,少数罕见的ECP没有拟合(比如,有一种把镧系元素开壳层f电子作为芯电子的ECP,这个库就没有)。4. Gaussian的EDF支持到103号元素,这个库支持到120号元素。
5. Gaussian的EDF代码不能移植。这个EDF数据库随便用。


评分

参与人数 2威望 +1 eV +4 收起 理由
leebo + 4 好物!
sobereva + 1

查看全部评分

2120

帖子

8

威望

2583

eV
积分
4863

Level 6 (一方通行)

首席卖萌官

发表于 2017-3-18 21:53:55 | 显示全部楼层
如果真如此,那么做波函数分析的时候估计也就可以不必一定要在小核赝势下分析了,只需要对内部的赝势进行合理的描述就可以了……我说得没错吧,Zork老师?

9953

帖子

16

威望

1万

eV
积分
22850

管理员

公社社长

发表于 2017-3-18 23:27:11 | 显示全部楼层
我抽空研究一下,若有可能会尽量加到multiwfn里
思想家公社的门口Blog:http://sobereva.com
北京科音自然科学研究中心:http://www.keinsci.com
Multiwfn主页:http://sobereva.com/multiwfn
计算化学公社:http://bbs.keinsci.com
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论内容相同,可加入任意其一但不可都加入,申请信息必须注明研究方向,否则一概不批。

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

9953

帖子

16

威望

1万

eV
积分
22850

管理员

公社社长

发表于 2017-3-18 23:28:27 | 显示全部楼层
我本是个娃娃 发表于 2017-3-18 21:53
如果真如此,那么做波函数分析的时候估计也就可以不必一定要在小核赝势下分析了,只需要对内部的赝势进行合 ...

还是应当用小核
思想家公社的门口Blog:http://sobereva.com
北京科音自然科学研究中心:http://www.keinsci.com
Multiwfn主页:http://sobereva.com/multiwfn
计算化学公社:http://bbs.keinsci.com
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论内容相同,可加入任意其一但不可都加入,申请信息必须注明研究方向,否则一概不批。

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

188

帖子

10

威望

1441

eV
积分
1829

Level 5 (御坂)

 楼主| 发表于 2017-3-18 23:36:38 | 显示全部楼层
我本是个娃娃 发表于 2017-3-18 21:53
如果真如此,那么做波函数分析的时候估计也就可以不必一定要在小核赝势下分析了,只需要对内部的赝势进行合 ...

是。尤其是large-core ecp,用edf可以得到合理结果。不过距离全电子结果还是有一定的差异

2120

帖子

8

威望

2583

eV
积分
4863

Level 6 (一方通行)

首席卖萌官

发表于 2017-3-19 09:27:10 | 显示全部楼层

哦吼吼,那就再做一次高精度单点,在做波函数分析咯

2120

帖子

8

威望

2583

eV
积分
4863

Level 6 (一方通行)

首席卖萌官

发表于 2017-3-19 09:28:08 | 显示全部楼层
beefly 发表于 2017-3-18 23:36
是。尤其是large-core ecp,用edf可以得到合理结果。不过距离全电子结果还是有一定的差异

那是自然,全电子对于每一个电子都有精确函数描述,而不是把一堆内核电子赝化成一个平均势场

9953

帖子

16

威望

1万

eV
积分
22850

管理员

公社社长

发表于 2017-3-20 16:02:47 | 显示全部楼层
我本是个娃娃 发表于 2017-3-19 09:27
哦吼吼,那就再做一次高精度单点,在做波函数分析咯


可以这样。
关键是,EDF只对于纯粹基于密度的量有贡献,比如AIM的电子密度拓扑分析,而诸如ELF这样基于波函数的量,EDF没有效果,大核时的结果肯定显著不如小核。
思想家公社的门口Blog:http://sobereva.com
北京科音自然科学研究中心:http://www.keinsci.com
Multiwfn主页:http://sobereva.com/multiwfn
计算化学公社:http://bbs.keinsci.com
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论内容相同,可加入任意其一但不可都加入,申请信息必须注明研究方向,否则一概不批。

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

2120

帖子

8

威望

2583

eV
积分
4863

Level 6 (一方通行)

首席卖萌官

发表于 2017-3-20 16:40:13 | 显示全部楼层
sobereva 发表于 2017-3-20 16:02
可以这样。
关键是,EDF只对于纯粹基于密度的量有贡献,比如AIM的电子密度拓扑分析,而诸如ELF这样基 ...

明白了大博士,赝势只是把内层电子的分布做了平均化处理,这样基于电子分布而得到的物理量是有效的,而对于要把每一个电子的行为,至少对价层和亚价层电子进行描述的物理量,就不适合了。
copy that!

9953

帖子

16

威望

1万

eV
积分
22850

管理员

公社社长

发表于 2017-5-24 10:49:38 | 显示全部楼层
我在刚更新的Multiwfn 3.4(dev)版中已经把这个库纳入了。默认情况下,载入任意含有波函数信息的文件时,只要程序发现有的原子用了赝势,就会根据原子序号与核电荷的差值得到被赝势代替的内核电子数,并在这个EDF库中读取相应的EDF信息。而对于高斯产生的.wfx文件,默认还是直接读取wfx里的EDF信息。
在Multiwfn手册附录4中做了更详细的说明。
思想家公社的门口Blog:http://sobereva.com
北京科音自然科学研究中心:http://www.keinsci.com
Multiwfn主页:http://sobereva.com/multiwfn
计算化学公社:http://bbs.keinsci.com
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论内容相同,可加入任意其一但不可都加入,申请信息必须注明研究方向,否则一概不批。

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

188

帖子

10

威望

1441

eV
积分
1829

Level 5 (御坂)

 楼主| 发表于 2017-5-24 11:56:03 | 显示全部楼层
本帖最后由 beefly 于 2017-5-24 11:57 编辑
sobereva 发表于 2017-5-24 10:49
我在刚更新的Multiwfn 3.4(dev)版中已经把这个库纳入了。默认情况下,载入任意含有波函数信息的文件时,只 ...

谢谢!

请更新到最新版本的库。最近做了两次升级,一是径向密度ρ(r)的拟合改为径向密度分布函数ξ(r)=4πr2ρ(r)的拟合(因为成键区受原子核附近的拟合精度影响不大,对大r处的误差反而更敏感),二是为了提高拟合精度,用了更柔和的Gaussian生成函数。

举个粒子:

EDF

EDF

188

帖子

10

威望

1441

eV
积分
1829

Level 5 (御坂)

 楼主| 发表于 2017-5-24 12:08:36 | 显示全部楼层
我看到已经是最新版了!

第二个例子
Graph1.png

9953

帖子

16

威望

1万

eV
积分
22850

管理员

公社社长

发表于 2017-5-25 02:15:27 | 显示全部楼层
beefly 发表于 2017-5-24 12:08
我看到已经是最新版了!

第二个例子

不知这个库之后还有什么进一步改进的计划?打算发表么?非常好的工作
思想家公社的门口Blog:http://sobereva.com
北京科音自然科学研究中心:http://www.keinsci.com
Multiwfn主页:http://sobereva.com/multiwfn
计算化学公社:http://bbs.keinsci.com
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论内容相同,可加入任意其一但不可都加入,申请信息必须注明研究方向,否则一概不批。

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

188

帖子

10

威望

1441

eV
积分
1829

Level 5 (御坂)

 楼主| 发表于 2017-5-26 10:35:48 | 显示全部楼层
sobereva 发表于 2017-5-25 02:15
不知这个库之后还有什么进一步改进的计划?打算发表么?非常好的工作

3x!正在准备论文,方法已经定型,估计不会再改了
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

Archiver|手机版|小黑屋|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2017-8-18 18:46 , Processed in 0.098966 second(s), 28 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.

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