24小时热门版块排行榜    

查看: 3195  |  回复: 8

xiaoxuan1111

铁杆木虫 (正式写手)

[求助] 求问教程里VASP算例中晶体Si的输入文件POSCAR是怎么得到的 已有3人参与

算例中的POSCAR文件是这样的
system Si
5.430
0.5 0.5 0.0
0.0 0.5 0.5
0.5 0.0 0.5
2
cart
0.00 0.00 0.00
0.25 0.25 0.25
请问为什么晶体硅的原胞中只有两个原子呢?我用MS建模得到的原胞中有4个原子啊
我搜到两种获得POSCAR的方法,一个是从MS中castep的.cell文件里复制
这是我用MS中castep得到的.cell中的文件
%BLOCK LATTICE_CART
      -0.000000000000000       2.713550000000001       2.713550000000000
       2.713550000000001       0.000000000000000       2.713550000000000
       2.713550000000000       2.713550000000000      -0.000000000000001
%ENDBLOCK LATTICE_CART

%BLOCK POSITIONS_FRAC
Si   0.1250000000000000   0.1250000000000000   0.1250000000000001
Si   1.1250000000000000   0.1249999999999999  -0.3749999999999998
Si   0.1249999999999999  -0.3750000000000000   1.1250000000000000
Si  -0.3749999999999999   1.1250000000000002   0.1250000000000001
第二中方法是用MS建模后保存成.cif格式,然后导入VESTA软件,再用这个软件保存成POSCAR格式后,得到
Si
1.0
        3.8375000954         0.0000000000         0.0000000000
        1.9187500477         3.3233725696         0.0000000000
        1.9187500477         1.1077908565         3.1333057072
   Si
    4
Cartesian
     0.959375024         0.553895428         0.391663213
     1.918750048         1.107790857         1.958316067
     1.918750048         2.215581713         0.391663213
     2.878125072         0.553895428         0.391663213
这两种方法得到的基矢和原子坐标是不一样的,难道说它们是等价的吗?
回复此楼

» 猜你喜欢

你想要什么,如何得到想要的,现在应该怎么做?
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

未岸亡灵

铁杆木虫 (著名写手)

暂时离开小木虫


【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
xiaoxuan1111: 金币+5, ★★★很有帮助, 非常感谢 2016-03-23 21:11:31
一样的,基矢不用说。原子坐标第一种目测是相对坐标,第二种是笛卡尔坐标。实质是一样的
2楼2016-03-23 18:43:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lq4866

铜虫 (初入文坛)

3楼2016-03-23 20:35:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoxuan1111

铁杆木虫 (正式写手)

引用回帖:
2楼: Originally posted by 未岸亡灵 at 2016-03-23 18:43:44
一样的,基矢不用说。原子坐标第一种目测是相对坐标,第二种是笛卡尔坐标。实质是一样的

非常感谢!
还有问题哈,为什么VASP算例中晶体硅的原胞中只有两个原子呢,而我用MS建模得到的原胞有4个原子,怎么才能得到跟算例中一样的POSCAR文件呢?
你想要什么,如何得到想要的,现在应该怎么做?
4楼2016-03-23 21:14:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

半尾残狐

金虫 (正式写手)

【答案】应助回帖

引用回帖:
4楼: Originally posted by xiaoxuan1111 at 2016-03-23 21:14:58
非常感谢!
还有问题哈,为什么VASP算例中晶体硅的原胞中只有两个原子呢,而我用MS建模得到的原胞有4个原子,怎么才能得到跟算例中一样的POSCAR文件呢?...

用vi 照实例自己写一个文件就可以了

发自小木虫Android客户端
小师妹啊啊啊
5楼2016-03-31 07:21:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzuer

铜虫 (小有名气)

【答案】应助回帖

个人认为你初始在MS中的建模就是不正确的
没有用正确的空间群信息,而是手动输入了4个分数坐标?

对于Si来说,你最简单的方法,是导入MS中自带的结构(~\share\Structures\semiconductors\Si.msi)

多说一句,Si的原胞是2个原子,单胞是8个原子,DFT计算中一般情况用原胞进行计算

» 本帖已获得的红花(最新10朵)

功能材料的DFT计算
6楼2016-03-31 09:26:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoxuan1111

铁杆木虫 (正式写手)

引用回帖:
5楼: Originally posted by 半尾残狐 at 2016-03-31 07:21:37
用vi 照实例自己写一个文件就可以了
...

我想知道教程中Si的POSCAR怎么来的,如果我能重复出来,其它没有实例的复杂晶体的POSCAR我也就有办法得到了
你想要什么,如何得到想要的,现在应该怎么做?
7楼2016-04-01 17:10:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoxuan1111

铁杆木虫 (正式写手)

送红花一朵
引用回帖:
6楼: Originally posted by fzuer at 2016-03-31 09:26:07
个人认为你初始在MS中的建模就是不正确的
没有用正确的空间群信息,而是手动输入了4个分数坐标?

对于Si来说,你最简单的方法,是导入MS中自带的结构(~\share\Structures\semiconductors\Si.msi)

