Dynamics of counterions in dendrimer polyelectrolyte solutions

The Role of Oxidation Pattern and Water Content in the Spatial Arrangement and Dynamics of Oxidized Graphene-Based Aqueous Dispersions

In this work, we employ fully atomistic molecular dynamics simulations to elucidate the effects of the oxidation pattern and of the water content on the organization of graphene sheets in aqueous dispersions and on the dynamic properties of the different moieties at neutral pH conditions. Analysis of the results reveals the role of the oxidation motif (peripherally or fully oxidized flakes) in the tendency of the flakes to self-assemble and in the control of key structural characteristics, such as the interlayer distance between the sheets and the average size and the distribution of the formed aggregates. In certain cases, the results are compared to a pertinent experimental system, validating further the relevant computational models. Examination of the diffusional motion of the oxidized flakes shows that different degrees of spatial restriction are imposed upon the decrease in the water content and elucidates the conditions under which a motional arrest of the flakes takes place. At constant water content, the structural differences between the formed aggregates appear to additionally impart distinct diffusional characteristics of a water molecule. A detailed examination of the counterion dynamics describes their interaction with the oxidized flakes and their dependence on the water content and on the oxidation pattern, offering new insight into the expected electrical properties of the dispersions. The detailed information provided by this work will be particularly useful in applications such as molecular sieving, nanofiltration, and in cases where conductive membranes based on oxidized forms of graphene are used.

Explicit characterization of counterion dynamics around a flexible polyelectrolyte

The article presents a comprehensive study of counterion dynamics around a generic linear polyelectrolyte chain with the help of coarse-grained computer simulations. The ion-chain coupling is discussed in the form of binding time, mean-square displacement (MSD) relative to the chain, local ion transport coefficient, and spatiotemporal correlations in the effective charge. We have shown that a counterion exhibits subdiffusive behavior 〈δR^{2}〉∼t^{δ}, δ≈0.9 w.r.t. chain's center of mass. The MSD of ions perpendicularly outward from the chain segment exhibits a smaller subdiffusive exponent compared to the one along the chain backbone. The effective diffusivity of ion is the lowest in chain's close proximity, extending up to the length-scale of radius of gyration R_{g}. Beyond R_{g} at larger distances, they attain diffusivity of free ion with a smooth cross-over from the adsorbed regime to the free ion regime. We have shown that the effective diffusivity drastically decreases for the multivalent ions, while the crossover length scale remains the same. Conversely, with increasing salt concentration the coupling-length scale reduces, while the diffusivity remains unaltered. The effective diffusivity of adsorbed-ion reveals an exponential reduction with electrostatic interaction strength. We further corroborate this from the binding time of ions on the chain, which also grows exponentially with the coupling strength of the ion-polymer duo. Moreover, the binding time of ions exhibits a weak dependence with salt concentration for the monovalent salt, while for multivalent salts the binding time decreases dramatically with concentration. Our work also elucidates fluctuations in the effective charge per site, where it exhibits strong negative correlations at short length-scales.

A microscopic view of graphene-oxide/poly(acrylic acid) physical hydrogels: effects of polymer charge and graphene oxide loading

In this work we have examined in detail by means of fully atomistic molecular dynamics simulations, physical hydrogels formed by a polymer electrolyte, poly(acrylic acid), and graphene oxide, at two different charging states of the polymer and two different graphene oxide concentrations. It was found that variations of these parameters incurred drastic changes in general morphological characteristics of the composite materials, the degree of physical adsorption of polyelectrolyte chains onto the graphene oxide surface, the polymer dynamic response at local and global length scales, in the charge distributions around the components, and in the mobility of the counterions. All these microscopic features are expected to significantly affect macroscopic physical properties of the hydrogels, such as their mechanical responses and their electrical behaviors.

Counter-ion distribution around flexible polyelectrolytes having different molecular architecture

We explore the monovalent counter-ion distribution around flexible highly-charged polyelectrolytes with different molecular architectures (linear chains, stars, and unknotted and trefoil rings) using molecular dynamics simulations that include an explicit solvent that interacts with the polyelectrolyte. In particular, we find that the molecular topology influences the fraction of counter-ions transiently associating with the polyelectrolyte on a scale of the order of the chain segments, forming a "condensed" counter-ion interfacial layer. As with the hydrogen bonding of water to proteins and other polymers, the persistence time of these interfacial "bound" counter-ions is relatively short, O(1 ps), and we characterize the fluctuations in the number of the counter-ions populating the interfacial layer. We also find that the counter-ions are distributed in a non-uniform fashion on the polyelectrolyte backbone, forming dynamical clusters whose form and average size is sensitive to molecular architecture. In addition, we find that the residual bound counter-ions, not located in either the interfacial layer or the bulk solution, form a diffuse ionic cloud around the polyelectrolyte due to the uncompensated polyelectrolyte charge along the backbone. Generally charge valence strongly influences the extent of the diffuse counter-ion cloud, but in the case of monovalent counter-ions, we find that the size of the diffuse counter-ion cloud nearly coincides with the polyelectrolyte radius of gyration, independent of molecular topology.

