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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 一次型NiAl结构优化，晶格矢量变化非常大，求助~
51/1返回列表
查看: 1625  |  回复: 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的回帖

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的回帖
查看全部 6 个回答

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埃，在顶点和体心分别有一种元素的 ...

不好意思想再问一下，如果我在此基础上把这个晶胞扩展成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的回帖
查看全部 6 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程专硕283求调剂 5+5 ，！？ 2026-03-02 7/350 2026-03-02 20:27 by hypershenger
[考研] 材料学硕318求调剂 +10 February_Feb 2026-03-01 10/500 2026-03-02 20:18 by hypershenger
[考研] 材料270求调剂 6+4 Eiiiio 2026-03-01 7/350 2026-03-02 20:17 by hypershenger
[考研] 085600材料工程一志愿中科大总分312求调剂 +9 吃宵夜1 2026-02-28 11/550 2026-03-02 20:14 by hypershenger
[考研] 材料复试调剂 +5 学材料的点 2026-03-01 6/300 2026-03-02 20:01 by hypershenger
[考研] 275求调剂 +7 明远求学 2026-03-01 7/350 2026-03-02 19:22 by zhukairuo
[考研] 化工京区271求调剂 +6 11ing 2026-03-02 6/300 2026-03-02 18:52 by caszguilin
[考研] 一志愿山东大学材料与化工325求调剂 +5 半截的诗0927 2026-03-02 5/250 2026-03-02 18:37 by 明亮9527
[考研] 清华大学 材料与化工 353分求调剂 +4 awaystay 2026-03-02 5/250 2026-03-02 17:33 by yeahyou
[考研] 一志愿华中科技大学，化学专业344分，求调剂 +3 邢xing1 2026-03-02 3/150 2026-03-02 17:32 by houyaoxu
[基金申请] 面上模板改不了页边距吧？ +6 ieewxg 2026-02-25 7/350 2026-03-02 12:44 by stidwellNK
[考研] 292求调剂 +7 yhk_819 2026-02-28 7/350 2026-03-02 12:43 by 无际的草原
[考研] 281求调剂 +5 2026计算机_诚心 2026-03-01 8/400 2026-03-02 11:05 by 汪！？！
[考研] 调剂 +3 13853210211 2026-03-02 4/200 2026-03-02 10:16 by 13853210211
[考研] 化工299分求调剂 一志愿985落榜 +5 嘻嘻(*^ω^*) 2026-03-01 5/250 2026-03-01 19:47 by 无际的草原
[考研] 298求调剂 +6 axyz3 2026-02-28 6/300 2026-03-01 19:00 by 18137688336
[考研] 328求调剂 +3 aaadim 2026-03-01 5/250 2026-03-01 17:29 by njzyff
[硕博家园] 2025届双非化工硕士毕业，申博 +3 更多的是 2026-02-27 4/200 2026-03-01 10:04 by ztg729
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
关闭
请填处理意见
关闭