[lammps-users] Is there heat current leak on fixed boundary, when set the boundary atoms velocity to zero.

It sounds like you have a non-perioidic system. Fix thermal/conductivity
is set up for a periodic system. You could alter the code to do the
swaps between different layers to allow for non-periodicity.