Spatial segregation of mixed-sized counterions in dendritic polyelectrolytes.

Langevin dynamics simulations are utilized to study the structure of a dendritic polyelectrolyte embedded in two component mixtures comprised of conventional (small) and bulky counterions. We vary two parameters that trigger conformational properties of the dendrimer: the reduced Bjerrum length, [Formula: see text], which controls the strength of electrostatic interactions and the number fraction of the bulky counterions, [Formula: see text], which impacts on their steric repulsion. We find that the interplay between the electrostatic and the counterion excluded volume interactions affects the swelling behavior of the molecule. As compared to its neutral counterpart, for weak electrostatic couplings the charged dendrimer exists in swollen conformations whose size remains unaffected by [Formula: see text]. For intermediate couplings, the absorption of counterions into the pervaded volume of the dendrimer starts to influence its conformation. Here, the swelling factor exhibits a maximum which can be shifted by increasing [Formula: see text]. For strong electrostatic couplings the dendrimer deswells correspondingly to [Formula: see text]. In this regime a spatial separation of the counterions into core-shell microstructures is observed. The core of the dendrimer cage is preferentially occupied by the conventional ions, whereas its periphery contains the bulky counterions.

Introduction

Dendritic polyelectorolytes such as PAMAM (polyamidoamine) are extensively studied from the experimental and the theoretical point of view. Under physiological conditions only the terminal groups of PAMAMs bear positive charges. As the solution pH is lowered, the branching groups become protonated as well[1]. Ionization of the monomers provides a supply of counterions in solution, which has a tremendous effect on the structure and dynamics of the molecules. In this respect, absorption of both monovalent and divalent counterions[2] into the pervaded volume of dendrimers and counterion condensation are of particular importance. These phenomena are challenging problems which are tackled through different methods including the Donnan, the Poisson-Boltzmann and the Manning–McGhee–von Hippel models[2,3] as well as computer simulations[4-16]. Despite intensive studies of charged dendrimers in various environments and the role of counterions, still little is known about the effect of ion specificities on these polymers. In particular, an important factor is the size of counterions which is typically assumed to be either identical to the size of monomers or negligible[3,5-16]. Such an assumption is in contrast to several investigations which clearly demonstrated that the counterion size is a crucial control parameter. For instance, the behavior of polyelectrolyte gels in reaction to the size of counterions was examined using the virial expansion approximation[17]. It was recognized that the ion size has a substantial effect on the swelling ratio. Counterion specificity was also incorporated into the theory of phase transitions in polyelectrolyte gels[18]. It was demonstrated that small counterions generate volume collapse at low salt concentrations. Moreover, dissipative dynamics simulations were performed to study the conformational properties of polyelectrolyte chains accompanied by small and bulky counterions with various possibilities of the charged bead location[19]. Here, the main focus was on polymer collapse and formation of ionic microstructures. Molecular dynamics simulations were also employed to investigate the influence of the counterion steric effects on the structure of strongly charged polyelectrolytes in a dilute solution[20]. The authors reported formation of a dense globule in the presence of counterions of the size of monomers. For bulky counterions with the size much larger as compared to the monomer size, the polymers adopt extended conformations. Finally, for solutions containing a mixture of counterions of different size, formation of a core–shell globule takes place with the smaller counterions located within the globular core and the bulky counterions forming a shell on the globule surface. Some insights into the phenomena of ion selectivity between different counterions in the mixture was also gained for charged gels[21-24], polyelectrolyte solutions[20,25], the electric double layer around DNA[26] and rodlike polyelectrolytes[27-30]. Last but not least, in our recent study we identified that the counterion size considerably affects conformations of charged dendrimers in a way that they swell to a larger extent in the presence of bulky counterions[31].

Inspired by these studies, in this paper using Langevin dynamics simulations we continue to explore the structural properties of dendritic polyelectrolytes with the emphasis on the role of mixed-size counterions. We present the equilibrium properties of a charged dendrimer in a wide range of Coulomb interaction strength characterized by the reduced Bjerrum length, . In addition, we focus on the impact of the counterion excluded volume induced by varied amounts of conventional and bulky counterions. For this purpose the composition of the ionic mixture is controlled by the number fraction of the bulky counterions, . We find that both parameters and strongly affect conformations of charged dendrimers.

