A Monte Carlo study of solutions of oppositely charged polyelectrolytes

Self-Cross-Linkable Chitosan-Alginate Complexes Inspired by Mussel Glue Chemistry

In this study, mussel-inspired chemistry, based on catechol-amine reactions, was adopted to develop self-cross-linkable chitosan-alginate (Chi-Alg) complexes. To do so, the biopolymers were each substituted with ∼20% catechol groups (ChiC and AlgC), and then four complex combinations (Chi-Alg, ChiC-Alg, Chi-AlgC, ChiC-AlgC) were prepared at the surface and in bulk solution. Based on QCM-D and lap shear adhesion tests, the complex with catechol only on Chi (ChiC-Alg) did not show a significant variation from the control complex (Chi-Alg). Conversely, the complexes with catechol on alginate (Chi-AlgC and ChiC-AlgC) rendered a self-cross-linking property and enhanced cohesive properties.

Polyelectrolyte-multivalent molecule complexes: physicochemical properties and applications

The complexation of polyelectrolytes with other oppositely charged structures gives rise to a great variety of functional materials with potential applications in a wide spectrum of technological fields. Depending on the assembly conditions, polyelectrolyte complexes can acquire different macroscopic configurations such as dense precipitates, nanosized colloids and liquid coacervates. In the past 50 years, much progress has been achieved to understand the principles behind the phase separation induced by the interaction of two oppositely charged polyelectrolytes in aqueous solutions, especially for symmetric systems (systems in which both polyions have similar molecular weight and concentration). However, in recent years, the complexation of polyelectrolytes with alternative building blocks such as small charged molecules (multivalent inorganic species, oligopeptides, and oligoamines, among others) has gained attention in different areas. In this review, we discuss the physicochemical characteristics of the complexes formed by polyelectrolytes and multivalent small molecules, putting a special emphasis on their similarities with the well-known polycation-polyanion complexes. In addition, we analyze the potential of these complexes to act as versatile functional platforms in various technological fields, such as biomedicine and advanced materials engineering.

Coarse-grained Simulations of the Impact of Chain Length and Stiffness on the Formation and Aggregation of Polyelectrolyte Complexes

Polyelectrolyte complexes formed from nucleic acids and synthetic polycations have been studied because of their potential in gene delivery. Coarse-grained molecular dynamics simulations are performed to examine the impact of chain length and polyanion stiffness on polyplex formation and aggregation. Polyplexes containing single polyanion chain fall into three structural regimes depending on polyanion stiffness: flexible polyanions form collapsed complexes, semiflexible polyanions form various morphologies including toroids and hairpins, and stiff polyanions form rod-like structures. Polyplex size generally decreases as polycation length increases. Aggregation (i.e., formation of complexes containing multiple polyanions) is observed in some simulations containing multiple polyanions and an excess of short polycations. Aggregation is observed to only occur for semiflexible and stiff polyanions and is promoted by shorter polycation lengths. Simulations of short, stiff polyanions condensed by long polycations are used as a model for siRNA gene delivery complexes. These simulations show multiple polyanions are spaced out along the polycation with polyanion-polyanion interactions, usually limited to overlapping chain ends. These structures differ from aggregates of longer polyanions in which the polyanions are packed together in parallel, forming bundles.

Coarse-grained simulation studies on the adsorption of polyelectrolyte complexes upon lipid membranes

Conformations of a polyelectrolyte complex irreversibly bound to a zwitterionic lipid bilayer.

Potential of mean force and transient states in polyelectrolyte pair complexation

The pair association between two polyelectrolytes (PEs) of the same size but opposite charge is systematically studied in terms of the potential of mean force (PMF) along their center-of-mass reaction coordinate via coarse-grained, implicit-solvent, explicit-salt computer simulations. The focus is set on the onset and the intermediate transient stages of complexation. At conditions above the counterion-condensation threshold, the PE association process exhibits a distinct sliding-rod-like behavior where the polymer chains approach each other by first stretching out at a critical distance close to their contour length, then "shaking hand" and sliding along each other in a parallel fashion, before eventually folding into a neutral complex. The essential part of the PMF for highly charged PEs can be very well described by a simple theory based on sliding charged "Debye-Hückel" rods with renormalized charges in addition to an explicit entropy contribution owing to the release of condensed counterions. Interestingly, at the onset of complex formation, the mean force between the PE chains is found to be discontinuous, reflecting a bimodal structural behavior that arises from the coexistence of interconnected-rod and isolated-coil states. These two microstates of the PE complex are balanced by subtle counterion release effects and separated by a free-energy barrier due to unfavorable stretching entropy.

