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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 为什么我进行Si结构优化的时候只进行了一次离子步的计算？
101/1返回列表
查看: 2243  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

发光的橘子

铜虫 (小有名气)

[求助] 为什么我进行Si结构优化的时候只进行了一次离子步的计算？

不知道为什么,我的离子步只进行了一次，求高手解答~
输出如下：
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1     0.494326194676E+03    0.49433E+03   -0.71069E+04  1088   0.107E+03
DAV:   2    -0.243342847152E+03   -0.73767E+03   -0.70315E+03  1312   0.172E+02
DAV:   3    -0.295807755585E+03   -0.52465E+02   -0.51726E+02  1336   0.409E+01
DAV:   4    -0.297111774949E+03   -0.13040E+01   -0.12964E+01  1288   0.794E+00
DAV:   5    -0.297154959370E+03   -0.43184E-01   -0.43072E-01  1352   0.131E+00    0.158E+01
RMM:   6    -0.293713652362E+03    0.34413E+01   -0.19656E+00  1088   0.236E+00    0.942E+00
RMM:   7    -0.292143136879E+03    0.15705E+01   -0.43861E+00  1088   0.360E+00    0.793E-01
RMM:   8    -0.292181064951E+03   -0.37928E-01   -0.13008E-01   821   0.706E-01
   1 F= -.29218106E+03 E0= -.29218286E+03  d E =-.292181E+03
curvature:   0.00 expect dE= 0.000E+00 dE for cont linesearch  0.000E+00
trial: gam= 0.00000 g(F)=  0.374E-03 g(S)=  0.543E+00 ort = 0.000E+00 (trialstep = 0.100E+01)
search vector abs. value=  0.543E+00
reached required accuracy - stopping structural energy minimisation
writing wavefunctions


我的INCAR如下：

SYSTEM = Silicon
  NPAR = 5
  LPLANE = .TRUE.
  LSCALU = .FALSE.
  NSIM = 4
  ALGO = FAST
  LREAL = .FALSE.
  ISTART = 0
  PREC = High
  EDOS = 900
  GGA = PE
  NSW = 100
  IBRION = 2
  ISIF = 2
  LORBIT = 11
  ENCUT = 400.0 eV
  EDIFF = 1.0E-5
  EDIFFG = -0.01
  NELM = 100
  NELMIN = 2
  ISMEAR = 1
  SIGMA-tag
  ISPIN = 1

POSCAR:

Si
5.38936
0.5 0.5 0.0
0.0 0.5 0.5
0.5 0.0 0.5
2
Cartesian
0.0000000 0.0000000 0.0000000
0.2500000 0.2500000 0.2500000
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
还记得那个凉凉的深秋
1楼 2012-03-02 19:11:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzx2008

荣誉版主 (著名写手)

优秀版主优秀版主

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
liliangfang(金币+2): 多谢指教 2012-03-03 14:59:35
因为Si的对称性很高，面心立方结构，其中原子位置都由对称性固定下来了。

这儿唯一的变量就是晶胞参数a， 再说ISIF = 2，故只有一个离子步~
赞一下 回复此楼
2楼2012-03-02 19:53:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

发光的橘子

铜虫 (小有名气)
引用回帖:
楼: Originally posted by fzx2008 at 2012-03-02 19:53:11:
因为Si的对称性很高，面心立方结构，其中原子位置都由对称性固定下来了。

这儿唯一的变量就是晶胞参数a， 再说ISIF = 2，故只有一个离子步~

用ISIF=3也没用，一样的效果
赞一下 回复此楼
还记得那个凉凉的深秋
3楼2012-03-02 20:05:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

发光的橘子

铜虫 (小有名气)
刚刚把http://muchong.com/bbs/viewthread.php?tid=1404901里面的内容复制了一遍，结果一样，求助~
赞一下 回复此楼
还记得那个凉凉的深秋
4楼2012-03-02 20:08:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzx2008

荣誉版主 (著名写手)

优秀版主优秀版主

【答案】应助回帖

★
uuv2010(金币+1): 多谢提示 2012-03-03 11:07:32
引用回帖:
3楼: Originally posted by 发光的橘子 at 2012-03-02 20:05:27:
用ISIF=3也没用，一样的效果

是的，这个晶格常数可能是优化后的结果了
赞一下 回复此楼
5楼2012-03-02 20:12:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

发光的橘子

铜虫 (小有名气)
引用回帖:
楼: Originally posted by fzx2008 at 2012-03-02 20:12:49:
是的，这个晶格常数可能是优化后的结果了

但是，晶格常数减了1A也一样
赞一下 回复此楼
还记得那个凉凉的深秋
6楼2012-03-02 20:30:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzx2008

荣誉版主 (著名写手)

优秀版主优秀版主

【答案】应助回帖

★ ★ ★ ★
uuv2010(金币+1): 欢迎参加讨论 2012-03-03 11:08:12
发光的橘子(金币+3): ★★★很有帮助 2012-03-06 01:23:01
引用回帖:
6楼: Originally posted by 发光的橘子 at 2012-03-02 20:30:23:
但是，晶格常数减了1A也一样

