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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 按文献值优化Fe3O4晶胞后偏小（附有详细参数）
211/3返回列表123下一页
查看: 2775  |  回复: 20
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

Tea_Candy

金虫 (小有名气)

[求助] 按文献值优化Fe3O4晶胞后偏小（附有详细参数）

按文献上的参数设置后，优化结果差别比较大：
文献上优化后的晶胞参数是 8.473，我的是8.166
我先后尝试过直接用ISIF=3和ISIF=2（然后用状态方程拟合）优化，以及不加U的优化，但结果都在8.16左右。
不知道是哪儿设置错了，恳请大家指教！
文献参数：


这是我设置的参数
INCAR：
ISTART = 0
   ICHARG = 2
   ISPIN = 2

   IBRION = 2
   ISIF   = 3
   NSW   = 100

   PREC   = Medium
   ENCUT  = 500
   NELM   = 100
   NELMIN = 2
   EDIFF  = 1.0e-3
   ISMEAR =  0
   SIGMA  = 0.1
   LWAVE = .FALSE.
   LCHARG = .FALSE.
   IALGO  = 48

   LREAL = AUTO
   GGA = 91
   VOSKOWN = 1
   LMAXMIX = 4
  
   LDAU = .TRUE.
   LDAUL = 2
   LDAUU = 4.5
   LDAUJ = 0.89

POSCAR：
Title
1.0
   8.3941000    0.0000000    0.0000000
   0.0000000    8.3941000    0.0000000
   0.0000000    0.0000000    8.3941000
32 24
Direct
0.37990 0.37990 0.37990
0.37990 0.87990 0.87990
0.87990 0.37990 0.87990
0.87990 0.87990 0.37990
0.62010 0.12010 0.87990
0.62010 0.62010 0.37990
0.12010 0.12010 0.37990
0.12010 0.87990 0.62010
0.12010 0.37990 0.12010
......

KPOINTS:
Auto
0
M
6 6 6
0 0 0

赝势文件用的是PAW_GGA
回复此楼

» 收录本帖的淘帖专辑推荐

DFT计算 建模

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
Life is more than Chemistry.
1楼 2012-05-17 10:16:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzx2008

荣誉版主 (著名写手)

优秀版主优秀版主

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
liliangfang: 金币+2, 谢谢指教 2012-05-17 17:16:40
1.EDIFF  = 1.0e-3有点小，至少应该-4次方
2.GGA+U的设置，这样写应该会报错，因为POSCAR中有两种类型原子
赞一下 回复此楼
2楼2012-05-17 15:39:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gougou531

禁虫 (正式写手)
★
感谢参与，应助指数 +1
liliangfang: 金币+1, 谢谢指教 2012-05-17 17:16:48
本帖内容被屏蔽
3楼2012-05-17 16:36:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Tea_Candy

金虫 (小有名气)
引用回帖:
2楼: Originally posted by fzx2008 at 2012-05-17 15:39:24:
1.EDIFF  = 1.0e-3有点小，至少应该-4次方
2.GGA+U的设置，这样写应该会报错，因为POSCAR中有两种类型原子

谢谢回复。
请问我想把U加在第二个原子上，LDA+U这样设置对么：

   LDAU = .TRUE.
   LDAUL = -1 2
   LDAUU = 0.00 4.5
   LDAUJ = 0.00 0.89
赞一下 回复此楼
Life is more than Chemistry.
4楼2012-05-17 20:20:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Tea_Candy

金虫 (小有名气)
引用回帖:
3楼: Originally posted by gougou531 at 2012-05-17 16:36:08:
能量设置太低了，可能收敛到不正确的基态，我觉得和这个有很大的关系，还有就是你设置u的方式不对，同意2楼的说法。如果可以的话可以考虑以力为收敛标准，设置EDIFFG这个参数，一般设置为-0.01

谢谢回复。
我试试看
赞一下 回复此楼
Life is more than Chemistry.
5楼2012-05-17 20:21:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Tea_Candy

金虫 (小有名气)
引用回帖:
2楼: Originally posted by fzx2008 at 2012-05-17 15:39:24:
1.EDIFF  = 1.0e-3有点小，至少应该-4次方
2.GGA+U的设置，这样写应该会报错，因为POSCAR中有两种类型原子

按照这个参数优化后， 结果又偏大了一点8.5333。实验值是~8.39，文献值GGA+U优化时8.473。请问我这个结果能接受吗？
赞一下 回复此楼
Life is more than Chemistry.
6楼2012-05-17 23:07:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liu108

木虫 (正式写手)

【答案】应助回帖

★
感谢参与，应助指数 +1
WDD880227: 金币+1, 感谢交流提示 2012-05-18 08:56:03
计算值和实验值只差不超过5%是可以接受的，另外，没见过91+U计算的，只见过PBE+U计算的。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

Tea_Candy
会当凌绝顶，一览众山小。
7楼2012-05-18 08:39:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenyuchen

木虫 (正式写手)

【答案】应助回帖

★
感谢参与，应助指数 +1
fzx2008: 金币+1, 鼓励交流 2012-05-18 12:08:20
引用回帖:
7楼: Originally posted by liu108 at 2012-05-18 08:39:52:
计算值和实验值只差不超过5%是可以接受的，另外，没见过91+U计算的，只见过PBE+U计算的。

同意这个。
lz计算的和文献报道的差别应该是交换关联势不同造成的。
赞一下 回复此楼
选择科学就要耐得住寂寞！
8楼2012-05-18 10:20:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenyuchen

木虫 (正式写手)

【答案】应助回帖

★
fzx2008: 金币+1, 鼓励交流 2012-05-18 12:08:14
有个疑问，文献说GGA用的是PW91，GGA+U没说GGA用的是哪一种，难道默认和前面一样吗？
赞一下 回复此楼
选择科学就要耐得住寂寞！
9楼2012-05-18 10:24:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Tea_Candy

金虫 (小有名气)
送鲜花一朵
引用回帖:
7楼: Originally posted by liu108 at 2012-05-18 08:39:52:
计算值和实验值只差不超过5%是可以接受的，另外，没见过91+U计算的，只见过PBE+U计算的。

谢谢您的回复，那我计算的误差为~1.6%就是可以的了。
关于那个交换关联势的选择，我看文献中就是用的GGA-PW91，没提到PBE，我把文献上传上来，一起研究研究吧
赞一下 回复此楼

» 本帖附件资源列表
Life is more than Chemistry.
10楼2012-05-18 13:21:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Tea_Candy 的主题更新
211/3返回列表123下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