如何在lammps里加入新的势函数形式?_lammps,如何自己建立势函数_MrZhengGang的博客

可以先参考已经有的势函数，弄清楚框架和写法，再写自己的。有的时候不仅要写 .cpp 文件，还要写一个 .h 头文件。

把写好的东东放在lammps/src目录，然后按照以前的方法make即可。

我的个人网站：传送门

官方参考官方手册，先看看：传送门。自己编程如果没有，就自己写吧。可以先参考已经有的势函数，弄清楚框架和写法，再写自己的。有的时候不仅要写 .cpp 文件，还要写一个 .h 头文件。编译把写好的东东放在lammps/src目录，然后按照以前的方法make即可。可以参考文章：lammps 编译安装教程（一）。我的个人网站：传送门...

上期介绍了LJ 势函数参数写法，本期主要介绍金属体系势函数的设置。金属原子之间没有键连接，因此，在 lammps 模拟中，金属体系的势函数类型为pair_style，而不是bond_style。模拟金属体系时，可以用LJ 势描述金属原子之间的受力，不过更精确的是嵌入原子势（EAM），eam 势函数公式为： eam 势由两部分组成，在原子对势（pair）的基础上添加了电子云密度相关项，比单纯的对势精确度更高。 eam 势函数写在一个以“.eam”为后缀的文件中， lammps 自带的势函数包含一部分eam 势文件，也可以到

本文继续介绍 lammps 模拟中的势函数设置方法。在前面的文章中，介绍过高熵合金的建模：高熵合金建模高熵合金包含的原子数较多，势函数的设置相对复杂，在平时一对一的咨询中，这种问题也相对较多。下面介绍三种高熵合金势函数设置方法。 1.下载专用势函数下面的网站包含了大多数原子的势函数： https://www.ctcms.nist.gov/potentials/ 如需下载包含Fe原子的势函数，直接点Fe：出现势函数列表：比如Al-Co-Cr-Fe-Ni就是高熵合金AlCoC 本文介绍 lammps 模拟中pcff 势设置方法。 cvff、pcff是ms文件转换为 lammps data文件最常用的两种势。相比于cvff 势，pcff 势参数更多，但在设置方式上和cvff 势过程是一样的。本文以沥青材料为例，介绍pcff 势设置方法。 (1) 建模并转换为data文件在ms中使用AC模块建立沥青模型，使用forcite模块设置pcff力场，导出为asphalt.car。使用msi2lmp转换为data文件： msi2lmp asphalt -class II -f

使用 lammps 模拟金属体系，最常用的势函数为eam、meam 势。对于大多数金属或者合金来说，可以从网上下载相对应的eam或meam 势力场文件。而对于一些新的合金材料，可能就找不到相对应的力场文件。解决的方法有两种，一是采用混合势（pair_style hybrid），第二种方法是自己拟合势函数文件。混合势的写法在前面文章中已经介绍过，本文介绍一种拟合EAM 势的方法。该方法是运行Xiaowang Zhou编写的一个拟合程序，该程序可自动生成所需要的合金势。这个程序能够拟合的合金元素有：本文介绍OPLS 势函数的设置方法。 lammps 力场大体分为两种情况：一种类似于eam这种势，参数是写到eam力场文件中的，在使用的时候需要下载力场文件，在设置时不需设置力场参数，只需映射原子类型即可。另一种力场则不需要力场文件，如lj/cut。这种力场大多使用公式描述原子间相互作用力，这些公式内嵌在 lammps 代码中，在设置力场时，不需要下载力场文件，只需要设定公式中的参数即可。本文所讲的OPLS属于第二种力场，不需要力场文件，但是需要设置力场参数。 OPLS对应的公式为：