24小时热门版块排行榜    

查看: 1982  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jRl3mE6ddZGq 2026-08-03 4/200 2026-08-04 04:54 by k7OM8YghWbkC
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jRl3mE6ddZGq 2026-08-03 7/350 2026-08-04 04:50 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:46 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 5/250 2026-08-04 04:44 by k7OM8YghWbkC
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:41 by k7OM8YghWbkC
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 jRl3mE6ddZGq 2026-08-03 7/350 2026-08-04 03:35 by k7OM8YghWbkC
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 12/600 2026-08-04 03:17 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 16/800 2026-08-03 12:23 by fuweiguochen
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 12/600 2026-08-03 08:22 by syl200707
[基金申请] 微信指数没变化,科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
请填处理意见