24小时热门版块排行榜    

查看: 3349  |  回复: 8

郭芳8888

金虫 (正式写手)

[求助] VASP计算时,TOTAL-FORCE这么高是什么原因,求大神 已有2人参与

POSITION                                                 TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      2.07327      0.00000      1.95850        -1.058770      0.000000      0.000000
      3.80631      0.00000      0.00000         0.756301      0.000000      0.000000
      0.03047      3.91700      0.00000        -0.213615      0.000000      0.000000
      3.89534      3.91700      0.00000         0.070406      0.000000      0.000000
      2.00281      1.91209      0.00000        -0.052952      0.044286      0.000000
      5.83582      1.99980      0.00000         0.037517     -0.110169      0.000000
      2.00281      5.92191      0.00000        -0.052952     -0.044286      0.000000
      5.83582      5.83420      0.00000         0.037517      0.110169      0.000000
      0.05807      0.00000      0.00000        -0.090957      0.000000      0.000000
      5.82279      0.00000      1.95850         0.162128      0.000000      0.000000
      0.03286      2.00928      1.95850        -0.093962     -0.067175      0.000000
      0.03286      5.82472      1.95850        -0.093962      0.067175      0.000000
      3.87511      1.90140      1.95850         0.270643      0.205806      0.000000
      1.99360      3.91700      1.95850         0.043551      0.000000      0.000000
      5.83095      3.91700      1.95850         0.008461      0.000000      0.000000
      3.87511      5.93260      1.95850         0.270643     -0.205806      0.000000
-----------------------------------------------------------------------------------
    total drift:                               -0.001301     -0.000000      0.000000
回复此楼

» 猜你喜欢

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

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

c_chi

木虫 (正式写手)

你这个是优化了多少步的结果?
2楼2014-05-05 20:56:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

郭芳8888

金虫 (正式写手)

引用回帖:
2楼: Originally posted by c_chi at 2014-05-05 20:56:16
你这个是优化了多少步的结果?

100步
3楼2014-05-05 21:10:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c_chi

木虫 (正式写手)

引用回帖:
3楼: Originally posted by 郭芳8888 at 2014-05-05 21:10:12
100步

能贴出你的INCAR和POSCAR吗?
4楼2014-05-05 21:26:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

郭芳8888

金虫 (正式写手)

引用回帖:
4楼: Originally posted by c_chi at 2014-05-05 21:26:18
能贴出你的INCAR和POSCAR吗?

好的,这是原始POSCAR
FePt
3.917
2 0 0
0 2 0
0 0 1
4 12
Direct
0.25 0 0.5
0.5 0 0
0 0.5 0
0.5 0.5 0
0.25 0.25 0
0.75 0.25 0
0.25 0.75 0
0.75 0.75 0
0 0 0
0.75 0 0.5
0 0.25 0.5
0 0.75 0.5
0.5 0.25 0.5
0.25 0.5 0.5
0.75 0.5 0.5
0.5 0.75 0.5
5楼2014-05-06 09:16:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

郭芳8888

金虫 (正式写手)

引用回帖:
4楼: Originally posted by c_chi at 2014-05-05 21:26:18
能贴出你的INCAR和POSCAR吗?

System =  FePt
NPAR = 1
VOSKOWN = 1
# ISPIN = 2 spin: 1 paramegnetic; 2 spin polarized
ISTART = 0
ICHARG = 2
LWAVE  = .FALSE.   
LCHARG = .FALSE.   
PREC   =  Accurate
MAGMOM = 3 3 3 3 0 0 0 0 0 0 0 0 0 0 0 0
RWIGS  =  1.111 1.226
Electronic diagonalization:     
ISMEAR = 1  
Parameters related to DOS:
NEDOS  =     901   
EMIN   =    -5.0   
EMAX   =    10.0   
NSW    = 100  
IBRION =  2  
ISIF   =  2         
ISYM   =  2   
LPETIM =  F
6楼2014-05-06 09:24:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c_chi

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
引用回帖:
6楼: Originally posted by 郭芳8888 at 2014-05-06 09:24:52
System =  FePt
NPAR = 1
VOSKOWN = 1
# ISPIN = 2 spin: 1 paramegnetic; 2 spin polarized
ISTART = 0
ICHARG = 2
LWAVE  = .FALSE.   
LCHARG = .FALSE.   
PREC   =  Accurate
MAGMOM = 3 3  ...

