"From what I understand, in your system, you initially had a non-rigid
immobile sheet, i.e. the sheet was fixed in space but the atoms
constituting the sheet were not immobile"
No. atoms were immobile as they were not time integrated. e.g. let us
say sheet is made up of 2 type of atoms type 3 and type 4. Then
following commands were given
group nonMobile type 3 4
group mobile subtract all nonMobile
fix 1 mobile npt temp 300 300 100 z 1.0 1.0 1.0 dilate mobile
This gives pressure up to 6000 atm. Just tried delete_bonds, that was
also of no help. After delete bonds, system does expand a little bit
but still far from required density/pressure.
However
group nonMobile type 3 4
group mobile subtract all nonMobile
fix sheet nonMobile rigid single
fix 1 mobile npt temp 300 300 100 z 1.0 1.0 1.0 dilate mobile
results in rapid expansion of system, and reaches 'closest' to what
should be actual pressure. Only few hundred atm off!
so my main worry is the physics of my simulation. fix rigid discards
internal degree of freedom for the sheet, but as sheet is immobile and
only interaction sheet has with system is pairwise lj and coulombic
interaction, should it matter so much?
Now i want to relax my system by increasing volume and taking to
desired pressure,should i do it with just non time integrated sheet,
like
fix nvt mobile temp 300 300 100
or do it with fix rigid?
fix sheet nonMobile rigid single
fix nvt mobile temp 300 300 100
Also I am not sure why internal degree of freedom of a rigid body
shall influence any external particle? As per my understanding, motion
of each particle is governed by the force it feels.
What difference it should make if a immobile particle is having
internal strain? So if my sheet is not moving, rigid or not, results
shall be similar (whether they are right or wrong is separate), should
they not be?