| 查看: 3327 | 回复: 6 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
舞昂琦VoVo木虫 (小有名气)
|
[求助]
如何在lammps里加入新的势函数形式？
|
|||
| 想在lammps里加入一种势函数形式，但看了一下lammps里已有的势函数源代码，感觉比较复杂，有点无从下手，不知道哪位有类似的经验，能够指教一下，万分感激~ |
回复此楼» 收录本帖的淘帖专辑推荐
 资源收集 | 模拟软件安装精华淘 |
» 猜你喜欢
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有5人回复
-
找导师
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有7人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有9人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- advanced materials不接受full paper形式的文章？
已经有4人回复
- 油水界面，加入表面活性剂后，界面会减弱或消失吗
已经有10人回复
- 各位高手，怎么在windows下编译LAMMPS？
已经有4人回复
- 【讨论】想学lammps是否应该先学习C++
已经有28人回复
- 【求助】tight-binding gulp或lammps里有紧束缚势嘛 ？
已经有5人回复
- 【讨论】lammps里的boundary的位置能够随意更改？
已经有14人回复
- 【求助】LAMMPS中非周期性边界条件和Born-Mayer-Huggins势函数不会设置
已经有4人回复
- 【在线答疑】关于NAMD/VMD
已经有210人回复
舞昂琦VoVo
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1723.2
- 帖子: 179
- 在线: 308.2小时
- 虫号: 1075737
- 注册: 2010-08-16
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
5楼2012-04-10 13:58:00
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zh1987hs: 金币+5, 谢谢 2012-04-10 08:06:43
舞昂琦VoVo: 金币+2, ★★★很有帮助 2012-04-10 13:51:55
感谢参与，应助指数 +1
zh1987hs: 金币+5, 谢谢 2012-04-10 08:06:43
舞昂琦VoVo: 金币+2, ★★★很有帮助 2012-04-10 13:51:55
|
lammps的势函数已经非常全面了，你仔细遇到以下lammps手册，说不定你想添加的，人家已经有了呢。http://lammps.sandia.gov/doc/Manual.html 如果没有的话，建议你首先把lammps里面的源程序里面的势函数看明白，打个比方说，你想添加一个新的bond style ，那么你把 源程序里面的bond_fene.cpp 和 bond_harmonic.cpp 两个程序看明白了，再去写你想要添加的应该就比较容易了。你还需要弄明白lammps的整个程序框架。 看看这个吧：http://lammps.sandia.gov/doc/Section_modify.html#mod_14 |
2楼2012-04-09 22:25:13
【答案】应助回帖
★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zh1987hs: 金币+4, 谢谢 2012-04-10 08:06:58
舞昂琦VoVo: 金币+2, ★★★很有帮助 2012-04-10 13:52:03
感谢参与，应助指数 +1
zh1987hs: 金币+4, 谢谢 2012-04-10 08:06:58
舞昂琦VoVo: 金币+2, ★★★很有帮助 2012-04-10 13:52:03
|
lammps现在已经很庞大了，研究其框架可以找比较早期的版本。等大致了解了之后然后在最新版本上进行代码修改。lammps是用c++写的，虽然不算标准的模块化设计，但相对gromacs的代码来说比较模块了，添加一个新的势基本不影响其他的代码。 读代码从main.cpp开始看，了解主要框架，然后重点参考你和你要加势的最像近倒那一类，比如bond_harmonic,pair_lj等。在这个基础上完成自己的修改。一般是2个文件，一个.h，一个.cpp。你需要完成的方法在bond.h,angle.h,pair.h..里面有。 然后还要修改1个文件： style.h，让lammps知道你的势。 然后就make，debug，。。。 |
3楼2012-04-10 00:20:22
舞昂琦VoVo
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1723.2
- 帖子: 179
- 在线: 308.2小时
- 虫号: 1075737
- 注册: 2010-08-16
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
4楼2012-04-10 13:54:53
