24小时热门版块排行榜    

CyRhmU.jpeg
查看: 5735  |  回复: 43
【奖励】 本帖被评价23次,作者zhangguangping增加金币 18.2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

[资源] 【分享】SIESTA中Bader analysis

Mulliken分布会因基组的不同,方法的不同结果差异很大,甚至有时候出来错误的结果。
除了Mulliken分布,现在还有其他的一些原子电荷布局分析方法:例如Natural bond orbital analysis(NBO),Lowdin analysis,Bader analysis方法。

现在主要简单介绍一下Bader analysis分析方法在SIESTA中的应用。

首先给出Bader analysis的网站(http://theory.cm.utexas.edu/bader/),从上面可以免费下载到分析的代码。其中咱们就用grid version(Version 0.27c)作分析,这个版本可以读取cube文件。
(当然做高斯计算的,也可以准备好空间电荷分布的cube文件,做Bader分析)

下面就要给出SIESTA中空间电荷分布的cube文件。

在SIESTA计算中如果设置SaveTotalCharge  T的话,就会在输出文件中出现.TOCH文件,这个文件记录的是空间总电荷的分布。下面应用SIESTA自带的Util\Contrib\APostnikov中的rho2xsf的修改版本rho2cube(见附件,或者将Util\Grid下面的grid2cube.f文件头部添加一行:
      if (task .eq. 'toch') then
        fnamein = paste(sysname,'.TOCH')
重新编译,也可以将TOCH转化为cube格式。)将.TOCH转化为cube格式。然后做Bader analysis。

以上摘编自http://voznyy.elinity.com/blog/2 ... alysis-with-siesta/

[ Last edited by zhangguangping on 2011-4-9 at 03:54 ]

下面是单独的rho2cube的源代码,编译生成目标文件rho2cube[ Last edited by zhangguangping on 2011-12-8 at 19:14 ]

现在的siesta-3.1-pl9版本(以及相应更高版本)已经添加了bader分析的功能,可以参照手册来进行bader分析。并且分析中已经对有H存在的体系对H添加了一个artificial core of 1 electron.
具体分析参照手册。

[ Last edited by zhangguangping on 2012-7-16 at 10:35 ]

注:rho2cube可以处理晶格矢量不正交的情况,而grid2cube只能处理晶格矢量正交的情况。

[ Last edited by zhangguangping on 2013-5-26 at 19:56 ][ Last edited by zhangguangping on 2013-5-26 at 20:09 ]
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理相关文档 第一性原理和电化学 siesta问题 计算-siesta
vasp的相关资料 atk&transiesta 我的小木虫 siesta & vasp
vasp+wannier90

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

785524519

金虫 (正式写手)


★★★★★ 五星级,优秀推荐

请问安装rho2cube时候出错 用的编译器是gfortran不是ifort怎么才能编译成功?提示的错误是/usr/bin/gfortran -w -xW -mp -O3  -c -o rho2cube.o rho2cube.f
gfortran: error: unrecognized command line option ?.mp?
make: *** [rho2cube.o] Error 1
41楼2016-08-24 18:42:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhangguangping 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复(可上传附件)
信息提示
请填处理意见