Scaling of Star Polymers with 1−80 Arms

Entropic exponents of grafted lattice stars

The surface entropic exponents of half-space lattice stars grafted at their central nodes in a hard wall are estimated numerically using the PERM algorithm. In the square half-lattice the exact values of the exponents are verified, including Barber's scaling relation and a generalization for 2-stars with one and two surface loops respectively. This is the relation γ_{211}=2γ_{21}-γ_{20}, where γ_{21} and γ_{211} are the surface entropic exponents of a grafted 2-star with one and two surface loops, respectively, and γ_{20} is the surface entropic exponent with no surface loops. This relation is also tested in the cubic half-lattice where surface entropic exponents are estimated up to 5-stars, including many with one or more surface loops. Barber's scaling relation and the relation γ_{3111}=γ_{30}-3γ_{31}+3γ_{311} are also tested, where the exponents {γ_{31},γ_{311},γ_{3111}} are of grafted 3-stars with one, two, or three surface loops, respectively, and γ_{30} is the surface exponent of grafted 3-stars.

Entropic force of cone-tethered polymers interacting with a planar surface

Computer simulations are used to characterize the entropic force of one or more polymers tethered to the tip of a hard conical object that interact with a nearby hard flat surface. Pruned-enriched Rosenbluth method Monte Carlo simulations are used to calculate the variation of the conformational free energy F of a hard-sphere polymer with respect to the cone-tip-to-surface distance h from which the variation of the entropic force f≡|dF/dh| with h is determined. We consider the following cases: a single freely jointed tethered chain, a single semiflexible tethered chain, and several freely jointed chains of equal length each tethered to the cone tip. The simulation results are used to test the validity of a prediction by Maghrebi et al. [Maghrebi et al., Europhys. Lett. 96, 66002 (2011)0295-507510.1209/0295-5075/96/66002; Phys. Rev. E 86, 061801 (2012)1539-375510.1103/PhysRevE.86.061801] that f∝(γ_{∞}-γ_{0})h^{-1}, where γ_{0} and γ_{∞} are universal scaling exponents for the partition function of the tethered polymer for h=0 and h=∞, respectively. The measured functions f(h) are generally consistent with the predictions, with small quantitative discrepancies arising from the approximations employed in the theory. In the case of multiple tethered polymers, the entropic force per polymer is roughly constant, which is qualitatively inconsistent with the predictions.

Star Polymers vs. Dendrimers: Studies of the Synthesis Based on Computer Simulations

A generic model was developed for studies of the polymerization process of regular branched macromolecules. Monte Carlo simulations were performed employing the Dynamic Lattice Liquid algorithm to study this process. A core-first methodology was used in a living polymerization of stars with up to 32 arms, and dendrimers consisted of 4-functional segments. The kinetics of the synthesis process for stars with different numbers of branches and dendrimers was compared. The size and structure of star-branched polymers and dendrimers during the synthesis were studied. The influence of the functionality of well-defined cores on the structure and on the dispersity of the system was also examined. The differences in the kinetics in the formation of both architectures, as well as changes to their structures, were described and discussed.

Numerical estimates of square lattice star vertex exponents

We implement parallel versions of the generalized atmospheric Rosenbluth methods and Wang-Landau algorithms for stars and for acyclic uniform branched networks in the square lattice. These are models of monodispersed branched polymers, and we estimate the star vertex exponents σ_{f} for f stars, and the entropic exponent γ_{G} for networks with comb and brush connectivity in two dimensions. Our results verify the predicted (but not rigorously proven) exact values of the vertex exponents and we test the scaling relation [B. Duplantier, J. Stat. Phys. 54, 581 (1989)JSTPBS0022-471510.1007/BF01019770]γ_{G}-1=[under ∑]f≥1m_{f}σ_{f}for several acyclic branched networks in two dimensions.

The influence of arm composition on the self-assembly of low-functionality telechelic star polymers in dilute solutions

We combine synthesis, physical experiments, and computer simulations to investigate self-assembly patterns of low-functionality telechelic star polymers (TSPs) in dilute solutions. In particular, in this work, we focus on the effect of the arm composition and length on the static and dynamic properties of TSPs, whose terminal blocks are subject to worsening solvent quality upon reducing the temperature. We find two populations, single stars and clusters, that emerge upon worsening the solvent quality of the outer block. For both types of populations, their spatial extent decreases with temperature, with the specific details (such as temperature at which the minimal size is reached) depending on the coupling between inter- and intra-molecular associations as well as their strength. The experimental results are in very good qualitative agreement with coarse-grained simulations, which offer insights into the mechanism of thermoresponsive behavior of this class of materials.

