24小时热门版块排行榜    

查看: 5438  |  回复: 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的回帖

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的回帖
查看全部 96 个回答

fanarsenal

木虫 (正式写手)


★★★ 三星级,支持鼓励

问一下楼主这个单位是什么的》
12楼2012-03-07 21:48:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
huangll992楼
2012-03-06 09:19   回复  
五星好评  顶一下,感谢分享!
2012-03-06 09:43   回复  
五星好评  顶一下,感谢分享!
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   回复  
五星好评  顶一下,感谢分享!
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见