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

计算化学公社

 找回密码
 现在注册!
查看: 352|回复: 2

[其它程序] 在UHF级别估算U-J参数的AIDFT+U

[复制链接]

323

帖子

16

威望

2364

eV
积分
3007

Level 5 (御坂)

发表于 2018-10-28 01:25:29 | 显示全部楼层 |阅读模式
本帖最后由 beefly 于 2018-10-28 07:28 编辑

LDA+U(以及GGA+U)方法可以处理强关联体系,其中的U、J参数一般来自实验值拟合,或者文献给出的类似体系参数。2007年,Mosey和Carter提出在UHF级别估算U、J参数,并公布了程序源代码。AIDFT-U程序(也叫GAMESS-UJ)以及使用手册的最新下载网址是:https://github.com/EACcodes/AIDFT-U/
由于AIDFT-U所基于的GAMESS 2006-Feb(R5)版已经很难找到,我把它移植到GAMESS 2009-May(R3)版。之所以未采用新一点的版本,是因为GAMESS从2010版开始,scflib.src文件中的UHF代码结构完全变了,修改比较麻烦。我做的另一处修改是把AIDFT-U的独立输入文件OrbsAA中的参数直接写到GAMESS输入文件里。所有修改的部分放在
  1. C<<<CMOSEY
  2. ...
  3. C>>>
复制代码
这样的结构里,很容易搜索。

简短的用法

1,到http://bbs.keinsci.com/thread-727-1-2.html 下载GAMESS 2009-May(R3)版源代码,以及下载附件中的AIDFT-U源代码(也就是针对该版本修改后的GAMESS源代码)。
2,解压缩,用修改后的GAMESS源代码覆盖gamess/source/中的文件。
3,编译GAMESS。
4,运行附件中的算例。AIDFT-U的关键词在GAMESS输入文件的namelist $OrbsAA中指定,共有三个,IDFrst,IDLast,PopThr(也就是原程序的OrbsAA文件中的三个参数,见手册第三页),分别表示给定原子的(收缩)d基函数起始编号和终止编号,以及布居分析的阈值。由于Mulliken布居分析的结果对基组比较敏感,IDFrst—IDLast尽量避开弥散函数。PopThr不输入的话默认是0;对于较大的模型体系,适当增加PopThr(如5.0e-3)可以节省计算时间。
5,在输出文件的末尾找到U、J的结果,可以用来做LDA+U计算。

输入文件的详细解释,见GAMESS的使用手册。更多细节(例如,如何构建带有点电荷的团簇模型),参见AIDFT-U程序手册以及其中的两篇参考文献。

已知问题

1,ENERGY COMPONENTS部分输出的各种能量都是错的,不要使用。(我个人不喜欢这种粗暴的改代码风格)
2,AIDFT-U仅支持d区过渡元素。不要用于主族元素和f区过渡元素。




gamess-modified.zip

204.72 KB, 下载次数: 2

AIDFT-U

test.zip

132.54 KB, 下载次数: 3

AIDFT-U test

gamess-uj_documentation.pdf

135.95 KB, 下载次数: 5

AIDFT-U doc

评分

参与人数 1威望 +1 收起 理由
sobereva + 1 精品内容

查看全部评分

154

帖子

0

威望

341

eV
积分
495

Level 3 能力者

发表于 2018-10-28 10:43:54 来自手机 | 显示全部楼层
弱弱的问一下,UHF级别可靠吗?

323

帖子

16

威望

2364

eV
积分
3007

Level 5 (御坂)

 楼主| 发表于 2018-10-30 01:12:04 | 显示全部楼层
granvia 发表于 2018-10-28 10:43
弱弱的问一下,UHF级别可靠吗?

对于u-j已经可以了
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

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

GMT+8, 2018-12-19 20:49 , Processed in 0.108965 second(s), 28 queries .

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