The paper is organized as follows. In “Model and method” section we outline the model and the simulation method. Our results are presented and discussed in “Results” section. We draw conclusions and remarks in “Summary” section.

Model and method

Model

We perform Langevin dynamics simulations of G4S4 charged dendrimer in a good, implicit solvent using the bead-spring model. In the above acronym, G represents dendrimer generation and S is the spacer length expressed in the number of bonds joining adjacent beads along a spacer chain. The tree-like structure of the molecule consists of the core of two bonded monomers and the trifunctional branching groups. The overall number of the monomers, N, and the branching groups (including the terminal groups), , in the dendrimer are given byand

We assign monovalent positive charges to the branching and the terminal groups so that the arrangement of charges corresponds to the architecture of PAMAM dendrimers at low pH[1]. In the following we interchangeably refer to the groups carrying charges as the charged monomers, the charged groups and the dendrimer ions. In order to maintain the overall system neutral, monovalent negative counterions are included in the model. To study the dendrimer with mixed-size counterions, and ions are considered as bulky and conventional, respectively, see Fig. 1. In the numerical scheme we incorporate three types of interactions. The excluded volume between pairs of particles of sizes, , , with their centers separated by a distance r, accounted for via the Lennard–Jones (LJ) 12-6 truncated and shifted potentialwhere is the interaction strength, and is the cutoff radius. The indices  = m, c, b, where m stands for monomers, c for conventional counterions and b for bulky ones, respectively. The choice of the cutoff results in purely repulsive LJ interactions mimicking good solvent conditions and, in Eq. (3), the mixing rule, , is used. The bonds joining any two adjacent beads which are described by the Finite Extensible Non-linear Elastic potential (FENE)where k is the spring constant and is the maximum extension of the bonds. The electrostatic interactions are introduced by pairwise Coulomb potentialwhere denotes the Boltzmann constant, T the absolute temperature, the Bjerrum length of the solvent, and r the distance between the centers of the ith and the jth particles with the charge valence , . With standing for the electric permittivity of the vacuum, for the relative permittivity of the solvent and e for elementary charge, the Bjerrum length is defined as . In particular, for water at room temperature ( K), , so that  Å.

Figure 1

2D scheme of G2S2 charged dendrimer at low pH. The charged and the uncharged monomers are marked respectively by the red and the blue circles. The conventional and the bulky counterions are depicted by the yellow and the green circles.

Method

The LAMMPS molecular dynamics package[32] was used to carry out the Langevin dynamics in the reduced LJ units. To set the unitless simulation parameters, the strength of the LJ potential, , and the monomer size, , were taken as the real units of energy and length respectively. The unitless simulation input parameters (denoted with an asterisk) were chosen as for the strength of the LJ potential, , for the size of the monomers and the conventional counterions, and for the mass of the monomers and the counterions. For the FENE potential the parameters were and [33]. According to the formula for reduced elementary charge, , and for the reduced temperature, , we set these parameters to the values and . Thus with the aforementioned energy unit, , our simulations were carried out at the reduced temperature which corresponds to room temperature. Coulomb interactions were calculated using the Particle–Particle–Particle Mesh (PPPM) method with the error tolerance for force and the real space cutoff radius [32]. The dumping parameter in the Langevin equation of motion was set to , where is the LJ time unit and m is the assumed real unit of mass. For example, for g/mol, ps. The calculations were performed with a time step . Note that the reduced Bjerrum length, , for a solvent at room temperature can be expressed as , which with our unit of length yields . This enabled us to vary by changing the relative electric permittivity of the solvent, . According to Eq. (5), is the parameter that determines Coulomb interaction energy between two monovalent charges at a distance apart related to the thermal energy at room temperature. Furthermore, a crucial parameter is the density of particles in the simulation box. In the paper we present the results obtained in the counterion-only limit, which are valid for systems of well-dispersed dendritic polyelectrolytes and counterions at finite density fixed through periodic boundary conditions of a cubic simulation box of length . The box size was large enough to prevent the dendrimer from intersecting with its periodic images, and allowed to perform simulations at experimentally accessible particle density: . Typically, about equilibration integration steps were performed and followed by integration steps of the production runs. During the production runs we collected the data on the scale up to hundreds of the longest relaxation times , where is defined as time needed for dendrimer to diffuse its own size. All G4 dendrimers considered in the current study had  . The snapshots of the monomers’ positions used in the analysis of the structural properties were saved every 1000 time steps. Throughout the paper our focus is on the conformational properties of the dendritic polyelectrolyte accompanied by binary mixtures of the conventional and the bulky counterions. The size of the bulky ions was , 3, and their number fraction was set to . The reduced Bjerrum length, , was varied between 0.125 and 16, covering the extremes of both weak and strong electrostatic interactions. In the following we refer to the counterions with a diameter as the conventional (counter)ions and with a diameter , 3 as the bulky ones.

