Dear Lammps users and developers,
I am simulating a temperature gradient in a polymer slab and seeing how some gold atoms inside diffuse. I have a 3D box with height of about 170 angstrons.
So, I’ve put a top region (100 - 170 angs), that I call active-layer, with fix langevin for 450 K and a bottom region (10-70 angs) with fix langevin for 300K and fix nve for (10-170 angs).
No problem until now.
BUT, after some MD steps (~3x 10^6 steps) I would like to use fix tfmc for the active_layer alone, keeping fix langevin 300 K for bottom region and using fix nve for all the atoms EXCEPT that of active-layer.
Can I do that? Of course things work with no errors, but I feel I am mixing MD with MC and maybe that is very wrong conceptually.
I thought about use fix tfmc for two regions, but with two different temperatures, they would have different MC “times”.
And about using tad instead tfmc? Would it be better?
Some words of help?
Thanks in advance,
Fernanda S Teixeira