Arben Jusufi

Counterion-induced entropic interactions in solutions of strongly stretched, osmotic polyelectrolyte stars

We examine the conformations and effective interactions of star-branched polyelectrolytes with and without added salt, by employing monomer-resolved molecular dynamics simulations and an analytical theory. The simulations take into account the excluded-volume and Coulomb interactions between the individual monomers, as well as the counter- and coions. The theory is based on a variational free energy that is written as a sum of electrostatic, polymer, and entropic contributions of the counter- and coions. For the conformations of isolated polyelectrolyte stars, we find strong stretching of the chains, resulting in a linear scaling of the star radius with the degree of polymerization, as well as trapping and condensation of a large fraction of counterions. The effective interactions at arbitrarily strong overlaps between the stars are shown to be dominated by the entropic contributions of the trapped counterions, with the electrostatic contribution playing only a minor role due to an almost complete neutral...

Self-assembly of coarse-grained ionic surfactants accelerated by graphics processing units

Due to the relatively long time scales inherent to ionic surfactant self-assembly (>μs), an aggressive computational approach is needed to obtain converged data on micellar solutions. This work presents a study of micellization using a coarse-grained (CG) model of aqueous ionic surfactants in replicated molecular dynamics (MD) simulations run on graphics processing unit hardware. The performance of our implementation of the CG model with electrostatics into the HOOMD-Blue GPU-accelerated MD software package is comparable to that of a modest sized cluster running a highly optimized parallel CPU code. From 0.36 ms of cumulative trajectory data, we are able to predict equilibrium thermodynamic and morphological properties of ionic surfactant micellar solutions. Estimating the critical micelle concentrations (CMC) from the free monomer (ρ1) and premicellar concentrations obtained from simulations of sodium hexyl sulfate (S6S, CMC of 460 ± 6 mM) at high (1 M) concentration, a value in good agreement with experimental results is obtained; however, the same method applied to simulations of sodium nonyl sulfate (S9S, ρ1 of 2.4 ± 0.01 mM) and sodium dodecyl sulfate (SDS, ρ1 of 0.02 ± 0.01 mM) at the same total concentration systematically underestimates the CMCs. An alternative method for calculating the CMC is presented, where the free monomer concentration computed from high concentration CG-MD data is used as the input to a simple theoretical model which can be used to extrapolate to a more accurate prediction of the CMC. Better agreement between the empirical and predicted CMC is obtained from this theory for S9S (28.7 ± 0.3 mM) and SDS (3.32 ± 0.04 mM), though the CMC for S6S is slightly underestimated (304 ± 3 mM). We also present statistically converged morphological data, including aggregation number distributions and the principal components of the gyration tensor. This data suggest a transition from spherical micelles to rod-like at a specific aggregation number, which increases with increasing hydrocarbon length.

Effective interactions between star polymers and colloidal particles

Using monomer-resolved molecular dynamics simulations and theoretical arguments based on the radial dependence of the osmotic pressure in the interior of a star, we systematically investigate the effective interactions between hard, colloidal particles and star polymers in a good solvent. The relevant parameters are the size ratio q between the stars and the colloids, as well as the number of polymeric arms f (functionality) attached to the common centre of the star. By covering a wide range of qs ranging from zero (star against a flat wall) up to about 0.75, we establish analytical forms for the star-colloid interaction which are in excellent agreement with simulation results. A modified expression for the star-star interaction for low functionalities, f10 is also introduced.

Effective potentials for 1:1 electrolyte solutions incorporating dielectric saturation and repulsive hydration

Implicit water potentials are developed for the study of thermodynamic and structural properties of solutions of NaCl, LiCl, and KCl. The interaction potential between cations and anions is parametrized from the ionic crystal potential. Two short-range corrections were added to the system to account for the water solvent. The first is due to dielectric saturation which reduces the dielectric permittivity in the vicinity of an ion. The second is a repulsive Gaussian potential which represents the first hydration shell around the ions. Grand canonical Monte Carlo simulations were performed to calculate the mean ionic activity coefficients. Molecular dynamics simulations were performed to calculate the radial distribution functions of 1.0 molal solutions at 298 K which were used to compare the structure of the explicit and implicit water simulations. The implementation of dielectric saturation and a repulsive hydration potential results in an excellent description of the mean activity coefficient and is able to capture structural features of contact ion pairs and solvent separated ions.

Electrostatic Screening and Charge Correlation Effects in Micellization of Ionic Surfactants

We have used atomistic simulations to study the role of electrostatic screening and charge correlation effects in self-assembly processes of ionic surfactants into micelles. Specifically, we employed grand canonical Monte Carlo simulations to investigate the critical micelle concentration (cmc), aggregation number, and micellar shape in the presence of explicit sodium chloride (NaCl). The two systems investigated are cationic dodecyltrimethylammonium chloride (DTAC) and anionic sodium dodecyl sulfate (SDS) surfactants. Our explicit-salt results, obtained from a previously developed potential model with no further adjustment of its parameters, are in good agreement with experimental data for structural and thermodynamic micellar properties. We illustrate the importance of ion correlation effects by comparing these results with a Yukawa-type surfactant model that incorporates electrostatic screening implicitly. While the effect of salt on the cmc is well-reproduced even with the implicit Yukawa model, the aggregate size predictions deviate significantly from experimental observations at low salt concentrations. We attribute this discrepancy to the neglect of ion correlations in the implicit-salt model. At higher salt concentrations, we find reasonable agreement of the Yukawa model with experimental data. The crossover from low to high salt concentrations is reached when the electrostatic screening length becomes comparable to the headgroup size.

