24小时热门版块排行榜

返回列表

donggaomin

金虫 (正式写手)

应助: 1 (幼儿园)
金币: 770.6
散金: 805
红花: 9
帖子: 564
在线: 242.7小时
虫号: 853446
注册: 2009-09-21
性别: GG
专业: 凝聚态物性 II ：电子结构

[求助] 请帮忙一下，把下面的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,所以麻烦了！

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

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人回复

好好学习

1楼 2012-08-20 16:03:55

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

donggaomin

金虫 (正式写手)

应助: 1 (幼儿园)
金币: 770.6
散金: 805
红花: 9
帖子: 564
在线: 242.7小时
虫号: 853446
注册: 2009-09-21
性别: GG
专业: 凝聚态物性 II ：电子结构

有人可以帮忙吗？自己顶一下先

赞一下

回复此楼

好好学习

2楼2012-08-20 20:02:53

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

xh512

金虫 (正式写手)

应助: 2 (幼儿园)
金币: 149.1
散金: 1339
红花: 11
帖子: 309
在线: 923.6小时
虫号: 461655
注册: 2007-11-18
专业: 金属材料的合金相、相变及

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
fzx2008: 金币+2, 谢谢指教！ 2012-08-21 01:01:44
donggaomin: 金币+20, ★★★★★最佳答案, good idea! 我现在装一个。不知道安装是否麻烦？ 2012-08-22 09:09:02

VASP里可以实现你想要的。需要修改源代码： constr_cell_relax.F。
细节你可以在本版搜索一下。

赞一下

回复此楼