The counterion condensation around the constituting chains of a PE hydrogel may be described by the counterion condensation theory for dilute PE solutions . The condensation model dictates that not all but some of the counterions are released to the surrounding solution while a fraction of the counterions condense near the chain backbone, effectively renormalizing the net of charge of the PE chain . Accordingly, if the line charge density of a chain is 1/ℓ 0 , the counterion condensation occurs if the ratio, also known as Manning parameter, obeys ξ = ℓ B /ℓ 0 > 1, where the Bjerrum length ℓ B characterizes the strength of the Coulomb coupling in a solvent with dielectric constant, ϵ. Nevertheless, the counterion condensation can be affected by the dialectic heterogeneities near the PE chains . Such heterogeneity can arise due to different ordering of water dipoles around the hydrophobic and hydrophilic interfaces of the PE chains as compared to the bulk . Consequently, in a deforming hydrogel, the individual chains can change their conformation, and in turn, such conformational changes of the constituting chains can affect the counterion distribution in a non-trivial way. Furthermore, as the hydrogel is deformed (i.e., compressed or extended), the individual chains can be collapsed or stretched to adjust the total volume of the hydrogel. As a result, the cumulative counterions condensation of the hydrogel can change depending on the ratio of chains undergoing various conformational transitions under deformation. Investigating such scenarios requires going beyond the continuum solvent models, and implicit solvent simulations since steric interactions and the dipole nature of water molecules surrounding the PE chains can strongly affect the condensation and, in turn, the electromechanical response of PE hydrogels under deformation.