LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
[lammps-users] Electric_field
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[lammps-users] Electric_field

From: Mohammad Izadi <izadi0511@...24...>
Date: Wed, 2 Aug 2017 21:44:23 +0430

Dear lammps users,


I’m using lammps/reax/c for simulating decomposition of an organic compound under a constant electric field. From outputs (I mean and specie) I recognized that the electric field is worked but I’m not sure that it’s exactly correct or not. Please look into my input file and make recommends if it needs. This is my input file:


echo            both

units                real

dimension 3

boundary  p p p

atom_style       charge


#read_restart    restart22

restart 500      restart11 restart22

pair_style         reax/c NULL

pair_coeff        * * ffield.reax.rdx C H O N

compute reax all pair reax/c

variable eb        equal c_reax[1]

variable ea        equal c_reax[2]

variable et         equal c_reax[9]

variable ew       equal c_reax[11]

variable efi       equal c_reax[13]

variable eqeq    equal c_reax[14]

#define a region for efield

region    elecf block 0.0 100.0 0.0 100.0 0.0 100.0 units box

group     elf region elecf

neighbor          2 bin

neigh_modify  every 10 delay 0 check no

velocity        all create 300 235485 mom yes rot yes

fix                    1 all nvt temp 300 300 25    ###########timestep*100=25

fix             2 all qeq/reax 1 0.0 10.0 1e-6 reax/c

fix             3 all reax/c/species 1 50 250 species.txt

fix             4 elf efield 2.0 0.0 0.0 region elecf

fix_modify      4 energy yes

thermo_style    custom   step    atoms   temp   press   vol   etotal   ke   pe   v_eb   v_ea  v_et   v_ew v_efi   v_eqeq

thermo          250

timestep           0.25

dump               1 all atom 30 dump.reax.rdx

dump               2 all xyz  30

run    30000

Thanks in advance  

Mohammad Ebrahim Izadi