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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 热点前沿问题讨论 » VASP计算时，TOTAL-FORCE这么高是什么原因，求大神
91/1返回列表
查看: 2766  |  回复: 8
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

郭芳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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2014-05-05 20:49:58
已阅   回复此楼   关注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 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