求助：Reaxff中模拟周期电场下热解甲苯问题

hbfy321

参考文献：10.1080/00268976.2013.860245
文献中施加周期电场后，体系的温度变化是和电场变化相关的。



我编写的in文件程序：
units                real
neighbor        2.0 bin
neigh_modify         every 5 delay 0 check no
boundary        p p p

atom_style        charge
read_data        toluene50.data

pair_style        reaxff NULL safezone 2.0 mincap 100 minhbonds 100
pair_coeff        * * ffield.reax.cho C C H
fix                1 all qeq/reaxff 1 0.0 10.0 1.0e-6 param.qeq

variable         E0 equal 4
variable         E1 equal 0.6
variable         f equal  100
variable        1e equal v_E0*sin(2*PI*v_f*1e-6*step*dt)
variable        2e equal v_E1*sin(2*PI*v_f*1e-6*step*dt)
thermo                100
thermo_style        custom step v_1e temp press vol pe ke etotal
fix                3 all nve
fix                4 all temp/berendsen 300 2500 200
fix                2 all efield 0 0 ${1e}
fix_modify        2 energy yes
dump                1 all custom 100 micro.lammpstrj type x y z
run                100000


问题：
1.我按照文献中的模拟条件模拟，电场下的温度变化是线性的，可以说是电场对温度没有起到作用，这是什么原因？
2.我查手册了解到reaxff施加电场得在非周期边界条件下，但我换成ppf后会导致原子丢失，如何解决这个问题？(弛豫是在ppp条件下进行的，加电场后转换为ppf)

mengxiangidea

1.ppp边界施加电场不行
2.原子丢失的问题原因有很多，归结到底，还是模型出问题了。或者考虑盒子变大一点，或者检查是否存在原子重叠的过程，延长弛豫时间，存在原子重叠，更改密度。

wildon

fix 2 all efield 0 0 ${1e} 的意思是 fix  2 all efield 0 0 0，${1e}这个不是变量，而是取函数的第一个值。
所以实际上是加了一个恒定电场。

ppf边界如果出现原子丢失，可以加一个反射墙试试。

hbfy321

wildon 发表于 2023-5-15 15:10
fix 2 all efield 0 0 ${1e} 的意思是 fix  2 all efield 0 0 0，${1e}这个不是变量，而是取函数的第一个值 ...

我按照您之前的帖子中使用lable命令模拟交变电场，有一个问题是温度控制的程序怎么编写？
进入lable后温度控制的步数好像是按照lable中的run命令来执行的，我模拟的结果就是升温会非常快。
可以加您的联系方式吗？论坛的交流效率有点低

wildon

hbfy321 发表于 2023-5-15 23:21
我按照您之前的帖子中使用lable命令模拟交变电场，有一个问题是温度控制的程序怎么编写？
进入lable后温 ...

加电场后，温度（粒子运动速度）是由电场控制的。
升温很快可能与场强太大有关，另外频率也有影响。
私聊我联系方式。

懒羊羊当大厨

请问楼主是否解决了问题，我也遇到了相似的问题，但我更希望电场对温度有影响