I am performing an atomistic simulations regarding the lubrication where the liquid crystals are confined between two mica plates. So far I have the dynamics of liquid crystal when the upper plate is moving with a constant velocity. And everything works well. Now, I want to incorporate the dynamics of molecules inside the plates so that the energy dissipation through the plates can also be incorporated. Is there any easy way to implement the idea?