24小时热门版块排行榜    

查看: 4484  |  回复: 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 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 延安大学回应“副教授留绝笔信”:确已离世,正核实举报内容 +3 yexuqing 2026-11-04 5/250 2026-11-07 13:22 by 苏东坡二世
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 23/1150 2026-11-07 06:49 by Equinoxhua
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +10 ncTVBM6AOJHF 2026-11-04 22/1100 2026-11-06 17:18 by 5sJFDKp2kwuI
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 16:43 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 16:23 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ncTVBM6AOJHF 2026-11-04 17/850 2026-11-06 13:47 by BkA9Bi5cMhkD
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +7 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 13:47 by BkA9Bi5cMhkD
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 13:47 by BkA9Bi5cMhkD
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +8 ncTVBM6AOJHF 2026-11-04 18/900 2026-11-06 13:36 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 13:36 by BkA9Bi5cMhkD
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-05 10/500 2026-11-06 12:32 by BkA9Bi5cMhkD
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 09:01 by 2lRNluKrJpwL
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 05:32 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 12/600 2026-11-06 05:32 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 11/550 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
[基金申请] 各位基础老师们,看完中央《从严治党》的决定后,对学术界更加有信心了吧 +5 iaeyuan 2026-11-01 5/250 2026-11-04 14:43 by jinchengzi
信息提示
请填处理意见