The simulation snapshots were rendered using the visual molecular dynamics (VMD)[34]. All plots were prepared using XMGRACE[35].

Results

Spatial distribution of counterions

In Fig. 2 we display the pair correlation function and for the conventional and the bulky counterions respectively, calculated around the dendrimer’s center-of-mass. The profiles are presented as functions of the radial distance, , from the dendrimer’s center-of-mass rescaled by the dendrimer’s radius of gyration, , for different Bjerrum lengths, , and a given number fraction of the bulky counterions, . The data demonstrate that both types of ionic species accumulate at distances from the dendrimer’s core. The absorption effect is enhanced with increasing , which is manifested by an increase in the pair correlation function in this volume. The latter effect is attributed to stronger Coulomb attraction between the counterions and the dendrimer ions. Our finding confirms the previous numerical results, which however, were solely restricted to the case of conventional counterions[5-16]. As shown in Fig. 2, the pair correlation functions reveal a strong dependence on the ionic type. For the -profiles peak in a neighborhood of the center of mass, flatten in the molecule’s domain, and drop on its periphery. The -profiles exhibit a broad maximum in the domain and decrease near the center and on the periphery. Thus for sufficiently large -values, the dendrimer’s core is rich in the conventional counterions due to their weak excluded volume and Coulomb attraction with the charged monomers. On the other hand, on the periphery of the dendrimer an opposite tendency occurs due to more pronounced steric repulsion of the bulky ions with the dense core of the molecule, see Fig. 14. For the bulky counterions stronger excluded volume interactions prevent them from occupying the dendrimer’s core, which results in their accumulation at a large distance from it. Note that a similar preference of the counterions to spatially separate was also observed for linear polyelectrolytes[20]. On this ground we conclude that an ionic microseparation caused by the presence of mixed-size counterions is a universal phenomenon which occurs for polyelectrolytes with various molecular topologies.

Figure 2

Log-lin plot of the pair correlation function for the conventional counterions, (solid lines), and the bulky counterions with a diameter , (dashed lines), calculated around the dendrimer’s center-of-mass at various Bjerrum lengths, , for fixed number fractions of the bulky counterions, , as indicated. In all the panels the maximum error bar is less than .

Figure 14

Log-lin plot of the radial density profiles, , of the monomers around the center-of-mass of the dendrimer plotted for various fractions, , of the bulky counterions with a diameter at . The maximum error bar for all the data sets is less than .

To analyze the uptake of the counterions quantitatively, in Figs. 3a and 4a we display respectively the fractions of the absorbed conventional, , and bulky, , counterions versus at constant values of . The fractions are calculated with respect to the overall number of the counterions, . Furthermore, and denote the mean numbers of the conventional and the bulky counterions at distances less than from the dendrimer’s center-of-mass. Our criterion for absorption is based on that both the radial densities of the monomers and the counterions are nearly zero at , see Figs. 2 and 14. Figure 3a demonstrates that increases monotonically from nearly zero up to with increasing . Thus, for strong Coulomb interactions all the conventional counterions present in the mixture penetrate the dendrimer, whereas for weak electrostatic couplings absorption of these ions is minor. A monotonic increase is also found for , see Fig. 4a. In particular, no absorption of the bulky counterions occurs for weak electrostatic interactions, whereas in the limit of large -values it is almost complete since . Therefore, an uptake of the counterions of both types takes place and is enhanced accordingly to as is increased. In Figs. 3b and 4b we display both fractions as functions of at constant . Our data indicate that increasing results in a monotonic reduction of from the maximum at corresponding to the dendrimer with only conventional counterions. In turn, is found to increase monotonically up to the maximum at corresponding to the polyelectrolyte accompanied by only the bulky counterions. Note that absorption of the bulky counterions with a diameter is slightly weaker than of those with . More interestingly, we observe that the variation of both fractions is almost linear with and steeper at larger -values. Last but not least, we analyze the overall fraction of the absorbed counterions, . Figure 5a demonstrates versus at constant . In accordance with the behavior of its components is a monotonically increasing function of from nearly zero up to one. Note that the -profiles are only weakly affected by the -values and tend to collapse onto a single curve. This observation is also confirmed in Fig. 5b, which demonstrates that is almost a constant function of at fixed . Thus, the absorption process of the counterions is basically determined by Coulomb interactions as long as no distinction between the ions is made. Variations in the composition of the ionic mixture at constant Bjerrum length have only a weak impact on the overall amount of the absorbed counterions.

