应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » origin中如何自定义Birch–Murnaghan三阶拟合
411/1返回列表
查看: 10285  |  回复: 40
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【奖励】 本帖被评价14次，作者tonglijia增加金币 10.4 个

tonglijia

木虫 (正式写手)

[资源] origin中如何自定义Birch–Murnaghan三阶拟合

有的虫友问我在origin中如何实现b-m三阶拟合的，所以贴出这个资源供大家参考
1.tools--fitting functions organizer e9
2.具体设置见图片
3.应该注意的地方：在编辑公式的时候极容易出错，所以个人建议在编辑的过程中先将functions form设置为 origin c格式，然后在编辑函数框中可以每输入一个字符，点击一下输入框后面的那个红色斜向上的按钮（函数形式为equation时不显示此按钮），点击后的图片为第二张图，此时点击compie，如果显示done!那么公式编辑是成功的，如果失败请返回查看。如果是每编辑一个字符就检查对错的话一定要注意parameters names中出现的需要拟合的量必须要与公式中当前包含的需要拟合的量一致！否则报错

b-m.jpg


-1
未命名.jpg

关于单位换算问题
1 bohr =5.29177210818 x 10^(-11) m = 1 au = 0.52917埃; 1 hartree = 2625.5 kj/mol =27.2114 ev;
1 hartree = 4.3597441775 x 10^(-18) j = 1 ha
1 ha/au^3 = [4.3597441775 x 10^(-18) j ] / [5.29177210818^3 x 10^(-11) m^3]
                  = 2.94210120348 x 10^(13) j/m^3
                  = 2.94210120348 x 10^(13) n/m^2
                  = 2.94210120348 x 10^4 gpa

[ last edited by tonglijia on 2015-10-16 at 17:27 ]

遇到问题的同学，请自行联系我422248685（QQ微信都是这个）

[ Last edited by tonglijia on 2021-11-4 at 07:29 ]
回复此楼

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

MS资源 培训、总结 长远发展 第一性原理

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-10-02 19:45:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

tonglijia

木虫 (正式写手)
体积能量拟合晶格参数的时候一定要注意量纲！需要了解量纲转换的请留言，有时间会写出详细的转换过程
赞一下(3人) 回复此楼
2楼2013-10-02 19:49:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

may-sucess

金虫 (正式写手)

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

很好，学习了，但是我的自编函数也编译通过了，但是拟合时，没有想要的参数结果。图片如下：
-2
aaaaa.png
赞一下(2人) 回复此楼
13楼2014-07-07 09:31:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

奈落之夜

木虫 (正式写手)

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

引用回帖:
2楼: Originally posted by tonglijia at 2013-10-02 19:49:37
体积能量拟合晶格参数的时候一定要注意量纲！需要了解量纲转换的请留言，有时间会写出详细的转换过程

我需要啊 ，楼主可否详细讲一下
赞一下(1人) 回复此楼
7楼2013-11-26 10:52:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tonglijia

木虫 (正式写手)
引用回帖:
7楼: Originally posted by 奈落之夜 at 2013-11-26 10:52:41
我需要啊 ，楼主可否详细讲一下...

最近比较忙，你先试着自己推导下，如果有不知道的地方再问我。也可以把你的推导过程贴出来，有时间我帮你看下
赞一下(1人) 回复此楼
8楼2013-11-26 16:03:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

k8wyy

铁虫 (小有名气)

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

你好，请问您这个EOS的拟合的输入数据是E和V吗？初始值E0和V0以及B0和B0'是怎么设置的呢？
赞一下(1人) 回复此楼
25楼2015-05-17 21:25:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

怒冲星河

木虫 (正式写手)

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

好东西，先顶一下
回复此楼
3楼2013-10-03 11:48:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

乐石||

木虫 (正式写手)

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

留名
回复此楼
6楼2013-10-04 11:02:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

超级月亮

木虫 (著名写手)

★★★ 三星级,支持鼓励

公式里边儿是B'，不是B，我现在要做二阶拟合，可不可以把E看做E[(V0/V)^2/3-1],这样，式子E=E0+9/8*B0*V0*[(V0/V)^2/3-1]^2用origin拟合出的二次项系数就是等于9/8*B0*V0,从而求得体积模量B0
赞一下 回复此楼
9楼2014-04-14 16:15:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hailan2005

银虫 (小有名气)

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

非常感谢楼主，按您的图示，已拟合成功，不过有一个小问题恳请您帮助看看。就是那个B‘的没有拟合出来，还是原来的初始值4。
-3
P-V.png
赞一下 回复此楼
10楼2014-04-30 11:23:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sungjen

银虫 (著名写手)
请问二阶拟合和三阶有什么区别？这种拟合方法是否可以拟合出体弹性模量 ？急求。
赞一下 回复此楼
11楼2014-06-03 10:35:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