Charge-dependent transport switching of single molecular ions in a weak polyelectrolyte multilayer

The tunable nature of weak polyelectrolyte multilayers makes them ideal candidates for drug loading and delivery, water filtration, and separations, yet the lateral transport of charged molecules in these systems remains largely unexplored at the single molecule level. We report the direct measurement of the charge-dependent, pH-tunable, multimodal interaction of single charged molecules with a weak polyelectrolyte multilayer thin film, a 10 bilayer film of poly(acrylic acid) and poly(allylamine hydrochloride) PAA/PAH. Using fluorescence microscopy and single-molecule tracking, two modes of interaction were detected: (1) adsorption, characterized by the molecule remaining immobilized in a subresolution region and (2) diffusion trajectories characteristic of hopping (D ∼ 10(-9) cm(2)/s). Radius of gyration evolution analysis and comparison with simulated trajectories confirmed the coexistence of the two transport modes in the same single molecule trajectories. A mechanistic explanation for the probe and condition mediated dynamics is proposed based on a combination of electrostatics and a reversible, pH-induced alteration of the nanoscopic structure of the film. Our results are in good agreement with ensemble studies conducted on similar films, confirm a previously-unobserved hopping mechanism for charged molecules in polyelectrolyte multilayers, and demonstrate that single molecule spectroscopy can offer mechanistic insight into the role of electrostatics and nanoscale tunability of transport in weak polyelectrolyte multilayers.

Evidence and characterization of dynamic heterogeneity in binary mixtures of phosphoric acid and benzimidazole

We report here anomalous diffusions of components in mixtures of monomer of polybenzimidazole, i.e., 2-phenyl-1H,1'H-5,5'-bibenzo[d]imidazole (BI) and phosphoric acid (PA) from molecular dynamics simulations. We have observed initial drop and further increase in self-diffusion constant for both monomer molecule (BI) and PA with gradual increase in PA concentration. The origin of such anomalous diffusion is identified in this work, which happens to be the presence of dynamic heterogeneity in each component of the binary mixture. We characterized microscopic picture of dynamical heterogeneity by finding correlation between dynamical heterogeneity and structural arrangement among the components of the binary system. Different types of H-bonding arrangements in the BI-PA systems at different concentration of PA are observed. The stability of the H-bonded network consisting of different types of H-bonds between BI and PA in the system has been studied by calculating the lifetime of various H-bonds. The results indicate that there are fast and slow moving PA molecules in the mixtures because of coexistence of different types of hydrogen bonds among the components of the mixture.

Theoretical and computational studies of dendrimers as delivery vectors

It is a great challenge for nanomedicine to develop novel dendrimers with maximum therapeutic potential and minimum side-effects for drug and gene delivery. As delivery vectors, dendrimers must overcome lots of barriers before delivering the bio-agents to the target in the cell. Extensive experimental investigations have been carried out to elucidate the physical and chemical properties of dendrimers and explore their behaviors when interacting with biomolecules, such as gene materials, proteins, and lipid membranes. As a supplement of the experimental techniques, it has been proved that computer simulations could facilitate the progress in understanding the delivery process of bioactive molecules. The structures of dendrimers in dilute solutions have been intensively investigated by monomer-resolved simulations, coarse-grained simulations, and atom-resolved simulations. Atomistic simulations have manifested that the hydrophobic interactions, hydrogen-bond interactions, and electrostatic attraction play critical roles in the formation of dendrimer-drug complexes. Multiscale simulations and statistical field theories have uncovered some physical mechanisms involved in the dendrimer-based gene delivery systems. This review will focus on the current status and perspective of theoretical and computational contributions in this field in recent years. (275 references).

Modeling the formation of ordered nano-assemblies comprised by dendrimers and linear polyelectrolytes: the role of Coulombic interactions

Models of mixtures of peripherally charged dendrimers with oppositely charged linear polyelectrolytes in the presence of explicit solvent are studied by means of molecular dynamics simulations. Under the influence of varying strength of electrostatic interactions, these systems appear to form dynamically arrested film-like interconnected structures in the polymer-rich phase. Acting like a pseudo-thermodynamic inverse temperature, the increase of the strength of the Coulombic interactions drive the polymeric constituents of the mixture to a gradual dynamic freezing-in. The timescale of the average density fluctuations of the formed complexes initially increases in the weak electrostatic regime reaching a finite limit as the strength of electrostatic interactions grow. Although the models are overall electrically neutral, during this process the dendrimer/linear complexes develop a polar character with an excess charge mainly close to the periphery of the dendrimers. The morphological characteristics of the resulted pattern are found to depend on the size of the polymer chains on account of the distinct conformational features assumed by the complexed linear polyelectrolytes of different length. In addition, the length of the polymer chain appears to affect the dynamics of the counterions, thus affecting the ionic transport properties of the system. It appears, therefore, that the strength of electrostatic interactions together with the length of the linear polyelectrolytes are parameters to which these systems are particularly responsive, offering thus the possibility for a better control of the resulted structure and the electric properties of these soft-colloidal systems.