Figure 3

Fraction, , of the conventional counterions absorbed into the dendrimer versus the Bjerrum length, , plotted for a fixed number fraction, , of the bulky counterions with a diameter ( in the inset) on the lin-log scale (a), and versus for fixed on the lin–lin scale (b). In (b) the solid (dashed) lines represent the data for (). The red lines indicate the number fraction, , of the conventional counterions. Here and in the remaining plots the error bars for all the data points are smaller than the symbol size.

Figure 4

Fraction, , of the bulky counterions absorbed into the dendrimer versus the Bjerrum length, , plotted for a fixed number fraction, , of the bulky counterions with a diameter ( in the inset) on the lin-log scale (a), and versus for fixed on the lin–lin scale (b). In (b) the solid (dashed) lines represent the data for (). The red lines indicate the number fraction, , of the bulky counterions.

Figure 5

Overall fraction, , of the counterions absorbed into the dendrimer versus the Bjerrum length, , plotted for a fixed number fraction, , of the bulky counterions with a diameter ( in the inset) on the lin-log scale (a), and versus at fixed on the lin–lin scale (b). In (b) the solid (dashed) lines represent the data for (). The red lines indicate .

Condensed and delocalized counterions

As demonstated in a number of studies a fraction of the absorbed counterions condense on the charged groups of dendritic polyelectrolytes[8-11,14]. Our data indicate that this is also the case for dendrimers with a mixture of conventional and bulky counterions. This tendency is shown in Fig. 6, in which we display the charged monomer-conventional counterion, , and the charged monomer-bulky counterion, , pair correlation functions, respectively. Pronounced maxima of the profiles at ( =  c, b) indicate that the counterions of both types condense on the dendrimer ions. This phenomenon is enhanced for larger values due to stronger Coulomb attraction between the oppositely charged species, which is manifested by an increase in the height of the maxima. Obviously, the fact that is related to the difference in the counterion size. Note that at given and the maxima of the -profiles are sharper as compared with these of the -profiles, which indicates that the conventional counterions are more likely to penetrate the vicinity of the dendrimer ions. The small size of the conventional ions enables them to approach the charged groups of the dendrimer more efficiently and to be more strongly electrostatically bound with them. On the other hand, the bulky counterions are prevented from gathering in close proximity of the charged monomers due to their stronger excluded volume.

Figure 6

Log-lin plot of the pair correlation function, (solid lines), (dashed lines), for the conventional and the bulky counterions with a diameter respectively, calculated around the charged groups of the dendrimer at various Bjerrum lengths, , for fixed number fractions of the bulky counterions, , as indicated. In all the panels the maximum error bar is less than .

