计算化学公社

标题: 求助：Reaxff中模拟周期电场下热解甲苯问题 [打印本页]

作者
Author: hbfy321 时间: 2023-5-15 11:27
标题: 求助：Reaxff中模拟周期电场下热解甲苯问题
参考文献：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)

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

作者
Author: wildon 时间: 2023-5-15 15:10
fix 2 all efield 0 0 ${1e} 的意思是 fix 2 all efield 0 0 0，${1e}这个不是变量，而是取函数的第一个值。
所以实际上是加了一个恒定电场。

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

作者
Author: hbfy321 时间: 2023-5-15 23:21

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命令来执行的，我模拟的结果就是升温会非常快。
可以加您的联系方式吗？论坛的交流效率有点低

作者
Author: wildon 时间: 2023-5-16 19:29

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

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

作者
Author: 懒羊羊当大厨 时间: 2024-11-11 23:04
请问楼主是否解决了问题，我也遇到了相似的问题，但我更希望电场对温度有影响

欢迎光临计算化学公社 (http://bbs.keinsci.com/)