Confined colloidal bilayers under shear: steady state and relaxation back to equilibrium

Crystalline bilayers of charged colloidal suspensions which are confined between two parallel plates and sheared via a relative motion of the two plates are studied by extensive Brownian dynamics computer simulations. The charge-stabilized suspension is modeled by a Yukawa pair potential. The unsheared equilibrium configuration is two crystalline layers with a nested quadratic in-plane structure. For increasing shear rates (.)gamma, we find the following steady states: First, up to a threshold of the shear rate, there is a static solid which is elastically sheared. Above the threshold, there are two crystalline layers sliding on top of each other with a registration procedure. Higher shear rates melt the crystalline bilayers and even higher shear rates lead to a reentrant solid stratified in the shear direction. This qualitative scenario is similar to that found in previous bulk simulations. We have then studied the relaxation of the sheared steady state back to equilibrium after an instantaneous cessation of shear and found a nonmonotonic behavior of the typical relaxation time as a function of the shear rate (.)gamma. In particular, application of high shear rates accelerates the relaxation back to equilibrium since shear-ordering facilitates the growth of the equilibrium crystal. This mechanism can be used to grow defect-free colloidal crystals from strongly sheared suspensions. Our theoretical predictions can be verified in real-space experiments of strongly confined charged suspensions.