lammps相关博文导航

为了访问者较为系统了解本博客所提供关于LAMMPS的信息,建立本导航页面,更新关于LAMMPS的相关博文。

1. LAMMPS基本

2. 如何计算系列

3. 命令翻译系列

4. 常见数据查询表

5. 相关软件及教程

添加新评论 (无需注册,可直接评论)

已有 31 条评论

  1. 博主,如果lammps提示warning:Too many common neighbors in CNA 2 times此类的话,该如何解决呢?是neighbor skin bin设置大了还是小了

  2. Aim

    博主,您好,我想请教您一个问题。做拉伸的时候,大多数是设置应变率进行控制,我想通过应力控制怎么做?

    1. 我的qq:2571858550

  3. 安生

    博主,我做级联碰撞的,想问一下,首先我要确定一个位置,然后寻找与这个位置最近的一个原子,有且只能选择一个原子,用命令要怎么实现!

  4. 请问下博主,怎么将分子动力学仿真输出的动能转化为温度,然后用fix ave/time 进行平均后输出?我用的
    compute KE Gui ke/atom
    variable k1 atom c_KE/(1.5*1.986857e-3)
    fix 5 Gui ave/time 100 5 1000 v_k1 file Ke.profile 但提示有错误,请大神不吝赐教

    1. 输出的信息中有温度啊?

  5. 何伟

    fix 2 all npt temp 700.0 700.0 100.0 aniso 0.0 0.0 10.0 drag 1.0
    博主,您好,感谢上次您对我的问题的解答,让我受益匪浅。想请教您一个新的问题,就是我用上面的指令来控温,但是一运行,发现温度刚开始会一下子变成1600K,然后才降下来,不用这条指令,温度就没有这种爆炸的情况,想请问您这是怎么回事,有什么办法可以替代这条指令吗?对了,我的材料是铜铝。

  6. 何伟

    博主,您好,您的文章对我帮助很大,在此表示感谢!但是对势函数的调用命令还是感到很困惑,就是我做的是铜的拉伸,但是有很多单向铜的势函数,我应该怎么选择呢?我目前选择的是Cu_u3.eam,我用的指令是pair_style eam pair_coeff **Cu_u3.eam Cu,但是一直显示有问题,希望博主能帮忙解答,万分感谢!

    1. Aim

      你好,我做的是Ni-Ti的拉伸卸载循环模拟,咱两的有相似之处,能不能加个好友,讨论下

      1. 你好,我也是研究NiTi合金的,能否加你个好友?

        1. Aim

          可以呀。你qq多少?

    2. 试试下面的命令:
      pair_style eam
      pair_coeff * * Cu_u3.eam
      参考:http://lammps.sandia.gov/doc/pair_eam.html

  7. 玉玉

    您好,我是LAMMPS新人一枚,想问一下没有找到自己需要相关的势函数,如果要自己搞势函数难度大吗?

    1. 建议你看看这篇文章:如何找到自己需要的原子间相互作用势函数?
      http://www.52souji.net/how-to-find-interatomic-potential.html

  8. 博士,您好,
    一直在看你的博文,对我这个初学者帮助很大,在此感谢!因为刚接触这个没有多久,看了您关于势函数部分的博文,但还是不清楚该怎么调用这些命令,比如我做切削单晶Si的模拟,采用tersoff势函数,但不太清楚命令里的数据对应的是什么意思,能不能详细讲下LAMMPS中势函数调用命令该怎么写。
    万分感谢了,相关资料太少,现在无从下手,希望能得博士的解答~

    1. pair_style tersoff
      pair_coeff * * Si.tersoff Si
      参考:http://lammps.sandia.gov/doc/pair_tersoff.html

      1. 好的,谢谢博士~

  9. 项俊锋

    您好,博主!
    本人做切削方加工,之前做有限元的,博士课题是跨尺度模拟,因此涉及分子动力学!请问下能否推荐几本您认为的关于分子动力学学习的经典书籍,当然最好是关于lammps,麻烦您了!

    1. 我爱搜集网博主

      关于分子动力学的书:

      英文
      Ercolessi, Furio:Computer simulation of liquids
      Allen, M P:A molecular dynamics primer

      中文
      陈舜麟:计算材料科学
      张跃、谷景华等:计算材料学基础

      关于lammps我还没看到过。

  10. andy

    你好,请问想用分子动力学方法做模拟,但又找不到合适的势函数该怎么办?自己做一个势函数有什么难度?

    1. 我爱搜集网博主

      请参考:【如何找到自己需要的原子间相互作用势函数?】
      http://cms.52souji.net/viewtopic.php?id=21

  11. 刘柱

    您好,请教一下,对lammps的dihedral style opls 有个困惑的地方,比如还是我的体系二氯乙烷分子,C2H4Cl2, 有9个三种不同的二面角 即三种二面角的力场参数 对这9个二面角的能量贡献都得考虑是吧,谢谢答复

    1. 我爱搜集网博主

      我个人没有做过有机物的模拟。我个人感觉都是要考虑的。不过也许这9个并不都是独立的。你自己考虑下。

  12. 兮兮鱼

    先生您好!我是一名大三的学生,目前在做大创,从您的博文中得到了很多指引,谢谢!
    然后现在我导师给我的一些关于LAMMPS的基础问题,我头绪不是很清晰,希望您看到以后可以给我一些指导!需要我做的都可以联系我的QQ583130935,我做的项目是“金属基底上石墨烯涂层力学性质研究”,以下是导师现阶段给我的问题,希望您可以推荐一些关于这个方面的文献和作者,也可以给我一些思路或者指引。也不一定要关于这些问题,有什么建议对LAMMPS新手说的都可以告诉我~先行鞠躬!
    金属基底的模型如何导入?
    金属与石墨烯之间的作用如何描述?
    金属-石墨烯涂层的稳定平衡结构是什么形态?
    作用面内剪切或拉伸荷载后涂层如何变形?
    温度对变形有何影响?

    1. 我爱搜集网博主

      您好!您的问题太个人化了,我不好回答。

      1. weikt06

        楼主,您好。看了您的博文也没有找到相关的势函数,楼主知识丰富,对这方面的研究很深入,麻烦您能否帮我找一下Mo-C体系的势函数,万分感谢!

    2. 兮兮鱼

      我在您博文中找到C-Fe势函数,可以用于第二个问题“金属与石墨烯之间的作用如何描述”吗?

      1. 我爱搜集网博主

        理论上是可以的。但你需要测试。

  13. 朱盼

    您好,
    我是学生,现在想要做Fe和PP的界面模型。不知道您能给我提供一些资料和指导吗?

    界面模型建立的基本流程 模型中需要用到的势能,Fe、PP、Fe-C、Fe-H的势能,希望能够附上相应的文献出处 如果没有现成的势能,而我又暂时做不了实验标定,那么有什么建议呢您?
    请大神指点,不胜感激。如果需要收费,请告知,在我力所能及的范围内将尽力支付。
    万分感谢。
    联系方式:QQ782515448 (PS:可能需要验证,麻烦输入一下lammps就好了O(∩_∩)O~)
    1. 我爱搜集网博主
      如果你会使用建模软件,你就知道该如何构建。 请自行搜索数据库网站。据我了解,Fe, Fe-C, Fe-H的势都是有的。PP不清楚是什么。 如果没有势函数,要么放弃,要么自己发展。
  14. 大飞

    学了lammps一个多月才看到楼主的文章真是损失 万分感谢啊