Implicit Solvent Models for Micellization of Ionic Surfactants

We propose a method for parametrization of implicit solvent models for the simulation of the self-assembly of ionic surfactants into micelles. The parametrization is carried out in two steps. The first step involves atomistic molecular dynamics simulations of headgroups and counterions with explicit solvent to determine structural properties. An implicit solvent model of the headgroup/counterion system is obtained by matching structural quantities between explicit solvent and implicit solvent systems. In the second step, we identify the solvophobic attractions between the tail beads. We determine the solvophobic parameters using grand canonical Monte Carlo simulations with histogram reweighting techniques. The matching objective for the identification of solvophobic attractions is the critical micelle concentration (cmc). We choose sodium dodecyl sulfate as the reference system. On the basis of hydrophobic parameters obtained from this particular model, we study specific ion effects (lithium and potassium instead of sodium) as well as the effect of cationic headgroups (dodecyltrimethylammonium bromide/chloride). Furthermore, the chain length dependence of micellization properties is investigated for sodium alkyl sulfate, with alkyl lengths between 6 and 14. All cases considered give results in broad agreement with experimental data, confirming the transferability of parameters and the generality of the approach.

Nanoscale carbon particles and the stability of lipid bilayers

The transfer of various nano-scale fullerenes into lipid bilayers has been studied using all-atom (AA) and coarse-grained (CG) molecular dynamics (MD) simulations. The free energy change, when C60, C180, and C540 fullerenes are transferred from water to the interior of a lipid [dioleoylphosphatidylcholine (DOPC)] bilayer, has been calculated. Upon entering the lipid bilayer, the largest (2.4 nm diameter) fullerene causes local distortions in the bilayer surface, which were previously observed in carbon nanotube simulations. These local distortions, however, do not lead to any free energy barriers to bilayer entry. The free energy profiles confirm spontaneous absorption of all three fullerenes. Qualitative agreement was observed when comparing fullerene partitioning in water/bilayer systems to water–hexane systems. In contrast to these nonspecific single fullerene properties, extensive CG-MD simulations of fullerene rich lipid bilayers reveal substantial impact of fullerene-size on the bilayer stability. While previous CG-MD simulations indicated that bilayer bound C60 aggregates have little effect on the bilayer structure, the present MD simulations indicate that C540 aggregation has drastic effects. Specifically, the observed destabilization likely has implications for understanding the cytotoxic mechanisms of nano-carbon particles upon uptake by biological cells.

Molecular design of responsive fluids: molecular dynamics studies of viscoelastic surfactant solutions

Understanding how macroscopic properties depend on intermolecular interactions for complex fluid systems is an enormous challenge in statistical mechanics. This issue is of particular importance for designing optimal industrial fluid formulations such as responsive oilfield fluids, based on viscoelastic surfactant solutions. We have carried out extensive molecular dynamics simulations, resolving the full chemical details in order to study how the structure of the lamellar phase of viscoelastic surfactant solutions depends on the head group (HG) chemistry of the surfactant. In particular, we consider anionic carboxylate and quaternary ammonium HGs with erucyl tails in aqueous solutions together with their sodium and chloride counterions at room temperature. We find a strong HG dependence of the lamellar structure as characterized by suitable pair correlation functions and density distributions. The depth of penetration of water into the bilayer membrane, the nature of counterion condensation on the HGs and even the order and correlation of the tails in the lamellae depend sensitively on the chemical details of the HG. We also determine the compressibility of the lamellar system as a first step to using atom-resolved molecular dynamics in order to link the molecular and macroscopic scales of length and time. The results give important insight into the links between molecular details and surfactant phase structure which is being exploited to develop more systematic procedures for the molecular design and formulation of industrial systems.

Microsurface Potential Measurements: Repulsive Forces between Polyelectrolyte Brushes in the Presence of Multivalent Counterions

We propose a new way to determine weak repulsive forces operative between colloidal particles by measuring the rate of slow coagulation. The rate of slow coagulation is directly related to the competition of the repulsion with thermal motion. Since the thermal forces are weak, measurements of the coagulation rate can lead to precise information on repulsive potentials having a magnitude of just a few kT. We demonstrate this novel way by studying colloidal spherical polyelectrolyte brush (SPB) particles in aqueous solution containing trivalent La3+ counterions. The particles consist of a monodisperse polystyrene core of 121 nm radius from which linear sodium poly(styrenesulfonate) (PSS) chains are densely grafted (contour length 48 nm). We determine the rate of coagulation by time-resolved simultaneous static and dynamic light scattering in the presence of LaCl3 (0.2 to 150 mM). Direct measurements of the repulsive force between macroscopic brush layers demonstrate that the potential is decaying exponentially with distance. This is in good agreement with a simple theoretical treatment that furthermore leads to the effective surface potential Psi0. The good agreement of data obtained by the novel microscopic method with direct macroscopic measurements underscores the general validity of our approach.

Triplet interactions in star polymer solutions

We analyze the effective triplet interactions between the centers of star polymers in a good solvent. Using an analytical short distance expansion inspired by scaling theory, we deduce that the triplet part of the three-star force is attractive but only 11% of the pairwise part even for a close approach of three star polymers. We have also performed extensive computer simulations for different arm numbers to extract the effective triplet force. The simulation data show good correspondence with the theoretical predictions. Our results justify the effective pair potential picture even beyond the star polymer overlap concentration.