24小时热门版块排行榜    

查看: 2029  |  回复: 10

mazuju028

木虫 (正式写手)

[交流] 【求助】请教分数原子该如何写poscar? 已有7人参与

在一个位置上,A原子和B原子各占0.5,不知道这种情形该如何写POSCAR中的Direct ?
比如说NaCl构型中,Na原子一半被Li原子取代这种情形。
谢谢。
回复此楼

» 猜你喜欢

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

锐气藏于胸,和气浮于脸,才气见于事,义气施与人
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guitar2033

至尊木虫 (职业作家)

★ ★
mazuju028(金币+1):谢谢交流,那对于这种情形一般会有什么处理方法呢?随机将部分A原子用B原子代替吗? 2010-11-03 15:51:31
cenwanglai(金币+1):谢谢回复! 2010-11-05 09:45:03
vasp中无法做到这一点
2楼2010-11-03 15:45:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

江山2010

木虫 (小有名气)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
mazuju028(金币+1):谢谢交流 2010-11-04 15:54:20
cenwanglai(金币+1):谢谢回复! 2010-11-05 09:46:03
引用回帖:
Originally posted by mazuju028 at 2010-11-03 15:40:25:
在一个位置上,A原子和B原子各占0.5,不知道这种情形该如何写POSCAR中的Direct ?
比如说NaCl构型中,Na原子一半被Li原子取代这种情形。
谢谢。

我不太明白你的意思,是不是NaCl晶体中一半数量的Na原子被Li原子取代?
那样的话,你可以用Materials-studio  作图,直接将数据导出即可,粘贴到POSCAR中
3楼2010-11-03 17:08:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mazuju028

木虫 (正式写手)

引用回帖:
Originally posted by 江山2010 at 2010-11-03 17:08:23:

我不太明白你的意思,是不是NaCl晶体中一半数量的Na原子被Li原子取代?
那样的话,你可以用Materials-studio  作图,直接将数据导出即可,粘贴到POSCAR中

是啊,就是一半的Na原子被Li原子取代,但是在MS中可以在一个位置上设置0.5个Na及0.5个Li。导出之后的数据部分如下:
Na   0.0000000000000000   0.0000000000000000   0.0000000000000000 MIXTURE   1  0.500000)
Li   0.0000000000000000   0.0000000000000000   0.0000000000000000 MIXTURE   1  0.500000)

这种情况好像不能直接贴到POSCAR中吧?
锐气藏于胸,和气浮于脸,才气见于事,义气施与人
4楼2010-11-03 17:18:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mazuju028

木虫 (正式写手)

cenwanglai:呵呵,mazu版别急,这是个老问题——老解决不了的问题。 2010-11-05 09:48:30
有没有人知道怎么处理啊,
锐气藏于胸,和气浮于脸,才气见于事,义气施与人
5楼2010-11-04 14:56:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
mazuju028(金币+2):谢谢提示 2010-11-04 15:54:07
cenwanglai(金币+1):谢谢回复! 2010-11-05 09:46:21
cenwanglai:有空写写你们搭建集群的事情吧。现场文字直播都可以。呵呵。 2010-11-05 09:47:03
bingmou:搭建集群?我喜欢! 2010-11-05 10:15:10
VASP不能处理这种百分比占位结构,多建几个超胞,按照这个比例选取几个不同的位置可以近似模拟这种情形。
心之所向无惧无悔
6楼2010-11-04 15:49:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主


小木虫(金币+0.5):给个红包,谢谢回帖交流
sunyang的方法有人提过,我没有尝试过,结果怎么样不清楚。希望有人提供实战经验或者参考文献。
7楼2010-11-05 09:52:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
bingmou(金币+1):感谢交流 2010-11-05 10:15:27
mazuju028(金币+2):嗯,谢谢,我就这样试着做一下看了。 2010-11-05 14:39:53
我觉得比较好的解决办法是给分数占位的位置单独定义一种元素种类,然后单独生成一个PAW赝势文件。

但是这个我不会。

分数占位的问题,在POSCAR中是不好非常直接体现的。sunyang的提议也有别人说过,即用超胞来模拟分数占位情况。比如用2unit的超胞,两个对应位置上一个是A原子,另外一个是B原子。这样就构造出AB各0.5占位的结构了。

没有实战经验。期待做个这个计算的分享经验或者提供文献。

[ Last edited by cenwanglai on 2010-11-5 at 12:57 ]
8楼2010-11-05 09:59:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+1):谢谢交流 2010-11-05 12:06:31
其实LS说的就是CASTEP干的事情,CASTEP的百分比就是这样处理的。

从某种意义上来说,这种无序的系统已经不能用能带来描述了。

我觉得最好的方法是通过第一性原理计算得到一个emperial potential,然后进行大规模的分子模拟
9楼2010-11-05 10:17:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qiqi2926

铜虫 (小有名气)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
mazuju028(金币+1):谢谢交流,能否具体解释一下虚晶近似是怎么做的? 2010-11-09 18:46:11
这个应该是原子占据数不为1。。。也曾遇到问题,但至今未找到解决方法,可以试一试虚晶近似。。。
10楼2010-11-09 17:33:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mazuju028 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 C79jjtjAKjEn 2026-09-15 4/200 2026-09-18 06:44 by manzzl
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 C79jjtjAKjEn 2026-09-15 3/150 2026-09-18 03:40 by cITaGg3p5edr
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 mibUvS8DDCwf 2026-09-15 3/150 2026-09-18 03:04 by cITaGg3p5edr
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 mibUvS8DDCwf 2026-09-15 3/150 2026-09-18 02:54 by cITaGg3p5edr
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-18 02:30 by cITaGg3p5edr
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 s3fFTmArrBt6 2026-09-13 6/300 2026-09-17 23:28 by cuPZTDXS3VOd
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +7 LwdutQ8HoqWP 2026-09-13 7/350 2026-09-17 22:53 by cuPZTDXS3VOd
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-17 22:39 by cuPZTDXS3VOd
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 0pYnUiPDfdnk 2026-09-14 4/200 2026-09-17 14:15 by ObVzqQQh1rrE
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-17 12:39 by ObVzqQQh1rrE
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +8 3n8v2C8RimXI 2026-09-13 8/400 2026-09-17 11:27 by UAIjTRMuUFoj
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 mibUvS8DDCwf 2026-09-15 4/200 2026-09-17 06:26 by k06BKNblrGRH
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-17 04:28 by k06BKNblrGRH
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-16 22:26 by Tql5LQhh5rLK
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 23jxep3nCNZb 2026-09-14 4/200 2026-09-16 21:03 by MHf40yA0Uhyb
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-16 02:14 by DgNGHc3h5tPl
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:46 by C79jjtjAKjEn
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见