may-sucess

金虫 (正式写手)
很好，学习了，但是我的自编函数也编译通过了，但是拟合时，没有想要的参数结果。
赞一下 回复此楼
12楼2014-07-04 22:20:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

may-sucess

金虫 (正式写手)
已解决！
回复此楼
14楼2014-07-07 15:54:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qianyi1101

新虫 (小有名气)
引用回帖:
14楼: Originally posted by may-sucess at 2014-07-07 15:54:07
已解决！

初始值怎么设置啊？我也没拟合出来结果呢
赞一下 回复此楼
15楼2014-08-05 15:26:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tonglijia

木虫 (正式写手)
引用回帖:
7楼: Originally posted by 奈落之夜 at 2013-11-26 10:52:41
我需要啊 ，楼主可否详细讲一下...

很久没看这个帖子了，没想到还有这么多回复，关于量纲的推导年前一定做出来，不好意思太旧没有回复了！抱歉
赞一下 回复此楼
16楼2014-12-17 21:28:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tonglijia

木虫 (正式写手)
引用回帖:
9楼: Originally posted by 超级月亮 at 2014-04-14 16:15:31
公式里边儿是B'，不是B，我现在要做二阶拟合，可不可以把E看做E,这样，式子E=E0+9/8*B0*V0*^2用origin拟合出的二次项系数就是等于9/8*B0*V0,从而求得体积模量B0

你的比较简洁，发个帖子吧，让大家也受益一下
赞一下 回复此楼
17楼2014-12-17 21:29:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tonglijia

木虫 (正式写手)
引用回帖:
10楼: Originally posted by hailan2005 at 2014-04-30 11:23:11
非常感谢楼主，按您的图示，已拟合成功，不过有一个小问题恳请您帮助看看。就是那个B‘的没有拟合出来，还是原来的初始值4。

P-V.png
...

很久没看这个帖子了，没想到还有这么多回复，不好意思太旧没有回复了！抱歉
从图像看你这个拟合是不对的， 你换一个其它的初始值试一下
赞一下 回复此楼
18楼2014-12-17 21:30:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tonglijia

木虫 (正式写手)
引用回帖:
11楼: Originally posted by sungjen at 2014-06-03 10:35:24
请问二阶拟合和三阶有什么区别？这种拟合方法是否可以拟合出体弹性模量 ？急求。

可以。二阶和三阶的精度不同
赞一下 回复此楼
19楼2014-12-17 21:32:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tonglijia

木虫 (正式写手)
引用回帖:
15楼: Originally posted by qianyi1101 at 2014-08-05 15:26:41
初始值怎么设置啊？我也没拟合出来结果呢...

初始值一般不用太在意，大概有一个值就行，然后用第一次拟合出来的参数初始值再此拟合
赞一下 回复此楼
20楼2014-12-17 21:33:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tsglss

铁杆木虫 (著名写手)

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

请问，非线性拟合过程中怎么输入导数或积分的函数？还有体积能量拟合晶格参数的时候量纲问题怎么解决？谢谢
赞一下 回复此楼
21楼2014-12-18 12:25:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小熙小7

银虫 (小有名气)
引用回帖:
14楼: Originally posted by may-sucess at 2014-07-07 15:54:07
已解决！

你好  请问你是怎么解决的呢 我也遇到了这个问题 谢谢啦
赞一下 回复此楼
23楼2015-03-03 17:51:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小熙小7

银虫 (小有名气)
引用回帖:
14楼: Originally posted by may-sucess at 2014-07-07 15:54:07
已解决！

你好  请问你是怎么解决的呢 我也遇到了这个问题 谢谢啦
赞一下 回复此楼
24楼2015-03-04 11:26:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hearuing

新虫 (初入文坛)

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

引用回帖:
24楼: Originally posted by 小熙小7 at 2015-03-04 11:26:18
你好  请问你是怎么解决的呢 我也遇到了这个问题 谢谢啦...

你好，我遇到的问题和你一样，想问一下你是怎么解决的，谢谢了！
赞一下 回复此楼
26楼2015-09-24 14:23:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tonglijia

木虫 (正式写手)
引用回帖:
7楼: Originally posted by 奈落之夜 at 2013-11-26 10:52:41
我需要啊 ，楼主可否详细讲一下...

http://muchong.com/bbs/viewthread.php?tid=9502795&target=1
赞一下 回复此楼
27楼2015-10-16 17:22:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tonglijia

木虫 (正式写手)
引用回帖:
21楼: Originally posted by tsglss at 2014-12-18 12:25:23
请问，非线性拟合过程中怎么输入导数或积分的函数？还有体积能量拟合晶格参数的时候量纲问题怎么解决？谢谢

http://muchong.com/bbs/viewthread.php?tid=9502795&target=1
赞一下 回复此楼
28楼2015-10-16 17:22:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tonglijia

