在lammps中添加外部非常量力场

Question

我想添加一个依赖于原子位置的力场。我试着计算原子的位置，使用它作为变量，然后设置力场，但我得到了：

错误:修复addforce的变量名不存在

我尝试过的代码是：

variable a equal c_X
variable extforce equal ${maxforce}*step/${eqrun}*a
fix externforce A  addforce v_extforce 0 0 

Answer 1

如果您希望它随原子位置而变化，则需要使用atom样式的变量。我不知道你们的其他变量是什么，举个例子，这是力的y分量，作为x位置的函数：

variable addfx atom x*0.2
fix externforce all addforce 0. v_addfx 0.

你得到的错误意味着你需要首先“使用”变量，即使用thermo_style命令输出它。