To quantify the degree of counterion condensation, in Figs. 7a and 8a we plot respectively the fractions, and , of the mean number of the condensed conventional and bulky counterions with respect to their overall number in the mixture as functions of at constant . Our criterion for condensation is based on the distance between the counterions and the dendrimer ions. A counterion is considered condensed if its distance from at least one charged monomer is less than , where c,b. In accordance with the previous investigations of charged dendrimers with solely conventional counterions[8-11,14], we find a monotonic increase in and with . Specifically, irrespective of the -values, for the thermal energy dominates over Coulomb attraction between the opposite charges and condensation cannot be achieved, i.e, and are nearly zero. For the electrostatic attraction dominates and promotes counterion condensation. In this regime is already a selective parameter which causes a split of both fractions. In particular, for small -values increases abruptly, whereas the increase in is minor. An opposite tendency is observed for large -values. Note that in the limit of strong Coulomb couplings almost all the conventional counterions present in the mixture are condensed, whereas condensation of the bulky ions still remains incomplete. The latter observation is attributed to weaker electrostatic binding of the bulky counterions in contact with the charged monomers and strong steric repulsions that prevent such contacts. The effect of variation of on condensation of the counterions of both types at constant is presented in Figs. 7b and 8b. Here, both ionic fractions follow the behavior previously discussed for the fractions of the absorbed counterions. Namely, decreases from the maximum at and increases up to the maximum at almost linearly. As compared with the intermediate cases (), the system containing only the conventional counterions () is characterized by the largest fraction of the condensed conventional counterions. Similarly, the fraction of the condensed bulky counterions is the largest in the other extreme of the dendrimer accompanied by the bulky ions only (). Furthermore, the change in the number of both condensed components is more abrupt at larger -values as well. Note that the effect of on condensation of the conventional counterions is rather minor, whereas for the bulky ones clearly visible. In the latter case, at given a drop in is observed for the bulky counterions with a diameter since these ions are more weakly bound with the charged groups. This also results in a smaller slope of the -profiles as functions of . To complement our analysis of condensation in Fig. 9 we display the overall fraction of the condensed counterions, . For the overall amount of the condensed counterions is negligible, whereas in the other regime of strong Coulomb interactions the overall condensation effect is significantly enhanced, see Fig. 9a. In particular, for the -profiles increase abruptly according to . Note that the most abrupt increase of condensation takes place at , i.e., for the dendrimer with only the conventional counterions. The slowest growth occurs in the other extreme of the polyelectrolyte accompanied by the bulky ions only. Figure 9b demonstrates as a function of at constant . It confirms that counterion condensation does not take place in the limit of weak electrostatic interactions. Above the crossover at a nearly linear decrease in is observed, which is steeper for the bulky counterions with a diameter . Thus, replacing the conventional counterions with the bulky ones weakens the overall condensation effect. The maximum in the overall fraction of the condensed counterions is obtained for the dendrimer with solely the conventional counterions and decays with increasing the amount of the bulky ions.

Figure 7

Fraction, , of the condensed conventional counterions versus the Bjerrum length, , plotted for a fixed number fraction, , of the bulky counterions with a diameter ( in the inset) on the lin-log scale (a), and versus at fixed values of on the lin–lin scale (b). In (b) the solid (dashed) lines represent the data for (). The red lines indicate .

Figure 8

Fraction, , of the condensed bulky counterions versus the Bjerrum length, , plotted for a fixed number fraction, , of the bulky counterions with a diameter ( in the inset) on the lin-log scale (a), and versus at fixed values of on the lin–lin scale (b). In (b) the solid (dashed) lines represent the data for (). The red lines indicate .

Figure 9

Overall fraction, , of the condensed counterions versus the Bjerrum length, , plotted for a fixed number fraction, , of the bulky counterions with a diameter ( in the inset) on the lin-log scale (a), and versus at fixed values of on the lin–lin scale (b). In (b) the solid (dashed) lines represent the data for (). The red lines indicate .

Besides condensed counterions, inside the dendrimer there are also delocalized ions which freely penetrate the molecule and exert the osmotic pressure[9]. Figure 10a displays the fraction of the delocalized conventional counterions, , whereas in Fig. 11a we show the fraction of the delocalized bulky counterions, . Both plots are presented as functions of at constant . Our data indicate that both ionic fractions are non-monotonic functions of and exhibit a broad maximum. Before the maximum it is absorption of the counterions of both ionic types that increases more rapidly with than their condensation. For strong Coulomb couplings an opposite tendency takes place, i.e., counterion condensation is dominant and results in a decrease in and . Note that the -maximum is the highest in the presence of solely the conventional ions and flattens with increasing . On the other hand, the -fraction is the largest in the presence of solely the bulky ions and flattens as is decreased. Moreover, at there are no delocalized conventional counterions due to their complete condensation, whereas the amount of the delocalized bulky ions is still finite. In Figs. 10b and 11b we display both ionic fractions versus at constant . A monotonic decrease in and a monotonic increase in is observed. Deviations from linear behavior of the corresponding ionic fractions are weak though more noticeable for systems with the bulky ions with a diameter . For both ionic fractions we do not find a clear, monotonic dependence of the slopes on . Finally, in Fig. 12a we display the overall fraction of the delocalized counterions, , versus . Similar to its components, is a non-monotonic function of exhibiting a broad maximum around which corresponds to the crossover between the two regimes of counterion absorption dominating over counterion condensation and vice versa. Note that in line with the -profiles, the effect of varying becomes noticeable above the crossover. The latter observation is also confirmed in Fig. 12b which demonstrates as a function of . Namely, for only a weak dependence of on is observed, whereas above the crossover an almost linear, sharp increase in the -fraction takes place that is more pronounced for . The latter conclusion is another manifestation of the fact that the larger bulky counterions are more weakly bound with the charged groups as compared to the counterions with a smaller diameter.

