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

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

已有 154 条评论

  1. 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,需要在离散化之前修改。

  2. crystal

    楼主,请问里面的势函数有的只有一列,有的四、五列,有什么不同吗?

    1. 我猜应该都可以吧,只是格式不同,你可以测试一下。

  3. grap

    楼主,你好,请问有没有Ga、N、C三元的Tersoff势函数,急求,谢谢了

    1. 抱歉,这个我也没有。你可以参考我另外写的一篇文章:
      http://www.52souji.net/how-to-find-interatomic-potential.html

  4. grap

    楼主,你好。我有一个问题想请教一下您。我做的是金属与石墨烯界面热阻的研究,因此需要金属与石墨烯片之间的势函数,我是这样写的
    pair_style hybrid lj/cut 10.0 rebo eam
    pair_coeff * * eam Fe_mm.eam.fs
    pair_coeff * * rebo CH.airebo NULL C
    pair_coeff 1 2 lj/cut 0.09463 2.05
    运行时报错,Invalid mass value,之前没有接触过金属的,问一下楼主用eam势函数的时候有哪些要求,谢谢。

    1. lammps小白

      用的是eam/fs势
      第一行改成 pair_style hybrid li/cut 10.0 rebo eam/fs

    2. eam势还是比较好用的啊,我估计你的这问题是处在hybrid上。

  5. 庞皓升

    请问楼主有木有玻璃的BKS函数

    1. 我爱搜集网博主

      这种势我也不了解,你可以自行Google

  6. 傅依苔

    楼主,有没有 Li-Sn之间的势函数?谢谢!!

    1. 我爱搜集网博主

      我没有。你可以Google

  7. happy369

    楼主您好 非常谢谢您无私的分享~
    想向楼主请教一下有关copper 和carbon原子之间的BOP势能的问题,我在网上和文献里没有没有找到有过copper 和carbon原子之间的BOP势能的参数和势能文件,我想问问楼主怎么可以得到这些参数和哪里可以下载到BOP的势能文件?

    1. 我爱搜集网博主

      如果没有势参数,那就需要自己拟合。有点复杂。

  8. nn

    请问楼主,势函数下载下来的这些内容是怎么得来的?
    R.R. Zope and Y. Mishin, Phys. Rev. B 68, 024102 (2003).
    LAMMPS setfl format. Conversion by C. A. Becker from Y. Mishin files.
    26 Sept. 2009. http://www.ctcms.nist.gov/potentials
    2 Ti Al
    10000 0.6967034321178504E-03 10000 0.6724884000000001E-03 0.6724884000000000E+01
    22 0.4786700000E+02 0.4148000000E+01 fcc
    -7.7941043376E-0010
    -4.3651571088E-0003
    -8.7142355246E-0003
    -1.3047261839E-0002
    -1.7364261862E-0002
    -2.1665261408E-0002
    -2.5950286292E-0002
    -3.0219362339E-0002
    -3.4472515430E-0002
    -3.8709771215E-0002
    -4.2931155566E-0002
    -4.7136694064E-0002
    -5.1326412416E-0002
    -5.5500336255E-0002
    -5.9658491180E-0002
    -6.3800902842E-0002
    如果我是引用另外一个人写的势函数文件,只要他的文章内容,但是Lammps运算需要的势函数文件eam/alloy文件怎么得到呢?

    1. 我爱搜集网博主

      根据参数和势函数的公式自己编程实现。

  9. nnnn

    博主你好,请问您提供的Ti-Al势函数文件,我用于Ti3Al是否也可以行?我之前γ-TiAl是用了这个,请问您这两种都能用这个势函数吗?

  10. MA HAIDONG

    有没有C-Zn 和C-ZnO 这两种势函数呢?

  11. wulianp

    楼主,势函数里面的数字都代表什么?有没有相关资料

  12. 陈思

    请问楼主有Cu-Fe的eam/fs势函数吗?还请问eam/fs势文件里第四行开始的5列数是什么意思?是怎么来的呢?谢谢!

    1. song xuan

      您好,我想下载airebos势函数,但是下不下来,可以麻烦发我邮箱吗,谢谢。

  13. yan

    急切需要Cu50Zr50和Cu48Zr52的实验值熔点~~~谢谢谢!!!

  14. 魏磊

    请问楼主,很多晶体结构的势函数从哪里下载?比如说钙钛矿结构的KTaO3 ,谢谢!

    1. 我爱搜集网博主

      首先,你看看上面给出的链接里有没有。但其实很多都是没有的。
      然后,也是最普遍的方式,就是在数据库中搜索相关文献。在关键词上,你可以用你的体系加上“interatomic potential”。
      如果经过大量的搜索,仍然没有找到,那很可能还没有人做过这方面的模拟。或者你自己开发势函数,或者你放弃模拟这个体系。

  15. 张帅

    楼主,有没有 Cu/diamond之间的势函数?楼主见多识广,帮帮忙

    1. 我爱搜集网博主

      我也是google到的,供你参考:
      http://www.sciencedirect.com/science/article/pii/S0038109803011554
      http://www.sciencedirect.com/science/article/pii/S0928493101003083

      1. 张帅

        谢谢楼主,我会认真看看文献的,愿我爱搜集网越来越好,越来越全!

        1. 张帅

          抱歉这么迟才回复您。

  16. 张帅

    有没有Cu/diamond之间的势函数,楼主你见多识广,帮帮忙

    1. MA HAIDONG

      铜和碳的势函数用morse

  17. wei.yang

    您好,初学lammps,关于势函数想请教一下,我模拟Ti/TiC结构,需要用meam势函数,这个函数需要library.meam 和补充文件 Ti/TiC.meam, 对于library文件中的十几个参数,单独的Ti和C的参数都可以找到,但是我有TiC复合结构,是否需要复合结构的参数?另外我还需要模拟TiN,N是个气体没有晶体结构,那如何准备library内的参数呢? 谢谢您,不胜感激。

  18. yehaibo

    我刚刚开始学lammps,我想问下了lammps计算的整个过程是啥样的,需要什么文件,

    1. 我爱搜集网博主

      参考:http://www.52souji.net/lammps123/

  19. 宋宋

    您好!我用的是Al-Ni:eam/alloy-1 势,查资料说eam势不用设置质量,但是报错All masses are not set,怎么解决呀?谢谢!

    1. 我爱搜集网博主

      eam势确实是不需要设置质量的。至于你遇到这个问题,我无法只根据你这些描述判断问题的出处。
      建议你在看看pair_style eam/alloy的使用,相信应该不是很难解决。

      1. 宋宋

        感谢!问题已解决,这个势文件中似乎没有质量信息,加上就解决了!

  20. yan

    2.一般而言,高元合金的势函数适用于低元合金或纯相,只需在设置pair_coeff时选择相应元素即可。对于2 Mo Mo 这样的eam.fs势函数 ,而我的原子只有type 1,应该怎么写势函数? pair_coeff * * Mo.eam.fs Mo可以吗?请各位赐教~

    1. 我爱搜集网博主

      对于Mo.eam.fs,其中只包括Mo的参数,所以应该不需要再设置元素名称。
      应该pair_coeff * * Mo.eam.fs就可以了吧。你可以试试。