Suppose a benzene ring is connected with copper nanorod on both side along the z axis and they compose the scattering region. On both sides of the scattering region along the z axis, the left and right leads are both iron nanorod.
I want to simulate the thermal conductance through the Non-Equilibrium Molecular Dynamics (NEMD) with LAMMPS.
Would anyone please tell me how to divide the whole structure into several boxes, which could contain the scattering region, for the 'fix heat' method calculation? How to set up the parameter?
Thank you in advance.