At the beginning of the simulations, energy minimization was performed to eliminate any particle overlaps. After that, the system pressure and temperature were equilibrated to 1 atmosphere pressure and 300 K temperature using the NPzT ensemble by employing the Nosé-Hoover thermostat and the Nosé-Hoover barostat, respectively. Subsequently, the system was further equilibrated in an NVT ensemble (300 K) with a Nosé-Hoover thermostat. The damping parameters for the thermostat and barostat were set as 0.1 ps and 1 ps, respectively. The PPPM (particle-particle particle-mesh) solver was employed to calculate the long-range Coulomb interactions with a relative accuracy of 10 -4 and with adjustments for the slab geometry. A time step of 2 fs was used, and integration was done using the Velocity Verlet algorithm. After the system reached equilibrium, the production run was started to understand the system's behavior when no electric field was applied (E=0 V/nm). The production run was carried out in an NVT ensemble for 20 ns using a Nosé-Hoover thermostat. A detailed description of the equilibration time as well as the condition to determine equilibration (i.e., by checking the variation of the brush height with time) have been discussed in section S2 of the SI. These descriptions are similar to those provided in our previous articles.