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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 为什么我进行Si结构优化的时候只进行了一次离子步的计算？
101/1返回列表
查看: 1943  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