Application of Monte Carlo simulation in addressing key issues of complex coacervation formed by polyelectrolytes and oppositely charged colloids

This paper reviews the recent advance of Monte Carlo (MC) simulation in addressing key issues of complex coacervation between polyelectrolytes and oppositely charged colloids. Readers were first supplied with a brief overview of current knowledge and experimental strategies in the study of complex coacervation. In the next section, the general MC simulation procedures as well as representative strategies applied in complex coacervation were summarized. The unique contributions of MC simulation in either capturing delicate features, easing the experimental trials or proving the concept were then elucidated through the following aspects: i) identify phase boundary and decouple interaction contributions; ii) clarify composition distribution and internal structure; iii) predict the influences of physicochemical conditions on complex coacervation; iv) delineate the mechanisms for "binding on the wrong side of the isoelectric point". Finally, current challenges as well as prospects of MC simulation in complex coacervation are also discussed. The ultimate goal of this review is to provide readers with basic guideline for synergistic design of experiments in combination with MC simulation, and deliver convincing interpretation and reliable prediction for the structure and behavior in polyelectrolyte-macroion complex coacervation.

Modeling competitive substitution in a polyelectrolyte complex

We have simulated the invasion of a polyelectrolyte complex made of a polycation chain and a polyanion chain, by another longer polyanion chain, using the coarse-grained united atom model for the chains and the Langevin dynamics methodology. Our simulations reveal many intricate details of the substitution reaction in terms of conformational changes of the chains and competition between the invading chain and the chain being displaced for the common complementary chain. We show that the invading chain is required to be sufficiently longer than the chain being displaced for effecting the substitution. Yet, having the invading chain to be longer than a certain threshold value does not reduce the substitution time much further. While most of the simulations were carried out in salt-free conditions, we show that presence of salt facilitates the substitution reaction and reduces the substitution time. Analysis of our data shows that the dominant driving force for the substitution process involving polyelectrolytes lies in the release of counterions during the substitution.

Branched-linear polyion complexes at variable charge densities

Structural behavior of complexes formed by a charged and branched copolymer and an oppositely charged and linear polyion was examined by Monte Carlo simulations employing a coarse-grained bead-spring model. The fractional bead charge and the branching density were systematically varied; the former between 0e and 1e and the latter such that both the comb-polymer and the bottle-brush limits were included. The number of beads of the main chain of the branched copolymer and of the linear polyion was always kept constant and equal, and a single side-chain length was used. Our analysis involved characterization of the complex as well as investigation of size, shape, and flexibility of the charged moieties. An interplay between Coulomb interaction and side-chain repulsion governed the structure of the polyion complex. At strong Coulomb interaction, the complexes underwent a gradual transition from a globular structure at low branching density to an extended one at high branching density. As the electrostatic coupling was decreased, the transition was smoothened and shifted to lower branching density, and, eventually, a behavior similar to that found for neutral branched polymer was observed. Structural analogies and dissimilarities with uncharged branched polymers in poor solutions are discussed.

Phase behavior of electrostatically complexed polyelectrolyte gels using an embedded fluctuation model

Nanostructured, responsive hydrogels formed due to electrostatic interactions have promise for applications such as drug delivery and tissue mimics. These physically cross-linked hydrogels are composed of an aqueous solution of oppositely charged triblocks with charged end-blocks and neutral, hydrophilic mid-blocks. Due to their electrostatic interactions, the end-blocks microphase separate and form physical cross-links that are bridged by the mid-blocks. The structure of this system was determined using a new, efficient embedded fluctuation (EF) model in conjunction with self-consistent field theory. The calculations using the EF model were validated against unapproximated field-theoretic simulations with complex Langevin sampling and were found consistent with small angle X-ray scattering (SAXS) measurements on an experimental system. Using both the EF model and SAXS, phase diagrams were generated as a function of end-block fraction and polymer concentration. Several structures were observed including a body-centered cubic sphere phase, a hexagonally packed cylinder phase, and a lamellar phase. Finally, the EF model was used to explore how parameters that directly relate to polymer chemistry can be tuned to modify the resulting phase diagram, which is of practical interest for the development of new hydrogels.

