求助：在lammps中引入开关函数模拟偶氮苯类化合物N=N双键的翻转，如何实现？

wangyi

本人是lammps小白，想用lammps计算偶氮苯类衍生物N=N双键的翻转，具体的计算思路参考文献：Journal of Chemical Physics, 2013, 139(1):224-5882，（1）弛豫一定时间；
（2）输出C-N=N-C二面角信息，产生与其数量对应的0-1间的随机数；
（3）计算开关函数P（φ）（C-N=N-C二面角），比较随机数和开关函数的大小，分别赋予顺式-反式的力场和反式-顺式的力场；
（4）重新计算，计算完返回第一步；
（5）执行这样n个循环。
文献中的流程如图所示，主要是获得特定的二面角信息、产生随机数与在执行循环和判断不太会处理，希望各位老师不吝赐教。