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

计算化学公社

 找回密码
 现在注册!
查看: 122|回复: 8

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

[复制链接]

133

帖子

8

威望

1075

eV
积分
1368

Level 4 (黑子)

发表于 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数据库随便用。


评分

参与人数 1威望 +1 收起 理由
sobereva + 1

查看全部评分

1932

帖子

8

威望

2091

eV
积分
4183

Level 6 (一方通行)

首席卖萌官

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

8443

帖子

16

威望

1万

eV
积分
19271

管理员

公社社长

发表于 2017-3-18 23:27:11 | 显示全部楼层
我抽空研究一下,若有可能会尽量加到multiwfn里
思想家公社的门口Blog:http://sobereva.com
北京科音自然科学研究中心:http://www.keinsci.com
Multiwfn主页:http://Multiwfn.codeplex.com
计算化学公社:http://bbs.keinsci.com
思想家公社QQ群1号:18616395(2000人已满),思想家公社QQ群2号:466017436(讨论计算化学为主,加入时必须注明研究方向)

Money and papers are rubbish, get a real life!

8443

帖子

16

威望

1万

eV
积分
19271

管理员

公社社长

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

还是应当用小核
思想家公社的门口Blog:http://sobereva.com
北京科音自然科学研究中心:http://www.keinsci.com
Multiwfn主页:http://Multiwfn.codeplex.com
计算化学公社:http://bbs.keinsci.com
思想家公社QQ群1号:18616395(2000人已满),思想家公社QQ群2号:466017436(讨论计算化学为主,加入时必须注明研究方向)

Money and papers are rubbish, get a real life!

133

帖子

8

威望

1075

eV
积分
1368

Level 4 (黑子)

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

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

1932

帖子

8

威望

2091

eV
积分
4183

Level 6 (一方通行)

首席卖萌官

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

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

1932

帖子

8

威望

2091

eV
积分
4183

Level 6 (一方通行)

首席卖萌官

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

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

8443

帖子

16

威望

1万

eV
积分
19271

管理员

公社社长

发表于 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://Multiwfn.codeplex.com
计算化学公社:http://bbs.keinsci.com
思想家公社QQ群1号:18616395(2000人已满),思想家公社QQ群2号:466017436(讨论计算化学为主,加入时必须注明研究方向)

Money and papers are rubbish, get a real life!

1932

帖子

8

威望

2091

eV
积分
4183

Level 6 (一方通行)

首席卖萌官

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

明白了大博士,赝势只是把内层电子的分布做了平均化处理,这样基于电子分布而得到的物理量是有效的,而对于要把每一个电子的行为,至少对价层和亚价层电子进行描述的物理量,就不适合了。
copy that!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

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

GMT+8, 2017-3-28 12:19 , Processed in 0.089169 second(s), 25 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.

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