24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
查看: 17917  |  回复: 128
【奖励】 本帖被评价101次,作者锐利的碎片增加金币 82.4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

锐利的碎片

木虫 (正式写手)


[资源] 【原创】elk中的状态方程拟合工具eos

目前相当好的拟合工具。使用方法可以看下里面的说明。
-----------------------------
                        +---------------------------+
                         |     EOS Version 1.4.0     |
                         +---------------------------+

Equation of state (EOS) program for fitting energy-volume data. The following
variables are set in the file eos.in:

cname               : name of crystal up to 256 characters
natoms              : number of atoms in unit cell
etype               : equation of state type (see below)
vplt1, vplt2, nvplt : volume interval over which to plot energy, pressure etc.
                       as well as the number of points in the plot
nevpt               : number of energy-volume points to be inputted
vpt(i) ept(i)       : energy-volume points (atomic units)

Note that the input units are atomic - Bohr and Hartree (NOT Rydbergs).

The equations of state currently implemented are:
1. Universal EOS (Vinet P et al., J. Phys.: Condens. Matter 1, p1941 (1989))
2. Murnaghan EOS (Murnaghan F D, Am. J. Math. 49, p235 (1937))
3. Birch-Murnaghan 3rd-order EOS (Birch F, Phys. Rev. 71, p809 (1947))
4. Birch-Murnaghan 4th-order EOS
5. Natural strain 3rd-order EOS (Poirier J-P and Tarantola A, Phys. Earth
    Planet Int. 109, p1 (1998))
6. Natural strain 4th-order EOS
7. Cubic polynomial in (V-V0)

--------------------------------------------------------------------------------
J. K. Dewhurst
August 2005
http://d.namipan.com/d/4822647a13c38712c26286e291d0cdb12e7987b8fa430400
10楼有gfortran编译的可执行文件和源码。elk虽然是GPL的,大家使用时还是注意下版权问题。

[ Last edited by 锐利的碎片 on 2010-5-21 at 13:10 ]
回复此楼

» 收录本帖的淘帖专辑推荐

研究生资料集 第一性原理计算辅助工具 分子动力学 good resources
模拟,理论 VASP blabla 第一性原理计算-tg

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

★★★★★ 五星级,优秀推荐

不知道大家一般选哪个状态方程?还有一个就是单位的问题,我一般直接选择E,V的输入单位就是eV,A,除了拟合出来的E-V单位应该还是这个外,我觉得压强单位应该变了,就是得到结果的压强虽然以GPa为单位,但是值肯定不一样的,不知道大家怎么换算?如下为得到的P-V,明显P很大啊。
PVPAT.OUT:
   20.8900000000000        14969.9645177969     
   21.8050000000000        10902.7771633449     
   22.7500000000000        7467.83064949064     
   23.7200000000000        4597.20824565116     
   24.7150000000000        2210.71355519284     
   25.7400000000000        230.226553695922     
   26.7900000000000       -1392.22471282881     
   27.8700000000000       -2714.94690449880     
   28.9800000000000       -3779.92212009348     
   30.1200000000000       -4624.07873273586
13楼2010-10-14 22:42:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 129 个回答

lfhuang

木虫 (著名写手)


★★★★★ 五星级,优秀推荐

好贴,支持!并已经转载到个人博客里面了。
4楼2010-05-20 10:16:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangdan9604

铁杆木虫 (正式写手)


error while loading shared libraries: libguide.so: cannot open shared object file: No such file or directory.
楼主能否把使用说的详细一点, 我 ./eos < eos.in 为什么会产生上面的错误啊,苦恼中!!
6楼2010-05-20 12:40:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

calos818

木虫 (著名写手)


★★★★★ 五星级,优秀推荐

帮楼上的顶下,希望楼主解决下楼上的问题呵呵,谢谢楼主分享!!!
7楼2010-05-20 13:43:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复(可上传附件)
信息提示
请填处理意见