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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 热点前沿问题讨论 » VASP计算时，TOTAL-FORCE这么高是什么原因，求大神
北京石油化工学院2026年研究生招生接收调剂公告
91/1返回列表
查看: 2953  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 找调剂 +10 楚乔乔 2026-04-01 10/500 2026-04-05 22:19 by syh9288
[考研] 一志愿北京交通大学材料工程总分358求调剂 +4 cs0106 2026-04-04 4/200 2026-04-05 18:46 by imissbao
[考研] 326求调剂 +3 顾若浮生 2026-04-05 3/150 2026-04-05 18:32 by 蓝云思雨
[考研] 275求调剂 +16 waltzh 2026-04-01 16/800 2026-04-05 17:14 by Hdyxbekcb
[考研] 0703化学 +12 妮妮ninicgb 2026-04-04 13/650 2026-04-05 10:46 by 啊俊！
[考研] 311分 22408 求调剂 +3 bing_bot 2026-04-03 3/150 2026-04-05 00:43 by chongya
[考研] 材料调剂 +11 一样YWY 2026-04-02 13/650 2026-04-04 23:10 by 无际的草原
[考研] 环境科学与工程334分求调剂 +9 王一一依依 2026-03-30 12/600 2026-04-04 20:55 by dongzh2009
[考研] 363求调剂 +7 zh096 2026-04-04 7/350 2026-04-04 17:13 by dongzh2009
[考研] 一志愿沪985，326分求调剂 +3 刘墨墨 2026-04-03 3/150 2026-04-04 11:16 by 悲伤的芋头
[考研] 294求调剂 +6 Grey_Ey 2026-04-03 6/300 2026-04-03 20:46 by 欣喜777
[考研] 335求调剂 +7 沈清璃 2026-04-03 7/350 2026-04-03 18:55 by lijunpoly
[考研] 266求调剂 +18 阳阳哇塞 2026-04-01 18/900 2026-04-03 18:38 by zllcz
[考研] 0856，269分求调剂 +15 有学上就行求求� 2026-03-30 18/900 2026-04-03 16:50 by melodiousnow
[考研] 求调剂机会 +5 意染ivy 2026-04-03 5/250 2026-04-03 15:13 by qoooooo614
[考研] 282求调剂 +5 呼吸都是减肥 2026-03-31 5/250 2026-04-03 12:03 by 1753564080
[考研] 一志愿北京科技材料科学与工程288分，求调剂 +14 是辰啊 2026-04-02 14/700 2026-04-02 21:10 by dongzh2009
[考研] 0817化工学硕调剂 +11 努力上岸中！ 2026-03-31 11/550 2026-04-01 20:30 by 赖春艳
[考研] 一志愿 南京航空航天大学 ，080500材料科学与工程学硕 +10 @taotao 2026-03-31 11/550 2026-04-01 09:43 by xiayizhi
[考研] 一志愿浙江大学工科动力工程370,数一121,专业课135，现在能去哪里 +3 080700调剂 2026-03-30 4/200 2026-03-31 12:00 by KLMY666
信息提示
关闭
请填处理意见
关闭