where Fi, E, and qi respectively represent the force acting on particle i, the electric field (E), and the partial charge on particle i. The EOS flow was probed for three cases, namely for E=0.1 V/nm, 0.5 V/nm, 1 V/nm. In MD simulations, the high electric field is always applied to suppress thermal fluctuation and get insight into the system's behavior in a reasonable computational time. However, the electric field in this range has been applied experimentally and also considered in previous simulations. During the flow simulations, we applied a thermostat to the water and ion molecules exclusively within one-fourth of the nanochannel, to remove the heat dissipated during flow [see Fig. ]. Also, we removed the flow velocity bias, which ensured that correct thermodynamic temperature was used by the thermostat. Furthermore, the polymers were thermostated using a Langevin thermostat throughout the nanochannel according to the thermal wall model. Such a step ensured that one could mimic the experimental conditions where the temperature of the channel was controlled through the surface. The steady-state condition was reached when both the velocity profile and the average brush height (one of the slowest responding parameters) became constant (see S2). In summary, both the equilibrium and non-equilibrium simulation procedures were like those employed in our previously published articles as well as those in studies done by other researchers. All simulations were conducted using the LAMMPS software package. OVITO was utilized for visualization and illustration purposes. The initial structure was created using the MOLTEMPLATE package.