Hello! I am having trouble understanding why I am receiving the error (ERROR on proc 0: Particle on or inside fix wall surface) in the code I am trying to run:
------------- Initialization
units lj
dimension 2
atom_style sphere
boundary f f p
pair_style lj/cut 33.67384588865670
pair_modify shift yes
------------- System definition
region testbox block -10 10 -10 10 -0.25 0.25
create_box 2 testbox
create_atoms 1 single 5.0 5.0 0
create_atoms 1 single 5.0 -5.0 0
create_atoms 2 single -5.0 -5.0 0
create_atoms 2 single -5.0 5.0 0
group 1 type 1
group 2 type 2
fix wallxhi all wall/lj93 xhi EDGE 1.0 1.0 0.1
fix wallxlo all wall/lj93 xlo EDGE 1.0 1.0 0.1
fix wallylo all wall/lj93 ylo EDGE 1.0 1.0 0.1
fix wallyhi all wall/lj93 yhi EDGE 1.0 1.0 0.1
fix 1 1 brownian 100 44830747 gamma_t 1.0 rng gaussian
fix 2 2 brownian 50000 34490088 gamma_t 1.0 rng gaussian
fix 3 1 viscous 0.0001884955592153876
fix 4 2 viscous 0.0001884955592153876
pair_coeff 1 1 0.00000000000431 30.00000000000000 33.67384588865670
pair_coeff 1 2 0.00000000107932 30.00000000000000 33.67384588865670
pair_coeff 2 2 0.00000000215433 30.00000000000000 33.67384588865670
dump mydmp all xyz 1 dump_(16_26_16).xyz
timestep 5e-06
run 100000
Output:
LAMMPS (29 Sep 2021 - Update 3)
Created orthogonal box = (-10.000000 -10.000000 -0.25000000) to (10.000000 10.000000 0.25000000)
1 by 1 by 1 MPI processor grid
Created 1 atoms
using lattice units in orthogonal box = (-10.000000 -10.000000 -0.25000000) to (10.000000 10.000000 0.25000000)
create_atoms CPU = 0.000 seconds
Created 1 atoms
using lattice units in orthogonal box = (-10.000000 -10.000000 -0.25000000) to (10.000000 10.000000 0.25000000)
create_atoms CPU = 0.000 seconds
Created 1 atoms
using lattice units in orthogonal box = (-10.000000 -10.000000 -0.25000000) to (10.000000 10.000000 0.25000000)
create_atoms CPU = 0.000 seconds
Created 1 atoms
using lattice units in orthogonal box = (-10.000000 -10.000000 -0.25000000) to (10.000000 10.000000 0.25000000)
create_atoms CPU = 0.000 seconds
2 atoms in group 1
2 atoms in group 2
Neighbor list info …
update every 1 steps, delay 10 steps, check yes
max neighbors/atom: 2000, page size: 100000
master list distance cutoff = 33.973846
ghost atom cutoff = 33.973846
binsize = 16.986923, bins = 2 2 1
1 neighbor lists, perpetual/occasional/extra = 1 0 0
(1) pair lj/cut, perpetual
attributes: half, newton on
pair build: half/bin/atomonly/newton
stencil: half/bin/2d
bin: standard
Setting up Verlet run …
Unit style : lj
Current step : 0
Time step : 5e-06
Per MPI rank memory allocation (min/avg/max) = 5.156 | 5.156 | 5.156 Mbytes
Step Temp E_pair E_mol TotEng Press
0 0 0.0023089621 0 0.0023089621 0.00013863795
ERROR on proc 0: Particle on or inside fix wall surface (…/fix_wall_lj93.cpp:95)
Last command: run 100000