Branched-linear polyion complexes investigated by Monte Carlo simulations

Complexes formed by one charged and branched copolymer with an oppositely charged and linear polyion have been investigated by Monte Carlo simulations. A coarse-grained description has been used, in which the main chain of the branched polyion and the linear polyion possess the same absolute charge and charge density. The spatial extension and other structural properties, such as bond-angle orientational correlation function, asphericity, and scaling analysis of formed complexes, at varying branching density and side-chain length of the branched polyion, have been explored. In particular, the balance between cohesive Coulomb attraction and side-chain repulsions resulted in two main structures of a polyion complex. These structures are (i) a globular polyion core surrounded by side chains appearing at low branching density and (ii) an extended polyion core with side chains still being expelled at high branching density. The globule-to-extended transition occurred at a crossover branching density being practically independent of the side chain length.

Chitosan-DNA complexes: charge inversion and DNA condensation

The design of biocompatible polyelectrolyte complexes is a promising strategy for in vivo delivery of biologically active macromolecules. Particularly, the condensation of DNA by polycations received considerable attention for its potential in gene delivery applications, where the development of safe and effective non-viral vectors remains a central challenge. Among polymeric polycations, Chitosan has recently emerged as a very interesting material for these applications. In this study, we compare the observed aggregation behavior of Chitosan-DNA complexes with the predictions of existing models for the complexation of oppositely charged polyelectrolytes. By using different and complementary microscopy approaches (AFM, FESEM and TEM), light scattering and electrophoretic mobility techniques, we characterized the structures of the complexes formed at different charge ratios and Chitosan molecular weight. In good agreement with theoretical predictions, a reentrant condensation, accompanied by charge inversion, is clearly observed as the polycation/DNA charge ratio is increased. In fact, the aggregates reach their maximum size in correspondence of a value of the charge ratio where their measured net charge inverts its sign. This value does not correspond to the stoichiometric 1:1 charge ratio, but is inversely correlated with the polycation length. Distinctive "tadpole-like" aggregates are observed in excess polycation, while only globular aggregates are found in excess DNA. Close to the isoelectric point, elongated fiber-like structures appear. Within the framework of the models discussed, different apparently uncorrelated observations reported in the literature find a systematic interpretation. These results suggest that these models are useful tools to guide the design of new and more efficient polycation-based vectors for a more effective delivery of genetic material.

Adsorption and solution properties of bottle-brush polyelectrolyte complexes: effect of molecular weight and stoichiometry

Polyelectrolyte complexes (PECs) self-assembled from bottle-brush polyelectrolytes, having a cationic main chain and uncharged side chains, and linear anionic sodium polystyrenesulfonate (NaPSS) have been investigated with emphasis on (i) the charge density and side chain density of the bottle-brush polyelectrolyte, (ii) the molecular weight of NaPSS, and (iii) the charge stoichiometry of the mixture. Light scattering and electrophoretic mobility data demonstrate that small molecular complexes are formed when the PEO45 side chain density is sufficiently high to provide steric stabilization and prevent PEC aggregation. The adsorption of PECs on negatively charged silicon oxynitride was investigated using dual polarization interferometry, and the time evolution of the adsorbed amount and thickness was determined. Cationic, uncharged, and negatively charged complexes all adsorb to negatively charged silicon oxynitride, and maximum adsorption is achieved for positively charged complexes containing small amounts of PSS. The adsorbed amount and the kinetics of adsorption are reduced with increasing PSS content, and for any given stoichiometry with increasing PSS molecular weight. These findings are discussed in terms of the PEC structure and the ability of anionic polyelectrolytes to leave the PECs during adsorption.

Formulation of chitosan-TPP-pDNA nanocapsules for gene therapy applications

The encapsulation of DNA inside nanoparticles meant for gene delivery applications is a challenging process where several parameters need to be modulated in order to design nanocapsules with specific tailored characteristics. The purpose of this study was to investigate and improve the formulation parameters of plasmid DNA (pDNA) loaded in chitosan nanocapsules using tripolyphosphate (TPP) as polyanionic crosslinker. Nanocapsule morphology and encapsulation efficiency were analyzed as a function of chitosan degree of deacetylation and chitosan-TPP ratio. The manipulation of these parameters influenced not only the particle size but also the encapsulation and release of pDNA. Consequently the transfection efficiency of the nanoparticulated systems was also enhanced with the optimization of the particle characteristics. Overall, the differently formulated nanoparticulated systems possess singular properties that can be employed according to the desired gene delivery application.