The Daoud and Cotton blob model and the interaction of star-shaped polymers

Since it was first proposed in 1982, the Daoud and Cotton (DC) model for star-shaped polymers was intensively used also for self-assembled copolymers and small colloids grafted with long polymers. We try to clarify the position of the DC model and focus on the star partition function which plays a central role in self-assembly and gives access to the star-star interaction. While the predicted star-star interaction agrees with scattering data by Likos et al. (Phys. Rev. Lett. 80, 4450 (1998)), an extensive simulation by Hsu et al. (Macromolecules, 37, 4658 (2004)) does not recover the prediction for the partition function. We try to reconcile this seemingly conflicting results. We discuss star-star interactions, star free energy in θ -solvents, mixing of A/B branches in copolymer stars, within or beyond the Daoud and Cotton blob model.

Thermodynamics of star polymer solutions: A coarse-grained study

We consider a coarse-grained (CG) model with pairwise interactions, suitable to describe low-density solutions of star-branched polymers of functionality f. Each macromolecule is represented by a CG molecule with (f + 1) interaction sites, which captures the star topology. Potentials are obtained by requiring the CG model to reproduce a set of distribution functions computed in the microscopic model in the zero-density limit. Explicit results are given for f = 6, 12, and 40. We use the CG model to compute the osmotic equation of state of the solution for concentrations c such that Φ_p=c∕c^*≲1, where c^* is the overlap concentration. We also investigate in detail the phase diagram for f = 40, identifying the boundaries of the solid intermediate phase. Finally, we investigate how the polymer size changes with c. For Φ_p≲0.3, polymers become harder as f increases at fixed reduced concentration c∕c^*. On the other hand, for Φ_p≳0.3, polymers show the opposite behavior: At fixed Φ_p, the larger the value of f, the larger their size reduction is.

Polymer models with optimal good-solvent behavior

We consider three different continuum polymer models, which all depend on a tunable parameter r that determines the strength of the excluded-volume interactions. In the first model, chains are obtained by concatenating hard spherocylinders of height b and diameter rb (we call them thick self-avoiding chains). The other two models are generalizations of the tangent hard-sphere and of the Kremer-Grest models. We show that for a specific value [Formula: see text], all models show optimal behavior: asymptotic long-chain behavior is observed for relatively short chains. For [Formula: see text], instead, the behavior can be parametrized by using the two-parameter model, which also describes the thermal crossover close to the θ point. The bonds of the thick self-avoiding chains cannot cross each other, and therefore the model is suited for the investigation of topological properties and for dynamical studies. Such a model also provides a coarse-grained description of double-stranded DNA, so that we can use our results to discuss under which conditions DNA can be considered as a model good-solvent polymer.

Loop statistics in polymers in crowded environment

We analyze the probability to find a single loop in a long flexible polymer chain in disordered environment in d dimensions. The structural defects are considered to be correlated on large distances r according to a power law ∼r(-a). Working within the frames of continuous chain model and applying the direct polymer renormalization scheme, we obtain the values of critical exponents governing the scaling of probabilities to find the loops of various positions along the chain as function of loops' length. Our results quantitatively reveal that the presence of structural defects in environment decreases the probability of loop formation in polymer macromolecules.

Monte Carlo simulations of lattice models for single polymer systems

Single linear polymer chains in dilute solutions under good solvent conditions are studied by Monte Carlo simulations with the pruned-enriched Rosenbluth method up to the chain length N~O(10(4)). Based on the standard simple cubic lattice model (SCLM) with fixed bond length and the bond fluctuation model (BFM) with bond lengths in a range between 2 and √10, we investigate the conformations of polymer chains described by self-avoiding walks on the simple cubic lattice, and by random walks and non-reversible random walks in the absence of excluded volume interactions. In addition to flexible chains, we also extend our study to semiflexible chains for different stiffness controlled by a bending potential. The persistence lengths of chains extracted from the orientational correlations are estimated for all cases. We show that chains based on the BFM are more flexible than those based on the SCLM for a fixed bending energy. The microscopic differences between these two lattice models are discussed and the theoretical predictions of scaling laws given in the literature are checked and verified. Our simulations clarify that a different mapping ratio between the coarse-grained models and the atomistically realistic description of polymers is required in a coarse-graining approach due to the different crossovers to the asymptotic behavior.

Combined collapse by bridging and self-adhesion in a prototypical polymer model inspired by the bacterial nucleoid