木虫 (正式写手)
引用回帖:
23楼: Originally posted by 小熙小7 at 2015-03-03 17:51:33
你好  请问你是怎么解决的呢 我也遇到了这个问题 谢谢啦...

http://muchong.com/bbs/viewthread.php?tid=9502795&target=1
赞一下 回复此楼
29楼2015-10-16 17:22:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tonglijia

木虫 (正式写手)
引用回帖:
26楼: Originally posted by hearuing at 2015-09-24 14:23:28
你好，我遇到的问题和你一样，想问一下你是怎么解决的，谢谢了！...

http://muchong.com/bbs/viewthread.php?tid=9502795&target=1
赞一下 回复此楼
30楼2015-10-16 17:23:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tonglijia

木虫 (正式写手)
引用回帖:
25楼: Originally posted by k8wyy at 2015-05-17 21:25:09
你好，请问您这个EOS的拟合的输入数据是E和V吗？初始值E0和V0以及B0和B0'是怎么设置的呢？

Ｅ０，Ｖ０根据计算数值输入，Ｂ０以一个大概的值为参考就可以了
赞一下 回复此楼
31楼2015-10-16 17:24:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

矢车菊00

银虫 (正式写手)

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

引用回帖:
10楼: Originally posted by hailan2005 at 2014-04-30 11:23:11
非常感谢楼主，按您的图示，已拟合成功，不过有一个小问题恳请您帮助看看。就是那个B‘的没有拟合出来，还是原来的初始值4。

P-V.png
...

学习了，自己试了一遍，怎么没有出现拟合的曲线，跪求答案
-4
QQ截图20151127182216.png
赞一下 回复此楼
32楼2015-11-27 18:25:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

采采芣苡

木虫 (小有名气)

★ 一星级,一般

引用回帖:
2楼: Originally posted by tonglijia at 2013-10-02 19:49:37
体积能量拟合晶格参数的时候一定要注意量纲！需要了解量纲转换的请留言，有时间会写出详细的转换过程

您好！能不能把详细的过程再发一次，再分享一个能量转换的详细过程呢？菜鸟，急需，谢谢～

发自小木虫Android客户端
赞一下 回复此楼
33楼2017-05-01 22:54:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

冰淇淋吖

铁虫 (初入文坛)
引用回帖:
10楼: Originally posted by hailan2005 at 2014-04-30 11:23:11
非常感谢楼主，按您的图示，已拟合成功，不过有一个小问题恳请您帮助看看。就是那个B‘的没有拟合出来，还是原来的初始值4。

P-V.png
...

请问B0V0初始值在哪里设置啊
赞一下 回复此楼
34楼2018-09-03 10:48:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tonglijia

木虫 (正式写手)
引用回帖:
34楼: Originally posted by 冰淇淋吖 at 2018-09-03 10:48:08
请问B0V0初始值在哪里设置啊...

这个很好解决，简单查一下晶体库中的晶体性质就可以了
赞一下 回复此楼
35楼2021-11-04 07:16:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tonglijia

木虫 (正式写手)
查看私信
回复此楼
36楼2021-11-04 07:25:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tonglijia

木虫 (正式写手)
查看私信
回复此楼
37楼2021-11-04 07:25:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tonglijia

木虫 (正式写手)
引用回帖:
26楼: originally posted by hearuing at 2015-09-24 14:23:28
你好，我遇到的问题和你一样，想问一下你是怎么解决的，谢谢了！...

查看私信
回复此楼
38楼2021-11-04 07:25:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaonahappy

金虫 (小有名气)
引用回帖:
34楼: Originally posted by 冰淇淋吖 at 2018-09-03 10:48:08
请问B0V0初始值在哪里设置啊...

请问这个B0VO在origin里面怎么设置进去？
赞一下 回复此楼
39楼2022-01-05 14:59:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaonahappy

金虫 (小有名气)
引用回帖:
12楼: Originally posted by may-sucess at 2014-07-04 22:20:39
很好，学习了，但是我的自编函数也编译通过了，但是拟合时，没有想要的参数结果。

您好，请问您的问题解决了吗？我拟合的时候也没有得到想要的参数结果。能否帮忙说一下如何解决？
赞一下 回复此楼
40楼2022-01-05 15:34:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaonahappy

金虫 (小有名气)
内容已删除
回复此楼
41楼2022-01-05 15:36:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
zuocuiping4楼
2013-10-03 13:22   回复  
三星好评  
loveting5楼
2013-10-03 15:38   回复  
五星好评  
黄河honour22楼
2014-12-19 07:24   回复  
五星好评  顶一下先 [ 发自小木虫客户端 ]
相关版块跳转 我要订阅楼主 tonglijia 的主题更新
411/1返回列表
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