看了下你的INCAR,我觉得在结构优化的时候有关磁性的参数可以不添加上去如VOSKOWN,另外PREC没必要取得这么高,只有在你计算力的有关性质的时候可以取高点;还有就是可以把POTIM参数减小一点,K点也不用取高。这里贴上我用的结构优化INCAR,也推荐你看一个帖子,希望对你有用:http://muchong.com/bbs/viewthread.php?tid=2512497&fpage=1

###################INCAR###############
System = Obj

### Start parameters
ISTART = 0     #0 start a new computation, no WAVECAR file -- default 1
ICHARG = 2     #initial charge densities from atoms

### Electronic Relaxation
ENCUT = 400
LREAL = Auto   #projection done in real space
EDIFF = 1.0E-4 #Default 1.0E-4
IALGO = 48     #RMM-DIIS algorithm for electrons
NELMIN = 4     #default 2; do a minimum of four electronic steps

### Ionic relaxation
EDIFFG = -0.02 #Default EDIFF*10
NSW = 300      #default 0
IBRION = 2     #use CG algorithm
ISIF = 2       #relax ions only
POTIM = 0.1    #default 0.5; time-step for ion-motion

### DOS related values
ISMEAR = -5      #method to determine partial occupancies;-5 for semiconductor and non-metal 1 for metal
SIGMA = 0.1      #broadening in eV

### BANDS related values
#NBANDS = 280     #number of occupied bands

### others
PREC = Medium
GGA = PE
NPAR = 8

### Output control
LWAVE = F
LCHARG = F
7楼2014-05-06 09:44:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c_chi

木虫 (正式写手)

引用回帖:
6楼: Originally posted by 郭芳8888 at 2014-05-06 09:24:52
System =  FePt
NPAR = 1
VOSKOWN = 1
# ISPIN = 2 spin: 1 paramegnetic; 2 spin polarized
ISTART = 0
ICHARG = 2
LWAVE  = .FALSE.   
LCHARG = .FALSE.   
PREC   =  Accurate
MAGMOM = 3 3  ...

还有就是在结构优化前最好优化晶格以及必要参数,结构优化前几步的受力你要看一下一般不超过及时,过大则模型有问题;还有就是关注受力的收敛趋势,是不是有减小的趋势或者是在震荡,要视情况采取措施;遇到结构不容易收敛的时候可以查看一下是哪一个原子哪个方向受力较大,然后再作调整
8楼2014-05-06 09:49:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

零下1℃

铁虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
郭芳8888: 金币+5, ★有帮助, 谢谢大神 2014-05-09 23:54:13
在计算磁性材料的时候,初始的几何优化阶段建议先不要添加磁性参数,先优化结构,在静态计算的时候再将磁性参数加上,这样更容易优化,楼主可以试试
9楼2014-05-06 14:20:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 郭芳8888 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2rKfVQg2jQww 2026-11-06 4/200 2026-11-07 10:03 by byEJ59EAM6hR
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 23/1150 2026-11-07 06:49 by Equinoxhua
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 ELBZ6WwEaRj6 2026-11-04 11/550 2026-11-06 16:32 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 15/750 2026-11-06 13:49 by BkA9Bi5cMhkD
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ncTVBM6AOJHF 2026-11-04 17/850 2026-11-06 13:47 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ncTVBM6AOJHF 2026-11-04 11/550 2026-11-06 13:47 by BkA9Bi5cMhkD
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +7 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 13:47 by BkA9Bi5cMhkD
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +8 ncTVBM6AOJHF 2026-11-04 18/900 2026-11-06 13:36 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 13:36 by BkA9Bi5cMhkD
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-05 10/500 2026-11-06 12:32 by BkA9Bi5cMhkD
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 09:01 by 2lRNluKrJpwL
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 05:32 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 12/600 2026-11-06 05:32 by qL1cU46YIdK0
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ttH3f1oA6ace 2026-11-04 14/700 2026-11-06 05:11 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 11/550 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-04 8/400 2026-11-06 02:24 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-06 01:51 by qL1cU46YIdK0
信息提示
请填处理意见