应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 一次型NiAl结构优化，晶格矢量变化非常大，求助~
61/1返回列表
查看: 1621  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

luluste

银虫 (小有名气)

[求助] 一次型NiAl结构优化，晶格矢量变化非常大，求助~

刚开始做VASP，终极目标是做掺杂，先建立了一个B2_NiAl晶胞，双原子的，想一次型结构优化后再建立超胞，继而算掺杂的种种参数。但是第一步就非常不顺利，得到的CONTCAR和POSCAR值差别太大了，师姐简单看了一下说应该没问题，实在没办法进行后续工作，特地发帖求助~~拜托各位了，金币数虽少，算是一点心意。附上输入文件。
INCAR如下：
SYSTEM = NiAl_undoped
ENCUT = 459.9 eV
ISTART = 0
ICHARG = 2
ISMEAR = 1
SIGMA =0.2
NSW = 100
IBRION =2
ISIF = 3
POTIM = 0.2
EDIFF = 1E-6
EDIFFG = -1E-3
PREC = Accurate
POSCAR如下：
Al_modify
2.887
-0.5  0.5 -0.5
0.5 -0.5 -0.5
-0.5 -0.5  0.5
1 1
direct
-0.0000000000000000  -0.0000000000000000  -0.0000000000000000
0.4999999999999999   0.5000000000000000   0.5000000000000000
KPOINTS如下：
Automatic generation
0
Monhkorst-Pack
25  25  25
0.0 0.0 0.0
得到的CONTCAR：
Al_modify                              
   2.88700000000000     
    -0.6336163515992900    0.3877578537921768   -0.6336163515992900
     0.3877578537921768   -0.6336163515992900   -0.6336163515992900
    -0.6336163515992901   -0.6336163515992900    0.3877578537921768
   1   1
Direct
  0.0000000000000000  0.0000000000000000  0.0000000000000000
  0.5000000000000000  0.5000000000000000  0.5000000000000000

  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-06-05 20:27:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
fzx2008: 金币+3, 专家考核, 谢谢指教 2012-06-05 21:15:13
luluste: 金币+2, ★★★很有帮助, 太感谢了~ 2012-06-06 08:22:27
-0.5  0.5 -0.5
0.5 -0.5 -0.5
-0.5 -0.5  0.5
这个结构是BCC的吧？并且是BCC的primitive cell，而BCC的primitive cell中只有000一个原子而已。
如果LZ想建立成晶格常数为2.887埃，在顶点和体心分别有一种元素的超胞，那应该是用：
1 0 0
0 1 0
0 0 1
来建立晶胞。
建议LZ建立模型后都用VESTA看一看有没有问题再计算~
赞一下 回复此楼
思想重于技巧，内涵重于表象
2楼2012-06-05 20:32:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luluste

银虫 (小有名气)
引用回帖:
2楼: Originally posted by souledge at 2012-06-05 20:32:39
-0.5  0.5 -0.5
0.5 -0.5 -0.5
-0.5 -0.5  0.5
这个结构是BCC的吧？并且是BCC的primitive cell，而BCC的primitive cell中只有000一个原子而已。
如果LZ想建立成晶格常数为2.887埃，在顶点和体心分别有一种元素的 ...

是BCC的~~顶点为Ni原子，体心为Al原子。我马上试一下~
赞一下 回复此楼
3楼2012-06-05 20:38:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luluste

银虫 (小有名气)
引用回帖:
2楼: Originally posted by souledge at 2012-06-05 20:32:39
-0.5  0.5 -0.5
0.5 -0.5 -0.5
-0.5 -0.5  0.5
这个结构是BCC的吧？并且是BCC的primitive cell，而BCC的primitive cell中只有000一个原子而已。
如果LZ想建立成晶格常数为2.887埃，在顶点和体心分别有一种元素的 ...

不好意思想再问一下，如果我在此基础上把这个晶胞扩展成3*3*3的54原子超胞，相应的还是设置
1 0 0
0 1 0
1 1 0么？
赞一下 回复此楼
4楼2012-06-05 20:46:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mousekingadv

