24小时热门版块排行榜    

查看: 2828  |  回复: 12

afujie

木虫 (正式写手)

[求助] Gulp 拟合势函数 初值设置的问题

一、看了大家的很多的帖子 都说gulp拟合势函数初值的设置 对结果很重要 需要一定的经验 大家可以谈谈自己初值设置的经验吗?
其次 每次拟合完之后 需要进行迭代 是什么意思?是把输出文件中推荐的值带到gin文件中进行再次拟合吗?可是这样的一组参数值(上一次fit成功之后的推荐值)程序默认是拟合好的 不会再进行拟合啊 各位能不能帮我分析一下 问题出在哪?
回复此楼

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

动力学强帖

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

老虎大王

木虫 (著名写手)

引用回帖:
9楼: Originally posted by afujie at 2012-04-25 20:29:18:
大王 想请教您一个问题
gulp拟合势函数 截断一般怎样选去呢?
比如说buck函数,我见例子中一般都是8,12,或者24的选,可是根据理论,一般都是算到第一到第三近邻,后面的作用就很小了,也就是5埃左右,为什么例 ...

离一个中心原子比较远的地方,涉及的配位原子也比较多。这里面就产生两个问题,一个问题是,虽然函数值在远的地方很小,但是因为涉及原子多,在总能量上积累起来的误差有可能还是很可观的,尤其是对于非晶和液体来说。另一个问题是,因为越远涉及原子越多,所以相互作用计算就越费时间。这是一对矛盾。但现在的计算机性能提高了,所以第二个问题的严重性就不是那么大了。选大一些,第一可以保证精确性,第二可以保证截断的平滑性。

对这个问题你要做一个测试,不是光测试势数值,而是要测试它算出来的整个系统的能量。
10楼2012-04-26 14:34:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
zhangguangping: 金币+2, 谢谢指教! 2012-04-09 16:57:33
afujie: 金币+5, 有帮助 2012-04-12 09:12:12
拟合初值的确定,一靠经验二靠运气,三建议试试1stOpt,采用的是全局优化算法,不依赖初值。

» 本帖已获得的红花(最新10朵)

2楼2012-04-09 09:12:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

afujie

木虫 (正式写手)

送鲜花一朵
引用回帖:
2楼: Originally posted by dingd at 2012-04-09 09:12:14:
拟合初值的确定,一靠经验二靠运气,三建议试试1stOpt,采用的是全局优化算法,不依赖初值。

在经验这方面 都应该注意什么呢?
进行全局优化 会给出势函数的参数吗?这个我没有用过 能说详细点吗
3楼2012-04-09 09:18:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
afujie: 金币+5, ★★★很有帮助 2012-04-15 22:27:24
zhangguangping: 金币+2, 谢谢指导! 2012-04-16 10:06:17
这个没有好的办法啊。主要靠运气。另外是靠大量的试验。经验嘛就是从大量试验中总结来的。要观察每个参数的变动所引起的结果的变化,以及某几个参数之间的协调关系。这个不是能说出来的,全靠自己灵活掌握。你试的次数多了,再加上运气好,就能拟合好;如果你试的次数少,运气又不好,就做不出来;你试的少,但运气好,也有可能出来;如果你试得多,但运气不好,还是出不来。
4楼2012-04-15 21:57:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)

【答案】应助回帖

★ ★
zhangguangping: 金币+2, 谢谢建议! 2012-04-16 10:06:26
另外我有一个很具体的建议,就是建议你多发动几个人和你一起做,大家分头做,然后讨论自己做的过程中有哪些发现。这个GULP,反正输入文件写好之后,谁都可以去做。你找几个本科生,专科生,高中生,都可以。有的人,运气还是相当好的。
5楼2012-04-15 22:00:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

afujie

木虫 (正式写手)

引用回帖:
5楼: Originally posted by 老虎大王 at 2012-04-15 22:00:38:
另外我有一个很具体的建议,就是建议你多发动几个人和你一起做,大家分头做,然后讨论自己做的过程中有哪些发现。这个GULP,反正输入文件写好之后,谁都可以去做。你找几个本科生,专科生,高中生,都可以。有的人 ...

大王 那你说是不是拟合出来的数值 跟实验值的误差不能超过多少呢?5%?C44 和剪切模量都能符合的很好吗?
6楼2012-04-15 22:29:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)

★ ★
dbb627: 金币+2, 谢谢应助! 2012-04-26 14:06:20
引用回帖:
6楼: Originally posted by afujie at 2012-04-15 22:29:06:
大王 那你说是不是拟合出来的数值 跟实验值的误差不能超过多少呢?5%?C44 和剪切模量都能符合的很好吗?

