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

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

已有 178 条评论

  1. 下一站

    博主,si的tersoff势函数无法打开

  2. Voidshatter

    请教下博主,我要研究graphene和SiC组成的复合体,也就是graphene铺在SiC表面上。SiC我采用的是tersoff势函数,graphene我采用的是rebo势函数。想问下graphene和SiC之间的势函数怎么选取更精确一些?因为我要外加力使材料压缩,可以采用tersoff势么?

    1. songxuan

      老师您好,airebo文件一直下载不下来,可以发邮箱吗 2727075497@qq.com

  3. Lammps 新手

    airebo势函数 下载不了啊。有没有前辈分享一下啊?邮箱:lixue_wangjie@163.com

    感谢感谢!

  4. lammps新手

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

    1. Lammps 新手

      airebo势函数 下载不了啊。有没有前辈分享一下啊?邮箱:lixue_wangjie@163.com

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

      1. lammps新手

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

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

  5. 嘻嘻哈哈乐乐

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

    1. 忽忽悠悠

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

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

  6. 似水流年

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

  7. Ethan

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

  8. dd

    下不了啊 楼主

  9. qiqibao

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

  10. 新城

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

  11. 如意

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

  12. YK

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

  13. 湖大赵同学

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

  14. 湖大赵同学

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

  15. 虫子

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

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

  17. taotao

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

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

  18. 西瓜啦啦啦撒

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

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

      1. 西瓜啦啦啦撒

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

  19. Bruce

    楼主您好!

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

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

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