| 查看: 658 | 回复: 0 | ||
161002019新虫 (初入文坛)
|
[求助]
金原子跟碳管的相互作用 TB-SMA势函数改进
|
|
各位大神,这是我编的势函数,模型是把金原子放在碳管里面,然后对碳管加热,来观察温度跟金原子能量之间变化的关系。文献上金原子之间的相互作用能是TB-SMA势,但是我不会改,所以用的lammps手册里面的,算出来的总能量特别大,被老师凶了,哪位大神能帮我解释一下?或者教教我TB-SMA势要怎么改到lammps当中去,请指教,非常感谢,在线等! 这是势函数 pair_coeff 1 2 lj/cut 0.01273 2.9943 10 pair_coeff 1 1 lj/cut 0.00288 3.5000 10 pair_coeff 2 2 eam Au_u3.eam 这是算出来的部分结果 Step tempau PotEng KinEng TotEng Press 0 0 4522288 0 4522288 5.1996866e+08 100 307.07543 4522271.5 17.86165 4522289.4 5.1997868e+08 200 360.5601 4522270.6 20.972691 4522291.5 5.1998058e+08 300 225.81346 4522280.5 13.134887 4522293.6 5.1997687e+08 400 222.56927 4522283.5 12.946182 4522296.4 5.1998249e+08 500 298.79333 4522280.8 17.379905 4522298.1 5.1997551e+08 600 352.69761 4522277.8 20.515354 4522298.3 5.1997417e+08 700 291.82476 4522280.2 16.974564 4522297.2 5.1997758e+08 800 275.6671 4522279 16.034722 4522295.1 5.1997591e+08 900 408.21137 4522270.9 23.744421 4522294.7 5.1997991e+08 1000 260.68251 4522276.8 15.163114 4522292 5.1998038e+08 1100 288.10594 4522275 16.758252 4522291.8 5.1997747e+08 |
» 猜你喜欢
有没有人能给点建议
已经有5人回复
假如你的研究生提出不合理要求
已经有12人回复
实验室接单子
已经有7人回复
全日制(定向)博士
已经有5人回复
萌生出自己或许不适合搞科研的想法,现在跑or等等看?
已经有4人回复
Materials Today Chemistry审稿周期
已经有4人回复
参与限项
已经有3人回复
对氯苯硼酸纯化
已经有3人回复
所感
已经有4人回复
要不要辞职读博?
已经有7人回复













回复此楼