Molecular Dynamics Simulations of Rhodamine B Zwitterion Diffusion in Polyelectrolyte Solutions

Polyelectrolytes continue to find wide interest and application in science and engineering, including areas such as water purification, drug delivery, and multilayer thin films. We have been interested in the dynamics of small molecules in a variety of polyelectrolyte (PE) environments; in this paper, we report simulations and analysis of the small dye molecule rhodamine B (RB) in several very simple polyelectrolyte solutions. Translational diffusion of the RB zwitterion has been measured in fully atomistic, 2 μs long molecular dynamics simulations in four different polyelectrolyte solutions. Two solutions contain the common polyanion sodium poly(styrene sulfonate) (PSS), one with a 30-mer chain and the other with 10 trimers. The other two solutions contain the common polycation poly(allyldimethylammonium) chloride (PDDA), one with two 15-mers and the other with 10 trimers. RB diffusion was also simulated in several polymer-free solutions to verify its known experimental value for the translational diffusion coefficient, D_RB, of 4.7 × 10^-6 cm²/s at 300 K. RB diffusion was slowed in all four simulated PE solutions, but to varying degrees. D_RB values of 3.07 × 10^-6 and 3.22 × 10^-6 cm²/s were found in PSS 30-mer and PSS trimer solutions, respectively, whereas PDDA 15-mer and trimer solutions yielded values of 2.19 × 10^-6 and 3.34 × 10^-6 cm²/s. Significant associations between RB and the PEs were analyzed and interpreted via a two-state diffusion model (bound and free diffusion) that describes the data well. Crowder size effects and anomalous diffusion were also analyzed. Finally, RB translation along the polyelectrolytes during association was characterized.

Molecular modeling of interfacial layer-by-layer assembly towards functionalized capsule materials

Encapsulated nanomaterials, such as polymer-coated nanoemulsions, have highly tunable properties leading to versatile applications. A current lack of understanding of the fundamentals governing the choice of "capsule" materials (polyelectrolyte + surfactant) and its ensuing performance effectively precludes their widespread use. Computational methods can start to redress this by discovering molecule-scale attributes that significantly control the design of capsule materials tuned to fit desired properties. We use molecular dynamics (MD) to carry out the layer-by-layer (LbL) assembly of six unique polyelectrolyte bilayer systems at a surfactant-mediated interface, modeling early-stage capsule synthesis. Monolayer thickness is related to layer density and polyelectrolyte/surfactant interaction energy through polyelectrolyte molecular weight and radius of gyration, respectively, yielding a simple relationship between absorption kinetics and layer structure. For the second monolayer, faster absorption kinetics are observed for pairings of polyelectrolytes with similarly sized functional groups. Surfactants with a more delocalized charge on the head-group catalyze the build-up of ions at the interface, resulting in faster absorption kinetics and greater confinement of the encapsulated material but leading to thicker, less uniform bilayers. These relationships between capsule building block molecules and nanomaterial capsule properties provide a foundation for property prediction and rational design of optimized multi-functional capsule materials.

Ion-Specific and Solvent Effects on PDADMA–PSS Complexation and Multilayer Formation

Among various parameters that influence the formation of polyelectrolyte complexes and multilayers, special emphasis should be placed on ion-specific and solvent effects. In our study, we systematically examined the above-mentioned effects on poly(diallyldimethylammonium chloride) (PDADMACl)-sodium poly(4-styrenesulfonate) (NaPSS) complexation in solution and at the surface by means of dynamic light scattering, ellipsometry and atomic force microscopy measurements. As solvents, we used water and water/ethanol mixture. The obtained results confirm the importance of ion-specific and solvent effects on complexes prepared in solution, as well as on multilayers built up on a silica surface. The experiments in mixed solvent solution showed that at a higher ethanol mole fraction, the decrease in monomer titrant to titrand ratio, at which the increase in the size of complexes is observed, takes place. The difference between chloride and bromide ions was more pronounced at a higher mole fraction of ethanol and in the case of positive complex formation, suggesting that the larger amount of bromide ions could be condensed to the polycation chain. These findings are in accordance with the results we obtained for polyelectrolyte multilayers and could be helpful for designing polyelectrolyte multilayers with tuned properties needed for various applications, primarily in the field of biomedicine.

PDADMAC/PSS Oligoelectrolyte Multilayers: Internal Structure and Hydration Properties at Early Growth Stages from Atomistic Simulations

