24小时热门版块排行榜    

查看: 4092  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +32 医学老男孩 2026-08-13 74/3700 2026-08-16 21:05 by 医学老男孩
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
请填处理意见