Coarse-grained molecular dynamics simulations of DNA condensation by block copolymer and formation of core-corona structures

Coarse-grained molecular dynamics simulations are used to study the condensation of single polyanion chains with block copolymers composed of cationic and neutral blocks. The simulations are an effort to model complexes formed with DNA and cationic copolymers such as polyethylenimine-g-polyethylene glycol which have been used in gene delivery. The simulations reveal that increases in the cationic block length of the copolymer result in greater condensation of the polyanion. The ability of the complexes to form core-corona structures, with the neutral blocks of the copolymers forming a corona around a dense core formed from the charged beads, is investigated. The core-corona structure is shown to be dependent on both condensation of the polyanion chain and the length of the neutral block of the copolymer. Increasing the length of the cationic and neutral blocks of the copolymer both result in improvement in the core-corona structure. The internal structure of the complex core is shown to be a function of the architecture of the copolymer. Complexes formed from linear diblock copolymers have homogeneous cores with similarly arranged cationic and anionic beads; however, complexes formed with star-shaped copolymers have a layered core structure, with anionic beads found in the center of the cores.

Polyelectrolyte condensation in bulk, at surfaces, and under confinement

In this review we discuss recent results from computer simulations based on coarse-grained polyion models representing aqueous solutions of polyelectrolytes. The focus will be directed to the conformation of the polyions and, in particular, their condensation in bulk, induced by multivalent ions and oppositely charged polyelectrolytes, at responsive surfaces and under confinement.

Molecular dynamics simulations of DNA-polycation complex formation

Complexes formed from DNA and polycations are of interest because of their potential use in gene therapy; however, there remains a lack of understanding of the structure and formation of DNA-polycation complexes at atomic scale. In this work, molecular dynamics simulations of the DNA duplex d(CGCGAATTCGCG) in the presence of polycation chains are carried out to shed light on the specific atomic interaction that result in complex formation. The structures of complexes formed from DNA with polyethylenimine, which is considered one of the most promising DNA vector candidates, and a second polycation, poly-L-lysine, are compared. After an initial separation of approximately 50 A, the DNA and polycation come together and form a stable complex within 10 ns. The DNA does not undergo any major structural changes on complexation and remains in the B-form. In the formed complex, the charged amine groups of the polycation mainly interact with DNA phosphate groups, with polycation intrusion into the major and minor grooves dependent on the identity and charge state of the polycation. The ability of the polycation to effectively neutralize the charge of the DNA phosphate groups and the resulting influence on the DNA helix interaction are discussed.

Conformational characteristics of single flexible polyelectrolyte chain

The behaviour of a flexible anionic chain of 150 univalent and negatively charged beads connected by a harmonic-like potential with each other in the presence of an equal number of positive and free counterions, is studied in molecular dynamics simulations with Langevin thermostat in a wide range of temperatures. Simulations were carried out for several values of the bending parameter, corresponding to fully flexible polyion, moderately and strongly stiff polyion as well as for the case when bend conformation is preferable to the straight one. We have found that in all cases three regimes can be distinguished, which can be characterized as "random coil", observed at high temperatures; "extended conformation" observed at moderate temperatures (of the order of 1 in reduced units), and compact "globular conformation" attained at low temperatures. While the transition between high-temperature random and extended conformations is gradual, the transition from the extended coil to the globular state, taking place at a temperature of about 0.2 in reduced units, is of abrupt character resembling a phase transition.

Interaction between colloids with grafted diblock polyampholytes

The interaction between composite colloidal particles composed of a spherical core and grafted AB-diblock polyampholytes (diblock copolymers with oppositely charged blocks) are investigated by using a coarse-grained model solved with Monte Carlo simulations. The B block is end-grafted onto the core of the colloid and its linear charge density is varied, whereas the linear charge density of the A block is fixed. The brush structure of a single colloid, the mean force between two colloids, and the structure of solutions of such colloids have been determined for different linear charge densities of the B blocks and block lengths. Many features of the present system are controlled by the charge of the B blocks. In the limit of uncharged B blocks, (i) the grafted chains are stretched and form an extended polyelectrolyte brush, (ii) a strong repulsive force is operating between two colloids, (iii) and the solution is thermodynamic stable and displays strong spatial correlation among the colloids. In the limit where the charges of the two types of blocks exactly compensate each other, (i) the chains are collapsed and form a polyelectrolyte complex surrounding the cores, (ii) an attractive force appears between two colloids, and (iii) strong colloid clustering appears in the solution. These features become more pronounced as the length of the polymer blocks is increased, and a phase instability occurs at sufficiently long chains. A comparison with properties for other related colloidal particles is also provided.

