A suitable model for molecular dynamics simulation studies of an acrylic acid (AA) polyelectrolyte gel considering explicit water molecules is presented. It consists of a molecule of crosslinker bound to four AA strands completely ionized. Radial distribution functions and nonbonding energy results calculated with both molecular and quantum mechanics methods show that a differentiated solvation occurs for monomer units located at different distances from the crosslinker molecule. The used calculation methodology and the found differentiated solvation could be very interesting for polyelectrolyte gel studies and their applications