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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 热点前沿问题讨论 » VASP计算时，TOTAL-FORCE这么高是什么原因，求大神
91/1返回列表
查看: 2916  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 伙伴们，祝我生日快乐吧 +15 myrtle 2026-03-10 24/1200 2026-03-15 21:16 by 苏州_逗号
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 290求调剂 +4 @将就将就看 2026-03-10 8/400 2026-03-14 14:23 by 千千运气
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 材料与化工 一志愿山大 321分 求调剂 +7 每天散步 2026-03-09 8/400 2026-03-14 02:18 by JourneyLucky
[考研] 308求调剂 +4 是Lupa啊 2026-03-09 4/200 2026-03-14 02:06 by tranquil_ya
[考研] 一志愿湖师大化学289求调剂 +6 XMCMM3.14159 2026-03-10 6/300 2026-03-14 00:28 by JourneyLucky
[考研] 2026考研调剂+本科延边大学+山东大学+生物化学与分子生物学+有项目经验 +3 ccdsscjy 2026-03-10 3/150 2026-03-14 00:12 by JourneyLucky
[考研] 311求调剂 +8 zchqwer 2026-03-10 8/400 2026-03-14 00:01 by JourneyLucky
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 材料调剂，307分 +13 张泳铭1 2026-03-09 17/850 2026-03-13 11:09 by 薛云鹏
[考研] 282分材料专业求调剂院校 +18 枫桥ZL 2026-03-09 25/1250 2026-03-13 10:47 by 白夜悠长
[考研] 085600 材料与化工 295 求调剂 +10 dream…… 2026-03-10 12/600 2026-03-12 13:46 by dream……
[考研] 083000环境科学与工程调剂 +8 mingmingry 2026-03-09 9/450 2026-03-11 10:23 by 沙漠之狐994
[考研] 0856材料与化工353分求调剂 +11 NIFFFfff 2026-03-09 11/550 2026-03-10 18:36 by suyuanhai
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
信息提示
关闭
请填处理意见
关闭