A snapshot of the system in contact with the confining wall is shown in Fig. REF .
The initial dimensions of the simulation box are \((17\sigma _0,17\sigma _0,35\sigma _0)\)  with 7424 solution particles.
The average concentration of solutes in the whole volume of \(\bar{c}_s =0.15\) .
The box is periodic in the \(x\)  and \(y\)  directions. In the \(z\)  direction, there is a solid wall at the bottom and a moving surface at the top, where particles undergo specular reflection.
As the tangential momentum of particles remains unchanged upon reflection, this wall imposes “slip” boundary conditions.
In contrast, lower surface consists of a layer of immobile solid atoms with the structure of the (100) surface of a face-centred cubic (FCC) lattice with lattice constant \(\sqrt{2}\sigma _0\) . The interaction parameters of the solutes with the wall are \((\varepsilon _{sw},\sigma ^{LJ}_{sw})=(1.5,1.5)\) .
We used a Nosé-Hoover thermostat [1]} to fix \(k_BT/\varepsilon _0=1.0\)  for all the simulations.
<FIGURE>