Tunable effective interactions between dendritic macromolecules

We employ extensive Monte Carlo and molecular-dynamics simulations to investigate the effective interactions between the centers of mass of dendritic macromolecules of variable flexibility and generation number. Two different models for the connectivity and steric interactions between the monomers are employed, the first one being purely entropic in nature and the second explicitly involving energetic interactions. We find that the effective potentials have a generic Gaussian shape, whose range and strength can be tuned via modifications in the generation number and flexibility of the spacers. We supplement our simulation analysis by a density-functional approach in which the connectivity between the monomers is approximated by an external confining potential that holds the monomer beads together. Using a simple density functional for the interactions between the monomers, we find semiquantitative agreement between theory and simulation. The implications of our findings for the interpretation of scattering data from concentrated dendrimer solutions are also discussed.     

Counterion distributions and effective interactions of spherical polyelectrolyte brushes

We consider the electrosteric repulsion of colloidal particles whose surface carries a dense layer of long polyelectrolyte chains (spherical polyelectrolyte brushes). The theory of electrosteric repulsion of star polyelectrolytes developed recently is augmented to include particles with a finite core radius. It is shown that most of the counterions are confined within the brush layer. The strong osmotic pressure thus created within the brush layer dominates the repulsive interaction between two such particles. Because of this the pair interaction potential between spherical polyelectrolyte brushes can be given in terms of an analytic expression. The theoretical predictions are compared with available experimental data and semi-quantitative agreement between the two is found.     

Tailoring the flow of soft glasses by soft additives

We examine the vitrification and melting of asymmetric star polymer mixtures by combining rheological measurements with mode coupling theory. We identify two types of glassy states, a single glass, in which the small component is fluid in the glassy matrix of the big one, and a double glass, in which both components are vitrified. Addition of small-star polymers leads to melting of both glasses, and the melting curve has a nonmonotonic dependence on the star-star size ratio. The phenomenon opens new ways for externally steering the rheological behavior of soft matter systems.     

Phase behavior of ionic microgels

We employ effective interaction potentials between spherical polyelectrolyte microgels in order to investigate theoretically the structure, thermodynamics, and phase behavior of ionic microgel solutions. Combining a genetic algorithm with accurate free energy calculations we are able to perform an unrestricted search of candidate crystal structures. Hexagonal, body-centered orthogonal, and trigonal crystals are found to be stable at high concentrations and charges of the microgels, accompanied by reentrant melting behavior and fluid-fcc-bcc transitions below the overlap concentration.     

Polyelectrolyte-compression forces between spherical DNA brushes

Optical tweezers are employed to measure the forces of interaction within a single pair of DNA-grafted colloids, dependent on the molecular weight of the DNA chains, and the concentration and valence of the surrounding ionic medium. The resulting forces are short range and set in as the surface-to-surface distance between the colloidal cores reaches the value of the brush height. The measured force-distance relation is analyzed by means of a theoretical treatment that quantitatively describes the effects of compression of the chains on the surface of the opposite-lying colloid. Quantitative agreement with the experiment is obtained for all parameter combinations.     

Polyelectrolyte stars in planar confinement

We employ monomer-resolved molecular dynamics simulations and theoretical considerations to analyze the conformations of multiarm polyelectrolyte stars close to planar, uncharged walls. We identify three mechanisms that contribute to the emergence of a repulsive star-wall force, namely, the confinement of the counterions that are trapped in the star interior, the increase in electrostatic energy due to confinement as well as a novel mechanism arising from the compression of the stiff polyelectrolyte rods approaching the wall. The latter is not present in the case of interaction between two polyelectrolyte stars and is a direct consequence of the impenetrable character of the planar wall.     

Formation of polymorphic cluster phases for a class of models of purely repulsive soft spheres

We present results from density functional theory and computer simulations that unambiguously predict the occurrence of first-order freezing transitions for a large class of ultrasoft model systems into cluster crystals. The clusters consist of fully overlapping particles and arise without the existence of attractive forces. The number of particles participating in a cluster scales linearly with density, therefore the crystals feature density-independent lattice constants. Clustering is accompanied by polymorphic bcc-fcc transitions, with fcc being the stable phase at high densities.     

Soft effective interactions between weakly charged polyelectrolyte chains

We apply extensive molecular dynamics simulations and analytical considerations in order to study the conformations and the effective interactions between weakly charged, flexible polyelectrolyte chains in salt-free conditions. We focus on charging fractions lying below 20%, for which case there is no Manning condensation of counterions and the latter can be thus partitioned in two states: those that are trapped within the region of the flexible chain and the ones that are free in the solution. We examine the partition of counterions in these two states, the chain sizes and the monomer distributions for various chain lengths, finding that the monomer density follows a Gaussian shape. We calculate the effective interaction between the centers of mass of two interacting chains, under the assumption that the chains can be modeled as two overlapping Gaussian charge profiles. The analytical calculations are compared with measurements from molecular dynamics simulations. Good quantitative agreement is found for charging fractions below 10%, where the chains assume coil-like configurations, whereas deviations develop for charge fraction of 20%, in which case a conformational transition of the chain towards a rodlike configuration starts to take place.