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

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

已有 140 条评论

  1. lammps新手

    楼主您好我是做氢同位素在钨中的扩散系数计算的,想问问您有没有W-H的势函数或相关的文献资料,我看上面纯相的合金相的都不太适用。。感谢!

    1. 精度要求不高可以使用W-H势,要求高可能需要自己开发了。

      1. lammps新手

        您给的网站中有一个 W-H-HE的势函数,可以用于W-H系统吗,因为我没要找到单独的W-H势函数

        1. 应该是可以的,三元体系的势一般是适用于二元的。

  2. 嘻嘻哈哈乐乐

    楼主您好,我想研究铜晶体里有He掺杂时的晶体变形。发现Cu-He之间好像只有ZBL势函数,这是用来做碰撞的,不知道我这种普通低速度情况是否适用。或者楼主有什么好的势推荐一下。非常感谢!

    1. 忽忽悠悠

      ZBL主要是针对短程力的修饰,一般EAM势中修改了嵌入势的部分。理论上来说ZBL修正后的势对一般模拟影响不是很大。我想知道层主的Cu-He势函数在哪里找到的?谢谢

    2. 这个真不好说~需要根据自己的需要进行测试

  3. 似水流年

    您好,可以介绍一下lammps中pair_coeff命令中参数各指代什么吗?

  4. Ethan

    怎么办,之前没下载,现在进不去了,早知道把势函数全部下载下来了。。。。。

  5. dd

    下不了啊 楼主

  6. qiqibao

    楼主您好,麻烦问下有没有Co-Nb-Ta的势函数呢,或者多元的关于它们的也行,急求,谢谢!

  7. 新城

    楼主您好,有没有Ca-O的势函数啊

  8. 如意

    您好,我是lammps新手,想用它模拟结晶,一种有机物溶解在另一种有机物熔融液中,想模拟冷却结晶这个过程,请问可以实现么?有没有类似的例子可以借鉴?谢谢您

  9. YK

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

  10. 湖大赵同学

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

  11. 湖大赵同学

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

  12. 虫子

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

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

  14. taotao

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

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

  15. 西瓜啦啦啦撒

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

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

      1. 西瓜啦啦啦撒

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

  16. Bruce

    楼主您好!

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

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

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

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

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

  19. qazwsx

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

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

  20. 东大哆哆

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

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