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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 为什么我进行Si结构优化的时候只进行了一次离子步的计算？
101/1返回列表
查看: 2589  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 u1xkze8ybe 2026-07-30 7/350 2026-07-31 12:10 by syhhffkrym
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 u1xkze8ybe 2026-07-30 3/150 2026-07-31 12:08 by syhhffkrym
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 u1xkze8ybe 2026-07-30 6/300 2026-07-31 12:05 by syhhffkrym
[教师之家] 基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊 +6 锦衣卫寒战 2026-07-28 6/300 2026-07-31 11:27 by syl200707
[高分子] HXDI做水性聚氨酯乳液，是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 小木虫看见有人已经知道结果了 +15 1234567wang 2026-07-25 19/950 2026-07-31 08:30 by ZJTJZ
[基金申请] 2026年国自然面上资助率 +8 布布和一二 2026-07-30 8/400 2026-07-31 01:51 by alongwaytogo
[基金申请] 系统今天又提示维护了，估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 微信指数没变化，科研之友没阅读 +14 wangze12014 2026-07-28 17/850 2026-07-30 20:10 by kissu88
[考博] 申请2027年材料/化学类博士 10+3 考研老狗? 2026-07-25 5/250 2026-07-30 19:33 by 白色木鱼
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[基金申请] 今年的WR进展到哪一步了？ +5 wsgjhwz 2026-07-25 9/450 2026-07-30 12:59 by 可淡不可忘
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 时间戳他又来了 +14 晓晓爱翠翠 2026-07-26 17/850 2026-07-29 16:30 by 呵呵二哥
[有机交流] 这个自发加氧反应的机理是什么？ 50+4 YUAN2273 2026-07-25 7/350 2026-07-29 12:27 by xiaobenpu
[基金申请] 没消息就是被刷了呗 +13 lambert2014 2026-07-24 14/700 2026-07-29 10:38 by msh145800
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了，我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
[分析] 水蒸气蒸馏法乳化现象 +3 月光娇娇 2026-07-26 3/150 2026-07-26 23:25 by caiyun
信息提示
关闭
请填处理意见
关闭