Studying synthesis confinement effects on the internal structure of nanogels in computer simulations
We study the effects of droplet finite size on the structure of nanogel particles synthesized by random crosslinking of molecular polymers diluted in nanoemulsions. For this we use a bead-spring computer model of polymer-like structures that mimics the confined random crosslinking process corresponding to irradiation- or electrochemically-induced crosslinking methods. Our results indicate that random crosslinking under strong confinement can lead to unusual nanogel internal structures, with a central region less dense than the external one, whereas under moderate confinement the resulting structure has a denser central region. We analyze the topology of the polymer networks forming nanogel particles with both types of architectures, their overall structural parameters, their response to the quality of the solvent and compare the cases of nonionic and ionic systems.
keywords:Nanogels, randomly crosslinked polymer networks, computer simulations
Nearly years have passed since the concept of microgel appeared for the first time as a new type of polymer systemBaker1949 (), in which a gel—i.e., a diluted network of permanently crosslinked polymers—is synthesized with the volume of a colloidal particle—i.e., particles with typical sizes ranging from few tens of nanometers to several micrometers 2011-fernandez-nieves-bk (); 2017-hamzah-jpr (). These soft particles have very interesting properties derived from their internal semiflexible network structure. For instance, they can be penetrated by other smaller particles or molecules and frequently exhibit a very strong response to external stimuli, being able to shrink or swell with a change of several times their average volume. This makes them promising building blocks of smart materials, which has greatly stimulated in recent years the research on their synthesis techniques, experimental characterization and theoretical modeling. Nowadays, it is possible to create colloidal gel particles responsive to a large variety of stimuli, including temperature, pH of the solvent or external fields Hoare2004 (); Gorelikov2004 (); Zhang2004a (); Das2006 (); Yin2008 (); Bayliss2011 (); 2011-malmsten-ch (); Holmqvist2012 (); Mohanty2015 (); Kobayashi2016 (); Backes2017a (); Backes2017b (); Colla2018 (). Prospective applications are as diverse as targeted drug delivery, oil spill recovery or sensing and smart coating technologies Raquois1995 (); Saatweber1996 (); RETAMA2003 (); Guo2005 (); LOPEZ2005 (); 2011-malmsten-ch (); Sivakumaran2011 (); 2011-thorne-cps (); Hu2012 (); Bonham2014 (); Son2016 (); Aliberti2017 (); Pepe2017 (); Schimka2017 (); Agrawal2018 (); Alhuraishawy2018 (); PU2019 (). Even the term ‘microgel’ is still applied oftenly to gel particles of any size up to roughly 100 micrometers, a conventional distinction between nanogels—referred to particles with sizes from 1 to 100 nanometers—and microgels—naming only particles from 0.1 to 100 micrometers—is becoming widespread 2007-aleman-pac (). This is favored by the fact that some important applications strictly require the use of nanometric gel particles. For instance, acting as drug delivery systems, only nanogels can overcome the blood-brain barrier Vinogradov2010 ().
To date, one of the most fruitful aspects of the research on colloidal gel particles is the large amount of techniques developed for their synthesis. The most common approaches are based on the polymerization of monomers in solution in presence of crosslinking agents. The polymerization and crosslinking may take place either in homogeneous solutions, when the newly formed colloidal particles are insoluble and precipitate, or inside droplets of solution emulsions, usually in presence of surfactant agents 2011-fernandez-nieves-bk (); 2017-hamzah-jpr (). In the case of polymerization/crosslinking in emulsions, the droplets act as a finite size confinement for the reactions, determining the final size of the colloidal gel particle. In this way, both nano- or microgel particles can be obtained from the same polymerization/crosslinking process by only selecting the size of the emulsion droplets. Analogously, inter- and intramolecular crosslinking into colloidal particles can be also performed on already formed polymer molecules diluted in homogeneous solutions or emulsions. In most cases molecular crosslinking is performed by means of photoinduced creation of radicals or electrochemical methods 2016-mavila-chrv (); 2016-galia-eccm (). These latter approaches have the advantange of leaving aside the use of free monomers and crosslinking agents, whose residual presence after the synthesis of the colloidal particles might prevent their use for some applications, particularly in biomedicine 2008-kwon-oh-pps (). The current development of novel synthesis approaches aims mainly at controlling the characteristic size and dispersity of the particles, their stability and the internal distribution of every type of functional groups. Important examples of such novel techniques are the ones based on microfluidic devices Zhang2007 (); Shah2008 (); Tumarkin2009 () or the use of microporous membranes for a highly controlled emulsification of the precursor solutions Charcosset2004 (). The scaling down of these techniques, with the use of nanofluidic devices and emulsification with nanoporous membranes, opens up the possibility of creating nanogels with a higher control on their structure and more sophisticated properties.
Despite the significant progress achieved in synthesis techniques, the size scales involved in micro- and nanogel structures still represent a challenge for the accurate characterization and modeling of their properties. In one hand, direct experimental measurements are mainly limited to overall particle properties, as for example their size, solvent content, internal distribution of contrast agents or mechanical response 2018-backes-pol (), whereas the detailed internal structure of the polymer network is still extremly challenging for direct measurements. Quantitative estimations of the degree of crosslinking of the network are usually deduced from the synthesis conditions, whereas direct observations based on conventional microscopy techniques only provide indications on the internal relative distribution of crosslinkers and other functional groups 2006-hoare-jpcb (). Among conventional measurement methods, neutron and X-ray small angle scattering (SANS/SAXS) are the only ones able to provide statistical information on the internal structure and overall shape of the particles, represented as form factors. Even by combining different microscopy and scattering methods, the level of characterization detail is limited 2019-witte-sm (). Only very recently it has been possible to obtain detailed topological information on the structure of polymer networks from labelling based super-resolution microscopy measurements2016-gelissen-nnl (); 2018-siemes-angw (); 2018-bergmann-pccp (); 2018-karanastasis-mh (). However, application of such novel technique to microgel characterization is still very scarce, whereas in the case of nanogels it still remains, to our best knowledge, unexplored. On the other hand, computer simulations are the most extensively used approach for the modeling of micro- and nanogels, being a powerful theoretical tool for the search of an accurate connection between their nanoscopic structure and their overall properties. However, these systems are complex enough to make their full explicit atomistic modeling unfeasible. Therefore, computer modeling relies on a coarse-grained representation of the polymer network, usually based on bead-spring models Grest86 (). Simplest bead-spring models of nano- and microgels are based on regular lattices and have been extensively used for studying swelling properties and thermosensitivity of microgels, especially of that ones based on charged polymers Kobayashi2014 (); Ghavami2016 (); Kobayashi2016 (); Ahualli2017 (); Ghavami2017 (); Kobayashi2017 (); Hofzumahaus2018 (). However, since the actual topology of the experimental polymer networks in nano- and microgels is expected to be rather complex, several efforts to move away from regular lattice architectures in order to capture such complexity have been made in last years. One of such approaches is the one presented by Gnan and co-workers Gnan2017 (). Their model uses the self-assembly properties of patchy colloids in order to mimic the polymerization/crosslinking process of monomers diluted in droplet emulsions Gnan2017 (), generating microgel disordered structures. This model has been used to study the internal structure and swelling behaviour of microgels depending on their degree of crosslinking and size of the confining droplet Rovigatti2018 (), as well as the effects of explicit solvent-polymer interactions 2018-camerin-scrp (). Another model, aimed at a realistic representation of the polymer network of microgels synthesized by crosslinking of polymer molecules under confinement, has been also introduced very recently by Moreno and Lo Verso Moreno2018 (). In this case, their model assumes the presence of macromolecular precursors that have been functionalized for their intermolecular crosslinking at spots with prescribed distributions, obtaining a network with a low fraction of crosslinks that has faster deswelling kinetics than equivalent microgels with regular lattice structure. Further details on the progress of numerical modeling of colloidal gel particles can be found in two very recent reviews 2019-rovigatti-sm-rev (); 2019-martin-molina-jml-rev ().
Independently from the approaches mentioned above, recently we introduced a model for colloidal gel particles that also intends to mimic the realistic crosslinking of polymer molecules confined in droplets Minina2018 (). In difference with the model of Moreno and Lo Verso, in our case no localized functionalization of the polymer precursors diluted inside the droplet is assumed and polymer monomers are crosslinked randomly according to their proximity after equilibration. This corresponds to an experimental photonically or electrochemically induced crosslinking process, which, to our best knowledge, has not been represented by any other computer model to date. In our first work we analyzed the effects of the presence of magnetic nanoparticle inclusions on the structure of the polymer network. Here, we use such model to study conventional, nonmagnetic systems.
In this work we focus on the computer simulation study of the internal structure of simple nanogel particles synthesized by random crosslinking of polymer molecules diluted in nanodroplets, analyzing how it is affected by the relative size of the confining droplet, the quality of the solvent and the use of either non charged polymers or polyelectrolytes as precursors. To our best knowledge, this combination of aspects has not been studied to date, neither theoretically nor experimentally. For instance, the effects of the size of the confining droplets on the synthesis of gel microparticles have been addressed in several experimental works, particularly for microfluidic synthesis approaches Zhang2007 (); Tumarkin2009 (); Crassous2015 (); DiLorenzo2016 (). In computer simulations, as pointed above, such aspect has been also addressed only for gel microparticles obtained from different synthesis routes Rovigatti2018 (); Moreno2018 (). However, it is reasonable to expect the effects of confinement on the crosslinking process to become more important as its size decreases, having a much higher impact on the internal structure of gel nanoparticles than in microgels. Our simulations support this hypothesis, providing indications of a non trivial dependence of the structure on the relative droplet size when it compares to the contour length of the molecular precursors. As a consequence of such effects we found that, while keeping the typical overall characteristics of colloidal gel particles, nanogels may exhibit internal structures different from the ones predicted in previous works for microgels.
The paper is organized as follows: in next Section we introduce the simulation model and protocol; next we discuss the simulation results, first considering the case of nonionic systems under good solvent conditions, followed by the analysis of the effects of solvent quality and electrostatic interactions for the case of ionic systems; finally, we conclude with a summary of results and outlook.
2 Nanogel model and simulation method
In order to model the structure of nanogel particles formed from a random crosslinking of polymer molecules diluted in a nanodroplet, we adopt a coarse-grained approach based on a bead-spring representation of the polymers Grest86 (), simple interaction potentials and a system of reduced units. As usual in this type of approach, the polymers are modelled as linear chains of monomers, that are simple spherical beads with reduced mass and reduced characteristic diameter . The excluded volume interactions of such monomers are represented by a truncated and shifted Lennard-Jones pair potential, also known as Weeks-Chandler-Andersen (WCA) potential 1971-weeks ():
where is the conventional Lennard-Jones potential and is the center-to-center distance between the interacting beads. Under good solvent conditions, the truncation distance is set to in order to make the interaction purely repulsive. We take the energy scale of this interaction, , as well as the energy of the thermal fluctuations in the system, , as unity, . When studying poor solvent conditions, the cut-off distance for potential (1) acting on any pair of beads is changed to in order to make it attractive, setting a varying attraction strength .
The bonds between adjacent monomers along the polymer chains are represented by finitely extensible nonlinear elastic (FENE) springs Grest86 (),
for what we take as the maximal bond extension and as the bond strength. The purpose of these choices is simply to keep an average distance between bonded nearest neighbors close to unity under all the conditions we sampled, preventing bond crossing while not requiring a very small integration timestep.
During the crosslinking process, the confining droplet is represented as a fixed sphere of radius that keeps the polymers inside its volume. This is achieved by means of the potential (1) acting also between the polymer beads and the confining sphere, being in this case the distance from its surface to the center of the bead. Crosslinkers are simply introduced as additional bonds established between pairs of beads that were not originally bonded neighbors along the polymer chains. This is done according to the protocol described below. As a trick to improve computational efficiency, we use a different type of bonding potential for crosslinking. This is a simple harmonic bond,
By taking , we checked that the mean length of the harmonic bonds after relaxation is around with a standard deviation of , independently of the fraction of crosslinked monomers and the size of the nanogel, being basically equivalent to the average length of the FENE bonds. Therefore, with our choice of parameters, FENE and harmonic bonds have the same structural effects after the crosslinking process. The reason for using a different potential for the crosslinkers will be explained below, when describing the crosslinking protocol.
The case of nanogels formed by polyelectrolytes is considered by assigning electric charges to the polymer beads and introducing the Debye-Hückel electrostatic pair potential between them 1923-debye-hueckel-pz (),
where is the center-to-center distance between a pair of beads 1 and 2, and their respective electric charges, the Debye screening wave vector, the cut-off distance of the interaction and the prefactor is defined, as is usual in this type of electrostatic calculation approach, in terms of the distance at which the electrostatic potential between two elementary charges compares to the thermal fluctuations, or Bjerrum length , so that . The Debye-Hückel model is a reasonable approximation for the calculation of electrostatic interactions in polyelectrolyte systems under moderate screening, i.e., with a low degree of ion condensation conditions, and has been used to study the swelling behavior of nanogel particles with regular lattice internal structure Kobayashi2016 ().
We choose to study nanogel particles by means of equilibrium molecular dynamics simulations in the NVT ensemble with a velocity Verlet integrator. In order to avoid the costly explicit simulation of solvent molecules, we use a Langevin thermostat. With this method, the effects of the friction and the thermal fluctuations produced by the solvent are represented implicitly by friction and stochastic terms introduced in the translational and rotational Newtonian equations of motion 1987-allen (); 2002-frenkel (). Since here we are only interested in equilibrium properties, we are free to choose any value for the friction constants, that we set to unity in reduced units, whereas the stochastic terms are set to satisfy the usual fluctuation-dissipation rules.
Our simulation protocol is the following. First, we fix the confining sphere representing the nanodroplet in a simulation box with open boundaries. We place inside the sphere equivalent polymer chains of beads each. The radius of the confining sphere is taken according to the characteristic size of an unconstrained polymer chain of such length,
where is the Flory exponent. We perform an initial relaxation of the polymer chains inside the sphere, letting them to interact only through the excluded volume and confinement repulsions defined from potential (1), by integrating for steps with a time step of , where is the time scale in our system of reduced units, . This is followed by the crosslinking process, that essentially consists in the selection of pairs of beads from different polymer chains to be bonded with potential (3), so that the final crosslinking fraction, defined as the ratio between the amount of crosslinks and the amount of beads, is . Such selection is made by randomly picking from a list of all pairs of particles belonging to different chains that are separated by a distance not larger than an arbitrary cutoff. Initially, we set such cutoff distance to but, whenever it is necessary, we slowly increase it until the desired amount of crosslinks is reached. At this point, the use of a harmonic potential instead of a FENE potential for the crosslinking becomes clear: it allows to establish bonds with lengths that, initially, can be slightly larger than the maximum extension of the FENE bonds but, after equilibration, tend to relax to the same equilibrium length. This trick largely simplifyes the mimicking of the random crosslinking process. To limit the crosslinking to interchain bonds is a simple way to ensure the connectivity of the structure and is a reasonable assumption at least for moderately dense dilutions of polymers, which keep a high degree of mixing under confinement Jun2007 (), and/or for radiation-induced polymer crosslinking processes at room temperature 2011-an-pol (). The case when good mixing conditions are not satisfied will be discussed below. Once the crosslinks are established, the system is let to relax for integration steps. Finally, the confining sphere is removed and the crosslinked object, with the full set of interactions, is equilibrated and its properties are measured in a run of integration steps. The results presented here were obtained by averaging over 50 independent runs, i.e., over 50 different network configurations obtained from the crosslinking process, subsequently equilibrated for each set of parameters. The model and simulation protocol were implemented using the ESPResSo 3.2 simulation package arnold13a ().
3 Results and Discussion
We sample a set of parameters specifically chosen to analyze whether confinement effects can be important in the structure and properties of nanogels made by random intermolecular crosslinking of polymer precursors. For the crosslinking process we consider two different solutions consisting of monodisperse polymer chains with either or beads each, so that the resulting crosslinked particle will consist of either or beads in total. According to the criterium described above, the radius of the sphere representing the confining nanodroplet is therefore taken as and , respectively, corresponding to volume fractions of polymer beads of for the smaller system and for the larger one. Since the key point we want to address here is the effect of the degree confinement during crosslinking and not of the size of the system, in the following we will refer to the case as the system crosslinked under relatively strong confinement conditions (SCC), and to the case as that one crosslinked under moderate confinement conditions (MCC). Regarding the equivalence to physical quantities in both cases, assuming for example that beads in our model represent Kuhn lengths of the polymer chains, these values could correspond to poly-vinylpyrrolidone (PVP) molecules with contour lengths of 300 and 600 nm 2010-knappe-pol () diluted in droplets with diameters of approximately 45 and 68 nm, respectively. Such droplet sizes are comparable to that ones of nanogel particles synthesized by pulse electron irradiation- or electrochemically-induced crosslinking of PVP molecules 2011-an-pol (); 2016-galia-eccm (). However, we should underline that here we aim at addressing a fundamental question rather than at the accurate modeling of a particular system, thus the chosen values could also be scaled to represent other systems of polymers with different Kuhn lenghts and droplet sizes. Therefore, only ratios of the sampled values are important in the following discussion.
Experimental random crosslinking techniques make more difficult the estimation of the degree of crosslinking present in the resulting nanogel particles than in the case of conventional chemical crosslinking methods. Typical values of the fraction of crosslinks with respect to the total amount of monomers provided by the latter techniques in micro- and nanogels are around , but values in a range as broad as to have been reported Bonham2014 (). Since we aim at a qualitative discussion, we chose to sample three different values from such range, from very low to strong: , and . We assume all these values to be reachable by radiation- and electrochemically-induced crosslinking methods.
In the next sections we first discuss the case of nonionic systems under good solvent conditions, analyzing the topology of the polymer networks and their structural parameters.
We start our discussion by characterizing the structural properties of the polymer networks obtained from our random interchain crosslinking procedure. Figure 1 shows some examples of typical equilibrium configurations of the networks obtained under good solvent conditions, corresponding to the non ionic case and each of the sampled crosslinking fractions. As expected, one can see that the overall structure becomes more compact and spherical as the fraction of crosslinks increases. For the lowest crosslink fraction, the structure is so loose that it can hardly be considered an actual soft particle, whereas for and the configurations look similar to the ones obtained from other computer models. Therefore, with this set of crosslinks fractions we can compare the structure of a loose, small polymer network with the ones of nanogel particles.
Unlike in other simulation models, and despite the apparent resemblance in the resulting overall structure, our polymer networks have no predefined linear segment distribution. Therefore, it is interesting to analyze such parameter in our networks. A linear segment with length is defined as a part of the polymer network consisting of a sequence of bonded but not crosslinked beads that are bounded by crosslinked beads and/or chain free ends. Since we have different bonds for the polymer backbones and the crosslinks, we can compute easily by simply counting the number of FENE bonds in each bounded segment. The relative probability distributions of are presented in Figure 2 for networks obtained under strongly confined crosslinking and in Figure 2 for that ones obtained under moderate confinement.
As one can expect, for each system type the distributions become more steep and narrow as the fraction of crosslinks increases, leading to a lower average segment length. By comparing both system types, we can see that the relative distributions of are broader for SCC systems, showing larger relative average values, . Regarding the precursor lengths used in each case, by doubling the length of the polymer backbones we obtain only an increase of the average segment length of approximately for the least crosslinked system, for the intermediately crosslinked system and less than for the most crosslinked one. This indicates a weak dependence of on the length of the polymer precursors, that tends to vanish as the fraction of crosslinks increases.
The reason for the weak dependence of on is not obvious and one may wonder whether it is the result of a confinement effect. In order to analyze this, in Figure 3 we present several reduced number density profiles of different equilibrated systems still under confinement, measured from the center of the confining sphere. Figure 3 shows the density profiles of polymer beads obtained for both sampled confinement conditions right before the crosslinking takes place. A qualitative difference can be observed: while the MCC system shows a rather uniform distribution, in the more strongly confined a maximum of density is found close to the boundary of the confining sphere, indicating that the polymers tend to occupy the external region of the droplet rather than the center. This is the signature of the system being in a highly concentrated regime, as defined for polymers under spherical confinement Jun2007 (). Under such conditions individual chains can segregate at no penalty in their free energy and concentrate close to the sphere boundaries. Having a significantly lower density, the moderately confined system does not display such behavior. This difference cannot but affect the distribution of crosslinks. This is evidenced by Figures 3 and 3, that show the number density profiles of crosslinks for each system size and fraction of crosslinks. Without regard of the fraction of crosslinks, we can observe that in the MCC system (Figure 3) most crosslinks are located in the center of the sphere, with a relatively flat middle region that decays when approaching the boundary. However, in the SCC system such preference for the central region is absent and a small maximum can be observed close to the boundary.
Figure 4 provides a further indication of the rather different structure adopted by the polymers under distinct confinement conditions. It shows two configuration snapshots corresponding to the point in our simulation protocol right before removing the spherical confinement. To ease the visualization, only three polymer chains with different colored backbones from each system are depicted. By comparing the SCC (Figure 4) and MCC (Figure 4) cases, we can see that the longer polymer chains of the latter system have a more uniform distribution insde the cavity than the shorter chains of the former one. In fact, shorter chains tend to be more aligned with the confinement wall.
At this point, the weak dependence of the linear segment length on can be explained by the different configurations adopted by the polymer chains when the crosslinking is performed. If chains in a solution of polymers of length are perfectly mixed, from simple probability considerations one can estimate that the average linear segment length after randomly assigning interchain crosslinks is approximately
The values provided by this expression for the MCC system are, for increasing fraction of crosslinks, and , which are very close to the ones measured from our configurations (see values in caption of Figure 2). However, for the SCC system, expression (6) significantly underestimates the values of actually obtained, indicating that in this case the polymer chains were not completely mixed during crosslinking. Chain segregation in a dense system, even when it is only partial, makes individual chains to form blob-like structures. Such structures necessarily decrease the amount of nearly close contacts between different chains with respect to a system of perfectly mixed chains. Therefore, the random interchain crosslinking procedure based on a distance criterium will involve only chain segments at the interfaces of the blobs, leaving that parts in the inner region of the blobs without crosslinks. Compared to a well mixed case, this makes the distribution of more disperse, with more long non-crosslinked segments, increasing its average value. In summary, we found two manifestations of a significant effect of strong confinement conditions on the internal structure of randomly crosslinked polymer networks and nanogel particles: the increase of the average linear segment length with respect to moderate confinement conditions and, more importantly, the adoption of a structure slightly more dense and crosslinked in the external region of the droplet than in the center.
The confinement effect discussed above is interesting and might be relevant for nanogel particles synthesized under analogous conditions. However, one may argue that this effect is only an artifact of our approximations, since here we are considering exclusively interchain crosslinking, whereas intrachain crosslinking should be significant at least for the more tightly confined system, in which polymer mixing is poor. However, intrachain crosslinking can not prevent this effect to manifest, at least in the peculiar profile of the strongly confined structure. The reasoning behind this affirmation is the following. First, in any case some degree of interchain crosslinking is necessary in order to obtain nanogel particles with larger molecular weight than each individual precursor polymer, being single molecule particles out of the scope of this work. Second, intrachain crosslinking will be still selected only by the close distance criterium, not being affected by the formation of segregated chain blobs. Therefore, the distribution of polymer beads and crosslinks will be essentially the same, with the only difference that the latter will be a mixture of both kinds of crosslinks, inter- and intrachain ones.
Another topological parameter of interest that may evidence the effect of confinement in the structure of the crosslinked networks is the average shortest path between crosslinks, , defined as
where is the shortest linear segment that connects the crosslinked beads and . Figure 5 shows , normalized by the polymer backbone length, as a function of the fraction of crosslinks. We can see that the average shortest path decays with , as it is expected. Interestingly, is approximately times smaller for the MCC case than for the SCC one. This result also supports the existence of a more homogeneously crosslinked network in the less confined system than in the more strongly confined one.
3.2 Gyration radius and asphericity
Once we examinated the differences in the topology of the polymer networks led by the confinement conditions during crosslinking, we analyze larger scale structural properties of these systems once the confinement is removed.
The radius of gyration, , and asphericity, , are standard scalar parameters used to characterize the overall structure of soft colloidal particles. The radius of gyration provides the characteristic size of the particle, whereas the asphericity is a non-negative parameter that tends to zero as the shape of the particle approaches a perfect sphere. Both quantities can be computed from the eigenvalues of the gyration tensor, (whith ), in the following way:
Figure 6 shows the average values of the reduced radius of gyration, , and the average asphericity normalized by the square radius of gyration, , obtained for both system types, as a function of the fraction of crosslinks. In Figure 6 we can observe that, as expected, the size of the particle decreases in all cases when the fraction of crosslinks increases: decays by approximately when the value of is increased by 10 times, making the polymer network more compact. Such compaction is also clearly accompanied by the adoption of a more spherical overall profile (see Figure 6). This general behavior is qualitatively consistent with theoretical predictions and experimental measurements of nano- and microgel particles. Importantly, while the case of lowest fraction of crosslinks remains arguable, the compactness and low asphericity of the structures corresponding to the intermediate and the largest fractions of crosslinks supports that, despite their small size, these crosslinked networks can be considered actual nanogel particles.
3.3 Density profiles
At this point we are interested in checking whether the characteristic density profiles observed under different confinement conditions (Figures 3 and 3) remain once the networks and nanogel particles are equilibrated without confinement. For this, we compute the corresponding reduced number density profiles, , measured in this case as a function of the distance to the center of mass of the system. Figure 7 shows the results obtained for this parameter, plotted for distances normalized with the values of radius of gyration discussed above.
We can observe that, for both system types, profiles of nanogel particles ( or ) exhibit a plateau-like region for distances and a relatively abrupt drop at larger distances. For , however, the drop starts at smaller distances and is significantly smoother. This is due to the relative lack of average internal structure of the latter, especially when compared to nanogel particles. Importantly, the dependence of the central region of the profile on the crosslinking confinement not only persists here for the whole final unconfined system but, in some cases, it is enhanced: for MCC nanogels there is a clear maximum of density at their centers of mass, whereas for SCC ones such central region tends to have a density similar (for high fraction of crosslinks) or lower (for intermediate fraction of crosslinks) than the main outer region. Even the statistical fluctuations at the central region are relatively large, the significant difference between the SCC and MCC cases is clear. For SCC with , one can also observe a weak indication of layering in the intermediate region of the particle. All these aspects could be important for the design of nanogel particles aimed at applications that require or could benefit from specific internal structuring.
3.4 Poor solvent conditions
Structural collapse of colloidal gels, from soft particles with fuzzy boundaries to compact solid-like particles with sharp edges, may happen as a response to certain external stimulus, as for instance changes in solvent quality 2018-keidel-sa ().
Here we analyze the effects of solvent quality in the structure of our crosslinked networks by making attractive the isotropic soft core interaction (1) between polymer beads. Figure 8 shows configuration snapshots that provide two examples of the structural evolution as the strength of the attraction, , is increased. They correspond to SCC systems, with the lowest and highest sampled fractions of crosslinks. In general, the change from a good solvent ( repulsive, left column in Figure 8) to a poor one (, middle column) leads to the formation of more compact structures, with smaller characteristic size. By further increasing to reach a very strong attraction—i.e., tending to the limit of no-solvent conditions Rubinstein03 ()—the beads are forced to minimize the volume of the structures they form, eventually collapsing into a compact arrangement (, right column in Figure 8). In our simulations, the transition to a collapsed structure has been observed for in the range from to . Figure 8 also shows the existence of a particular case in the way the networks tend to collapse: whereas in all other cases we observed that the structures tend to become compact spherical particles (lower row in Figure 8) as reported in numerous studies, in the case of the polymer network formed with the lowest fraction of crosslinks (upper row in Figure!8) the collapsed structure is very anisometric, approaching a cylindrical shape, as a consequence of the low amount of crosslinks and their highly irregular distribution. This is another indication that such loose networks do not fulfill the common characteristics of nanogel particles.
The qualitative structural evolution observed in the snapshots can be better characterized by means of the radius of gyration and the asphericity. The results obtained for these parameters as a function of the attraction strength are presented in Figures 9 and 10, respectively. Specifically, Figure 9 shows the average reduced radius of gyration, , relative to the value corresponding to good solvent conditions, , for different fractions of crosslinks. From these plots one can conclude that decays abruptly as soon as solvent quality worsens. This decrease tends to saturate once the attraction strength becomes more important than the thermal fluctuations, . In general, this relative decay is larger for MCC system, reflecting its initial lower overall density. In addition, the strongest decay is observed for the most loosely crosslinked system, as a consequence of the strong nonuniform shrinking from a rather expanded network.
The effect of the internal structure also manifests itself in the normalized average asphericity, , shown in Figure 10. In agreement with the discussion above on the snapshots of Figure 8, in all cases except for the most loosely crosslinked SCC system the asphericity drops around to become very small for , indicating the rather spherical shape of the particles under such conditions. The exceptional case corresponding to a strongly anisometric collapsed configuration is reflected in a large growth of asphericity with (Figure ).
3.5 Ionic systems
An important part of nano- and microgels created to date consist of charged polymers that introduce large structural responses to variations in the ionic strength of the solvent, particularly strong swelling and deswelling behaviors. These ionic systems have been the subject of numerous theoretical studies, especially by means of computer simulations based on regular lattices Claudio2009 (); Kobayashi2014 (); Ghavami2016 (); Kobayashi2016 (); Ahualli2017 (); Ghavami2017 (); Kobayashi2017 (); Hofzumahaus2018 (). Here we also perform a preliminary study of the properties of our model nanogel particles and networks when synthesized from polyelectrolyte precursors, using the Debye-Hückel approximation (4). We stress that here we are interested only in a qualitative analysis of the structural changes led by electrostatic interactions and their dependence on the different network topologies produced by distinct confinement conditions during synthesis. Therefore, we sample reduced parameters aimed at underlining such changes rather than at the accurate representation of the system. Specifically, here we assign to each polymer bead a charge of and sample values of the electrostatic interaction strength, , in the range for two values of the inverse Debye screening length, and .
Figure 11 includes configuration snapshots of systems analogous to the ones presented in Figure 1, that show a clear swelling due to a strong, slightly screened electrostatic repulsion between the beads. This leads to a stiffening of the linear segments in the polymer networks, causing the overall expansion of the structure. By comparing both sets of snapshots, we can see that loosely crosslinked networks () show a very strong expansion, loosing any eventual particle-like shape. However, despite still experiencing a considerable swelling, nanogel particles () keep their nearly spherical shape.
Figure 12 illustrates in more detail, also by means of the relative average radius of gyration, the swelling experienced by systems of each type and fraction of crosslinks as the strength of the electrostatic interactions increases. As expected, the largest relative expansion (above within the sampled range) happens for the system that was proven to be structurally the weakest—i.e., the SCC with lowest —under the lower screening conditions, (see Figure 12, top row). The smallest expansion corresponds to MCC systems with largest and screening (around , see lower row in Figure 12). Also important is that, confirming the observation above, the average asphericity of the nanogel particles (not shown) does not change significantly with their electrostatically induced swelling.
Finally, we check the changes in the internal structure due to the electrostatic swelling by comparing the density profiles of the non ionic systems, already introduced in Figure 7, with the ionic ones. An example of this comparison, corresponding to the case of strong electrostatic interactions and screening, is shown in Figure 13. We can see that, besides the significant broadening of the profiles, the main qualitative features observed for non ionic systems remain in their ionic counterparts: a decay in the external region of the profiles, more or less smooth depending on the fraction of crosslinks; a plateau-like intermediate region and a central region where the differences in the internal topology led by the original confinement conditions manifest. Interestingly, the drop in density in the central region that corresponds to SCC nanogels with intermediate fraction of crosslinks is significantly enhanced by the electrostatic swelling, showing that the internal structure is basically preserved despite the general expansion of the interstitial regions. In turn, this type of expansion also makes the maximum in density near the center of mass of the larger nanogel particles to become much less prominent. Therefore, we can conclude that electrostatic interactions basically tend to enhance the particular features of the internal structure of nanogels obtained under strongly confined moderate crosslinking.
3.6 Form factor measurements
As a last result, we check whether the differences in the internal structure we found in our crosslinked networks and nanogel particles depending on the confinement conditions might be detected in scattering experiments. Note that here we do not aim at a formal characterization of small angle scattering data profiles, but only at providing a qualitative hint on the experimental feasibility of such analysis.
The component that contains the information corresponding to single particles in scattering measurements data is known as form factor. SANS measurements are frequently employed to deduce the form factors of micro- and nanogels particles, usually by means of fitting the experimental data to a given predefined model of the internal structure.2014-schneider-lm (); 2016-felberg-mph (); 2019-witte-sm () In simulation data one can calculate form factors directly from the coordinates of the elements of the particle. In our case, the following expression can be applied to each pair of polymer beads in the network, and :
where is the center-to-center distance between them, is the length of the scattering wave vector and the angular brackets stand for an averaging over different configurations. Figure 14 shows a selection of form factors corresponding to ionic systems under electrostatic swelling conditions ( and ), as these are the ones that exhibit the largest quantitative differences in their internal structure. In these plots, the length of the wave vector is rescaled by the corresponding radius of gyration to check that the Guinier regime is established at .Guinier1955 () This corresponds to the plateau observed for . The difference in the structure of the polymer networks has to manifest in the large region. For such region we set the limit , so that we do not reach shorter length scales than that of a single polymer bead diameter, . By comparing the form factors corresponding to moderately and strongly confined crosslinking conditions, in the region of large one can observe quantitative differences in their profiles, particularly for systems of loosely crosslinked networks. The latter show a rather smooth profile that is frequently observed in experimental measurements as a consequence of the smearing produced by a large configurational entropy in the system. In our case, this smearing comes from the averaging over independent realizations. The profiles corresponding to the intermediate and the largest fraction of crosslinks, however, show few clear maxima. This indicates that the internal structure of these nanogel particles, even being formed in a random crosslinking process, keeps a rather consistent distribution among independent realizations. This relatively low internal configurational variation is also favored by the small size of the particles. Therefore, these results suggest that the effects of the confinement conditions during random crosslinking of loose polymer networks and nanogel particles actually might be detected in small angle scattering experimental measurements.
By means of a coarse-grained, bead-spring computer model of polymers, we studied the effects of strong and moderate confinement conditions on the internal structure and properties of nanogel particles synthesized by random crosslinking of polymer molecules diluted in nanodroplet emulsions. Random crosslinking, that corresponds to current experimental irradiation- or electrochemically-induced crosslinking methods, makes the resulting polymer networks to not have an architecture directly predefined by their chemical inital components. Instead, our results indicate that such internal structure can be influenced by the relative size of the confining droplet.
Our simple random crosslinking model is able to produce polymer networks that, except for a very low fraction of crosslinks, show the characteristic behavior of actual nanogel particles, including structurally stable nearly spherical shapes. In order to understand the role of confinement before and during crosslinking, we compared two scenarios: first, a moderate density of polymer beads inside the confining droplet, or moderate confinement conditions, that makes polymer molecules tend to mix well and have a typical density profile, with a higher density in the center of the droplet, a relatively homogeneous intermediate region and a strong decay at the edge; second, a strong confinement due to a high density of polymer beads, that makes polymers tend to segregate and occupy preferently the region near the confinement boundary, leaving a region of slightly lower density in the center. We observed that random moderate crosslinking preserves such internal structural features in the resulting nanogel particles, which possess different network topologies as revealed by their distributions of linear segments and average shortest paths between crosslinked spots. In particular, a more inhomogeneous distribution of crosslinks is obtained for the crosslinking under stronger confinement conditions.
A qualitative similar compaction and final collapse into compact spherical structures as the quality of the solvent worsens has been also observed in both types of model nanogel particles. Only a collapse into very anisometric structures has been found in the extreme case of polymer networks obtained by the combination of a low fraction of crosslinks and a relatively high inhomogeneity in their distribution due to strongly confined crosslinking conditions.
The usual changes in overall size associated to the swelling/deswelling behavior due to electrostatic interactions have been observed when we consider the nanogel particles to consist of charged polymers. However, the impact on the internal structure is different depending on the topology of the polymer networks. Whereas the density profiles of nanogels synthesized under moderate confinement conditions tend to flatten with the electrostatic swelling, in the ones corresponding to nanogels obtained under strong confinement and moderate degree of crosslinking the central region of lower density tends to be enhanced.
In all aspects analyzed here, the most interesting features have been observed for a moderate value of fraction of crosslinks. This is due to the fact that such parameter controls the tradeoff between the structural stability, that makes the network to have a persistent shape and behave as a soft spherical particle, and the existence of internal structural inhomogeneities, that can be useful or even essential for given practical applications.
In summary, these results suggest the control of confinement conditions, for instance by means of membrane emulsification techniques, as a simple approach to tune qualitatively the structure of nanogel particles using random crosslinking methods, whereas the fraction of crosslinks can be used to tune quantitatively their properties.
As a final outlook comment, in our ongoing research we extend the current model to consider other sources of internal heterogeneities, able to produce well differentiated core-shell nanogel structures.
E. S. M. and S. S. K. acknowledge the financial support of the Austrian Science Fund (FWF): START-Projekt Y627-N27. S. S. K. and C. N. L. acknowledge the support of ETN-COLLDENSE (H2020-MSCA-ITN-2014, Grant No. 642774). P. A. S. is supported by the Act 211 of the Government of the Russian Federation, contract No. 02.A03.21.0006. Computer simulations were performed at the Vienna Scientific Cluster (VSC). We thank Eugenia Kumacheva for fruitful discussions concerning the experimental relevance of our model.
- (1) W. O. Baker, Microgel, a new macromolecule, Ind. Eng. Chem. 41 (3) (1949) 511–520. doi:10.1021/ie50471a016.
- (2) A. Fernandez-Nieves, H. M. Wyss, J. Mattsson, D. A. Weitz (Eds.), Microgel Suspensions: Fundamentals and Applications, WILEY-VCH Verlag & Co. KGaA, Weinheim, Germany, 2011. doi:10.1002/9783527632992.
- (3) Y. B. Hamzah, S. Hashim, W. A. W. A. Rahman, Synthesis of polymeric nano/microgels: a review, J. Polym. Res. 24 (9) (2017) 134. doi:10.1007/s10965-017-1281-9.
- (4) T. Hoare, R. Pelton, Highly pH and temperature responsive microgels functionalized with vinylacetic acid, Macromolecules 37 (7) (2004) 2544–2550. doi:10.1021/ma035658m.
- (5) I. Gorelikov, L. M. Field, E. Kumacheva, Hybrid microgels photoresponsive in the near-infrared spectral range, J. Am. Chem. Soc. 126 (49) (2004) 15938–15939. doi:10.1021/ja0448869.
- (6) J. Zhang, S. Xu, E. Kumacheva, Polymer microgels: Reactors for semiconductor, metal, and magnetic nanoparticles, J. Am. Chem. Soc. 126 (25) (2004) 7908–7914. doi:10.1021/ja031523k.
- (7) M. Das, H. Zhang, E. Kumacheva, Microgels: Old materials with new applications, Annu. Rev. Mater. Res. 36 (1) (2006) 117–142. doi:10.1146/annurev.matsci.36.011205.123513.
- (8) J. Yin, D. Dupin, J. Li, S. P. Armes, S. Liu, pH-induced deswelling kinetics of sterically stabilized poly(2-vinylpyridine) microgels probed by stopped-flow light scattering, Langmuir 24 (17) (2008) 9334–9340. doi:10.1021/la8014282.
- (9) K. Bayliss, J. S. van Duijneveldt, M. A. Faers, A. W. P. Vermeer, Comparing colloidal phase separation induced by linear polymer and by microgel particles, Soft Matter 7 (2011) 10345–10352. doi:10.1039/C1SM05917C.
- (10) M. Malmsten, Microgel Suspensions: Fundamentals and Applications, WILEY-VCH Verlag & Co. KGaA, 2011, Ch. 15, Microgels in Drug Delivery, pp. 375–405. doi:10.1002/9783527632992.ch15.
- (11) P. Holmqvist, P. S. Mohanty, G. Nägele, P. Schurtenberger, M. Heinen, Structure and dynamics of loosely cross-linked ionic microgel dispersions in the fluid regime, Phys. Rev. Lett. 109 (2012) 048302. doi:10.1103/PhysRevLett.109.048302.
- (12) P. S. Mohanty, P. Bagheri, S. Nöjd, A. Yethiraj, P. Schurtenberger, Multiple path-dependent routes for phase-transition kinetics in thermoresponsive and field-responsive ultrasoft colloids, Phys. Rev. X 5 (2015) 011030. doi:10.1103/PhysRevX.5.011030.
- (13) H. Kobayashi, R. G. Winkler, Universal conformational properties of polymers in ionic nanogels, Sci. Rep. 6 (2016) 19836. doi:10.1038/srep19836.
- (14) S. Backes, P. Krause, W. Tabaka, M. U. Witt, D. Mukherji, K. Kremer, R. von Klitzing, Poly(n-isopropylacrylamide) microgels under alcoholic intoxication: When a LCST polymer shows swelling with increasing temperature, ACS Macro Lett. 6 (10) (2017) 1042–1046. doi:10.1021/acsmacrolett.7b00557.
- (15) S. Backes, P. Krause, W. Tabaka, M. U. Witt, R. von Klitzing, Combined cononsolvency and temperature effects on adsorbed PNIPAM microgels, Langmuir 33 (50) (2017) 14269–14277. doi:10.1021/acs.langmuir.7b02903.
- (16) T. Colla, P. S. Mohanty, S. Nöjd, E. Bialik, A. Riede, P. Schurtenberger, C. N. Likos, Self-assembly of ionic microgels driven by an alternating electric field: Theory, simulations, and experiments, ACS Nano 12 (5) (2018) 4321–4337. doi:10.1021/acsnano.7b08843.
- (17) C. Raquois, J. Tassin, S. Rezaiguia, A. Gindre, Microgels in coating technology: structure and rheological properties, Prog. Org. Coat. 26 (2) (1995) 239–250. doi:10.1016/0300-9440(95)00573-0.
- (18) D. Saatweber, B. Vogt-Birnbrich, Microgels in organic coatings, Prog. Org. Coat. 28 (1) (1996) 33–41. doi:10.1016/0300-9440(95)00582-X.
- (19) J. R. Retama, B. Lopez-Ruiz, E. Lopez-Cabarcos, Microstructural modifications induced by the entrapped glucose oxidase in cross-linked polyacrylamide microgels used as glucose sensors, Biomaterials 24 (17) (2003) 2965–2973. doi:10.1016/S0142-9612(03)00095-4.
- (20) Z. Guo, H. Sautereau, D. E. Kranbuehl, Structural evolution and heterogeneities studied by frequency-dependent dielectric sensing in a styrene/dimethacrylate network, Macromolecules 38 (19) (2005) 7992–7999. doi:10.1021/ma050788b.
- (21) V. C. Lopez, J. Hadgraft, M. Snowden, The use of colloidal microgels as a (trans)dermal drug delivery system, Int. J. Pharm. 292 (1) (2005) 137–147. doi:10.1016/j.ijpharm.2004.11.040.
- (22) D. Sivakumaran, D. Maitland, T. Hoare, Injectable microgel-hydrogel composites for prolonged small-molecule drug delivery, Biomacromolecules 12 (11) (2011) 4112–4120. doi:10.1021/bm201170h.
- (23) J. B. Thorne, G. J. Vine, M. J. Snowden, Microgel applications and commercial considerations, Colloid Polym. Sci. 289 (5) (2011) 625. doi:10.1007/s00396-010-2369-5.
- (24) Y. Hu, Q. Wang, J. Wang, J. Zhu, H. Wang, Y. Yang, Shape controllable microgel particles prepared by microfluidic combining external ionic crosslinking, Biomicrofluidics 6 (2) (2012) 026502. doi:10.1063/1.4720396.
- (25) J. A. Bonham, M. A. Faers, J. S. van Duijneveldt, Non-aqueous microgel particles: synthesis, properties and applications, Soft Matter 10 (2014) 9384–9398. doi:10.1039/C4SM01834F.
- (26) H. A. Son, S. K. Choi, E. S. Jeong, B. Kim, H. T. Kim, W. M. Sung, J. W. Kim, Microbial activation of bacillus subtilis-immobilized microgel particles for enhanced oil recovery, Langmuir 32 (35) (2016) 8909–8915. doi:10.1021/acs.langmuir.6b02010.
- (27) A. Aliberti, A. Ricciardi, M. Giaquinto, A. Micco, E. Bobeico, V. La Ferrara, M. Ruvo, A. Cutolo, A. Cusano, Microgel assisted lab-on-fiber optrode, Sci. Rep. 7 (1) (2017) 14459. doi:10.1038/s41598-017-14852-5.
- (28) A. Pepe, P. Podesva, G. Simone, Tunable uptake/release mechanism of protein microgel particles in biomimicking environment, Sci. Rep. 7 (1) (2017) 6014. doi:10.1038/s41598-017-06512-5.
- (29) S. Schimka, N. Lomadze, M. Rabe, A. Kopyshev, M. Lehmann, R. von Klitzing, A. M. Rumyantsev, E. Y. Kramarenko, S. Santer, Photosensitive microgels containing azobenzene surfactants of different charges, Phys. Chem. Chem. Phys. 19 (2017) 108–117. doi:10.1039/C6CP04555C.
- (30) G. Agrawal, R. Agrawal, Stimuli-responsive microgels and microgel-based systems: Advances in the exploitation of microgel colloidal properties and their interfacial activity, Polymers 10 (4) (2018) 418. doi:10.3390/polym10040418.
- (31) A. K. Alhuraishawy, B. Bai, M. Wei, Combined ionically modified seawater and microgels to improve oil recovery in fractured carbonate reservoirs, J. Petrol. Sci. Eng. 162 (2018) 434–445. doi:10.1016/j.petrol.2017.12.052.
- (32) J. Pu, J. Geng, P. Han, B. Bai, Preparation and salt-insensitive behavior study of swellable, Cr3+-embedded microgels for water management, J. Mol. Liq. 273 (2019) 551–558. doi:10.1016/j.molliq.2018.10.070.
- (33) J. V. Alemán, A. V. Chadwick, J. He, M. Hess, K. Horie, R. G. Jones, P. Kratochvíl, I. Meisel, I. Mita, G. Moad, S. Penczek, R. F. T. Stepto, Definitions of terms relating to the structure and processing of sols, gels, networks, and inorganic-organic hybrid materials (IUPAC Recommendations 2007), Pure Appl. Chem. 79 (2007) 1801–1829. doi:10.1351/pac200779101801.
- (34) S. V. Vinogradov, Nanogels in the race for drug delivery, Nanomedicine 5 (2) (2010) 165–168. doi:10.2217/nnm.09.103.
- (35) S. Mavila, O. Eivgi, I. Berkovich, N. G. Lemcoff, Intramolecular cross-linking methodologies for the synthesis of polymer nanoparticles, Chem. Rev. 116 (3) (2016) 878–961. doi:10.1021/acs.chemrev.5b00290.
- (36) A. Galia, S. Lanzalaco, M. A. Sabatino, C. Dispenza, O. Scialdone, I. Sirés, Crosslinking of poly(vinylpyrrolidone) activated by electrogenerated hydroxyl radicals: A first step towards a simple and cheap synthetic route of nanogel vectors, Electrochem. Commun. 62 (2016) 64–68. doi:10.1016/j.elecom.2015.12.005.
- (37) J. K. Oh, R. Drumright, D. J. Siegwart, K. Matyjaszewski, The development of microgels/nanogels for drug delivery applications, Prog. Polym. Sci. 33 (4) (2008) 448–477. doi:10.1016/j.progpolymsci.2008.01.002.
- (38) H. Zhang, E. Tumarkin, R. M. A. Sullan, G. C. Walker, E. Kumacheva, Exploring microfluidic routes to microgels of biological polymers, Macromol. Rapid Commun. 28 (5) (2007) 527–538. doi:10.1002/marc.200600776.
- (39) R. K. Shah, J.-W. Kim, J. J. Agresti, D. A. Weitz, L.-Y. Chu, Fabrication of monodisperse thermosensitive microgels and gel capsules in microfluidic devices, Soft Matter 4 (2008) 2303–2309. doi:10.1039/B808653M.
- (40) E. Tumarkin, E. Kumacheva, Microfluidic generation of microgels from synthetic and natural polymers, Chem. Soc. Rev. 38 (2009) 2161–2168. doi:10.1039/B809915B.
- (41) C. Charcosset, I. Limayem, H. Fessi, The membrane emulsification process—a review, J. Chem. Technol. Biotechnol. 79 (3) (2004) 209–218. doi:10.1002/jctb.969.
- (42) S. Backes, R. Von Klitzing, Nanomechanics and nanorheology of microgels at interfaces, Polymers 10 (9) (2018) 978. doi:10.3390/polym10090978.
- (43) T. Hoare, D. McLean, Kinetic prediction of functional group distributions in thermosensitive microgels, J. Phys. Chem. B 110 (41) (2006) 20327–20336. doi:10.1021/jp0643451.
- (44) J. Witte, T. Kyrey, J. Lutzki, A. M. Dahl, J. Houston, A. Radulescu, V. Pipich, L. Stingaciu, M. Kühnhammer, M. U. Witt, R. von Klitzing, O. Holderer, S. Wellert, A comparison of the network structure and inner dynamics of homogeneously and heterogeneously crosslinked PNIPAM microgels with high crosslinker content, Soft Matter 15 (2019) 1053–1064. doi:10.1039/C8SM02141D.
- (45) A. P. H. Gelissen, A. Oppermann, T. Caumanns, P. Hebbeker, S. K. Turnhoff, R. Tiwari, S. Eisold, U. Simon, Y. Lu, J. Mayer, W. Richtering, A. Walther, D. Wöll, 3d structures of responsive nanocompartmentalized microgels, Nano Lett. 16 (11) (2016) 7295–7301. doi:10.1021/acs.nanolett.6b03940.
- (46) E. Siemes, O. Nevskyi, D. Sysoiev, S. K. Turnhoff, A. Oppermann, T. Huhn, W. Richtering, D. Wöll, Nanoscopic visualization of cross-linking density in polymer networks with diarylethene photoswitches, Angew. Chem. Int. Ed. 57 (38) (2018) 12280–12284. doi:10.1002/anie.201807741.
- (47) S. Bergmann, O. Wrede, T. Huser, T. Hellweg, Super-resolution optical microscopy resolves network morphology of smart colloidal microgels, Phys. Chem. Chem. Phys. 20 (2018) 5074–5083. doi:10.1039/C7CP07648G.
- (48) A. A. Karanastasis, Y. Zhang, G. S. Kenath, M. D. Lessard, J. Bewersdorf, C. K. Ullal, 3D mapping of nanoscale crosslink heterogeneities in microgels, Mater. Horiz. 5 (2018) 1130–1136. doi:10.1039/C8MH00644J.
- (49) G. S. Grest, K. Kremer, Molecular dynamics simulation for polymers in the presence of a heat bath, Phys. Rev. A 33 (1986) 3628–3631. doi:10.1103/PhysRevA.33.3628.
- (50) H. Kobayashi, R. G. Winkler, Structure of microgels with Debye-Hückel interactions, Polymers 6 (5) (2014) 1602–1617. doi:10.3390/polym6051602.
- (51) A. Ghavami, H. Kobayashi, R. G. Winkler, Internal dynamics of microgels: A mesoscale hydrodynamic simulation study, J. Chem. Phys. 145 (24) (2016) 244902. doi:10.1063/1.4972893.
- (52) S. Ahualli, A. Martín-Molina, J. A. Maroto-Centeno, M. Quesada-Pérez, Interaction between ideal neutral nanogels: A Monte Carlo simulation study, Macromolecules 50 (5) (2017) 2229–2238. doi:10.1021/acs.macromol.6b02333.
- (53) A. Ghavami, R. G. Winkler, Solvent induced inversion of core–shell microgels, ACS Macro Lett. 6 (7) (2017) 721–725. doi:10.1021/acsmacrolett.7b00318.
- (54) H. Kobayashi, R. Halver, G. Sutmann, R. G. Winkler, Polymer conformations in ionic microgels in the presence of salt: Theoretical and mesoscale simulation results, Polymers 9 (1) (2017) 15. doi:10.3390/polym9010015.
- (55) C. Hofzumahaus, P. Hebbeker, S. Schneider, Monte Carlo simulations of weak polyelectrolyte microgels: pH-dependence of conformation and ionization, Soft Matter 14 (2018) 4087–4100. doi:10.1039/C7SM02528A.
- (56) N. Gnan, L. Rovigatti, M. Bergman, E. Zaccarelli, In silico synthesis of microgel particles, Macromolecules 50 (21) (2017) 8777–8786. doi:10.1021/acs.macromol.7b01600.
- (57) L. Rovigatti, N. Gnan, E. Zaccarelli, Internal structure and swelling behaviour of in silico microgel particles, J. Phys.: Condens. Matter 30 (4) (2018) 044001. doi:10.1088/1361-648X/aaa0f4.
- (58) F. Camerin, N. Gnan, L. Rovigatti, E. Zaccarelli, Modelling realistic microgels in an explicit solvent, Sci. Rep. 8 (1) (2018) 14426. doi:10.1038/s41598-018-32642-5.
- (59) A. J. Moreno, F. Lo Verso, Computational investigation of microgels: synthesis and effect of the microstructure on the deswelling behavior, Soft Matter 14 (2018) 7083–7096. doi:10.1039/C8SM01407H.
- (60) L. Rovigatti, N. Gnan, L. Tavagnacco, A. J. Moreno, E. Zaccarelli, Numerical modelling of non-ionic microgels: an overview, Soft Matter 15 (2019) 1108–1119. doi:10.1039/C8SM02089B.
- (61) A. Martín-Molina, M. Quesada-Pérez, A review of coarse-grained simulations of nanogel and microgel particles, J. Mol. Liq. 280 (2019) 374–381. doi:10.1016/j.molliq.2019.02.030.
- (62) E. S. Minina, P. A. Sánchez, C. N. Likos, S. S. Kantorovich, The influence of the magnetic filler concentration on the properties of a microgel particle: Zero-field case, J. Magn. Magn. Mater. 459 (2018) 226–230. doi:10.1016/j.jmmm.2017.10.107.
- (63) J. J. Crassous, A. M. Mihut, L. K. Månsson, P. Schurtenberger, Anisotropic responsive microgels with tuneable shape and interactions, Nanoscale 7 (2015) 15971–15982. doi:10.1039/C5NR03827H.
- (64) F. Di Lorenzo, S. Seiffert, Effect of droplet size in acrylamide-based microgel formation by microfluidics, Macromol. React. Eng. 10 (3) (2016) 201–205. doi:10.1002/mren.201500061.
- (65) J. D. Weeks, D. Chandler, H. C. Andersen, Role of repulsive forces in determining the equilibrium structure of simple liquids, J. Chem. Phys. 54 (1971) 5237–5247. doi:10.1063/1.1674820.
- (66) P. Debye, E. Hückel, Zur theorie der elektrolyte. i. gefrierpunktserniedrigung und verwandte erscheinungen, Phys. Z. 24 (1923) 185–206.
- (67) M. P. Allen, D. J. Tildesley, Computer Simulation of Liquids, 1st Edition, Oxford Science Publications, Clarendon Press, Oxford, 1987.
- (68) D. Frenkel, B. Smit, Understanding molecular simulation, Academic Press, 2002.
- (69) S. Jun, A. Arnold, B.-Y. Ha, Confined space and effective interactions of multiple self-avoiding chains, Phys. Rev. Lett. 98 (2007) 128303. doi:10.1103/PhysRevLett.98.128303.
- (70) J.-C. An, A. Weaver, B. Kim, A. Barkatt, D. Poster, W. N. Vreeland, J. Silverman, M. Al-Sheikhly, Radiation-induced synthesis of poly(vinylpyrrolidone) nanogel, Polymer 52 (25) (2011) 5746–5755. doi:10.1016/j.polymer.2011.09.056.
- (71) A. Arnold, O. Lenz, S. Kesselheim, R. Weeber, F. Fahrenberger, D. Roehm, P. Košovan, C. Holm, ESPResSo 3.1 — molecular dynamics software for coarse-grained models, in: M. Griebel, M. A. Schweitzer (Eds.), Meshfree Methods for Partial Differential Equations VI, Vol. 89 of Lecture Notes in Computational Science and Engineering, Springer, 2013, pp. 1–23. doi:10.1007/978-3-642-32979-1_1.
- (72) P. Knappe, R. Bienert, S. Weidner, A. F. Thünemann, Characterization of poly(n-vinyl-2-pyrrolidone)s with broad size distributions, Polymer 51 (8) (2010) 1723–1727. doi:10.1016/j.polymer.2010.02.039.
- (73) R. Keidel, A. Ghavami, D. M. Lugo, G. Lotze, O. Virtanen, P. Beumers, J. S. Pedersen, A. Bardow, R. G. Winkler, W. Richtering, Time-resolved structural evolution during the collapse of responsive hydrogels: The microgel-to-particle transition, Sci. Adv. 4 (4) (2018) eaao7086. doi:10.1126/sciadv.aao7086.
- (74) M. Rubinstein, R. H. Colby, Polymer Physics, Oxford University Press, 2003.
- (75) G. C. Claudio, K. Kremer, C. Holm, Comparison of a hydrogel model to the Poisson–Boltzmann cell model, J. Chem. Phys. 131 (9) (2009) 094903. doi:10.1063/1.3207275.
- (76) F. Schneider, A. Balaceanu, A. Feoktystov, V. Pipich, Y. Wu, J. Allgaier, W. Pyckhout-Hintzen, A. Pich, G. J. Schneider, Monitoring the internal structure of poly(n-vinylcaprolactam) microgels with variable cross-link concentration, Langmuir 30 (50) (2014) 15317–15326. doi:10.1021/la503830w.
- (77) L. E. Felberg, A. Doshi, G. L. Hura, J. Sly, V. A. Piunova, W. C. Swope, J. E. Rice, R. Miller, T. Head-Gordon, Structural transition of nanogel star polymers with ph by controlling pegma interactions with acid or base copolymers, Mol. Phys. 114 (21) (2016) 3221–3231. doi:10.1080/00268976.2016.1224942.
- (78) A. Guinier, G. Fournet, Small angle scattering of X-rays, Wiley, New York, 1955. doi:10.1002/pol.1956.120199326.