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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 一次型NiAl结构优化，晶格矢量变化非常大，求助~
61/1返回列表
查看: 1622  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料化工调剂 +12 今夏不夏 2026-03-01 14/700 2026-03-02 16:09 by 今夏不夏
[考研] 求调剂 +9 yunziaaaaa 2026-03-01 10/500 2026-03-02 15:18 by 半截的诗0927
[考研] 0856化工专硕求调剂 +15 董boxing 2026-03-01 15/750 2026-03-02 15:06 by 晃晃不许晃
[考研] 290分材料工程085601求调剂 数二英一 +3 llx0610 2026-03-02 3/150 2026-03-02 14:15 by yc258
[考研] 材料与化工328求调剂 +3 。，。，。，。i 2026-03-02 3/150 2026-03-02 13:09 by houyaoxu
[考研] 求调剂 +3 熬夜的猫头鹰 2026-03-02 3/150 2026-03-02 11:45 by 刘兵
[考研] 材料学硕318求调剂 +14 February_Feb 2026-03-01 16/800 2026-03-02 11:17 by yuchj
[考研] 一志愿郑大材料学硕298分，求调剂 +6 wsl111 2026-03-01 6/300 2026-03-02 11:00 by ydudjddnd
[考研] 274求调剂 +3 cgyzqwn 2026-03-01 7/350 2026-03-02 10:38 by lature00
[考研] 材料类求调剂 +11 wana_kiko 2026-02-28 14/700 2026-03-02 08:46 by 聪明的大松鼠
[基金申请] 成果系统访问量大，请一小时后再尝试。---NSFC啥时候好哦，已经两天这样了 +4 NSFC2026我来了 2026-02-28 4/200 2026-03-01 22:37 by 铁门栓
[考研] 0856求调剂285 +10 吕仔龙 2026-02-28 10/500 2026-03-01 21:37 by 公瑾逍遥
[考研] 299求调剂 +3 Y墨明棋妙Y 2026-02-28 5/250 2026-03-01 21:01 by tangxiaotian
[考研] 298求调剂 +6 axyz3 2026-02-28 6/300 2026-03-01 19:00 by 18137688336
[考研] 0856材料求调剂 +11 hyf hyf hyf 2026-02-28 12/600 2026-03-01 18:57 by 18137688336
[考研] 课题组接收材料类调剂研究生 +3 gaoxiaoniuma 2026-02-28 4/200 2026-03-01 14:30 by jjj三跨
[考研] 302材料工程求调剂 +4 Doleres 2026-03-01 5/250 2026-03-01 11:52 by liqiongjy
[硕博家园] 2025届双非化工硕士毕业，申博 +3 更多的是 2026-02-27 4/200 2026-03-01 10:04 by ztg729
[考研] 307求调剂 +4 73372112 2026-02-28 6/300 2026-03-01 00:04 by ll247
[考研] 304求调剂 +3 52hz~~ 2026-02-28 5/250 2026-03-01 00:00 by 52hz~~
信息提示
关闭
请填处理意见
关闭