金虫 (小有名气)
引用回帖:
4楼: Originally posted by luluste at 2012-06-05 20:46:15
不好意思想再问一下，如果我在此基础上把这个晶胞扩展成3*3*3的54原子超胞，相应的还是设置
1 0 0
0 1 0
1 1 0么？...

是的。但在下面要写出54个原子每个原子的位置。
另：正格矢部分可以写成（打个比方）
1.0000000000000
10 0 0
0 10 0
0 0 10
这样优化后的晶格常数看起来比较方便。
晶胞变大后KPOINTS可以相应减小（有资源的请略过）；即使是两个原子的晶胞，25*25*25的撒点也是非常（非常！）多的。
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

luluste
5楼2012-06-05 22:07:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luluste

银虫 (小有名气)
送鲜花一朵
引用回帖:
5楼: Originally posted by mousekingadv at 2012-06-05 22:07:41
是的。但在下面要写出54个原子每个原子的位置。
另：正格矢部分可以写成（打个比方）
1.0000000000000
10 0 0
0 10 0
0 0 10
这样优化后的晶格常数看起来比较方便。
晶胞变大后KPOINTS可以相应减小（有资源 ...

真是非常感谢
回复此楼
6楼2012-06-06 08:23:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 luluste 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 289求调剂 +4 yang婷 2026-03-02 4/200 2026-03-02 14:31 by 向上的胖东
[考研] 290分材料工程085601求调剂 数二英一 +3 llx0610 2026-03-02 3/150 2026-03-02 14:15 by yc258
[考研] 282求调剂 +4 2103240126 2026-03-02 5/250 2026-03-02 13:45 by littlehu66
[考研] 291 求调剂 +3 化工2026届毕业� 2026-03-02 3/150 2026-03-02 12:55 by houyaoxu
[考研] 26考研报考西工大材料308分求调剂 +4 weizhong123 2026-03-01 4/200 2026-03-02 12:46 by 无际的草原
[考研] 292求调剂 +7 yhk_819 2026-02-28 7/350 2026-03-02 12:43 by 无际的草原
[考研] 265分求调剂不调专业和学校有行学上就 +6 礼堂丁真258 2026-02-28 9/450 2026-03-02 12:04 by 52hz~~
[考研] 268求调剂 +4 简单点0 2026-03-02 5/250 2026-03-02 11:54 by ms629
[考研] 295求调剂 +8 19171856320 2026-02-28 8/400 2026-03-02 11:19 by yuchj
[考研] 281求调剂 +5 2026计算机_诚心 2026-03-01 8/400 2026-03-02 11:05 by 汪！？！
[考研] 284求调剂 +10 天下熯 2026-02-28 11/550 2026-03-02 11:03 by 无际的草原
[考研] 0854复试调剂 276 +4 wmm9 2026-03-01 6/300 2026-03-02 09:28 by 热情沙漠
[考研] 材料工程269求调剂 +3 白刺玫 2026-03-02 3/150 2026-03-02 09:25 by 一休哥FU
[考研] 306分材料调剂 +4 chuanzhu川烛 2026-03-01 5/250 2026-03-01 19:48 by 无际的草原
[考研] 化工299分求调剂 一志愿985落榜 +5 嘻嘻(*^ω^*) 2026-03-01 5/250 2026-03-01 19:47 by 无际的草原
[考研] 0856材料求调剂 +11 hyf hyf hyf 2026-02-28 12/600 2026-03-01 18:57 by 18137688336
[考研] 285求调剂 +8 满头大汗的学生 2026-02-28 8/400 2026-03-01 16:47 by caszguilin
[考研] 课题组接收材料类调剂研究生 +3 gaoxiaoniuma 2026-02-28 4/200 2026-03-01 14:30 by jjj三跨
[考研] 317一志愿华南理工电气工程求调剂 +6 Soliloquy_Q 2026-02-28 11/550 2026-03-01 11:14 by 歌liekkas
[考研] 304求调剂 +3 52hz~~ 2026-02-28 5/250 2026-03-01 00:00 by 52hz~~
信息提示
关闭
请填处理意见
关闭