| 查看: 1071 | 回复: 4 | ||
donggaomin金虫 (正式写手)
|
[求助]
请帮忙一下，把下面的VASP输入文件转成ABINIT的输入文件,送金币
|
|
近来需要补计算一点工作，就是固定一条轴的晶格（C轴），只驰豫另两条轴（这两条轴a,b 与固定轴是垂直的）。这个在VASP里难以实现，而在ABINIT里可以容易地实现，所以想麻烦一下您： INCAR: SYSTEM = BFO GGA = PS #赝势用PBEsol LDAU = .TRUE. # +U LDAUTYPE = 2 #U的类型为Dudarev提出的，只有一个Ueff=4 LDAUL = -1 2 -1 #对第二种原子使用+U LDAUU = 0.0 4.0 0.0 #U=4 LDAUJ = 0.0 0.0 0.0 #J＝0 LMAXMIX = 4 #说明是对d电子使用U ISPIN = 2 #打开自旋 MAGMOM = 18*0 6 6 6 -6 -6 -6 6*0 #每个原子的初始磁矩 ISTART = 0 ICHARG = 2 ENCUT = 500 ISMEAR = 0; SIGMA = 0.2 NSW= 100; IBRION=2 #驰豫步数 ISIF= 3 #原想是全优化，但现在是需要在[001]方向加单轴压力（即固定C轴），其它两轴全优化 ISYM = 0 #去掉对称性 EDIFF =1E-5; EDIFFG = -1E-3 #收敛精度设置，单位eV PREC = Accurate KPOINTS: BFO 0 Gamma 4 4 2 0 0 0 POSCAR: #就是想固定C轴而优化a,b 轴 BiFeO3 1.000000000000000 4.9844546183992420 -2.8777762156962770 0.0000000000000000 0.0000000000000000 5.7555524313925510 0.0000000000000000 0.0000000000000000 0.0000000000000000 11.7002132207541400 O Fe Bi 18 6 6 Direct 0.3140625285079395 0.2445237791076238 0.1278824980496154 0.9807291951746109 0.5778571124409595 0.4612158313829440 0.6473958618412823 0.9111904457742881 0.7945491647162797 0.7554762208923762 0.0695387494003228 0.1278824980496154 0.4221428875590405 0.4028720827336585 0.4612158313829440 0.0888095542257119 0.7362054160669871 0.7945491647162797 0.9304612505996772 0.6859374714920534 0.1278824980496154 0.5971279172663415 0.0192708048253962 0.4612158313829440 0.2637945839330129 0.3526041381587248 0.7945491647162797 0.7554762208923762 0.6859374714920534 0.6278824980496154 0.4221428875590405 0.0192708048253962 0.9612158313829511 0.0888095542257119 0.3526041381587248 0.2945491647162797 0.9304612505996772 0.2445237791076238 0.6278824980496154 0.5971279172663415 0.5778571124409595 0.9612158313829511 0.2637945839330129 0.9111904457742881 0.2945491647162797 0.3140625285079395 0.0695387494003228 0.6278824980496154 0.9807291951746109 0.4028720827336585 0.9612158313829511 0.6473958618412823 0.7362054160669871 0.2945491647162797 0.0000000000000000 0.0000000000000000 0.7289617914542390 0.6666666666666643 0.3333333333333357 0.0622951247875747 0.3333333333333357 0.6666666666666643 0.3956284581209033 0.0000000000000000 0.0000000000000000 0.2289617914542390 0.6666666666666643 0.3333333333333357 0.5622951247875747 0.3333333333333357 0.6666666666666643 0.8956284581209033 0.0000000000000000 0.0000000000000000 0.9823906908530873 0.6666666666666643 0.3333333333333357 0.3157240241864230 0.3333333333333357 0.6666666666666643 0.6490573575197516 0.0000000000000000 0.0000000000000000 0.4823906908530873 0.6666666666666643 0.3333333333333357 0.8157240241864230 0.3333333333333357 0.6666666666666643 0.1490573575197516 我不懂ABINIT,所以麻烦了！ |
回复此楼» 猜你喜欢
能否退出参与的面上项目解除限项
已经有18人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有3人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
估计是周四
已经有4人回复
-
我面上完蛋了
已经有9人回复
-
2026国自然函评费到账
已经有18人回复
-
建议基金发布提前给出明确的时间点
已经有15人回复
-
范进中举一文的中心思想
已经有7人回复
-
朋友圈看到的
已经有9人回复
-
让我中一个面上吧！
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- ABINIT和VASP相比，计算GW，哪个更快啊？
已经有13人回复
- 求助VASP计算有效质量的输入文件示例
已经有3人回复
- 使用PGI编译VASP中无法使用长数组
已经有5人回复
- 求VASP简介。。。。。想学习下VASP
已经有21人回复
- VASP计算能带时是否读取WAVECAR
已经有16人回复
- vasp自洽计算K点不匹配问题，急！！
已经有15人回复
- Vasp中删除一个键 输入文件如何表现出来？
已经有8人回复
- 【求助】关于输入参数MD.MaxForceTol
已经有5人回复
- 【求助】求助abinit和vasp计算的总能的差异
已经有6人回复
- 【求助】求将vasp的输出文件转化成xyz坐标的脚本
已经有11人回复
- 【转帖】VASP和ABINIT优缺点比较
已经有18人回复
- 【求助】请帮下决心：买VASP（经费不多）还是免费的（如abinit）？
已经有11人回复
- 【求助】vasp磁矩计算
已经有8人回复
donggaomin
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 651.8
- 散金: 838
- 红花: 9
- 帖子: 567
- 在线: 244.4小时
- 虫号: 853446
- 注册: 2009-09-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2012-08-20 20:02:53
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
fzx2008: 金币+2, 谢谢指教！ 2012-08-21 01:01:44
donggaomin: 金币+20, ★★★★★最佳答案, good idea! 我现在装一个。不知道安装是否麻烦？ 2012-08-22 09:09:02
感谢参与，应助指数 +1
fzx2008: 金币+2, 谢谢指教！ 2012-08-21 01:01:44
donggaomin: 金币+20, ★★★★★最佳答案, good idea! 我现在装一个。不知道安装是否麻烦？ 2012-08-22 09:09:02
|
VASP里可以实现你想要的。需要修改源代码： constr_cell_relax.F。 细节你可以在本版搜索一下。 |
» 本帖已获得的红花（最新10朵）
donggaomin3楼2012-08-20 23:10:07
donggaomin
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 651.8
- 散金: 838
- 红花: 9
- 帖子: 567
- 在线: 244.4小时
- 虫号: 853446
- 注册: 2009-09-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2012-08-21 09:55:00
5楼2012-08-21 23:22:45
