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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助成功】晶格优化中，murn.f的使用方法
51/1返回列表
查看: 2378  |  回复: 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的回帖

锐利的碎片

木虫 (正式写手)

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 个回答

锐利的碎片

木虫 (正式写手)

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的回帖

gleerat

木虫 (正式写手)
引用回帖:
Originally posted by 锐利的碎片 at 2010-08-25 18:18:46:
只拟合参数的话，两个单位统一就行了，要算别的就要统一成原子单位。
另外注意下最后列数据的时候列的是体积不是晶格参数。

这是输入文件：
引用回帖:
CaFe2As afm-opt
20
2
11.654 11.674 50
11
349.7400 -138.628584
349.8000 -138.625185
349.8600 -138.621861
349.9200 -138.618444
349.9800 -138.615030
350.0400 -138.611600
350.1000 -138.608153
350.1600 -138.604723
350.2200 -138.601295
350.2800 -138.597855
350.3400 -138.594417

这是输出文件EVPAI.OUT：
引用回帖:
  0.58272000000000002       -7.7674340427867836     
  0.58274000000000004       -7.7674332047055623     
  0.58276000000000006       -7.7674323666222618     
  0.58277999999999996       -7.7674315285369193     
  0.58279999999999998       -7.7674306904495172     
  0.58282000000000000       -7.7674298523600553     
  0.58284000000000002       -7.7674290142685329     
  0.58285999999999993       -7.7674281761749517     
  0.58288000000000006       -7.7674273380793109     
  0.58289999999999997       -7.7674264999816112     
  0.58291999999999999       -7.7674256618818518     
  0.58294000000000001       -7.7674248237800345     
  0.58296000000000003       -7.7674239856761584     
  0.58297999999999994       -7.7674231475702236     
  0.58299999999999996       -7.7674223094622512     
  0.58301999999999998       -7.7674214713522005     
  0.58304000000000000       -7.7674206332400928     
  0.58305999999999991       -7.7674197951259272     
  0.58308000000000004       -7.7674189570097045     
  0.58309999999999995       -7.7674181188914249     
  0.58311999999999997       -7.7674172807710882     
  0.58313999999999999       -7.7674164426486954     
  0.58316000000000001       -7.7674156045242464     
  0.58317999999999992       -7.7674147663977209     
  0.58319999999999994       -7.7674139282691597     
  0.58322000000000007       -7.7674130901385423     
  0.58323999999999998       -7.7674122520058502     
  0.58326000000000000       -7.7674114138711214     
  0.58328000000000002       -7.7674105757343188     
  0.58330000000000004       -7.7674097375954609     
  0.58331999999999995       -7.7674088994545487     
  0.58333999999999997       -7.7674080613115821     
  0.58335999999999999       -7.7674072231665612     
  0.58338000000000001       -7.7674063850194859     
  0.58339999999999992       -7.7674055468703376     
  0.58342000000000005       -7.7674047087191358     
  0.58343999999999996       -7.7674038705659001     
  0.58345999999999998       -7.7674030324105914     
  0.58348000000000000       -7.7674021942532301     
  0.58350000000000002       -7.7674013560937958     
  0.58351999999999993       -7.7674005179323293     
  0.58353999999999995       -7.7673996797687899     
  0.58355999999999997       -7.7673988416031996     
  0.58357999999999999       -7.7673980034355568     
  0.58359999999999990       -7.7673971652658427     
  0.58362000000000003       -7.7673963270940964     
  0.58363999999999994       -7.7673954889202799     
  0.58365999999999996       -7.7673946507444116     
  0.58367999999999998       -7.7673938125664934     
  0.58370000000000000       -7.7673929743865040

这是输出文件EVPAP.OUT：
引用回帖:
   17.487000000000002       -6.9314291999999993     
   17.490000000000002       -6.9312592499999992     
   17.493000000000002       -6.9310930499999994     
   17.496000000000002       -6.9309222000000004     
   17.499000000000002       -6.9307514999999995     
   17.502000000000002       -6.9305800000000009     
   17.505000000000003       -6.9304076499999994     
   17.508000000000003       -6.9302361500000007     
   17.511000000000003       -6.9300647499999997     
   17.513999999999999       -6.9298927500000005     
   17.516999999999999       -6.9297208499999998     

我的问题是：如何从EVPAI.out中得到优化结果的晶格参数？
赞一下 回复此楼
明察、慎思、笃行
5楼2010-08-25 20:19:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 明天应该可查了！？ +4 chengyan1220 2026-08-23 4/200 2026-08-24 08:47 by gloomy6159
[基金申请] 能否退出参与的面上项目解除限项 +6 koalala 2026-08-24 6/300 2026-08-24 08:21 by 正能量1斤
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode，4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 时间戳变了，能看出什么问题？ +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