Monte Carlo simulations of a charged dendrimer with explicit counterions and salt ions

Static properties of a dendrimer with generation g = 5 with positively charged terminal groups in an athermal solvent are studied by lattice Monte Carlo simulations using the cooperative motion algorithm as the tossing scheme. The calculations are performed both for a salt-free system with neutralizing counterions and for a small amount of added monovalent and divalent salt. The full Coulomb potential and the excluded volume interactions between ions and beads are taken explicitly into account with the reduced temperature tau, the number of salt cations (anions) n(s), and salt valence z(s) as the simulation parameters. The bahaviour of the systems is analyzed by the mean effective charge per end-bead <Q>, Coulomb mean energy <E>, mean-square radius of gyration <R(g)(2)>, pair correlation functions g(alphabeta), and charge density rho(ch). The simulations show that for n(s)> or = 0 and decreasing tau: (a) there is encapsulation in the dendrimer and condensation onto the terminal groups of anions accompanied by a monotonic decrease in <Q> and <E> and by subsequent swelling and shrinking of the molecule; (b) encapsulation, condensation and shrinking are the most significant and swelling weaker for |z(s)| = 2; (c) penetration of salt cations into the dendrimer is minor when compared to that of anions; (d) rho(ch) is reduced and becomes negative close to the center of mass of the dendrimer and on its periphery; (e) for the considered n(s) > 0, unlike divalent salt ions the monovalent ones cause slight effects when compared to the salt-free case.

Structure of polyelectrolyte complexes by Brownian dynamics simulation: Effects of the bond length asymmetry of the polyelectrolytes

Brownian dynamics simulations were performed to study the structure of polyelectrolyte complexes formed by two flexible, oppositely charged polyelectrolyte chains. The distribution of monomers in the complex as well as the radius of gyration and structure factor of complexes and individual polyelectrolytes are reported. These structural properties were calculated for polyelectrolyte chains with equal number of monomers, keeping constant the bond length of the negative chain and increasing the bond length of the positive chain. This introduces an asymmetry in the length of the chains that modulates the final structure of the complexes. In the symmetric case the distribution of positive and negative monomers in the complex are identical, producing clusters that are locally and globally neutral. Deviations from the symmetric case lead to nonuniform, asymmetric monomer distributions, producing net charge oscillations inside the complex and large changes in the radius of gyration of the complex and individual chains. From the radius of gyration of the polyelectrolyte chains it is shown that the positive chain is much more folded than the negative chain when the chains are asymmetric, which is also confirmed through the scaling behavior of the structure factors.

Coil-globule coexistence and compaction of DNA chains

In this work we discuss different factors governing coil-globule coexistence in the compaction process of DNA. We initially analyse the role played by fluctuations in the degree of binding of an external compacting agent in the conformational behavior of the chain backbone. The analysis relies both on Monte Carlo simulation results and simple statistical approaches. Compacting agents of various binding characteristics are taken into consideration and the degree of charge neutralization upon the chain is related to conformational indicators. Selected model systems comprising stiff chains in the presence of multivalent ions are employed to assess intrinsic single-chain conformational fluctuation, in the presence of external agents but not resulting from differences in binding. It is shown that trends found for a variety of compacting agents, including the extension of the coil-globule coexistence regions, can be rationalised on the basis of this analysis.

Complexation of semiflexible chains with oppositely charged cylinder

We study the complexation of long thin semiflexible polymer chains with an oppositely charged cylinder. Starting from the linear Poisson-Boltzmann equation, we calculate the electrostatic potential and the energy of such a charge distribution. We find that sufficiently flexible chains prefer to wrap around the cylinder in a helical manner, when their charge density is smaller than that of the cylinder. The optimal value of the helical pitch is found by minimization of the sum of electrostatic and bending energies. The dependence of the pitch on the number of chains, their rigidity, and salt concentration in solution is analyzed. We discuss our results in the light of recent experiments on DNA complexation with cylindrical dendronized polymers.