多说一句, ...

谢谢回复!
我空间群倒没错,我是从spring material上查到的Si晶体结构信息
http://materials.springer.com/is ... ic/docs/sd_1622696#

按它给的信息,Si的原子坐标是(1/8,1/8,1/8),用MS建模后原胞确实有4个原子,单胞有16个。
我试着改了一下原子坐标,改成(0,0,0),用MS建模后就正常了,原胞2个,单胞8个,跟MS自带的结构是一样的。

很奇怪,是我建模方法错了,还是这个数据库给的原子坐标错了呢?
你想要什么,如何得到想要的,现在应该怎么做?
8楼2016-04-01 17:20:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoxuan1111

铁杆木虫 (正式写手)

引用回帖:
6楼: Originally posted by fzuer at 2016-03-31 09:26:07
个人认为你初始在MS中的建模就是不正确的
没有用正确的空间群信息,而是手动输入了4个分数坐标?

对于Si来说,你最简单的方法,是导入MS中自带的结构(~\share\Structures\semiconductors\Si.msi)

多说一句, ...

我今天又试了试,用MS建模时,227空间群的option中有两个选项,Origin-1和Origin-2
用origin-1时,Si原子坐标使用(0,0,0)
用origin-2时,Si原子坐标使用(1/8,1/8,1/8),这样都能得到8个原子的单胞和2个原子的原胞,生成的POSCAR文件分别是
Si_2
1.0
        3.8375000954         0.0000000000         0.0000000000
        1.9187500477         3.3233725696         0.0000000000
        1.9187500477         1.1077908565         3.1333057072
   Si
    2
Direct
     0.000000000         0.000000000         0.000000000
     0.250000000         0.250000000         0.250000000

Si_3
1.0
        3.8375000954         0.0000000000         0.0000000000
        1.9187500477         3.3233725696         0.0000000000
        1.9187500477         1.1077908565         3.1333057072
   Si
    2
Direct
     0.125000000         0.125000000         0.125000000
     0.875000000         0.875000000         0.875000000
晶格基矢相图,原子坐标不同,用VASP算了一下,总能和费米能级是一样的,所以它们应该是等价的。
这样新问题又来了,怎么根据晶体结构信息http://materials.springer.com/is ... ic/docs/sd_1622696#,判断该选用origin-1还是origin-2呢?
你想要什么,如何得到想要的,现在应该怎么做?
9楼2016-04-02 10:14:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaoxuan1111 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕找调剂 +4 哈哈哈吼吼吼哈 2026-03-23 4/200 2026-03-24 17:45 by dick_runner
[考研] 0805 316求调剂 +4 大雪深藏 2026-03-18 4/200 2026-03-24 17:34 by dick_runner
[考研] 070300化学求调剂 +9 苑豆豆 2026-03-20 9/450 2026-03-24 17:15 by licg0208
[考研] 289求调剂 +5 硕星赴 2026-03-23 5/250 2026-03-24 17:09 by hello七七
[考研] 一志愿北京化工大学材料与化工 264分各科过A区国家线 +3 哈哈157349 2026-03-21 3/150 2026-03-24 14:11 by zhyzzh
[考研] 305分求调剂(食品工程) +5 Sxy112 2026-03-21 7/350 2026-03-24 12:27 by 544594351
[考研] 086003食品工程求调剂 +4 淼淼111 2026-03-24 4/200 2026-03-24 11:53 by 544594351
[考研] 07化学280分求调剂 +4 722865 2026-03-23 4/200 2026-03-24 00:01 by chixmc
[考研] 一志愿武理材料工程348求调剂 +6  ̄^ ̄゜汗 2026-03-19 9/450 2026-03-23 19:53 by pswait
[考研] 上海电力大学材料防护与新材料重点实验室招收调剂研究生(材料、化学、电化学,环境) +3 我爱学电池 2026-03-23 3/150 2026-03-23 17:16 by AZMK
[考研] 接收2026硕士调剂(学硕+专硕) +4 allen-yin 2026-03-23 6/300 2026-03-23 15:04 by 汪!?!
[考研] 298求调剂 +8 上岸6666@ 2026-03-20 8/400 2026-03-23 11:02 by laoshidan
[考研] 333求调剂 +6 87639 2026-03-21 10/500 2026-03-23 10:41 by Iveryant
[考研] 324求调剂 +6 lucky呀呀呀鸭 2026-03-20 6/300 2026-03-22 16:01 by ColorlessPI
[考研] 085600材料与化工306 +4 z1z2z3879 2026-03-21 4/200 2026-03-21 23:44 by ms629
[考研] 材料与化工(0856)304求B区调剂 +3 邱gl 2026-03-20 7/350 2026-03-21 19:05 by 15709483992
[考研] 332求调剂 +4 ydfyh 2026-03-17 4/200 2026-03-21 02:20 by JourneyLucky
[考研] 华东师范大学-071000生物学-293分-求调剂 +3 研究生何瑶明 2026-03-18 3/150 2026-03-21 01:30 by JourneyLucky
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
信息提示
请填处理意见