24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3271  |  回复: 14

qiqi2926

铜虫 (小有名气)

[交流] 【求助】拟合晶格常数已有10人参与

以下是我优化结构得到的EtVo.dat的数据:(第一行是晶格常数,第二行是体积,第三行是能量):
4.076 95.2900 -48.176301
4.126 98.8400 -48.233218
4.176 102.4800 -48.280047
4.226 106.2000 -48.314084
4.276 110.0200 -48.331987
4.326 113.9200 -48.336034

4.376 117.9200 -48.326447
4.426 122.0100 -48.298049
4.476 126.1900 -48.272207
4.526 130.4700 -48.231560
可以看出能量最低点出现在晶格常数为4.276~4.326之间。而通过拟合得到的晶格常数却远离这个范围。。。。这又是为什么?
顺便求拟合验证。。。!!!谢谢。。。
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

bingmou

金虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+1):谢谢指教 2010-08-09 21:45:41
我这里有一个LINUX下面的程序,是用Murn...状态方程拟合的。自己编译以下,很容易,而且拟合很好。
这里是下载地址:http://good.gd/647372.htm

[ Last edited by bingmou on 2010-8-9 at 11:29 ]
9楼2010-08-09 11:28:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

fanchen021

铁杆木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
cenwanglai(金币+1):hehe 2010-08-01 20:05:45
注意单位。。
2楼2010-08-01 19:19:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qiqi2926

铜虫 (小有名气)

引用回帖:
Originally posted by fanchen021 at 2010-08-01 19:19:13:
注意单位。。

单位是统一的
3楼2010-08-04 11:26:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

渔火江枫

金虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
我也遇到过相似的情况,你用的是Primitive cell吗?
IdeasWorthSpreading!
4楼2010-08-04 12:47:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

condensed

木虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
aylayl08(金币+1):谢谢 2010-08-04 16:39:16
拟合出来的显然是体积.
5楼2010-08-04 13:37:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小康康

铁杆木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
cenwanglai(金币+1):谢谢回复! 2010-08-06 11:35:54
你可以用pwscf去拟合晶格常数啊
6楼2010-08-05 22:55:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snail594

金虫 (小有名气)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zxzj05(金币+1):谢谢回帖交流 2010-08-06 15:26:23
同样用vasp计算拟合没有出现问题,你可以统一采用单位原子体积和单位原子能量,避免一些不一致错误,还有单位,最后可能是公式写错了,我是自己写代码拟合的
Inlearningwetrust
7楼2010-08-06 15:09:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qiqi2926

铜虫 (小有名气)

引用回帖:
Originally posted by snail594 at 2010-08-06 15:09:50:
同样用vasp计算拟合没有出现问题,你可以统一采用单位原子体积和单位原子能量,避免一些不一致错误,还有单位,最后可能是公式写错了,我是自己写代码拟合的

我是把状态方程利用关系式化简了,然后利用MATLAB进行拟合的。。。能不能把你的代码传给我一份。。。我自己试试。。。谢谢。。。qiqi2926@163.com
8楼2010-08-08 21:32:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qiqi2926

铜虫 (小有名气)

引用回帖:
Originally posted by bingmou at 2010-08-09 11:28:08:
我这里有一个LINUX下面的程序,是用Murn...状态方程拟合的。自己编译以下,很容易,而且拟合很好。
这里是下载地址:http://good.gd/647372.htm

[ Last edited by bingmou on 2010-8-9 at 11:29  ...

貌似对立方晶系很好,但是对六角的不太好。。。
10楼2010-08-09 18:35:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 qiqi2926 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见