LAMMPS常用原子间势函数下载

本页面提供LAMMPS下常用体系的原子间势函数的下载。持续更新。

对于各种不同势的介绍可以先参考LAMMPS文档。右键另存为下载。

一. 纯相的势函数

二. 二元合金相的势函数

三. 三元合金相的势函数

四. 几点说明

  1. 调换元素顺序,不影响势函数的选用。比如Ag-Cu与Cu-Ag势函数相同。
  2. 一般而言,高元合金的势函数适用于低元合金或纯相,只需在设置pair_coeff时选择相应元素即可。比如Cu-Fe-Ni的势函数适用于Cu-Fe或Fe-Ni等。
  3. 因为来源问题,可能有些势函数原本是一样的,请自行鉴别。

五. 势函数来源

  1. Interatomic Potentials Repository Project
  2. 链接: http://redmine.scorec.rpi.edu/anonsvn/lammps-cuda/potentials/
  3. 链接:https://cmse.postech.ac.kr/html/research/2nnmeam.htm

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

已有 126 条评论

  1. YK

    楼主您好!我是lammps新手,想请教一下势函数文件格式应该怎样弄啊,看到lammps官网上面写的是setfl文件,但具体不会弄。谢谢!

  2. 湖大赵同学

    难混溶的Cu-W也行,我主要是想做难混溶的Cu基合金模拟,谢谢你!

  3. 湖大赵同学

    你好,能麻烦你帮我找Cu-Co的势函数吗?我急着要用,我自己实在找不到,非常感谢!!!

  4. 虫子

    楼主你好,想请教一个问题。我用fix npt temp 300 300 0.5 iso 0 0 5,(模型时间步长为0.005) 进行水分子压力调控,但是总是达不到目标压力。npt弛豫期间压力会在10bar到-20bar之间波动,进入nve弛豫期间,压机波动变大会在-200到200之间波动,请求您的帮助,谢谢

  5. HIPPO

    楼主您好,想请教一下势函数文件里面第一行设置的参数是什么含义呀?然后下面五列数据分别是什么呢?比如Fe的EAM势函数:
    1 Fe
    10000 3.00000000000000E-0002 10000 5.30000000000000E-0004 5.30000000000000E+0000
    26 5.58500000000000E+0001 2.85531200000000E+0000 bcc
    0 -1.73205399240757E-0001 -2.44950248213795E-0001 -3.00002866354823E-0001 -3.46415257255683E-0001

    1. 参考:http://lammps.sandia.gov/doc/pair_eam.html
      lines 1,2,3 = comments (ignored)
      line 4: Nelements Element1 Element2 … ElementN
      line 5: Nrho, drho, Nr, dr, cutoff

      line 1 = atomic number, mass, lattice constant, lattice type (e.g. FCC)
      embedding function F(rho) (Nrho values)
      density function rho(r) (Nr values)

  6. taotao

    楼主您好,我是研一lammps新手,请问能加一下您的QQ或者微信吗?求您的QQ号或者微信号。。。。。。

    1. 你好!有问题可以在这里提问,我会尽力回答。

  7. 西瓜啦啦啦撒

    楼主您好,我想下个势函数文件,怎么下不下来呢

    1. 多尝试几次,应该是没问题的。

      1. 西瓜啦啦啦撒

        那个C-Si-Ge 的tersoff函数 还是打不开

  8. Bruce

    楼主您好!

    能不能介绍一下如何用pair_style eff/cut 势函数模拟一个元素在不同温度压力下,求出平衡常数的方法?

  9. 博主你好,我想请教关于升温速率的设置,如果我用fix npt temp Tstart Tstop 这时候温度是波动上升的,而且不一定会完完全全到达指定温度。除了用循环恒温的方法,能不能按照一定的速率升温呢

    1. 你的这个升温过程是很正常的。按我理解,在分子动力学里面,很难做到你期望的“一定速率”,总会有波动。

  10. 楼主你好,感谢你分享了这么多势函数。
    我想问下你对Al-Zn的势函数了解吗?请问你能帮忙找到吗?

    1. 我也没有~自己找找看~参考:
      http://www.52souji.net/how-to-find-interatomic-potential.html

  11. qazwsx

    楼主你好,感谢你分享了这么多有用的势函数信息 。
    我现在非常需要 Fe-Mo-C,或Ti-Mo-C的势函数,请问你能帮忙找到 吗?

    1. 我也没有~自己找找看~参考:
      http://www.52souji.net/how-to-find-interatomic-potential.html

  12. 东大哆哆

    楼主你好,最近想做氦离子轰击硅衬底的分子动力学模拟,用的是lammps。请问He-He和He-Si之间选择哪一种势函数呢?

    1. 这个你要看看文献中一般用什么类型。就我了解,对于He-He一般的二体势就可以了。

  13. 111

    想知道铜和水分子之间的势能参数,以及二氧化硅和水之间的势能参数

    1. 这个势不好找,你在谷歌学术中找找吧。

  14. 一棵树

    楼主你好,请问有SI和H的eam函数吗?

    1. eam一般用于金属和合金的体系。

  15. 自由

    楼主你好,还是lammps入门,想模拟三元合金相Al-Cu-Fe,但现成势函数有Al-Cu,Al-Fe,Cu-Fe三个二元合金相的势函数,在考虑可以结合几个二元的势函数来模拟三元的么?

  16. 小豆闯天下

    您好,有Ni-Ti的eam 的势函数吗?

    1. Aim

      您好,你也是做镍钛的吗?可以加个好友相互交流吗?qq:215810645

    2. 参考这个:https://www.ctcms.nist.gov/potentials/Ni.html

  17. 123456

    有没有Co-C相关的势函数???

    1. 我也没有~自己找找看~参考:
      http://www.52souji.net/how-to-find-interatomic-potential.html

  18. 王浩

    博主您好,求助一个MO-W合金的势函数,俺不会拟合,能给一个么,谢谢您!

    1. 我也没有~自己找找看~参考:
      http://www.52souji.net/how-to-find-interatomic-potential.html

  19. 李晨

    楼主您好。 我是一名刚刚入手lammps的新人。我的研究方向是锂离子电池。我想问一下
    Li-Si 的势函数您有么? 因为我看过一些文献都是做分子模拟的 但是我找不到势函数。 还有一个问题就是随着Li数量的增加,Li-Si 的势函数 会不会发生变化? 谢谢

    1. 据我所知,Li-Si势是有的,应该是2NN MEAM,具体你可以google scholar一下。
      至于第二个问题,在使用势函数的时候一般不需要考虑成分的变化。

      1. 李晨

        非常感谢博主的回复。 我找到一片关于Li-Si 势函数的文献: A second nearest-neighbor embedded atom method interatomic potential for Li–Si alloys.
        里面有一章 2.1. 2NN MEAM potential 。 但是里面都是一些公式,请问怎么才能应用成MEAM的文件格式。

        1. lammps里有相关的命令,具体参阅:
          http://lammps.sandia.gov/doc/pair_meam.html

  20. 李强

    想请问博主是否有Cu-Zn的势函数呢,找了好久都没有找到,能不能帮帮忙呢,十分感谢!

    1. 你可以看看这篇文章,兴许对你有帮助。
      http://www.52souji.net/how-to-find-interatomic-potential.html

      1. 李强

        十分感谢博主的回复,另外我还想问个问题,就是博主在几点说明里提到“调换元素顺序,不影响势函数的选用。比如Ag-Cu与Cu-Ag势函数相同”,我有些不明白,比如说我想找Cu-Al的势函数,但是查到的几乎都是Al-Cu,那我要如何才能得到不同含量Al的Cu-Al的势函数呢,谢谢!

        1. 一般来说,开发势函数的目的就是为了应用到各种不同的场景中。所以,对于不同含量Al的Cu-Al体系,可以使用同样一个势函数。