atoms lost even i used fix wall/reflact command

Thanks for quick reply,

I equilibrated crystal at 300K with periodic boundary condition. now I used last dump configuration of equilibrated crystal and want to heat up to around 1500K. but during this, I used fix wall/refelct command.

boundary f f p

fix 11 all wall/reflect xlo EDGE xhi EDGE ylo EDGE yhi EDGE units lattice

I can use that command if and only if I have fix boundary condition in a certain direction on which I want to apply that command. The problem is when I start my simulation, it shows that my atoms are lost.

sorry i reply personally.

From: Mathur Nil <[email protected]...>
Date: Mon, Jul 17, 2017 at 7:05 PM
Subject: Re: [lammps-users] atoms lost even i used fix wall/reflact command
To: Axel Kohlmeyer <[email protected]>

Thanks for quick reply,

I equilibrated crystal at 300K with periodic boundary condition. now I used
last dump configuration of equilibrated crystal and want to heat up to
around 1500K. but during this, I used fix wall/refelct command.

boundary f f p

fix 11 all wall/reflect xlo EDGE xhi EDGE ylo EDGE yhi EDGE units lattice

I can use that command if and only if I have fix boundary condition in a
certain direction on which I want to apply that command. The problem is when
I start my simulation, it shows that my atoms are lost.

if you want to have fixed boundaries, why don't you start with them right away?
you are making your life needlessly complicated. when changing
boundary conditions, you will need to re-equilibrate, so your previous
periodic equilibration is mostly a wasted effort.

also, in your original e-mail you said that you wanted to keep x- and
y- dimensions constant and let the system expand in z-direction.
this is _not at all_ what you are inducing by making the changes you
did. your changes set up a completely different system, i.e. rod that
is replicated in z-directions (but still with fixed box dimensions).

why don't you just keep all periodic boundaries and apply fix npt and
couple it only to the z direction?

regardless of whether you have fixed boundary conditions or periodic
boundary conditions, the box dimensions are unchanged unless you issue
a command or fix that can change them. perhaps you need to read up
some more about what boundary conditions mean in LAMMPS.

axel.