Entropy and enthalpy of polyelectrolyte complexation: Langevin dynamics simulations

We report a systematic study by Langevin dynamics simulation on the energetics of complexation between two oppositely charged polyelectrolytes of same charge density in dilute solutions of a good solvent with counterions and salt ions explicitly included. The enthalpy of polyelectrolyte complexation is quantified by comparisons of the Coulomb energy before and after complexation. The entropy of polyelectrolyte complexation is determined directly from simulations and compared with that from a mean-field lattice model explicitly accounting for counterion adsorption. At weak Coulomb interaction strengths, e.g., in solvents of high dielectric constant or with weakly charged polyelectrolytes, complexation is driven by a negative enthalpy due to electrostatic attraction between two oppositely charged chains, with counterion release entropy playing only a subsidiary role. In the strong interaction regime, complexation is driven by a large counterion release entropy and opposed by a positive enthalpy change. The addition of salt reduces the enthalpy of polyelectrolyte complexation by screening electrostatic interaction at all Coulomb interaction strengths. The counterion release entropy also decreases in the presence of salt, but the reduction only becomes significant at higher Coulomb interaction strengths. More significantly, in the range of Coulomb interaction strengths appropriate for highly charged polymers in aqueous solutions, complexation enthalpy depends weakly on salt concentration and counterion release entropy exhibits a large variation as a function of salt concentration. Our study quantitatively establishes that polyelectrolyte complexation in highly charged Coulomb systems is of entropic origin.

Monte Carlo simulations of oppositely charged macroions in solution

The structure and phase behavior of oppositely charged macroions in solution have been studied with Monte Carlo simulations using the primitive model where the macroions and small ions are described as charged hard spheres. Size and charge symmetric, size asymmetric, and charge asymmetric macroions at different electrostatic coupling strengths are considered, and the properties of the solutions have been examined using cluster size distribution functions, structure factors, and radial distribution functions. At increasing electrostatic coupling, the macroions form clusters and eventually the system displays a phase instability, in analogy to that of simple electrolyte solutions. The relation to the similar cluster formation and phase instability occurring in solutions containing oppositely charged polymers is also discussed.

Diblock polyampholytes grafted onto spherical particles: Monte Carlo simulation and lattice mean-field theory

Spherical brushes composed of diblock polyampholytes (diblock copolymers with oppositely charged blocks) grafted onto solid spherical particles in aqueous solution are investigated by using the primitive model solved with Monte Carlo simulations and by lattice mean-field theory. Polyampholyte chains of two compositions are considered: a copolymer with a long and a short block, A100B10, and a copolymer with two blocks of equal length, A50B50. The B block is end-grafted onto the surface, and its charge is varied, whereas the charge of the A block is fixed. Single-chain properties, radial and lateral spatial distributions of different types, and structure factors are analyzed. The brush structure strongly depends on the charge of the B block. In the limit of an uncharged B block, the chains are stretched and form an extended polyelectrolyte brush. In the other limit with the charges of the blocks compensating each other, the chains are collapsed and form a polyelectrolyte complex surrounding the particles. At intermediate charge conditions, a polyelectrolyte brush and a polyelectrolyte complex coexist and constitute two substructures of the spherical brush. The differences of the brush structures formed by the A100B10 and A50B50 polyampholytes are also analyzed. Finally, a comparison of the predictions of the two theoretical approaches is made.

Lattice Monte Carlo simulations of three-dimensional charged polymer chains

The configurational properties of strongly charged polyelectrolytes accompanied by neutralizing counterions in dilute solutions are simulated using the cooperative motion algorithm on the face-centered-cubic lattice. The full Coulomb potential and the excluded volume condition between different ions/beads are taken into account and the reduced temperature T* is considered the main, variable parameter. The calculations that have been carried out for solutions of both single and several chains indicate a few regions of their behavior: (1) for T*--> infinity, it corresponds to that of neutral, self-avoiding polymers under good solvent conditions; (2) for T* approximately 1, due to the electrostatic interactions being effectively stronger, the chains are more outstretched compared to their size at other temperatures; (3) for T* well below one, the counterion condensation becomes more and more dominant, which gradually leads to strongly collapsed chains; and (4) at the lowest temperatures the chains and counterions assume low-energy configurations in the form of neutral, compact aggregates.