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

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

已有 185 条评论

  1. 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)

  2. taotao

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

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

  3. 西瓜啦啦啦撒

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

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

      1. 西瓜啦啦啦撒

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

  4. Bruce

    楼主您好!

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

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

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

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

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

  7. qazwsx

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

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

  8. 东大哆哆

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

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

  9. 111

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

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

  10. 一棵树

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

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

  11. 自由

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

  12. 小豆闯天下

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

    1. Aim

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

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

  13. 123456

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

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

  14. 王浩

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

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

  15. 李晨

    楼主您好。 我是一名刚刚入手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

  16. 李强

    想请问博主是否有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体系,可以使用同样一个势函数。

  17. 大大

    请问博主,您有Li-Fe的势函数吗?我通过你给的网址目前没有,还有其他的方向可以找吗?谢谢回答

    1. 可以找找文献,感觉应该不难找到。

  18. zbkxbk

    请问楼主,我想找个Rebo-CHO的势,参考文献已找到,这个势文件应该从哪里找呢?

    1. LAMMPS支持这种势类型吗?

      1. zbkxbk

        我想跑一个包含CHO元素的有机分子和C纳米管的结构,纳米管的势用rebo势,但是不知道和这个分子之间的势该怎么办?求指教

        1. 抱歉,本人没用过这种势。你自己摸索看吧~

  19. 一包花生米

    请问有没有SiC.sw势函数,找不到

    1. 我看lammps安装包中提供了SiC的tersoff势,sw势并没有。

      1. 一包花生米

        好的~谢谢~谢谢~

  20. yuan

    我要用金的势能函数,我在哪里修改里面的cut_off呀
    Au functions (universal 3)
    79 196.97 4.0800 FCC
    500 5.0100200400801306e-04 500 1.1212121212121229e-02 5.5500000000000114e+00

    -4.8957152905617285e-01 -7.9715640025473178e-01 -1.0578883960846852e+00 -1.2926127947875230e+00

    -1.5095569422429875e+00 -1.7122648936144031e+00 -1.9048413422057351e+00 -2.0891201727137485e+00

    1. 你这个应该是从势函数文件中复制的吧?这个是已经被离散化的,如果你要修改cutoff,需要在离散化之前修改。