the behavior of the microgel close to a surface, we model the latter
as two layers of wall particles that are initially located on a compact
square lattice of size σ at the bottom of the simulation box;
the separation between layers is 0.7σ. To avoid crystallization
of the monomers close to the surface, the wall particles are randomly
displaced from the lattice sites, including the direction perpendicular
to the plane, following a Gaussian distribution with standard deviation
σsd = 0.2. The obtained layers are then subsequently
fixed throughout the whole simulation runs.
Microgel monomers
interact with wall particles via the WCA potential (
the monomer–monomer attraction αmm with the
monomer–surface attraction, αms, now encoding
the surface hydrophobicity.
To mimic experimental conditions
where the microgel is physically
anchored to the wall, we also consider the case where permanent bonds
between a few monomers and the surface particles are formed. This
is obtained by the following procedure (illustrated in
αms = 0): (i) a swollen microgel (equilibrated in
bulk at αmm = 0) is pushed toward the wall; (ii)
when it comes in contact with the surface, monomers with distance
less than dz = 1.5σ from the upper layer of
the wall are considered, and among them; (iii) b monomers
are randomly chosen and anchored to their closest wall-particle via
the harmonic potential V(r) = K(r – r0)2 with K = 15 and r0 = 21/6σ; finally (iv) the microgel
is left to relax to its equilibrium state. The procedure is then repeated
for different surface αms conditions. We did it for
both the hydrophilic αms = 0 and hydrophobic αms = 0.9 conditions, yielding different bonding patterns. Density
profiles calculated with microgels anchored to a hydrophilic surface
are comparable to experiments in all cases, except for the measurements
at T = 35 °C close to a hydrophobic surface,
where we found that the extension of the tail is better captured by
simulations of microgels initially anchored to a hydrophobic surface
(see
to take into account is the number of bonds b that
we should consider in the simulations. As mentioned in
to experimental data. For the hydrophobic case, we expect in experiments
a much larger number of bonds due to the additional attraction to
the surface in the anchoring procedure. For this reason, we performed
all simulations (that take much longer, also due to the long aging
regime) with a fixed value b = 200, roughly 1 order
of magnitude difference with respect to the hydrophilic case. This
value was then found to be in rather good agreement with experiments
in terms of the tail of the density profiles.