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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】用murn.f来拟合晶格常数
51/1返回列表
查看: 1409  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lugang911

[交流] 【求助】用murn.f来拟合晶格常数 已有3人参与

各位虫友，你们好！我用网上的murn.f来拟合晶格常数，得到的out.m文件如下，请问那个是晶格常数？另外侯老师给出的g/alat=是什么意思？
# unit of length =       0.529177  angstroem
# energies input in ELECTRONVOLTS
# unit cell volume = a**3 *       0.250000
# a(min)=           3.5000000
# a(max)=           4.3000000
# points:  50
#  Zero point lattice energy by Moruzzi-Janak-Williams:
#  Grueneisen constant =   3.500
#  T(Debye) =  -5.73
#  Volume/at for T(Debye) =     3.60

#  Tested only for fcc lattice (Al)

# number of data points: n =    3
#----------
#    i     alatt         volume              E
#     1       3.800000      13.718000      -7.7110664890        -7.7110663000
#     2       3.900000      14.829750      -7.8155879439        -7.8155878000
#     3       4.000000      16.000000      -7.7645597103        -7.7645596000
#----------
#     alatt0   vol0 (ula**3)  b0 (mbar)     b0prime        e0
#    3.9000     14.82975      1.00000      4.00000      0.00000
# min:   1.9500      1.85372      0.01000      0.01000
# max:   5.8500     50.05041     10.00000     10.00000
#----------
# fit by murnaghan equation equation of state
#----------
#iter, alatt0, vol0 (ula**3),  b0 (mbar),  b0prime,  e0,  sq sum, ierr
#  20   3.9006    14.83620    17.31752    14.43339    -7.8155902     0.00000 0
#----------
# a (a.u.),   E
#----------
   3.50000          1.262708  9.078298
   3.51633         -0.372188  7.443402
   3.53265         -1.718709  6.096881
   3.54898         -2.828148  4.987442
   3.56531         -3.742495  4.073095
   3.58163         -4.496152  3.319438
   3.59796         -5.117333  2.698257
   3.61429         -5.629198  2.186392
   3.63061         -6.050774  1.764816
   3.64694         -6.397705  1.417885
   3.66327         -6.682869  1.132721
   3.67959         -6.916873  0.898717
   3.69592         -7.108463  0.707127
   3.71224         -7.264853  0.550737
   3.72857         -7.392004  0.423586
   3.74490         -7.494842  0.320749
   3.76122         -7.577440  0.238151
   3.77755         -7.643172  0.172418
   3.79388         -7.694834  0.120756
   3.81020         -7.734744  0.080846
   3.82653         -7.764825  0.050766
   3.84286         -7.786674  0.028916
   3.85918         -7.801621  0.013969
   3.87551         -7.810769  0.004821
   3.89184         -7.815038  0.000552
   3.90816         -7.815195  0.000395
   3.92449         -7.811878  0.003712
   3.94082         -7.805620  0.009970
   3.95714         -7.796866  0.018725
   3.97347         -7.785987  0.029603
   3.98980         -7.773296  0.042295
   4.00612         -7.759050  0.056540
   4.02245         -7.743469  0.072121
   4.03878         -7.726735  0.088855
   4.05510         -7.709000  0.106591
   4.07143         -7.690392  0.125198
   4.08776         -7.671019  0.144571
   4.10408         -7.650972  0.164618
   4.12041         -7.630326  0.185264
   4.13673         -7.609145  0.206446
   4.15306         -7.587482  0.228108
   4.16939         -7.565382  0.250208
   4.18571         -7.542884  0.272706
   4.20204         -7.520019  0.295571
   4.21837         -7.496814  0.318776
   4.23469         -7.473291  0.342299
   4.25102         -7.449470  0.366120
   4.26735         -7.425366  0.390224
   4.28367         -7.400993  0.414598
   4.30000         -7.376361  0.439229
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-06-23 09:45:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lugang911

谢谢！
还有一个问题我一直想不明白就是：
在优化晶格常数时默认情况下是NSW=0，IBRION=-1（原子不动），ISIF=0（离子位置驰豫），其实这个时候就已经是在优化了，ISIF=2仅比0多了项stress tensor,这个有很大区别吗？
还有必要再用ISIF=2吗？另外，也搞不懂的是IBRION=-1时已经表明原子是不动的，而紧跟着ISIF=0又表示离子位置驰豫，这不相互矛盾吗？向各位虫友求证。
赞一下 回复此楼
5楼2010-06-23 15:04:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

condensed

木虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
cenwanglai(金币+1):谢谢回复 2010-06-23 13:06:08
#iter, alatt0, vol0 (ula**3),  b0 (mbar),  b0prime,  e0,  sq sum, ierr
#  20   3.9006    14.83620    17.31752    14.43339    -7.8155902     0.00000 0

alatto 3.9006

做了20步fit后的结果
赞一下(2人) 回复此楼
2楼2010-06-23 09:58:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lugang911

谢谢
是不是要把这个晶格常数写到POSCAR中，并把NSW=15；IBRION=2；ISIF=3加入到INCAR中再做结构优化，从而得到一个稳定的基态结构?
赞一下 回复此楼
3楼2010-06-23 11:15:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
hedaors(金币+1):Thanks, 其实用 isif=3 还是更保险点 2010-06-23 13:20:18
引用回帖:
Originally posted by lugang911 at 2010-06-23 11:15:33:
谢谢
是不是要把这个晶格常数写到POSCAR中，并把NSW=15；IBRION=2；ISIF=3加入到INCAR中再做结构优化，从而得到一个稳定的基态结构?

晶格常数已经优化出来了,就没有必要用ISIF=3俩,2就好.

其它随便,收敛就好.

呵呵:认真想了想,结构优化我还是没掌握到位.

对于原胞,我直接优化晶格参数,然后就用保持形状和体积来优化原子位置了。这个对于简单的晶体结构在能量和结合结构上也能满足一定精度要求.如果是复杂的结构,可能有必要做ISIF=3的细致优化.

这种全面优化在能量和几何参数上能在精度上提高多少(比如受力从 0.02->0.01),我没有概念.希望有人可以提供参考.

[ Last edited by cenwanglai on 2010-6-23 at 13:48 ]
赞一下(2人) 回复此楼
4楼2010-06-23 13:07:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