We analyze the internal structure and hydration properties of poly(diallyl dimethyl ammonium chloride)/poly(styrene sulfonate sodium salt) oligoelectrolyte multilayers at early stages of their layer-by-layer growth process. Our study is based on large-scale molecular dynamics simulations with atomistic resolution that we presented recently [Sánchez et al., Soft Matter2019, 15, 9437], in which we produced the first four deposition cycles of a multilayer obtained by alternate exposure of a flat silica substrate to aqueous electrolyte solutions of such polymers at 0.1M of NaCl. In contrast to any previous work, here we perform a local structural analysis that allows us to determine the dependence of the multilayer properties on the distance to the substrate. We prove that the large accumulation of water and ions next to the substrate observed in previous overall measurements actually decreases the degree of intrinsic charge compensation, but this remains as the main mechanism within the interface region. We show that the range of influence of the substrate reaches approximately 3 nm, whereas the structure of the outer region is rather independent from the position. This detailed characterization is essential for the development of accurate mesoscale models able to reach length and time scales of technological interest.

Atomistic simulation of PDADMAC/PSS oligoelectrolyte multilayers: overall comparison of tri- and tetra-layer systems

By employing large-scale molecular dynamics simulations of atomistically resolved oligoelectrolytes in aqueous solutions, we study in detail the first four layer-by-layer deposition cycles of an oligoelectrolyte multilayer made of poly(diallyl dimethyl ammonium chloride)/poly(styrene sulfonate sodium salt) (PDADMAC/PSS). The multilayers are grown on a silica substrate in 0.1 M NaCl electrolyte solutions and the swollen structures are then subsequently exposed to varying added salt concentration. We investigated the microscopic properties of the films, analyzing in detail the differences between three- and four-layer systems. Our simulations provide insights into the early stages of growth of a multilayer, which are particularly challenging for experimental observations. We found rather strong complexation of the oligoelectrolytes, with fuzzy layering of the film structure. The main charge compensation mechanism is for all cases intrinsic, whereas extrinsic compensation is relatively enhanced for the layer of the last deposition cycle. In addition, we quantified other fundamental observables of these systems, such as the film thickness, water uptake, and overcharge fractions for each deposition layer.

Coarse-grained simulations of polyelectrolyte complexes: MARTINI models for poly(styrene sulfonate) and poly(diallyldimethylammonium)

We present simulations of aqueous polyelectrolyte complexes with new MARTINI models for the charged polymers poly(styrene sulfonate) and poly(diallyldimethylammonium). Our coarse-grained polyelectrolyte models allow us to study large length and long time scales with regard to chemical details and thermodynamic properties. The results are compared to the outcomes of previous atomistic molecular dynamics simulations and verify that electrostatic properties are reproduced by our MARTINI coarse-grained approach with reasonable accuracy. Structural similarity between the atomistic and the coarse-grained results is indicated by a comparison between the pair radial distribution functions and the cumulative number of surrounding particles. Our coarse-grained models are able to quantitatively reproduce previous findings like the correct charge compensation mechanism and a reduced dielectric constant of water. These results can be interpreted as the underlying reason for the stability of polyelectrolyte multilayers and complexes and validate the robustness of the proposed models.

Properties of water in the interfacial region of a polyelectrolyte bilayer adsorbed onto a substrate studied by computer simulations

We study the static and dynamic properties of water near a poly(styrene sulfonate)/poly(diallyldimethylammonium) (PSS/PDADMA) bilayer adsorbed onto a substrate by atomistic molecular dynamics simulations. Qualitative changes in the dynamics of water in the proximity of the adsorbed bilayer are observed - such as in the lateral diffusion, residence time and hydrogen-bonding lifetime - as compared with water in the presence of the bare substrate. Static properties of water are similarly influenced, and a high polarization of water molecules is found to be present surprisingly far from the adsorbed bilayer.

Aqueous NaCl and CsCl solutions confined in crystalline slit-shaped silica nanopores of varying degree of protonation

All-atom molecular dynamics simulations were conducted to study the dynamics of aqueous electrolyte solutions confined in slit-shaped silica nanopores of various degrees of protonation. Five degrees of protonation were prepared by randomly removing surface hydrogen atoms from fully protonated crystalline silica surfaces. Aqueous electrolyte solutions containing NaCl or CsCl salt were simulated at ambient conditions. In all cases, the ionic concentration was 1 M. The results were quantified in terms of atomic density distributions within the pores, and the self-diffusion coefficient along the direction parallel to the pore surface. We found evidence for ion-specific properties that depend on ion-surface, water-ion, and only in some cases ion-ion correlations. The degree of protonation strongly affects the structure, distribution, and the dynamic behavior of confined water and electrolytes. Cl(-) ions adsorb on the surface at large degrees of protonation, and their behavior does not depend significantly on the cation type (either Na(+) or Cs(+) ions are present in the systems considered). The cations show significant ion-specific behavior. Na(+) ions occupy different positions within the pore as the degree of protonation changes, while Cs(+) ions mainly remain near the pore center at all conditions considered. For a given degree of protonation, the planar self-diffusion coefficient of Cs(+) is always greater than that of Na(+) ions. The results are useful for better understanding transport under confinement, including brine behavior in the subsurface, with important applications such as environmental remediation.