求助：使用REAXFF和周期性电场出现错误

wildon

大家好，我在使用REAXFF和周期性电场出现错误，“ERROR: Cannot (yet) use fix qeq/reaxff/omp with variable efield (src/REAXFF/fix_qeq_reaxff.cpp:408)
”我的边界条件为 p p f（reaxff下电场只能加在非周期性边界）
电场部分的代码为：
    variable E0 equal 1
    variable f equal  1e2
    variable e equal v_E0*sin(2*PI*v_f*time*1e-15)

    fix             fxef all efield 0 0 v_e
运行后报错“ERROR: Cannot (yet) use fix qeq/reaxff/omp with variable efield (src/REAXFF/fix_qeq_reaxff.cpp:408)。

但我修改命令  fix  fxef all efield 0 0 v_e 为 fix  fxef all efield 0 0 ${e} , 代码可以正常运行，v_e和${e}的区别参考variable命令两种不同的使用方式“v_"和"$"的区别 - 知乎 (zhihu.com)
然而fix  fxef all efield 0 0 ${e}，似乎并没有起任何的效果，我对比了 有/无（fix  fxef all efield 0 0 ${e}）这条命令的区别，温度变化上没有任何区别。

请问上述问题应该如何解决，谢谢！

wildon

想了个笨方法解决这个问题。写一个循环，在一定步长内将电场视作常数，根据步长改变电场值。
从大尺度看来，只要步长足够小，可近似认为是周期性变化电场。
循环代码如下，可根据自己需要修改：
“    variable a loop 10
    label loop_start
    fix   xef all efield 0 0 ${a}
    run  100
    unfix fxef
    next a
    jump SELF loop_start
    ”

hbfy321

我在reaxff中模拟交变电场，使用ppf边界会直接产生原子丢失问题。换成ppp也出现你问题里的fix efield命令没有作用的现象。
请问如何在ppf边界条件中添加电场，不会使原子丢失呢？

wildon

hbfy321 发表于 2023-5-14 16:09
我在reaxff中模拟交变电场，使用ppf边界会直接产生原子丢失问题。换成ppp也出现你问题里的fix efield命令没 ...

可以添加反射墙。
ppf边界可以加交变电场吗？