24小时热门版块排行榜    

查看: 5431  |  回复: 95
【奖励】 本帖被评价85次,作者freewain增加金币 66.8

freewain

铜虫 (初入文坛)


[资源] 提供一个 Birch-Murnaghan 拟合的算法(更新1)

提供一个计算 Birch-Murnaghan EOS 的思路和 MatLab 算法,请大家指正。更新:
最近发现,提供的这个算法并不太好,虽然原理上没有问题,但是在数值算法中有时会产生误差。造成这个的原因在于中间涉及的一些数太小,或者说需要的精度太高。当我试图用 C 或者 Fortran 实现这种算法时,发现连 double 精度都无法保证正确结果。
现在看来,最好的方法还是用 least squares 方法直接拟合 Birch-Murnaghan 方程。我最近会试着用 Fortran 写一下这个功能,如果成功,会更新这个帖子放到这里。
很抱歉之前的不成熟的算法耽误大家的时间了。
另外,这里有一个拟合 Murnaghan 方程的 Fortran 代码,作者是 SIESTA 的维护者之一 Postnikov。http://www.home.uni-osnabrueck.de/apostnik/Software/mur_fit.f

[ Last edited by freewain on 2013-4-22 at 21:21 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : Birch-Murnaghan Regression.pdf
  • 2012-03-06 05:36:43, 94.58 K

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

@个人收集@ 第一性原理计算-tg Photochemistry 分子动力学
表界面计算

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有2个 )

nkhjh

金虫 (初入文坛)


lz,我的运行结果总是出现错误,要么是精度不够,要么是函数出现问题,输出结果就只有B‘=4,其余的没有值。是能量太大了?还是什么原因?谢谢您!
16楼2012-03-13 09:23:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xmqlzf

银虫 (小有名气)


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

不错,受教了!感谢分享
20楼2012-04-17 21:06:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanarsenal

木虫 (正式写手)


★★★ 三星级,支持鼓励

问一下楼主这个单位是什么的》
12楼2012-03-07 21:48:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

freewain

铜虫 (初入文坛)


引用回帖:
12楼: Originally posted by fanarsenal at 2012-03-07 21:48:00:
问一下楼主这个单位是什么的》

能量是eV; 晶格常熟是Angstrom.
13楼2012-03-08 03:21:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kuiziyang

银虫 (小有名气)


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

顶一下,感谢分享!虽然还没有明确怎么用
25楼2012-06-17 10:34:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

harry_zhou00

银虫 (初入文坛)


楼主 请教你一下, 如果这个是Hcp的情况,应该怎么拟合呢?
39楼2012-12-05 18:36:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

haitian123

银虫 (小有名气)


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

不错哦
46楼2013-03-15 22:59:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

claudshone

新虫 (初入文坛)


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

琢磨了一个下午,实在无解,感谢分享!
47楼2013-04-20 16:45:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
huangll992楼
2012-03-06 09:19   回复  
五星好评  顶一下,感谢分享!
2012-03-06 09:43   回复  
五星好评  顶一下,感谢分享!
2012-03-06 09:47   回复  
五星好评  
cheerian5楼
2012-03-06 11:00   回复  
五星好评  
乐石||6楼
2012-03-06 16:45   回复  
五星好评  顶一下,感谢分享!
nkhjh7楼
2012-03-07 09:25   回复  
五星好评  顶一下,感谢分享!
zhxd10288楼
2012-03-07 14:08   回复  
五星好评  顶一下,感谢分享!
thinkllh9楼
2012-03-07 15:32   回复  
五星好评  顶一下,感谢分享!
luxia021010楼
2012-03-07 15:50   回复  
五星好评  顶一下,感谢分享!
tonyngch11楼
2012-03-07 16:14   回复  
三星好评  顶一下,感谢分享!
2012-03-12 16:53   回复  
五星好评  顶一下,感谢分享!
jswoo15楼
2012-03-12 19:21   回复  
五星好评  顶一下,感谢分享!
若人生17楼
2012-04-01 21:51   回复  
五星好评  
zhuzhudgyg18楼
2012-04-12 15:08   回复  
五星好评  顶一下,感谢分享!
liuziran19楼
2012-04-12 15:34   回复  
五星好评  
xzyme21楼
2012-04-17 21:46   回复  
五星好评  顶一下,感谢分享!
2012-05-16 08:44   回复  
五星好评  顶一下,感谢分享!
bjfang23楼
2012-05-18 17:04   回复  
五星好评  顶一下,感谢分享!
enola24楼
2012-05-19 09:30   回复  
三星好评  顶一下,感谢分享!
allendean26楼
2012-07-27 19:38   回复  
三星好评  顶一下,感谢分享!
2012-07-28 16:07   回复  
五星好评  顶一下,感谢分享!
albertmarx28楼
2012-07-31 15:18   回复  
五星好评  顶一下,感谢分享!
caesarmu29楼
2012-08-05 15:40   回复  
三星好评  顶一下,感谢分享!
ljyang32030楼
2012-09-01 16:56   回复  
五星好评  顶一下,感谢分享!
lvfei11131楼
2012-10-07 21:17   回复  
五星好评  顶一下,感谢分享!
cs12329232楼
2012-10-12 12:06   回复  
五星好评  顶一下,感谢分享!
losenq33楼
2012-10-15 12:40   回复  
五星好评  顶一下,感谢分享!
thinkllh34楼
2012-10-15 14:58   回复  
顶一下,感谢分享!
2012-10-15 22:59   回复  
五星好评  顶一下,感谢分享!
dragon506136楼
2012-11-16 16:12   回复  
五星好评  顶一下,感谢分享!
nancy_sun37楼
2012-12-05 10:59   回复  
顶一下,感谢分享!
2012-12-05 18:33   回复  
五星好评  顶一下,感谢分享!
ZHC_Y40楼
2012-12-25 09:33   回复  
五星好评  顶一下,感谢分享!
guodonglin41楼
2012-12-31 20:26   回复  
五星好评  顶一下,感谢分享!
liangpei42楼
2013-01-02 14:22   回复  
五星好评  顶一下,感谢分享!
shuamfs43楼
2013-02-03 17:50   回复  
五星好评  顶一下,感谢分享!
fsfirst44楼
2013-02-04 20:05   回复  
五星好评  顶一下,感谢分享!
tzl72245楼
2013-03-15 22:38   回复  
五星好评  顶一下,感谢分享!
huichli48楼
2013-05-17 22:21   回复  
五星好评  顶一下,感谢分享!
guodonglin49楼
2013-05-18 12:41   回复  
顶一下,感谢分享!
kervinzhao50楼
2013-05-28 22:59   回复  
五星好评  顶一下,感谢分享!
相关版块跳转 我要订阅楼主 freewain 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见