24小时热门版块排行榜    

查看: 3637  |  回复: 10
本帖产生 1 个 1ST强帖 ,点击这里进行查看

ljjhb1

银虫 (正式写手)


[交流] abinit能否处理一个位置被两种元素占据的固溶体的能带和态密度

现在遇到一个固溶体结构,iscd的数据显示某一个位置被两种元素占据,不知abinit软件能否处理这种情况下的能带和态密度计算,谢谢
回复此楼

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

Abinit 相关资料

» 猜你喜欢

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

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )
★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
ljjhb1(金币+1): 谢谢 2011-05-07 18:40:16
mazuju028(金币+2, 1ST强帖+1): 谢谢提示。 2011-05-08 22:43:13
一个位置被多重元素占据的固溶体,计算能带和态密度 可以采用 超胞方法,SQS,VCA,CPA 等途径实现

Abinit 可以通过构造 多种原子的混合赝势(VCA方法)来实现计算固溶体的性质

可以参考文献
http://iopscience.iop.org/0965-0393/13/7/010/pdf/msmse5_7_010.pdf
http://people.gl.ciw.edu/cohen/m ... ceedings/ghosez.pdf
http://iopscience.iop.org/0953-8 ... f/cm7_30_306203.pdf

VCA
http://ftp.abinit.org/ws07/Liege_Mikami.PDF
http://forum.abinit.org/viewtopic.php?f=7&t=703#p2086
http://forum.abinit.org/viewtopic.php?f=10&t=1027
2楼2011-05-07 17:52:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖
★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
ljjhb1(金币+1): 谢谢 2011-05-07 18:40:28
mazuju028(金币+3): 谢谢交流 2011-05-08 22:43:52
generating ALCHemical pseudopotentials


相关参数 ntypalch
             npspalch
             algalch
             mixalch  ------- 这个关键字中含有例子

比如Ba(0.25) Sr(0.75) Ti O3
npsp 4
znucl 8 40 56 38
ntypat 3
ntypalch 1
mixalch 0.25 0.75


Note that alchemical mixing cannot be used with PAW.
3楼2011-05-07 18:25:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
mazuju028(金币+2): 谢谢交流 2011-05-08 22:44:13
引用回帖:
Originally posted by fzx2008 at 2011-05-07 18:25:17:
generating ALCHemical pseudopotentials


相关参数 ntypalch
             npspalch
             algalch
             mixalch  ------- 这个关键字中含有例子

比如Ba(0.25) Sr(0.75) Ti O3
npsp 4 ...

对 PAW暂时还不支持

具体参数:
http://www.abinit.org/documentat ... /vargs.html#algalch
4楼2011-05-07 18:28:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljjhb1

银虫 (正式写手)


引用回帖:
Originally posted by sars518 at 2011-05-07 17:52:34:
一个位置被多重元素占据的固溶体,计算能带和态密度 可以采用 超胞方法,SQS,VCA,CPA 等途径实现

Abinit 可以通过构造 多种原子的混合赝势(VCA方法)来实现计算固溶体的性质

可以参考文献
[url]http:// ...

您好,我现在是遇到这样一个问题:icsd给出的是模型是同一位置被两种元素占据,我现在不想用VCA虚晶方法计算,想使用超胞方法计算,也就是想把模型转变为一个位置被一种元素占据,不知这种方法可行否,如何建呢,能否详细说下,有重谢。
比如您可以以ZrTiO4为例,我给出其坐标数据(Zr和Ti位置是重合的,分别占50%)
*data for   ICSD #69619
Coll Code   69619
Rec  Date   1992/08/20
Mod  Date   1999/01/19
Chem Name   Zirconium Titanium Oxide
Structured  (Zr Ti) O4
Sum         O4 Ti1 Zr1
ANX         AX2
Min Name    Srilankite
D(calc)     5.12
Title       Zr Sn0.5 Ti0.5 O4 als Wirtsstruktur fuer keramische Farbkoerper.
            Strukturuntersuchungen an Zr Ti O4 und Zr Sn0.5 Ti0.5 O4
Author(s)   Siggel, A.;Jansen, M.
Reference   Zeitschrift fuer Anorganische und Allgemeine Chemie
            (1990), 582, 93-102
Unit Cell   4.7855(13) 5.4755(12) 5.0277(14) 90. 90. 90.
Vol         131.74
Z           2
Space Group P b c n
SG Number   60
Cryst Sys   orthorhombic
Pearson     oP12
Wyckoff     d c
R Value     .0274
Red Cell    P  4.785 5.027 5.475 89.999 89.999 90 131.741
Trans Red   1.000 0.000 0.000 / 0.000 0.000 1.000 / 0.000 -1.000 0.000
Comments    Compound with mineral name: Srilankite
            The structure has been assigned a PDF number: 35-584
            Calculated density unusual but tolerable.
Atom  #   OX   SITE      x           y           z           SOF      H
Zr   1  +4    4 c   0           0.2025(1)   0.25           0.5        0   
Ti   1  +4    4 c   0           0.2025(1)   0.25           0.5        0   

