24小时热门版块排行榜    

查看: 6308  |  回复: 26
【奖励】 本帖被评价13次,作者identation增加金币 9.8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

[资源] 【原创】ABINIT的一种晶体结构建模方法

以前,在论坛上看到不少帖子,ABINIT的建模方法都是采用natom+rprim+xred+acell来进行建模的,这包括abinit的说明中的例子也是这样。因此,这需要预先用MS建立晶体结构,然后make P1,再导出坐标,然后才能计算。
这里,介绍一种用natom+natrd+spgroup+brvltt+acell+angdeg+xred进行建模的方法
这种建模方法和MS的利用空间群和等效的晶体学位置进行建模的方式是一致的,个人感觉挺方便的。
具体参数的意义可在www.abinit.org中找到。
以PRL 100,186402的一篇文章为例,La2NiMnO6的晶体学数据如下表所示



我们用Monoclinic型为例子,其晶胞有20个原子,因此建模部分如下
spgroup 14
ntypat 57 28 25 8
natom 20
brvltt -1
natrd 6 #读入的原子数目,就是晶体学数据给出的有坐标的原子数目
typat 1 2 3 4 4 4 #表中给出的Monoclinic型La2NiMnO6晶体的数据分别是
                             # La Ni  Mn  O1 O2 O3
acell 5.467 5.510 7.751 Angstrom     #用埃做单位后面要加上Angstrom
                                                          #也可用Bohr做单位
angdeg 90.000 90.119 90.000
xred


[ Last edited by identation on 2010-4-8 at 23:57 ]
回复此楼

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

Abinit 相关资料 有意思的 Abinit及相关小程序

» 猜你喜欢

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

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

zxling

木虫 (小有名气)


我在建群号为12的模型的时候
当读到ntypat 3 7
时出错,说是我的赝势只有两个 读不到三个赝势
而3 7表示的是 Li 和 N原子呢?
其他的输入都都套用了lz的
14楼2010-08-25 21:55:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 27 个回答

wenjlu

木虫 (正式写手)


★★★ 三星级,支持鼓励

★ ★ ★
wuli8(金币+1):感谢交流 2010-04-10 19:21
ym23(金币+2):欢迎常来交流 2010-09-07 22:24:28
(1)这种方法是采用 crystal builder 方法。 必须保证 spgroup不等于0。也就是说必须自己设定值,abinit默认为0。我看你的input里面没有输入

(2)这种方法的缺点是:对于掺杂(特别是无序),这种输入方式无法建模。 只能用unit cell 或者primitive cell (对于个别有序情况) 全部原子输入方式

[ Last edited by wenjlu on 2010-4-8 at 22:12 ]
5楼2010-04-08 22:10:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

非天之翼蓝羽

金虫 (初入文坛)


★★★ 三星级,支持鼓励

同样的建立结构方式。。
6楼2010-04-09 09:14:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ym23(金币+1):欢迎常来交流 2010-09-07 22:24:37
这个方法有的时候很好使,直接敲很少的原子坐标输入就行了。
但是对No.14群用的时候有问题。
希望能一起讨论
8楼2010-04-24 13:45:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天放榜吗? +12 布布和一二 2026-08-19 13/650 2026-08-19 14:33 by 陈晨晨陈啊
[基金申请] filecode,4个jtjc了 +8 ziyangfang 2026-08-19 9/450 2026-08-19 13:08 by blacksandy
[基金申请] 今天基金会出结果吗?20260819 +8 kkkl_v 2026-08-19 8/400 2026-08-19 12:34 by Vivilian
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +52 医学老男孩 2026-08-13 115/5750 2026-08-19 12:09 by 医学老男孩
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 时间戳变了,能看出什么问题? +10 基诺咪客 2026-08-17 14/700 2026-08-19 10:34 by 基诺咪客
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +8 启萌科技 2026-08-12 30/1500 2026-08-19 10:22 by 启萌科技
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +4 Tide man 2026-08-12 5/250 2026-08-19 07:41 by 201120518022
[基金申请] 重要消息,中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +3 瞬息宇宙 2026-08-15 3/150 2026-08-18 21:40 by xiaohuixin
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
信息提示
请填处理意见