Figure 10

Fraction, , of the delocalized conventional counterions versus the Bjerrum length, , plotted for a fixed number fraction, , of the bulky counterions with a diameter ( in the inset) on the lin-log scale (a), and versus at fixed values of on the lin–lin scale (b). In (b) the solid (dashed) lines represent the data for ().

Figure 11

Fraction, , of the delocalized bulky counterions versus the Bjerrum length, , plotted for a fixed number fraction, , of the bulky counterions with a diameter ( in the inset) on the lin-log scale (a), and versus at fixed values of on the lin–lin scale (b). In (b) the solid (dashed) lines represent the data for ().

Figure 12

Overall fraction, , of the delocalized counterions versus the Bjerrum length, , plotted for a fixed number fraction, , of the bulky counterions with a diameter ( in the inset) on the lin-log scale (a), and versus at fixed values of on the lin–lin scale (b). In (b) the solid (dashed) lines represent the data for ().

Radius of gyration

In Fig. 13 we display the swelling factor of the dendritic polyelectrolyte, , defined as the ratio between the radius of gyration of the charged dendrimer and its neutral analogue versus and . The data indicate that as a function of the dendrimer’s size changes non-monotonically in the whole range of the -values, see Fig. 13a. Increasing causes swelling of the polymer up to a broad maximum followed by deswelling observed for strong Coulomb couplings. The dependence of is qualitatively similar to the outcomes of the previous simulations of dendritic polyelectrolytes carried out for solely conventional ions[8,10,11,14]. The equilibrium size of the dendrimer results from a balance between the elasticity of the molecule and intramolecular interactions of different origin. Namely, for absorption of the counterions is minor and is unaffected by both and . In this regime swelling is mostly due to Coulomb repulsion between the charged groups of the dendrimer. Ionic specificity becomes essential for where penetration of the dendrimer’s pervaded volume by the counterions is significant and leads to a split in depending on the value of . The -profiles exhibit maxima at which correspond with the maximum of the overall fraction of the delocalized counterions, . Here, the high osmotic pressure due to the delocalized counterions of both types and the unscreened electrostatic repulsion between the dendrimer ions lead to the pronounced swelling effect[8,9]. Furthermore, the maximum of increases with varying from zero to one. Given the fact that in the considered -range the drop in the fractions and is accompanied by the increase in the fractions and such that their sums, and , only weakly depend on , the -induced enhancement of swelling is attributed to a larger contribution of the excluded volume interactions to the osmotic pressure caused by the absorbed bulky counterions. Subsequently, at counterion condensation gradually becomes dominant. It screens the intramolecular Coulomb repulsion and suppresses the impact of the delocalized ions, which effectively results in polymer deswelling. Note that for the dendrimer accompanied by ionic mixtures with there exist -values at which the molecule recovers its size in the neutral state, and a further increase in promotes its collapse. In the latter case the radius of gyration is mostly determined by excluded volume repulsions and multiple attractions between ion pairs formed by the condensed counterions and the charged groups[10,36,37]. For counterion mixtures with a reduction of the dendrimer size occurs as well, though the molecule still remains swollen. The effect of on is also displayed in Fig. 13b. This figure clearly shows that increases monotonically with at constant . In particular, at increasing from zero to one leads to a structural transition of the dendrimer from collapsed to swollen conformations. This kind of transiton is due to weaker counterion condensation which reduces the effect of the multiple attraction promoting shrinking and enhances the impact of the aforementioned sources of intramolecular repulsions. The conformational changes of the dendrimer are also reflected by the spatial distribution of the monomers shown in Fig. 14. Here, we display the density profiles of the monomers, , as functions of the rescaled radial distance from the dendrimer’s center-of-mass at . The collapse of the polymer with decreasing is signaled by an increase in , especially within the dendrimer’s pervaded volume. Moreover, our simulations confirm that the dendrimer exists in dense-core conformations[8]. The monomer densities are characterized by a sharp peak in a neighborhood of the center-of-mass, a broad plateau in the polymer’s domain and a drop on the periphery. Note that the shape of the -profiles to some extent follows the density profiles of the monomers. This is another manifestation of pronounced condensation of the conventional ions which predominantly occupy a close proximity of the charged groups at large Bjerrum lengths. From this we also conclude that the emergence of the ionic core–shell structures inside the molecule is mostly driven by selective condensation of the counterions. We visualize counterion absorption into the dendrimer in Fig. 15 by displaying a series of snapshots from our MD simulations at for various values of .

