[lammps-users] question about fix rigid


I want to do simulation about liquid crystal phase transition with LAMMPS. Object is composed of atoms, which are treated as an independent rigid body with fix rigid command. But I want to use NPT ensemble, How can I do? System only includes thousands of rigid body, not other atoms. Is there other methods to

Thanks in advance.

Fix npt or fix nph work with rigid bodies. The harder
part is defining what temperature is. I've had more
luck thermostatting with fix langevin for rigid bodies.
So you may want to do that and use fix nph.