误差多少算比较好呢?一般来说,能量和结构应该尽可能接近,我看文献上能量的优化值的误差基本上都是1%以内的,结构(晶格常数)也差不多是这样吧,也有误差稍微大一些的。这个不光要看相对误差,也要看绝对误差,总之感觉上要很接近,人家才会接受。各种模量,误差稍微大一些没有很大关系,毕竟不同的实验值之间也有不小的差别,你选定一组实验值,你的计算值和它差5%,或者10%以内我觉得都可以接受。比如100Gpa,你优化后出来的是90GPa,这也算是比较接近了。

» 本帖已获得的红花(最新10朵)

7楼2012-04-16 10:08:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

afujie

木虫 (正式写手)

送鲜花一朵
引用回帖:
7楼: Originally posted by 老虎大王 at 2012-04-16 10:08:14:
误差多少算比较好呢?一般来说,能量和结构应该尽可能接近,我看文献上能量的优化值的误差基本上都是1%以内的,结构(晶格常数)也差不多是这样吧,也有误差稍微大一些的。这个不光要看相对误差,也要看绝对误差 ...

非常感谢 大王的指导
8楼2012-04-16 10:50:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

afujie

木虫 (正式写手)

引用回帖:
4楼: Originally posted by 老虎大王 at 2012-04-15 21:57:02:
这个没有好的办法啊。主要靠运气。另外是靠大量的试验。经验嘛就是从大量试验中总结来的。要观察每个参数的变动所引起的结果的变化,以及某几个参数之间的协调关系。这个不是能说出来的,全靠自己灵活掌握。你试的 ...

大王 想请教您一个问题
gulp拟合势函数 截断一般怎样选去呢?
比如说buck函数,我见例子中一般都是8,12,或者24的选,可是根据理论,一般都是算到第一到第三近邻,后面的作用就很小了,也就是5埃左右,为什么例子中要选那么大呢?
9楼2012-04-25 20:29:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 afujie 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 8rmuugja8q 2026-02-22 8/400 2026-02-23 12:22 by alian_214
[基金申请] 体制内长辈说体制内绝大部分一辈子在底层,如同你们一样大部分普通教师忙且收入低 +10 瞬息宇宙 2026-02-20 13/650 2026-02-23 11:23 by holypower
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 khieu8v8m0 2026-02-22 8/400 2026-02-23 09:35 by w4l55oybr1
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 8/400 2026-02-23 09:29 by w4l55oybr1
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +5 usprnugpzw 2026-02-21 11/550 2026-02-23 09:24 by w4l55oybr1
[教师之家] 为什么中国大学工科教授们水了那么多所谓的顶会顶刊,但还是做不出宇树机器人? +5 欢乐颂叶蓁 2026-02-21 8/400 2026-02-23 09:19 by 欢乐颂叶蓁
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 w89i99eaeh 2026-02-22 5/250 2026-02-23 08:04 by w4l55oybr1
[博后之家] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 6/300 2026-02-23 07:59 by w4l55oybr1
[博后之家] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +6 3dfhjxgsh7 2026-02-22 9/450 2026-02-23 07:49 by w4l55oybr1
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 4/200 2026-02-23 06:46 by jsjzfl
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 khieu8v8m0 2026-02-22 5/250 2026-02-23 06:29 by w4l55oybr1
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 8/400 2026-02-23 06:24 by w4l55oybr1
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +5 3dfhjxgsh7 2026-02-22 6/300 2026-02-23 02:04 by 5jlh3qtdvx
[教师之家] 版面费该交吗 +7 苹果在哪里 2026-02-22 8/400 2026-02-22 22:37 by otani
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 usprnugpzw 2026-02-21 6/300 2026-02-22 19:48 by w89i99eaeh
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 usprnugpzw 2026-02-22 3/150 2026-02-22 16:37 by khieu8v8m0
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 usprnugpzw 2026-02-21 4/200 2026-02-22 16:27 by khieu8v8m0
[基金申请] “人文社科而论,许多学术研究还没有达到民国时期的水平” +4 苏东坡二世 2026-02-18 5/250 2026-02-22 16:07 by liangep1573
[基金申请] 什么是人一生最重要的? +4 瞬息宇宙 2026-02-21 4/200 2026-02-22 11:44 by huagongfeihu
[基金申请] 今年春晚有几个节目很不错,点赞! +11 瞬息宇宙 2026-02-16 12/600 2026-02-21 21:14 by lq493392203
信息提示
请填处理意见