24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1556  |  回复: 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 的主题更新
信息提示
请填处理意见