| 查看: 18437 | 回复: 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 |
» 猜你喜欢
085600材料与化工
已经有6人回复
-
281求调剂（0805）
已经有24人回复
-
265求调剂
已经有8人回复
-
0817 化学工程 299分求调剂 有科研经历 有二区文章
已经有21人回复
-
295材料求调剂，一志愿武汉理工085601专硕
已经有4人回复
-
0856调剂，是学校就去
已经有9人回复
-
298-一志愿中国农业大学-求调剂
已经有9人回复
-
求调剂
已经有4人回复
-
学校已经提交到NSFC，还能修改吗？
已经有7人回复
-
材料学硕297已过四六级求调剂推荐
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Murnaghan和Birch–Murnaghan状态方程（EOS）的拟合结果
已经有28人回复
- origin拟合Birch-Murnaghan状态方程求B0，B0'和v0
已经有38人回复
- 急切求助EOS拟合
已经有21人回复
- 有关拟合方程方面的问题，求大神帮忙
已经有11人回复
- 请教如何用VASP计算氧化物的晶格常数？以alpha-Al2O3为例
已经有7人回复
- 求助 用Oringin拟合转状态方程
已经有6人回复
- GCA-EOS是什么状态方程
已经有6人回复
- Elk与QE联合使用，是否可行？
已经有10人回复
- 在Oringin里面怎么拟合状态方程EOS
已经有7人回复
- ELK的电子局域函数
已经有21人回复
- PRSV方程中的k
已经有4人回复
- 求助 ELK软件的学习！
已经有4人回复
- 大家看看我的excell拟合方程，怎么能合理的解释相关性呢？
已经有8人回复
- 提供一个 Birch-Murnaghan 拟合的算法（更新1）
已经有95人回复
- ELK 中的原子位置问题
已经有3人回复
- ELK发布1.4.5了
已经有17人回复
- 这种z=ax+by方程的拟合
已经有4人回复
- Brich-Murnaghan状态方程去拟合平衡时的晶格体积出问题！
已经有7人回复
- 第五类等温线用什么方程拟合
已经有3人回复
- 【原创】处理elk能带的程序
已经有16人回复
- 【求助】E-V曲线晶体状态方程拟合问题
已经有8人回复
- 【求助】f90子程序怎么编译错误？请大家帮忙，谢谢
已经有3人回复
- 【讨论】我收集到的状态方程编程
已经有7人回复
- 【其他】vasp相关技术有奖征集！
已经有35人回复
14楼2010-10-14 23:03:29
4楼2010-05-20 10:16:54
5楼2010-05-20 12:17:11
6楼2010-05-20 12:40:38
7楼2010-05-20 13:43:37
8楼2010-05-20 23:02:36
9楼2010-05-21 11:40:01
10楼2010-05-21 13:08:17
11楼2010-05-27 09:40:17
12楼2010-05-27 11:47:43
★★★★★ 五星级,优秀推荐
|
不知道大家一般选哪个状态方程？还有一个就是单位的问题，我一般直接选择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
15楼2010-10-15 09:13:49
16楼2010-10-15 10:10:41
|
请问下载自带的eos.in中参数是什么意思 "Silicon" : cname 2 : natoms 1 : etype 140.0 450.0 1000 : vplt1, vplt2, nvplt 7 : nevpt 165.8207473 -578.0660968 196.8383062 -578.1728409 231.5010189 -578.2305598 270.0113940 -578.2548962 312.5719400 -578.2566194 359.3851654 -578.2453281 410.6535788 -578.2253154 466.5796888 -578.2028836 第四行是什么意思，README里说volume interval over which to plot energy, pressure etc. as well as the number of points in the plot 看不懂请指教 还有就是从第六行开始的数据左边是体积右边是能量吗，如果是都用了什么单位呢。为什么体积变化这么大，右边能量变化却很小，谢谢了 |
17楼2010-10-15 10:17:48
18楼2010-10-18 10:38:39
19楼2010-11-25 21:19:55
20楼2010-11-25 22:54:55
21楼2010-11-26 10:06:10
22楼2011-01-11 16:30:01
23楼2011-05-25 17:18:03
24楼2011-05-25 17:20:14
25楼2011-08-15 14:31:50
27楼2011-09-27 09:54:46
liancsh28楼2011-09-27 10:32:11
29楼2011-10-05 13:19:12
★★★★★ 五星级,优秀推荐
|
以下两种结果与实验相差都很大，高手指点一下! [whz@localhost test]$ more PARAM.OUT Silicon Birch-Murnaghan 3rd-order EOS Birch F, Phys. Rev. 71, p809 (1947) (Default units are atomic: Hartree, Bohr etc.) V0 = 296.7512868 E0 = -578.2575138 B0 = 0.2486647064E-02 B0' = 4.117362248 B0 (GPa) = 73.15967018 [whz@localhost test]$ more eos.in "Silicon" : cname 2 : natoms 3 : etype 140.0 450.0 1000 : vplt1, vplt2, nvplt 8 : nevpt 165.8207473 -578.0660968 196.8383062 -578.1728409 231.5010189 -578.2305598 270.0113940 -578.2548962 312.5719400 -578.2566194 359.3851654 -578.2453281 410.6535788 -578.2253154 466.5796888 -578.2028836 [whz@localhost test]$ [whz@localhost test]$ more PARAM.OUT Silicon Universal EOS Vinet P et al., J. Phys.: Condens. Matter 1, p1941 (1989) (Default units are atomic: Hartree, Bohr etc.) V0 = 295.4576079 E0 = -578.2578353 B0 = 0.2591304888E-02 B0' = 4.291344871 B0 (GPa) = 76.23880915 [whz@localhost test]$ more eos.in "Silicon" : cname 2 : natoms 1 : etype 140.0 450.0 1000 : vplt1, vplt2, nvplt 8 : nevpt 165.8207473 -578.0660968 196.8383062 -578.1728409 231.5010189 -578.2305598 270.0113940 -578.2548962 312.5719400 -578.2566194 359.3851654 -578.2453281 410.6535788 -578.2253154 466.5796888 -578.2028836 [whz@localhost test]$ |
30楼2011-10-27 21:56:05
34楼2012-02-15 04:41:56
37楼2012-03-25 00:46:13
39楼2012-05-08 21:24:59
40楼2012-05-12 11:18:31
41楼2012-10-16 20:41:52
42楼2012-10-27 18:57:39
简单回复
2010-05-16 19:27
回复
感谢分享
kgdu3楼
2010-05-16 20:11
回复
谢谢分享
uuv201026楼
2011-08-15 19:54
回复
五星好评 
liancsh31楼
2011-10-27 23:40
回复
五星好评 顶一下，感谢分享！
kgdu32楼
2011-12-04 13:31
回复
顶一下，感谢分享！
qqpie33楼
2011-12-17 03:03
回复
五星好评 顶一下，感谢分享！
fanarsenal35楼
2012-02-15 14:34
回复
五星好评 顶一下，感谢分享！
tider36楼
2012-03-05 11:39
回复
五星好评 顶一下，感谢分享！
seventrees38楼
2012-03-30 09:21
回复
五星好评 
bigcontinent43楼
2013-06-24 10:03
回复
五星好评 顶一下，感谢分享！
heju44楼
2013-08-28 15:24
回复
五星好评 顶一下，感谢分享！
tonglijia45楼
2013-09-29 19:41
回复
五星好评 顶一下，感谢分享！
swuhaizhilan46楼
2013-10-16 19:53
回复
五星好评 顶一下，感谢分享！
tsin1847楼
2013-10-24 18:48
回复
五星好评 顶一下，感谢分享！
superstat5448楼
2013-10-30 22:24
回复
五星好评 顶一下，感谢分享！
白菜胡萝卜49楼
2013-11-01 17:20
回复
五星好评 顶一下，感谢分享！
excellencer50楼
2013-11-05 23:11
回复
五星好评 顶一下，感谢分享！
10