- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
ABINIT-5.8.4软件的Windows版本实战守则
(by H.J.Zhao)
前言 第一性原理计算发表文章的要点
1. 你算的是什么,
2. 你为什么要算这个东西,
3. 你用了怎样的算法,
4. 你算出来了什么,
5. 你怎样解释你算出来的东西,有什么结论。
第一章 ABINIT软件的运行
/user1/11542/archives/2005/363031.shtml
以BaTiO3为例,需建立输入文件*.in,在abinit网站下载Ba、Ti、O的赝势文件,之后执行abinis程序,依次输入:
1. bto.in
2. bto.out
3. btoi
4. btoo
5. bto.temp
6. 把三个赝势文件拖入程序执行框中
运行即可!
Abinit使用篇abinis的输入文件分为三类:定义了所有输入文件名的文件(比如in.files,这个文件就是用来告诉abinit哪些文件是abinit读入参数的文件名,以及输出主要结果的文件的名称),定义了计算的控制参数的文件(比如取名为INP),赝势的文件。下面举例如何做好输入文件计算金刚石结构的Si的状态方程,由此也得到Si的晶格常数和体弹性模量。大致的步骤是在Si的晶格常数的实验值附近取11个数据点,也就是说取11个晶格常数或体积,然后计算在这些晶格常数下的总能。在计算得到总能后,采用状态方程拟合得到状态方程、平衡态时的体积(或晶格常数)和体弹性模量。本例子中采用的是LDA-HGH赝势。赝势的文件名为:14si.4.hgh。in.files的内容为:#####BeginINP??????????????????? #设置关键词的文件名为INPOUT????????????????? #主要的输出文件为OUT,该文件将被写入计算最重要的结果sii???????????????????? sisi14si.4.hgh?????? #赝势的文件名######ENDINP文件的内容为:# Crystalline cubic Si #ndtset? 11????????????????????????? #说明下面将有11组数据acell: 3*9.8112 ? ? ? ? ? ? ? ??? # 晶格常数a=b=c,将从9.8112.... a.u.开始增加acell+ 3*0.09???????????????? ? ? #晶格常数将以0.09 a.u.的间隔进行增加#Ground state calculation?kptopt?? 1??????????????????????? #在 k点网格取样时根据对称性来取样,并由下面的?????????????????????????????????????? # ngkpt和kptrlatt, 或者nshiftk和shiftk来确定k点的数目?iscf?? 5?????????????????????????? #采用CG方法对能量进行优化,用在基态计算中。########################################################################Definition of the unit cell? rprim? ?????? 0.0? 0.5? 0.5????? #下面三行定义了原胞的基矢,本例子中Si是fcc结构?????? 0.5? 0.0? 0.5??? ?????? 0.5? 0.5? 0.0?? #Definition of the atom typesntypat 1??????????????? #定义原胞中原子的类别的数目,本例子中只有1类原子znucl 14?????????????? #定义原胞中原子的核电荷数#Definition of the atomsnatom 2????????????? #定义原胞中原子的总个数,本例子中有2个原子typat? 2*1??????????? #定义每类原子的个数,本例子中第一类原子有2个xred??????????????????? #下面定义了原胞中原子的坐标?? 0.0? 0.0? 0.0?????? ?? 0.25 0.25 0.25???? #Gives the number of band, explicitely (do not take the default)nband? 16??????? #定义了要计算的能带的数目,最好按这样来设置: ???????????????????????? # nband 原胞中总的价电子数目/2 + 10#Exchange-correlation functiona
文档评论(0)