Recent experimental results suggest that the E. coli chromosome feels a self-attracting interaction of osmotic origin, and is condensed in foci by bridging interactions. Motivated by these findings, we explore a generic modeling framework combining solely these two ingredients, in order to characterize their joint effects. Specifically, we study a simple polymer physics computational model with weak ubiquitous short-ranged self attraction and stronger sparse bridging interactions. Combining theoretical arguments and simulations, we study the general phenomenology of polymer collapse induced by these dual contributions, in the case of regularly spaced bridging. Our results distinguish a regime of classical Flory-like coil-globule collapse dictated by the interplay of excluded volume and attractive energy and a switch-like collapse where bridging interactions compete with entropy loss terms from the looped arms of a star-like rosette. Additionally, we show that bridging can induce stable compartmentalized domains. In these configurations, different "cores" of bridging proteins are kept separated by star-like polymer loops in an entropically favorable multi-domain configuration, with a mechanism that parallels micellar polysoaps. Such compartmentalized domains are stable, and do not need any intra-specific interactions driving their segregation. Domains can be stable also in the presence of uniform attraction, as long as the uniform collapse is above its theta point.

Degradation versus self-assembly of block co-polymer micelles

The stability of micelles self-assembled from block co-polymers can be altered by the degradation of the blocks. Slow degradation shifts the equilibrium size distribution of block co-polymer micelles and changes their properties. The quasi-equilibrium scaling theory shows that the degradation of hydrophobic blocks in the core of micelles destabilizes the micelles, reducing their size, while the degradation of hydrophilic blocks forming coronas of micelles favors larger micelles and may, at certain conditions, induce the formation of micelles from individual chains.

N-Body interactions in soft-sphere coarse-grained models of star polymers

We consider the coarse-grained approach in which star polymers are mapped onto atoms located at the centers of the stars interacting by means of ultrasoft potentials. We generalize the Witten-Pincus formula to all n-body potentials: in the good-solvent regime all potentials show an ultrasoft logarithmic divergence when the relative distance of all n stars goes to zero, with coefficients that can be related to the partition-function exponents γ(f).

Looped star polymers show conformational transition from spherical to flat toroidal shapes

Inspired by the topological organization of the circular Escherichia coli chromosome, which is compacted by separate domains, we study a polymer architecture consisting of a central ring to which either looped or linear side chains are grafted. A shape change from a spherical to a toroidal organization takes place as soon as the inner ring becomes large enough for the attached arms to fit within its circumference. Building up a torus, the system flattens, depending on the effective bending rigidity of the chain induced by entropic repulsion of the attached loops and, to a lesser extent, linear arms. Our results suggest that the natural formation of a toroidal structure with a decreased amount of writhe induced by a specific underlying topology could be one driving force, among others, that nature exploits to ensure proper packaging of the genetic material within a rod-shaped, bacterial envelope.

Aggregation of amphiphilic polymers in the presence of adhesive small colloidal particles

The interaction of amphiphilic polymers with small colloids, capable to reversibly stick onto the chains, is studied. Adhesive small colloids in solution are able to dynamically bind two polymer segments. This association leads to topological changes in the polymer network configurations, such as looping and cross-linking, although the reversible adhesion permits the colloid to slide along the chain backbone. Previous analyses only consider static topologies in the chain network. We show that the sliding degree of freedom ensures the dominance of small loops, over other structures, giving rise to a new perspective in the analysis of the problem. The results are applied to the analysis of the equilibrium between colloidal particles and star polymers, as well as to block copolymer micelles. The results are relevant for the reversible adsorption of silica particles onto hydrophilic polymers, used in the process of formation of mesoporous materials of the type SBA or MCM, cross-linked cyclodextrin molecules threading on the polymers and forming the structures known as polyrotaxanes. Adhesion of colloids on the corona of the latter induce micellization and growth of larger micelles as the number of colloids increase, in agreement with experimental data.

Star copolymers in porous environments: scaling and its manifestations

We consider star polymers, consisting of two different polymer species, in a solvent subject to quenched correlated structural obstacles. We assume that the disorder is correlated with a power-law decay of the pair-correlation function g(x)~x(-a). Applying the field-theoretical renormalization group approach in d dimensions, we analyze different scenarios of scaling behavior working to first order of a double ɛ=4-d, δ=4-a expansion. We discuss the influence of the correlated disorder on the resulting scaling laws and possible manifestations such as diffusion-controlled reactions in the vicinity of absorbing traps placed on polymers as well as the effective short-distance interaction between star copolymers.

Elastic lattice polymers