Figure 13

Swelling factor, , versus the Bjerrum length, , plotted for a fixed number fraction, , of the bulky counterions with a diameter ( in the inset) on the lin-log scale (a), and versus at fixed values of on the lin–lin scale (b). In (b) the solid (dashed) lines represent the data for (). The red lines indicate .

Figure 15

Snapshots of conformations of the dendrimer with mixed-size counterions at for fixed fractions of the bulky counterions with a diameter : (a), (b), (c), (d), (e), (f). The monomers (and the conventional/bulky counterions) are displayed by red (yellow/blue) spheres.

Summary

In this work using Langevin dynamics simulations we studied G4S4 dendritic polyelectrolyte with two component mixtures of conventional and bulky counterions. The main focus of our investigation was on the conformational behavior of the molecule modified by the reduced Bjerrum length, , and the number fraction of the bulky counterions, . Our simulations indicate that irrespective of the amount of the bulky counterions in the ionic mixture the dendrimer swells non-monotonically with increasing . At small Bjerrum lengths, , the counterions are distributed uniformly in solution and their density within the dendrimer’s pervaded volume is negligible. As a consequence, the radius of gyration of the polyelectrolyte exceeds that of the neutral dendrimer due to Coulomb repulsion between the like-charged monomers and remains unaffected by the ions. At intermediate Bjerrum lengths, , the electrostatic attraction between the counterions and the dendrimer ions promotes absorption of the former into the dendrimer’s pervaded volume and their condensation. Inside the dendrimer the densities of the counterions of both types become significant and starts to affect the molecule’s conformations. In this regime, for a given value of , the -dependent radius of gyration of the dendrimer exhibits a broad maximum, which corresponds with the maximum of the overall fraction of the absorbed, delocalized counterions. Here, swelling is primarily attributed to the high osmotic pressure exerted by the delocalized counterions on the polymer and unscreened intramolecular Coulomb repulsion between the dendrimer ions. Moreover, swelling increases monotonically with due to the enhanced contribution of the excluded volume interactions between the particles penetrating the dendrimer’s interior to the osmotic pressure. Subsequently, for large Bjerrum lengths, , condensation of the conventional and the bulky counterions is significant and determines the conformational behavior of the polymer. In this regime deswelling and eventually collapse of the dendrimer into a globule-like state takes place, the latter being attributed to formation of condensed-ion/charged-monomer pairs possessing a net attraction. Increasing weakens counterion condensation as well as the impact of the intramolecular attractive forces, and a departure of the polyelectrolyte from collapsed to swollen conformations is observed. Our simulations show that the absorbed counterions tend to separate into core–shell microstructures. The core of the dendrimer is rich in the conventional ions, whereas the volume around it contains the bulky ions. A similar observation was made for strongly charged linear chains[20]. Our study indicates that conformations of highly branched macromolecules tend to reveal an ionic separation facilitated by selective condensation of conventional and bulky counterions. It seems likely that this phenomenon is generic in nature and does not depend on specific arrangements of branching nodes. Thus a broad range of polymeric architectures including dendrimers of higher generations and with different spacer lengths, dendronized, pom-pom-shaped, star-burst and bottlebrush polymers[38-40] as well as microgels[41] are expected to demonstrate similar microseparation effects. These systems are clearly a motivation for further studies. Another interesting subject for the future work is the influence of pH conditions and charge arrangements on the counterion segregation.