24小时热门版块排行榜    

查看: 2133  |  回复: 5

csu_lxc

金虫 (正式写手)

[求助] 求助,高手帮帮忙,关于势参数的问题

我最近看到一篇,其中有个位置不明白,就是关于TB-SMA势函数。
势函数是这样



我也明白如何利用势函数计算晶格常数、弹性模量等等
然后,文献里面写的是“首先,给定势参数的初始值;然后,利用该多体势计算出系统的物理性质,如晶格常数、弹性模量等等,并与实验值进行比较;接着,再调整势参数的值,重复上述过程直至得到满意的结果。多体势的拟合是一个非常耗时的实验校正过程。”
我想得到势函数里面的那些势参数,不知道怎么做,怎么入手
谁能告诉我具体步骤是怎么样的?要用到什么软件吗?
谢谢了,非常感谢,纠结了很久,我很菜,刚学,第一次接触这种问题。谢谢指教!!

[ Last edited by csu_lxc on 2011-12-26 at 21:30 ]
回复此楼

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

动力学强帖

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chaizhm

木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
csu_lxc(金币+2): ★有帮助 谢谢回复 2011-12-27 08:37:37
我说说我的理解啊

你看的那篇文献已经说了,是设一个初始值,然后进行迭代,直到满足结果。显然你这种问题属于非线性迭代的问题,非线性的问题不一定有唯一解,而且很可能不收敛,所以选取初始值很重要,这个一般是得要经验。如果用初始值求得的参数与已知不符合,就要用些算法来调整初始值,最常用的非线性问题算法有最小梯度法,共轭梯度法等等,都可以尝试一下。得到下一组值做为初始值,知道迭代收敛为止
壁立万仞,无欲则刚
2楼2011-12-26 21:43:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chaizhm

木虫 (著名写手)

【答案】应助回帖

我说说我的理解啊

你看的那篇文献已经说了,是设一个初始值,然后进行迭代,直到满足结果。显然你这种问题属于非线性迭代的问题,非线性的问题不一定有唯一解,而且很可能不收敛,所以选取初始值很重要,这个一般是得要经验。如果用初始值求得的参数与已知不符合,就要用些算法来调整初始值,最常用的非线性问题算法有最小梯度法,共轭梯度法等等,都可以尝试一下。得到下一组值做为初始值,知道迭代收敛为止
壁立万仞,无欲则刚
3楼2011-12-26 21:43:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ussy06

银虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
csu_lxc(金币+8): ★★★很有帮助 2011-12-27 14:44:39
势函数的拟合是一个比较复杂的过程,简单说,可以用matlab对势函数进行拟合,势参数一般通过优化方法拟合系统的一些物理性质得到。程序中首先要给出系统这些基本物理量的实验值(实验上查到的或者第一性原理计算值),在拟合过程中代入一组势参数的初值,此时势函数会计算出一组得到这些物理量的计算值会进行比对,并与实验值进行比对,如果相差比较多,则调整势参数的初值,再次拟合.....直到某一组势参数计算出来的物理量的值与实验值很接近,并且满足结构稳定条件,这样就能初步得到一组比较好的势参数。
在使用之前还要对势函数进行一步验证,根据应用情况,一般需检验结构稳定性、缺陷和非平衡太性质,如果要解决晶格动力学方面的问题,还需要检验声子谱等晶格动力学性能。
4楼2011-12-27 11:39:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csu_lxc

金虫 (正式写手)

引用回帖:
楼: Originally posted by ussy06 at 2011-12-27 11:39:20:
势函数的拟合是一个比较复杂的过程,简单说,可以用matlab对势函数进行拟合,势参数一般通过优化方法拟合系统的一些物理性质得到。程序中首先要给出系统这些基本物理量的实验值(实验上查到的或者第一性原理计算值 ...

谢谢
5楼2011-12-27 14:44:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

161002019

新虫 (初入文坛)

不知道楼主把这个势函数的问题解决了没有?我现在遇到的问题是怎么把TB-SMA势函数改到lammps中去,请指教。非常感谢
fighting!
6楼2017-09-27 10:17:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 csu_lxc 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 12/600 2026-10-13 20:40 by Merryya
[硕博家园] 极限拉扯,如何抉择 +9 Isak0320 2026-10-12 13/650 2026-10-13 11:57 by Isak0320
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见