We study a model of "elastic" lattice polymer in which a fixed number of monomers m is hosted by a self-avoiding walk with fluctuating length l . We show that the stored length density ρm≡1-<l>/m scales asymptotically for large m as ρm=ρ∞(1-θ/m+…) , where θ is the polymer entropic exponent, so that θ can be determined from the analysis of ρm. We perform simulations for elastic lattice polymer loops with various sizes and knots, in which we measure ρm. The resulting estimates support the hypothesis that the exponent θ is determined only by the number of prime knots and not by their type. However, if knots are present, we observe strong corrections to scaling, which help to understand how an entropic competition between knots is affected by the finite length of the chain.

Simulation of dilute solutions of linear and star-branched polymers by dissipative particle dynamics

A most promising off-lattice technique in order to simulate not only static but in addition dynamic behavior of linear and star-branched chains is the dissipative particle dynamics (DPD) method. In this model the atomistic representation of polymer molecules is replaced by a (coarse-grained) equivalent chain consisting of beads which are repulsive for each other in order to mimic the excluded volume effect (successive beads in addition are linked by springs). Likewise solvent molecules are combined to beads which in turn are repulsive for each other as well as for the polymer segments. The system is relaxed by molecular dynamics solving Newton's laws under the influence of short ranged conservative forces (i.e., repulsion between nonbonded beads and a proper balance of repulsion and attraction between bonded segments) and dissipative forces due to friction between particles, the latter representing the thermostat in conjunction with proper random forces. A variation of the strength of the repulsion between different types of beads allows the simulation of any desired thermodynamic situation. Static and dynamic properties of isolated linear and star-branched chains embedded in athermal, exothermal, and endothermal solvent are presented and theta conditions are examined. The generally accepted scaling concept for athermal systems is fairly well reproduced by linear and star-branched DPD chains and theta conditions appear for a unique parameter independent of functionality as in the case of Monte Carlo simulations. Furthermore, the correspondence between DPD and Monte Carlo data referring to the shape of chains and stars is fairly well, too. For dilute solutions the Zimm behavior is expected for dynamic properties which is indeed realized in DPD systems.

Coarse-grained models for semi-dilute polymer solutions under good-solvent conditions

We determine coarse-grained models with a relatively small number of units which reproduce the universal behavior of polymer solutions in the semi-dilute regime under good-solvent conditions. We check both the thermodynamical behavior (osmotic pressure and chemical potential) and structural properties (polymer size).

Single DNA denaturation and bubble dynamics

While the Watson-Crick double-strand is the thermodynamically stable state of DNA in a wide range of temperature and salt conditions, even at physiological conditions local denaturation bubbles may open up spontaneously due to thermal activation. By raising the ambient temperature, titration, or by external forces in single molecule setups bubbles proliferate until full denaturation of the DNA occurs. Based on the Poland-Scheraga model we investigate both the equilibrium transition of DNA denaturation and the dynamics of the denaturation bubbles with respect to recent single DNA chain experiments for situations below, at, and above the denaturation transition. We also propose a new single molecule setup based on DNA constructs with two bubble zones to measure the bubble coalescence and extract the physical parameters relevant to DNA breathing. Finally we consider the interplay between denaturation bubbles and selectively single-stranded DNA binding proteins.

Escape transition of a polymer chain from a nanotube: how to avoid spurious results by use of the force-biased pruned-enriched Rosenbluth algorithm

A polymer chain containing N monomers confined in a finite cylindrical tube of diameter D grafted at a distance L from the open end of the tube may undergo a rather abrupt transition, where part of the chain escapes from the tube to form a "crownlike" coil outside of the tube. When this problem is studied by Monte Carlo simulation of self-avoiding walks on the simple cubic lattice applying a cylindrical confinement and using the standard pruned-enriched Rosenbluth method (PERM), one obtains spurious results, however, with increasing chain length the transition gets weaker and weaker, due to insufficient sampling of the "escaped" states, as a detailed analysis shows. In order to solve this problem, a new variant of a biased sequential sampling algorithm with resampling is proposed, force-biased PERM: the difficulty of sampling both phases in the region of the first order transition with the correct weights is treated by applying a force at the free end pulling it out of the tube. Different strengths of this force need to be used and reweighting techniques are applied. Using rather long chains (up to N=18000 ) and wide tubes (up to D=29 lattice spacings), the free energy of the chain, its end-to-end distance, the number of "imprisoned" monomers can be estimated, as well as the order parameter and its distribution. It is suggested that this algorithm should be useful for other problems involving state changes of polymers, where the different states belong to rather disjunct "valleys" in the phase space of the system.

Ring polymer simulations with global radius of curvature

