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

24小时热门版块排行榜    

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

gleerat

木虫 (正式写手)

[交流] 【求助成功】晶格优化中，murn.f的使用方法

在晶格参数优化中，拟合Murn状态方程需要用到murn.f编译生成的可执行文件，使用的输入文件为inp.m，其内容和格式为：
引用回帖:
2
0.25
6.00 7.45 50
8
6.0849 -.62120891E+01
6.2739 -.64553428E+01
6.4629 -.65246916E+01
6.6518 -.64694519E+01
6.8408 -.63284020E+01
7.0298 -.61271095E+01
7.2188 -.58843994E+01
7.4077 -.56183090E+01
-----------------------------------------

第二行是原胞与晶胞的的体积之比。
这里的晶胞是否就是在优化过程中使用的POSCAR所构建的晶胞，即如果使用原胞构建POSCAR的话，这里的比值就应该是1？

对于VASP来说，使用惯用晶胞优化得出的数据再次构建原胞而进行计算的话，对计算有什么影响么？也就是说使用原胞还是使用惯用晶胞优化对结果有无影响？

在静态自洽计算中使用的kmesh的大小与非自洽中的密集kmesh的大小有什么关系，也就是说满足条件的非自洽计算的kmesh至少要比静态计算的kmesh大多少？

[ Last edited by gleerat on 2010-11-13 at 13:01 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
明察、慎思、笃行
1楼 2010-08-24 17:45:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)
引用回帖:
Originally posted by 锐利的碎片 at 2010-08-25 20:25:30:

最低点体积在param.out里面

我用Murn拟合的结果：
引用回帖:
FIT BY MURNAGHAN EQUATION EQUATION OF STATE
-------------------------------------------------------------------------------

ITER ALATT0   VOL0 (ULA**3)  B0 (MBAR)     B0PRIME      E0 (EV)    SUM OF SQUARES
                                                                           IERR

   20  11.6070    390.93557     44.22016    236.91881   -138.68376      0.00000
                                                                             10
   40  11.6049    390.71578     51.35014    279.16539   -138.69018      0.00000
                                                                             10
   60  11.6074    390.96808     54.84225    298.57001   -138.68727      0.00000
                                                                             10
   80  11.6097    391.20231     57.16045    311.12902   -138.68418      0.00000
                                                                             10
  100  11.6115    391.38848     60.24171    328.23583   -138.68202      0.00000
                                                                             10
  120  11.6140    391.63791     63.73295    347.38332   -138.67885      0.00000
                                                                             10
  140  11.6139    391.63337     63.96723    348.76331   -138.67898      0.00000
                                                                              1
  160  11.6154    391.78090     66.11759    360.53890   -138.67708      0.00000
                                                                             10
  180  11.6176    392.00145     70.13895    382.68500   -138.67434      0.00000
                                                                             10
  200  11.6195    392.19658     73.85296    403.11390   -138.67187      0.00000
                                                                             10
  220  11.6195    392.19099     74.58457    407.31816   -138.67209      0.00000
                                                                              1
  240  11.6219    392.43936     79.45171    433.91890   -138.66886      0.00000
                                                                             10
  260  11.6239    392.63890     85.90905    469.79775   -138.66655      0.00000
                                                                             10
  280  11.6258    392.83065     90.77614    496.44831   -138.66405      0.00000
                                                                             10
  300  11.6257    392.82255     90.96167    497.60759   -138.66420      0.00000
                                                                              1
  320  11.6281    393.06961     98.66533    539.95532   -138.66103      0.00000
                                                                             10
  340  11.6300    393.25841    105.75179    578.99435   -138.65864      0.00000
                                                                             10
  360  11.6312    393.38182    111.21954    609.14271   -138.65709      0.00000
                                                                             10
  380  11.6330    393.56618    120.91667    662.70032   -138.65482      0.00000
                                                                             10
  400  11.6344    393.70253    128.29361    703.35551   -138.65307      0.00000
                                                                             10
  420  11.6356    393.83222    136.90084    750.87243   -138.65145      0.00000
                                                                             10
  440  11.6370    393.97242    147.20622    807.70666   -138.64969      0.00000
                                                                             10
  460  11.6377    394.04404    153.29677    841.35978   -138.64880      0.00000
                                                                             10
  480  11.6377    394.04397    153.37512    841.79671   -138.64880      0.00000
                                                                              1
  500  11.6377    394.04407    153.37242    841.77987   -138.64880      0.00000
                                                                              1
  520  11.6377    394.04441    153.38161    841.82848   -138.64879      0.00000
                                                                              1
  540  11.6377    394.04446    153.41227    842.00057   -138.64879      0.00000
                                                                              1
  560  11.6377    394.04474    153.41587    842.01819   -138.64879      0.00000
                                                                              1
  580  11.6377    394.04466    153.42263    842.05719   -138.64879      0.00000
                                                                              1
  600  11.6380    394.07095    155.59950    854.05452   -138.64846      0.00000
                                                                              1
  620  11.6380    394.07058    155.62739    854.21573   -138.64846      0.00000
                                                                              0
  620  11.6380    394.07058    155.62739    854.21573   -138.64846      0.00000
                                                                              0
-------------------------------------------------------------------------------

      ALATT0   VOL0 (ULA**3)  B0 (MBAR)     B0PRIME       E0 (EV)      E0 (RY)

