| 查看: 899 | 回复: 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,所以麻烦了！ |
回复此楼» 猜你喜欢
308求调剂
已经有3人回复
-
321求调剂
已经有11人回复
-
一志愿中海洋材料工程专硕330分求调剂
已经有8人回复
-
330求调剂
已经有4人回复
-
304求调剂
已经有7人回复
-
一志愿南京理工大学085701资源与环境302分求调剂
已经有6人回复
-
南京大学化学376求调剂
已经有3人回复
-
考研调剂求学校推荐
已经有5人回复
-
323求调剂
已经有3人回复
-
317求调剂
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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 (幼儿园)
- 金币: 673.8
- 散金: 819
- 红花: 9
- 帖子: 565
- 在线: 244.3小时
- 虫号: 853446
- 注册: 2009-09-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
送鲜花一朵 
4楼2012-08-21 09:55:00
donggaomin
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 673.8
- 散金: 819
- 红花: 9
- 帖子: 565
- 在线: 244.3小时
- 虫号: 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朵）
3楼2012-08-20 23:10:07
5楼2012-08-21 23:22:45
