We have also assessed the nonlinear model through implicit-solvent, coarse-grained simulations in a wide range of pH values and salt concentrations, parameters that are of great relevance when preparing these types of suspensions in experiments. The simulations consisted of a single diamond-lattice microgel containing weak acid sites in a simulation box whose size was determined by the suspension concentration. We constructed the microgel model using the Kremer-Grest bead-spring model for the polymer and the restrictive primitive model for the electrostatic interactions. Conformational degrees of freedom were sampled using MD with the Langevin thermostat. The chemical reactions and the coupling to a reservoir of constant pH and constant salt concentration were achieved using the grandreaction method. In contrast to previous works, this method allowed us to accurately simulate microgels even for extreme pH-values.