alat=  11.63799 b0=  155.6274 b0p=   854.216 E0=  -5.095122
       11.6380    394.07058    155.62739    854.21573   -138.64846    -10.19037
        6.1586
-------------------------------------------------------------------------------

        alat(A)          energy(Ryd)            alat(bohr)          ener gy(H)
-------------------------------------------------------------------------------

        6.16703         -10.1889123       11.65400          -5.0944562
        6.16725         -10.1888623       11.65441          -5.0944312
        6.16746         -10.1888121       11.65482          -5.0944061
        6.16768         -10.1887618       11.65522          -5.0943809
        6.16789         -10.1887115       11.65563          -5.0943557
        6.16811         -10.1886610       11.65604          -5.0943305
        6.16833         -10.1886104       11.65645          -5.0943052
        6.16854         -10.1885597       11.65686          -5.0942799
        6.16876         -10.1885090       11.65727          -5.0942545
        6.16897         -10.1884581       11.65767          -5.0942291
        6.16919         -10.1884073       11.65808          -5.0942036
        6.16941         -10.1883563       11.65849          -5.0941782
        6.16962         -10.1883053       11.65890          -5.0941527
        6.16984         -10.1882543       11.65931          -5.0941271
        6.17005         -10.1882032       11.65971          -5.0941016
        6.17027         -10.1881521       11.66012          -5.0940760
        6.17049         -10.1881009       11.66053          -5.0940504
        6.17070         -10.1880497       11.66094          -5.0940248
        6.17092         -10.1879985       11.66135          -5.0939992
        6.17113         -10.1879472       11.66176          -5.0939736
        6.17135         -10.1878959       11.66216          -5.0939479
        6.17157         -10.1878446       11.66257          -5.0939223
        6.17178         -10.1877932       11.66298          -5.0938966
        6.17200         -10.1877418       11.66339          -5.0938709
        6.17221         -10.1876904       11.66380          -5.0938452
        6.17243         -10.1876390       11.66420          -5.0938195
        6.17265         -10.1875876       11.66461          -5.0937938
        6.17286         -10.1875361       11.66502          -5.0937681
        6.17308         -10.1874847       11.66543          -5.0937423
        6.17329         -10.1874332       11.66584          -5.0937166
        6.17351         -10.1873817       11.66624          -5.0936909
        6.17373         -10.1873302       11.66665          -5.0936651
        6.17394         -10.1872787       11.66706          -5.0936393
        6.17416         -10.1872272       11.66747          -5.0936136
        6.17437         -10.1871756       11.66788          -5.0935878
        6.17459         -10.1871241       11.66829          -5.0935620
        6.17481         -10.1870725       11.66869          -5.0935363
        6.17502         -10.1870209       11.66910          -5.0935105
        6.17524         -10.1869694       11.66951          -5.0934847
        6.17545         -10.1869178       11.66992          -5.0934589
        6.17567         -10.1868662       11.67033          -5.0934331
        6.17589         -10.1868146       11.67073          -5.0934073
        6.17610         -10.1867630       11.67114          -5.0933815
        6.17632         -10.1867114       11.67155          -5.0933557
        6.17653         -10.1866598       11.67196          -5.0933299
        6.17675         -10.1866081       11.67237          -5.0933041
        6.17697         -10.1865565       11.67278          -5.0932782
        6.17718         -10.1865049       11.67318          -5.0932524
        6.17740         -10.1864532       11.67359          -5.0932266
        6.17761         -10.1864016       11.67400          -5.0932008

其中有
引用回帖:
      ALATT0   VOL0 (ULA**3)  B0 (MBAR)     B0PRIME       E0 (EV)      E0 (RY)

alat=  11.63799 b0=  155.6274 b0p=   854.216 E0=  -5.095122
       11.6380    394.07058    155.62739    854.21573   -138.64846    -10.19037
        6.1586

是不是说alat=  11.63799是我要的平衡状态下的晶胞参数a的结果？
赞一下 回复此楼
明察、慎思、笃行
8楼2010-08-25 22:43:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

锐利的碎片

木虫 (正式写手)

star watcher
★
zzy870720z(金币+1):谢谢提示 2010-08-25 18:44:45
gleerat(金币+5): 2010-11-13 13:00:35
自我推荐下
http://muchong.com/bbs/viewthread.php?tid=2056902
赞一下(1人) 回复此楼
2楼2010-08-24 20:40:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)
引用回帖:
Originally posted by 锐利的碎片 at 2010-08-24 20:40:09:
自我推荐下
http://muchong.com/bbs/viewthread.php?tid=2056902

谢谢。
我还想知道另外的两个问题的答案，望不吝赐教。

我试着使用了一下eos，它是否不支持以埃和电子伏特为单位？
赞一下 回复此楼
明察、慎思、笃行
3楼2010-08-25 18:07:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

锐利的碎片

木虫 (正式写手)

star watcher
★
zzy870720z(金币+1):谢谢补充 2010-08-25 18:45:01
gleerat(金币+5): 2010-11-13 13:00:42
只拟合参数的话，两个单位统一就行了，要算别的就要统一成原子单位。
另外注意下最后列数据的时候列的是体积不是晶格参数。
赞一下(1人) 回复此楼
4楼2010-08-25 18:18:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