Dear Developers and Users,
I would like to incorporate the grand canonical Monte Carlo method in order to maintain a fixed chemical potential in a system consisting of small rigid bodies influenced by langevin dynamics. To demonstrate what I am trying to accomplish, I am attaching an input file which constitutes a simple diatomic system. This particular system is able to exchange molecules, but once molecules have been inserted into the box the relevant dynamics do not appear to take place. I would appreciate insight into what I can do to get this to work.
Thanks!
Carlos
dimer.mol (84 Bytes)
test.in (1.45 KB)