We simulate three-dimensional flexible off-lattice ring polymers of length L up to L=4000 for various values of the global radius of curvature Rgrc=0.25 , 0.48, and 1.0 and Rgrc=2.0 . We utilize two different ensembles: one with a delta -function constraint on the radius, and the other with a theta -function. For both cases the global radius of curvature provides a valid regularization of polymers with thickness D=2Rgrc . The Flory-type critical exponent nu SAW of self-avoiding rings at D=2 is found to be nu SAW=0.5869(5) from the radii of gyration chain length scaling, while other D values produce consistent results. For our current implementation, the numerical effort of chain thickness calculations is bounded by a number O(LlnL) per single update. We also study low-temperature configurations of spatially dense Lennard-Jones homopolymers on a ring and identify some conformational building blocks.

Discontinuous molecular dynamics simulation study of polymer collapse

Discontinuous molecular dynamics simulations were used to study the coil-globule transition of a polymer in an explicit solvent. Two different versions of the model were employed, which are differentiated by the nature of monomer-solvent, solvent-solvent, and nonbonded monomer-monomer interactions. For each case, a model parameter lambda determines the degree of hydrophobicity of the monomers by controlling the degree of energy mismatch between the monomers and solvent particles. We consider a lambda-driven coil-globule transition at constant temperature. The simulations are used to calculate average static structure factors, which are then used to determine the scaling exponents of the system in order to determine the theta-point values lambda(theta) separating the coil from the globule state. For each model we construct coil-globule phase diagrams in terms of lambda and the particle density rho. Additionally, we explore for each model the effects of varying the range of the attractive interactions on the phase boundary separating the coil and globule phases. The results are analyzed in terms of a simple Flory-type theory of the collapse transition.

Virial coefficients and osmotic pressure in polymer solutions in good-solvent conditions

We determine the second, third, and fourth virial coefficients appearing in the density expansion of the osmotic pressure Pi of a monodisperse polymer solution in good-solvent conditions. Using the expected large-concentration behavior, we extrapolate the low-density expansion outside the dilute regime, obtaining the osmotic pressure for any concentration in the semidilute region. Comparison with field-theoretical predictions and experimental data shows that the obtained expression is quite accurate. The error is approximately 1%-2% below the overlap concentration and rises at most to 5%-10% in the limit of very large polymer concentrations.

Polymer size in dilute solutions in the good-solvent regime

We determine the density expansion of the radius of gyration, of the hydrodynamic radius, and of the end-to-end distance for a monodisperse polymer solution in good-solvent conditions. We consider the scaling limit (large degree of polymerization), including the leading scaling corrections. Using the expected large-concentration behavior, we extrapolate these low-density expansions outside the dilute regime, obtaining a prediction for the radii for any concentration in the semidilute region. For the radius of gyration, comparison with field-theoretical predictions shows that the relative error should be at most 5% in the limit of very large polymer concentrations.

Entropy-induced separation of star polymers in porous media

We present a quantitative picture of the separation of star polymers in a solution where part of the volume is influenced by a porous medium. To this end, we study the impact of long-range-correlated quenched disorder on the entropy and scaling properties of f-arm star polymers in a good solvent. We assume that the disorder is correlated on the polymer length scale with a power-law decay of the pair correlation function g(r) approximately r-a. Applying the field-theoretical renormalization group approach we show in a double expansion in epsilon=4-d and delta=4-a that there is a range of correlation strengths delta for which the disorder changes the scaling behavior of star polymers. In a second approach we calculate for fixed space dimension d=3 and different values of the correlation parameter a the corresponding scaling exponents gammaf that govern entropic effects. We find that gammaf-1, the deviation of gammaf from its mean field value is amplified by the disorder once we increase delta beyond a threshold. The consequences for a solution of diluted chain and star polymers of equal molecular weight inside a porous medium are that star polymers exert a higher osmotic pressure than chain polymers and in general higher branched star polymers are expelled more strongly from the correlated porous medium. Surprisingly, polymer chains will prefer a stronger correlated medium to a less or uncorrelated medium of the same density while the opposite is the case for star polymers.

Corrections to scaling in multicomponent polymer solutions

We calculate the correction-to-scaling exponent omegaT that characterizes the approach to the scaling limit in multicomponent polymer solutions. A direct Monte Carlo determination of omegaT in a system of interacting self-avoiding walks gives omegaT=0.415+/-0.020. A field-theory analysis based on five- and six-loop perturbative series leads to omegaT=0.41+/-0.04. We also verify the renormalization-group predictions for the scaling behavior close to the ideal-mixing point.