O    1  -2    8 d   0.2689(6)   0.4024(5)   0.4337(6)      1.         0   
Lbl  Type      U11         U22         U33         U12         U13         U23
Zr1  Zr4+ 0.0285(4)   0.0314(4)   0.0339(4)   0           -.0130(3)   0           
Ti1  Ti4+ 0.0285(4)   0.0314(4)   0.0339(4)   0           -.0130(3)   0           
O1   O2-  0.0307(15)  0.0268(14)  0.0563(19)  -.0020(18)  0.0058(15)  0.0154(13)  
*end for    ICSD #69619



[ Last edited by ljjhb1 on 2011-5-7 at 19:16 ]
5楼2011-05-07 19:06:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljjhb1

银虫 (正式写手)


引用回帖:
Originally posted by fzx2008 at 2011-05-07 18:25:17:
generating ALCHemical pseudopotentials


相关参数 ntypalch
             npspalch
             algalch
             mixalch  ------- 这个关键字中含有例子

比如Ba(0.25) Sr(0.75) Ti O3
npsp 4 ...

谢谢,能否看看我的回帖,帮忙解决超胞建模问题,有重谢
6楼2011-05-07 19:10:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫(金币+0.5):给个红包,谢谢回帖
ljjhb1(金币+4): 感谢 2011-05-08 08:55:47
引用回帖:
Originally posted by ljjhb1 at 2011-05-07 19:06:54:
您好,我现在是遇到这样一个问题:icsd给出的是模型是同一位置被两种元素占据,我现在不想用VCA虚晶方法计算,想使用超胞方法计算,也就是想把模型转变为一个位置被一种元素占据,不知这种方法可行否,如何 ...

你可以考虑构建2*2*2 或者 n*m*l的超胞

按照原始比例称份,分别用Zr和Ti原子 替换(ZrTi)位置的原子,最后Zr:Ti保持为1:1.由于替换方案有很多种,所以,你可以计算下这几种方案的能量,最后用能量低的做进一步计算。
7楼2011-05-07 19:40:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫(金币+0.5):给个红包,谢谢回帖
ljjhb1(金币+4): 特别感谢 2011-05-08 08:56:30
写在前面:超单胞方法,不论何种组合都只是一种近似模拟,都不等价实际意义上的无序排列。


建模的实际操作中,首先4c格位全部用Zr,然后去掉对称性,单胞中应有4个Zr,任取其中两个替换成Ti即可。这样单胞Zr2Ti2O8。

若考虑2x2x2的超单胞,排列组合数目更多!

下面给出单胞中替换示意图,仅供参考

8楼2011-05-07 20:08:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljjhb1

银虫 (正式写手)


引用回帖:
Originally posted by sars518 at 2011-05-07 19:40:08:
你可以考虑构建2*2*2 或者 n*m*l的超胞

按照原始比例称份,分别用Zr和Ti原子 替换(ZrTi)位置的原子,最后Zr:Ti保持为1:1.由于替换方案有很多种,所以,你可以计算下这几种方案的能量,最后用能量低的 ...

您好,感谢你的帮助。看来超胞计算计算量挺大加上每种晶格优化的话。
还有一个疑问,如果使用VCA方法计算,abinit能否接着算光学性质呢(我知道castep不能计算的,所以想建超胞),再次感谢。
9楼2011-05-08 09:08:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljjhb1

银虫 (正式写手)


引用回帖:
Originally posted by fzx2008 at 2011-05-07 20:08:08:
写在前面:超单胞方法,不论何种组合都只是一种近似模拟,都不等价实际意义上的无序排列。


建模的实际操作中,首先4c格位全部用Zr,然后去掉对称性,单胞中应有4个Zr,任取其中两个替换成 ...

非常感谢您的帮助。
知道你对abinit熟悉,不知道如果我使用VCA方法建模,abinit能否计算VCA模型的光学性质,谢谢
10楼2011-05-08 09:10:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫(金币+0.5):给个红包,谢谢回帖
ljjhb1你好,Zr和Ti属于chemically similar atoms,VCA应该可以算。
建议调研一下文献,有不少。
11楼2011-05-08 10:35:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ljjhb1 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 93BebMhtakh前后11位开头都是大写 +3 且听虎啸 2026-08-17 4/200 2026-08-17 22:31 by ccec101
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +36 医学老男孩 2026-08-13 85/4250 2026-08-17 21:23 by lcmdlut
[基金申请] 感觉是下周放榜了 +5 angus9576 2026-08-17 9/450 2026-08-17 20:57 by ssxclkj
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 时间戳变了,能看出什么问题? +4 基诺咪客 2026-08-17 4/200 2026-08-17 16:10 by Vivilian
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +12 Ldrop2023 2026-08-13 15/750 2026-08-17 15:02 by 小豌豆_发芽
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +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
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 重要来源:本周末出结果 +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
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
信息提示
请填处理意见