1.如果晶格常数减了1A，ISIF=3，优化时还只用1个离子步就比较奇怪了。

2.想这种对称性特高的，只有一个自由度的，很多人采用能量拟合的方法。
在实验晶格常数左右割取一些值进行单点能量计算。再拟合E-V得到优化的结果。
赞一下 回复此楼
7楼2012-03-02 20:35:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

发光的橘子

铜虫 (小有名气)
引用回帖:
楼: Originally posted by fzx2008 at 2012-03-02 20:35:30:
1.如果晶格常数减了1A，ISIF=3，优化时还只用1个离子步就比较奇怪了。

2.想这种对称性特高的，只有一个自由度的，很多人采用能量拟合的方法。
在实验晶格常数左右割取一些值进行单点能量计算。再拟合E-V得到 ...

可以确定的是，VASP程序应该没什么问题，因为算其他体系都很正常
赞一下 回复此楼
还记得那个凉凉的深秋
8楼2012-03-02 20:42:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

发光的橘子

铜虫 (小有名气)
引用回帖:
楼: Originally posted by fzx2008 at 2012-03-02 20:35:30:
1.如果晶格常数减了1A，ISIF=3，优化时还只用1个离子步就比较奇怪了。

2.想这种对称性特高的，只有一个自由度的，很多人采用能量拟合的方法。
在实验晶格常数左右割取一些值进行单点能量计算。再拟合E-V得到 ...

不过你的看法很有可能是正确的，极有可能是对称性的问题，我试着把一个Si换成了N原子，结果就正常运行了。但是纯Si也不至于才运行一步啊，搞不大懂……
赞一下 回复此楼
还记得那个凉凉的深秋
9楼2012-03-02 21:38:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rex881026

木虫 (小有名气)

【答案】应助回帖

★
感谢参与，应助指数 +1
liliangfang(金币+1): 嘻嘻交流 2012-03-03 15:01:43
IBRION = 2的话
那你POTIM取个0.5试试
赞一下 回复此楼
10楼2012-03-03 14:27:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 发光的橘子 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 291 求调剂 +29 化工2026届毕业� 2026-04-09 29/1450 2026-04-10 22:55 by dick_runner
[考研] 08工学 309分求调剂 +6 Yin DY 2026-04-08 6/300 2026-04-10 09:18 by Delta2012
[考研] 312求调剂 +3 李鸿飞飞 2026-04-06 3/150 2026-04-09 17:32 by wp06
[考研] 337求调剂 +4 Gky09300550， 2026-04-09 4/200 2026-04-09 17:18 by 帕尔马拉特
[考研] 复试调剂，一志愿郑州大学材料与化工289分 +31 硕星赴 2026-04-08 31/1550 2026-04-09 16:54 by Delta2012
[考研] 262求调剂 +10 天下第一文 2026-04-04 13/650 2026-04-09 15:16 by 探123
[考研] 材料调剂 +14 一样YWY 2026-04-06 14/700 2026-04-08 23:00 by 猪会飞
[考研] 生物学学硕，初试351分，求调剂 +4 …～、王…～ 2026-04-08 5/250 2026-04-08 21:49 by limeifeng
[考研] 一志愿南昌大学，085600，344分求调剂 +11 调剂上岸玘 2026-04-05 12/600 2026-04-08 16:17 by luoyongfeng
[考研] 315求调剂 +17 欣喜777 2026-04-04 18/900 2026-04-08 13:54 by hangsimei
[考研] 085602调剂 初试总分335 +3 19123253302 2026-04-06 3/150 2026-04-07 18:00 by jp9609
[考研] 325 调剂 +6 QQ小虾 2026-04-07 6/300 2026-04-07 15:17 by Ccclqqq
[考研] 一志愿南航，数一英一学硕317求调剂！！ +6 Acaciad 2026-04-04 6/300 2026-04-06 12:13 by 考研学校招点人
[考研] 332求调剂 +17 小小孟... 2026-04-05 18/900 2026-04-06 09:51 by 蓝云思雨
[考研] 081200-11408-276学硕求调剂 +4 崔wj 2026-04-04 5/250 2026-04-05 14:06 by imissbao
[考研] 358求调剂 +7 秋gk 2026-04-04 7/350 2026-04-05 13:29 by huangmoli
[考研] 313求调剂 +3 海日海日 2026-04-04 3/150 2026-04-05 07:48 by 544594351
[考研] 调剂 +8 熊二想上岸 2026-04-04 8/400 2026-04-05 05:27 by houyaoxu
[考研] 320求调剂 +3 一样圆 2026-04-04 3/150 2026-04-04 22:29 by 啵啵啵0119
[考研] 325求调剂 +4 春风不借意 2026-04-04 4/200 2026-04-04 22:08 by 啵啵啵0119
信息提示
关闭
请填处理意见
关闭