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, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _B^*$$\end{document}λB∗, which controls the strength of electrostatic interactions and the number fraction of the bulky counterions, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_b$$\end{document}fb, 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 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_b$$\end{document}fb. 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 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_b$$\end{document}fb. For strong electrostatic couplings the dendrimer deswells correspondingly to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_b$$\end{document}fb. 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.

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][5][6][7][8][9][10][11][12][13][14][15][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][6][7][8][9][10][11][12][13][14][15][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][22][23][24] , polyelectrolyte solutions 20,25 , the electric double layer around DNA 26 and rodlike polyelectrolytes [27][28][29][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, * B . 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, f b . We find that both parameters * B and f b 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), N bg , in the dendrimer are given by and 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, N bg monovalent negative counterions are included in the model. To study the dendrimer with mixed-size counterions, N b and N c = N bg − N b 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 potential where ǫ is the interaction strength, and R c αβ = 2 1/6 σ αβ 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, σ αβ = σ α + σ β /2 , is used. The bonds joining any two adjacent beads which are described by the Finite Extensible Non-linear Elastic potential (FENE)  www.nature.com/scientificreports/ where k is the spring constant and R 0 is the maximum extension of the bonds. The electrostatic interactions are introduced by pairwise Coulomb potential where k B denotes the Boltzmann constant, T the absolute temperature, B the Bjerrum length of the solvent, and r the distance between the centers of the ith and the jth particles with the charge valence z i , z j = ±1 . With ε 0 standing for the electric permittivity of the vacuum, ε r for the relative permittivity of the solvent and e for elementary charge, the Bjerrum length is defined as B (ε r , T) = e 2 /(4πε 0 ε r k B T) . In particular, for water at room temperature ( T r ≈ 298 K), ε rw ≈ 81 , so that B (ε rw , T r ) ≈ 7 Å.
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, ǫ = 1k B T r , and the monomer size, σ m = B (ε rw , T r ) , were taken as the real units of energy and length respectively. The unitless simulation input parameters (denoted with an asterisk) were chosen as ǫ * = 1 for the strength of the LJ potential, σ * m , σ * c = 1 for the size of the monomers and the conventional counterions, and m * = 1 for the mass of the monomers and the counterions. For the FENE potential the parameters were k * = kσ 2 m /ǫ = 30 and R * 0 = R 0 /σ m = 1.5 33 . According to the formula for reduced elementary charge, q * = e/(4πε 0 ǫσ m ) 1/2 , and for the reduced temperature, T * = k B T/ǫ , we set these parameters to the values q * = 9 and T * = 1 . 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 10 −4 and the real space cutoff radius R c * = 10 32 . The dumping parameter in the Langevin equation of motion was set to γ * = γ /τ = 1 , where τ = σ m (m/ǫ) 1/2 is the LJ time unit and m is the assumed real unit of mass. For example, for m ≈ 30 g/mol, τ ≈ 2.4 ps. The calculations were performed with a time step �t * = �t/τ = 0.005 . Note that the reduced Bjerrum length, * B (ε r , T r ) = B (ε r , T r )/σ m , for a solvent at room temperature can be expressed as * B (ε r , T r ) = * B (ε rw , T r )ε rw /ε r , which with our unit of length yields * B (ε r , T r ) = * B = ε rw /ε r . This enabled us to vary * B by changing the relative electric permittivity of the solvent, ε r . According to Eq. (5), * B is the parameter that determines Coulomb interaction energy between two monovalent charges at a distance σ m 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 L * = 100 . 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 10 6 equilibration integration steps were performed and followed by 10 7 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 τ D , where τ D is defined as time needed for dendrimer to diffuse its own size. All G4 dendrimers considered in the current study had τ D < 1000 τ . 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 σ * b = σ b /σ m = 2 , 3, and their number fraction was set to f b = N b /N bg = 0, 0.16, 0.32, 0.48, 0.64, 0.8, 1 . The reduced Bjerrum length, * B , 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 σ * c = 1 as the conventional (counter)ions and with a diameter σ * b = 2 , 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 g cm,c and g cm,b 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, r * /R * g , from the dendrimer's center-of-mass rescaled by the dendrimer's radius of gyration, R * g , for different Bjerrum lengths, * B , and a given number fraction of the bulky counterions, f b . The data demonstrate that both types of ionic species accumulate at distances r * /R * g 2 from the dendrimer's core. The absorption effect is enhanced with increasing * B , 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][6][7][8][9][10][11][12][13][14][15][16] . As shown in Fig. 2, the pair correlation functions reveal a strong dependence on the ionic type. For * B 4 the g cm,c -profiles peak in a neighborhood of the center of mass, flatten in the molecule's domain, and drop on its periphery. The g cm,b -profiles exhibit a broad maximum in the domain and decrease near the center and on the periphery. Thus for sufficiently large * B -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 www.nature.com/scientificreports/    www.nature.com/scientificreports/ 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.
To analyze the uptake of the counterions quantitatively, in Figs. 3a and 4a we display respectively the fractions of the absorbed conventional, f in = N in /N bg , and bulky, f bin = N bin /N bg , counterions versus * B at constant values of f b . The fractions are calculated with respect to the overall number of the counterions, N bg . Furthermore, N in and N bin denote the mean numbers of the conventional and the bulky counterions at distances less than 2R * g 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 r * /R * g ≈ 2 , see Figs. 2 and 14. Figure 3a demonstrates that f in increases monotonically from nearly zero up to 1 − f b with increasing * B . 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 f bin , see Fig. 4a. In particular, no absorption of the bulky counterions occurs for weak electrostatic interactions, whereas in the limit of large * B -values it is almost complete since f bin ≈ f b . Therefore, an uptake of the counterions of both types takes place and is enhanced accordingly to f b as * B is increased. In Figs. 3b and 4b we display both fractions as functions of f b at constant * B . Our data indicate that increasing f b results in a monotonic reduction of f in from the maximum at f b = 0 corresponding to the dendrimer with only conventional counterions. In turn, f bin is found to increase monotonically up to the maximum at f b = 1 corresponding to the polyelectrolyte accompanied by only the bulky counterions. Note that absorption of the bulky counterions with a diameter σ * b = 3 is slightly weaker than of those with σ * b = 2 . More interestingly, we observe that the variation of both fractions is almost linear with f b and steeper at larger * B -values. Last but not least, we analyze the overall fraction of the absorbed counterions, f oin = f in + f bin . Figure 5a demonstrates f oin versus * B at constant f b . In accordance with the behavior of its components f oin is a monotonically increasing function of * B from nearly zero up to one. Note that the f oin -profiles are only weakly affected by the f b -values and tend to collapse onto a single curve. This observation is also confirmed in Fig. 5b, which demonstrates that f oin is almost a constant function of f b at fixed * B . 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.
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][9][10][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, g chm,c , and the charged monomer-bulky counterion, g chm,b , pair correlation functions, respectively. Pronounced maxima of the profiles at r * α ≈ 1 + σ * α /2 ( α = c, b) indicate that the counterions of both types condense on the dendrimer ions. This phenomenon is enhanced for larger * B 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 r * b > r * c is related to the difference in the counterion size. Note that at given f b and * B the maxima of the g chm,c -profiles are sharper as compared with these of the g chm,b -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.
To quantify the degree of counterion condensation, in Figs. 7a and 8a we plot respectively the fractions, f cond = N cond /N bg and f bcond = N bcond /N bg , of the mean number of the condensed conventional and bulky counterions with respect to their overall number in the mixture as functions of * B at constant f b . 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 2 + σ * α /2 , where α = c,b.   1 ). Furthermore, the change in the number of both condensed components is more abrupt at larger * B -values as well. Note that the effect of σ * b on condensation of the conventional counterions is rather minor, whereas for the bulky ones clearly visible. In the latter case, at given f b a drop in f bcond is observed for the bulky counterions with a diameter σ * b = 3 since these ions are more weakly bound with the charged groups. This also results in a smaller slope of the f bcond -profiles as functions of f b . To complement our analysis of condensation in Fig. 9 we display the overall fraction of the condensed counterions, f ocond = f cond + f bcond . For * B 1 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 * B 4 the f ocond -profiles increase abruptly according to f b . Note that the most abrupt increase of condensation takes place at f b = 0 , 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 f ocond as a function of f b at constant * B . It confirms that counterion condensation does not take place in the limit of weak electrostatic interactions. Above the crossover at * B ≈ 1 a nearly linear decrease in f ocond is observed, which is steeper for the bulky counterions with a diameter σ * b = 3 . 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.
Besides condensed counterions, inside the dendrimer there are also delocalized ions which freely penetrate the molecule and exert the osmotic pressure 9 . Figure 10a   www.nature.com/scientificreports/ delocalized counterions, f od = f d + f bd , versus * B . Similar to its components, f od is a non-monotonic function of * B exhibiting a broad maximum around * B ≈ 4 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 f ocond -profiles, the effect of varying f b becomes noticeable above the crossover. The latter observation is also confirmed in Fig. 12b which demonstrates f od as a function of f b . Namely, for * B 4 only a weak dependence of f ocond on f b is observed, whereas above the crossover an almost linear, sharp increase in the f od -fraction takes place that is more pronounced for σ * b = 3 . 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. Fig. 13 we display the swelling factor of the dendritic polyelectrolyte, α = R * g /R * g0 , defined as the ratio between the radius of gyration of the charged dendrimer and its neutral analogue versus * B and f b . The data indicate that as a function of * B the dendrimer's size changes non-monotonically in the whole range of the f b -values, see Fig. 13a. Increasing * B 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 * B 1 absorption of the counterions is minor and α is unaffected by both f b and σ * b . In this regime swelling is mostly due to Coulomb repulsion between the charged groups of the dendrimer. Ionic specificity becomes essential for * B 1 where penetration of the dendrimer's pervaded volume by the counterions is significant and leads to a split in α depending on the value of f b . The α-profiles exhibit maxima at 1 * B 4 which correspond with the maximum of the overall fraction of the delocalized counterions, f od . 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 f b from zero to one. Given the fact that in the considered * B -range the drop in the fractions f in and f d is accompanied by the increase in the fractions f bin and f bd such that their sums, f oin and f od , only weakly depend on f b , the f b -induced enhancement of swelling is attributed to a larger contribution of the excluded volume interactions to the osmotic pressure caused by the absorbed bulky  -values at which the molecule recovers its size in the neutral state, and a further increase in * B 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 f b 0.32 a reduction of the dendrimer size occurs as well, though the molecule still remains swollen. The effect of f b on α is also displayed in Fig. 13b. This figure clearly shows that α increases monotonically with f b at constant * B 1 . In particular, at * B ≈ 16 increasing f b 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, ρ * cm,m , as functions of the rescaled radial distance from the dendrimer's center-of-mass at * B = 16 . The collapse of the polymer with decreasing f b is signaled by an increase in ρ * cm,m , 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-ofmass, a broad plateau in the polymer's domain and a drop on the periphery. Note that the shape of the g cm,c -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 selectivecondensation of the counterions. We visualize counterion absorption into the dendrimer in Fig. 15 by displaying a series of snapshots from our MD simulations at * B = 16 for various values of f b .

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, * B , and the number fraction of the bulky www.nature.com/scientificreports/ counterions, f b . Our simulations indicate that irrespective of the amount of the bulky counterions in the ionic mixture the dendrimer swells non-monotonically with increasing * B . At small Bjerrum lengths, * B 1 , 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, 1 * B 4 , 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 f b starts to affect the molecule's conformations. In this regime, for a given value of f b , the * B -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 f b 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, * B 4 , 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 f b 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, pompom-shaped, star-burst and bottlebrush polymers [38][39][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.