Simulation of surfactant adsorption at liquid-liquid interface: What we may expect from soft-core models?

The present work attempts to systematically explore the surfactant sorption at liquid-liquid interfaces with coarse-grained models targeting thermodynamic properties of reference liquid solutions. We employ dissipative particle dynamics with soft-core forcefield tested against experimental data on micellization of surfactants in water, and the previous results are reproduced in this work. We consider three different nonionic surfactants: hexaethylene glycol monododecyl ether (C12E6), 2-[4-(2,4,4-trimethylpentan-2-yl)phenoxy]ethanol) known as Triton X-100 (TX-100), and two alkyl glucoside surfactants (CnG1) with n-alkane tail fragments and a saccharide hydrophilic head at decane-water and toluene-water interfaces. For TX-100, we composed a model based on the literature forcefield and found good agreement with the experimental critical micelle concentrations (CMCs). The head-head interactions are of different origins for different surfactant groups: entropic repulsion between ethylene oxide chains of C12E6 and TX-100, and more chemically specific and complex interactions between the maltose heads of alkyl glucosides. We interpret our results with the Redlich-Peterson equation of monolayer adsorption in order to relate the adsorption to the bulk concentration of the surfactant and the interfacial tension. The densities of the adsorbed monolayer at CMC mostly agree with the experimental data, and a reasonable agreement was obtained for the interfacial tension at CMC. At the same time, we found significant discrepancies between the simulated and experimental adsorption isotherms. We explain them by the oversimplified forcefield: when the parameters are fitted to the free energies of bulk solutions, they may not correctly reproduce the interfacial free energies.

Adhesion and Separation of Nanoparticles on Polymer-Grafted Porous Substrates.

This work explores interactions of functionalized nanoparticles (NP) with polymer brushes (PB) in a binary mixture of good and poor solvents. NP-PB systems are used in multiple applications, and we are particularly interested in the problem of chromatographic separation of NPs on polymer-grafted porous columns. This process involves NP flow through the pore channels with walls covered by PBs. NP-PB adhesion is governed by adsorption of polymer chains to NP surface and entropic repulsion caused by the polymer chain confinement between NP and the channel wall. Both factors depend on the solvent composition, variation of which causes contraction or expansion of PB. Using dissipative particle dynamics simulations in conjunction with the ghost tweezers free energy calculation technique, we examine the free energy landscapes of functionalized NPs within PB-grafted channels depending on the solvent composition at different PB grafting densities and polymer-solvent affinities. The free energy landscape determines the probability of NP location at a given distance to the surface, positions of equilibrium adhesion states, and the Henry constant that characterizes adsorption equilibrium and NP partitioning between the stationary phase of PB and mobile phase of flowing solvent. We analyze NP transport through a polymer-grafted channel and calculate the mean velocity and retention time of NP depending on the NP size and solvent composition. We find that, with the increase of the bad (poor) solvent fraction and respective PB contraction, NP separation exhibits a transition from the hydrodynamic size exclusion regime with larger NPs having shorter retention time to the adsorption regime with smaller NPs having shorter retention time. The observed reversal of the sequence of elution is reminiscent of the critical condition in polymer chromatography at which the retention time is molecular weight independent. This finding suggests the possibility of the existence of an analogous special regime in nanoparticle chromatography at which NPs with like surface properties elute together regardless of their size. The latter has important practical implications: NPs can be separated by surface chemistry rather than by their size employing the gradient mode of elution with controlled variation of solvent composition.

Coarse-grained model of nanoscale segregation, water diffusion, and proton transport in Nafion membranes.

We present a coarse-grained model of the acid form of Nafion membrane that explicitly includes proton transport. This model is based on a soft-core bead representation of the polymer implemented into the dissipative particle dynamics (DPD) simulation framework. The proton is introduced as a separate charged bead that forms dissociable Morse bonds with water beads. Morse bond formation and breakup artificially mimics the Grotthuss hopping mechanism of proton transport. The proposed DPD model is parameterized to account for the specifics of the conformations and flexibility of the Nafion backbone and sidechains; it treats electrostatic interactions in the smeared charge approximation. The simulation results qualitatively, and in many respects quantitatively, predict the specifics of nanoscale segregation in the hydrated Nafion membrane into hydrophobic and hydrophilic subphases, water diffusion, and proton mobility. As the hydration level increases, the hydrophilic subphase exhibits a percolation transition from a collection of isolated water clusters to a 3D network of pores filled with water embedded in the hydrophobic matrix. The segregated morphology is characterized in terms of the pore size distribution with the average size growing with hydration from ∼1 to ∼4 nm. Comparison of the predicted water diffusivity with the experimental data taken from different sources shows good agreement at high and moderate hydration and substantial deviation at low hydration, around and below the percolation threshold. This discrepancy is attributed to the dynamic percolation effects of formation and rupture of merging bridges between the water clusters, which become progressively important at low hydration, when the coarse-grained model is unable to mimic the fine structure of water network that includes singe molecule bridges. Selected simulations of water diffusion are performed for the alkali metal substituted membrane which demonstrate the effects of the counter-ions on membrane self-assembly and transport. The hydration dependence of the proton diffusivity reproduces semi-qualitatively the trend of the diverse experimental data, showing a sharp decrease around the percolation threshold. Overall, the proposed model opens up an opportunity to study self-assembly and water and proton transport in polyelectrolytes using computationally efficient DPD simulations, and, with further refinement, it may become a practical tool for theory informed design and optimization of perm-selective and ion-conducting membranes with improved properties.

Adhesion of Phospholipid Bilayers to Hydroxylated Silica: Existence of Nanometer-Thick Water Interlayers.

Lipid bilayers attached to solid surfaces play an important role in bioinspired materials and devices and serve as model systems for studies of interactions of cell membranes with particles and biomolecules. Despite active experimental and theoretical studies, the interactions of lipid membranes with solid substrates are still poorly understood. In this work, we explore, using atomistic molecular dynamics simulations, the equilibrium and stability of a phospholipid 1,2-dimyristoyl-sn-glycero-3-phosphocholine membrane supported on hydroxylated amorphous silica. We reveal two distinct types of thermodynamically stable states, characterized by different widths of the water layer between the membrane and the substrate. In α-states, the membrane is closely attached with the lipid head groups interacting directly with surface hydroxyls; however, because of the molecular level roughness of the amorphous silica surface, there exists an inhomogeneous water layer trapped between the substrate and the membrane. In ß-states, the membrane is separated from the silica surface by a water film of ∼2.5 nm in thickness. The thermodynamic equilibrium is quantified in terms of the disjoining pressure isotherm as a function of membrane-substrate separation, which has a double sigmoidal shape with two minima and one maximum, which correspond to the limits of stability of α- and ß-states. The thermodynamic properties and bilayer structure are compared with experimental findings and simulation results for relevant systems.

Coarse-grained model of water diffusion and proton conductivity in hydrated polyelectrolyte membrane.

Using dissipative particle dynamics (DPD), we simulate nanoscale segregation, water diffusion, and proton conductivity in hydrated sulfonated polystyrene (sPS). We employ a novel model [Lee et al. J. Chem. Theory Comput. 11(9), 4395-4403 (2015)] that incorporates protonation/deprotonation equilibria into DPD simulations. The polymer and water are modeled by coarse-grained beads interacting via short-range soft repulsion and smeared charge electrostatic potentials. The proton is introduced as a separate charged bead that forms dissociable Morse bonds with the base beads representing water and sulfonate anions. Morse bond formation and breakup artificially mimics the Grotthuss mechanism of proton hopping between the bases. The DPD model is parameterized by matching the proton mobility in bulk water, dissociation constant of benzenesulfonic acid, and liquid-liquid equilibrium of water-ethylbenzene solutions. The DPD simulations semi-quantitatively predict nanoscale segregation in the hydrated sPS into hydrophobic and hydrophilic subphases, water self-diffusion, and proton mobility. As the hydration level increases, the hydrophilic subphase exhibits a percolation transition from isolated water clusters to a 3D network. The analysis of hydrophilic subphase connectivity and water diffusion demonstrates the importance of the dynamic percolation effect of formation and breakup of temporary junctions between water clusters. The proposed DPD model qualitatively predicts the ratio of proton to water self-diffusion and its dependence on the hydration level that is in reasonable agreement with experiments.

In Situ Growth and Characterization of Metal Oxide Nanoparticles within Polyelectrolyte Membranes.

This study describes a novel approach for the in situ synthesis of metal oxide-polyelectrolyte nanocomposites formed via impregnation of hydrated polyelectrolyte films with binary water/alcohol solutions of metal salts and consecutive reactions that convert metal cations into oxide nanoparticles embedded within the polymer matrix. The method is demonstrated drawing on the example of Nafion membranes and a variety of metal oxides with an emphasis placed on zinc oxide. The in situ formation of nanoparticles is controlled by changing the solvent composition and conditions of synthesis that for the first time allows one to tailor not only the size, but also the nanoparticle shape, giving a preference to growth of a particular crystal facet. The high-resolution TEM, SEM/EDX, UV-vis and XRD studies confirmed the homogeneous distribution of crystalline nanoparticles of circa 4 nm and their aggregates of 10-20 nm. The produced nanocomposite films are flexible, mechanically robust and have a potential to be employed in sensing, optoelectronics and catalysis.

Adhesion of nanoparticles to polymer brushes studied with the ghost tweezers method.

Mechanisms of interactions between nanoparticles (NPs) and polymer brushes (PBs) are explored using dissipative particle dynamics simulations and an original "ghost tweezers" method that emulates lab experiments performed with optical or magnetic tweezers. The ghost tweezers method is employed to calculate the free energy of adhesion. Ghost tweezers represents a virtual harmonic potential, which tethers NP with a spring to a given anchor point. The average spring force represents the effective force of NP-PB interaction as a function of the NP coordinate. The free energy landscape of NP-PB interactions is calculated as the mechanical work needed to transfer NP from the solvent bulk to a particular distance from the substrate surface. With this technique, we explore the adhesion of bare and ligand-functionalized spherical NPs to polyisoprene natural rubber brush in acetone-benzene binary solvent. We examine two basic mechanisms of NP-PB interactions, NP adhesion at PB exterior and NP immersion into PB, which are governed by interplay between entropic repulsive forces and enthalpic attractive forces caused by polymer adsorption at the NP surface and ligand adsorption at the substrate. The relative free energies of the equilibrium adhesion states and the potential barriers separating these states are calculated at varying grafting density, NP size, and solvent composition.

Morphological transformations in polymer brushes in binary mixtures: DPD study.

Morphological transformations in polymer brushes in a binary mixture of good and bad solvents are studied using dissipative particle dynamics simulations drawing on a characteristic example of polyisoprene natural rubber in an acetone-benzene mixture. A coarse-grained DPD model of this system is built based on the experimental data in the literature. We focus on the transformation of dense, collapsed brush in bad solvent (acetone) to expanded brush solvated in good solvent (benzene) as the concentration of benzene increases. Compared to a sharp globule-to-coil transition observed in individual tethered chains, the collapsed-to-expanded transformation in brushes is found to be gradual without a prominent transition point. The transformation becomes more leveled as the brush density increases. At low densities, the collapsed brush is highly inhomogeneous and patterned into bunches composed of neighboring chains due to favorable polymer-polymer interaction. At high densities, the brush is expanded even in bad solvent due to steric restrictions. In addition, we considered a model system similar to the PINR-acetone-benzene system, but with the interactions between the solvent components worsened to the limit of miscibility. Enhanced contrast between good and bad solvents facilitates absorption of the good solvent by the brush, shifting the collapsed-to-expanded transformation to lower concentrations of good solvent. This effect is especially pronounced for higher brush densities.

Translocation dynamics of freely jointed Lennard-Jones chains into adsorbing pores.

Polymer translocation into adsorbing nanopores is studied by using the Fokker-Planck equation of chain diffusion along the energy landscape calculated with Monte Carlo simulations using the incremental gauge cell method. The free energy profile of a translocating chain was found by combining two independent sub-chains, one free but tethered to a hard wall, and the other tethered inside an adsorbing pore. Translocation dynamics were revealed by application of the Fokker-Planck equation for normal diffusion. Adsorption of polymer chains into nanopores involves a competition of attractive adsorption and repulsive steric hindrance contributions to the free energy. Translocation times fell into two regimes depending on the strength of the adsorbing pore. In addition, we found a non-monotonic dependence of translocation times with increasing adsorption strength, with sharp peak associated with local free energy minima along the translocation coordinate.

Molecular modeling of organophosphorous agents and their aqueous solutions.

Using molecular dynamics simulations, we modeled solvation and diffusion in aqueous solutions of organophosphorous compounds, including nerve G-agents sarin and soman (methylphosphonofluoridates) and their common simulants DMMP (dimethyl methylphosphonate) and DIFP (diisopropyl fluorophosphate). The aqueous solutions of the organophosphorous compounds were found to display complex molecular scale structures and dynamic properties due to competing interactions between strongly hydrophobic and hydrophilic groups. The mixing of agents with water was proved to be exothermic with negative excess mixing volume, indicating a strongly hydrophilic solvation. This effect was confirmed in a specially performed experiment. We discuss to what extent DMMP and DIFP are suitable simulants for G-agents in experimental studies, as far as their interactions with water are concerned. We also focus on the relevance of the structural features and mobilities of agents in water to their interactions with permselective polyelectrolyte membranes that may be employed as protective barriers against chemical warfare agents.

Calculation of chemical potentials of chain molecules by the incremental gauge cell method.

The gauge cell Monte Carlo method is extended to calculations of the incremental chemical potentials and free energies of linear chain molecules. The method was applied to chains of Lennard-Jones beads with stiff harmonic bonds up to 500 monomers in length. We show that the suggested method quantitatively reproduces the modified Widom particle insertion method of Kumar et al. [S. K. Kumar, I. Szleifer, and A. Z. Panagiotopoulos, Phys. Rev. Lett. 66(22), 2935 (1991)], and is by an order of magnitude more efficient for long chains in terms of the computational time required for the same accuracy of chemical potential calculations. The chain increment ansatz, which suggests that the incremental chemical potential is independent of the chain length, was tested at different temperatures. We confirmed that the ansatz holds only for coils above the θ temperature. Special attention is paid to the effects of the magnitude of adsorption potential and temperature on the behavior of single chains in confinements that are comparable in size with the free chain radius of gyration. At sufficiently low temperatures, the dependence of the incremental chemical potential on the chain length in wetting pores is superficially similar to a capillary condensation isotherm, reflecting monolayer formation following by pore volume filling, as the chain length increases. We find that the incremental gauge cell method is an accurate and efficient technique for calculations of the free energies of chain molecules in bulk systems and nanoconfinements alike. The suggested method may find practical applications, such as modeling polymer partitioning on porous substrates and dynamics of chain translocation into nanopores.

Self-assembly in block polyelectrolytes.

The self-consistent field theory (SCFT) complemented with the Poisson-Boltzmann equation is employed to explore self-assembly of polyelectrolyte copolymers composed of charged blocks A and neutral blocks B. We have extended SCFT to dissociating triblock copolymers and demonstrated our approach on three characteristic examples: (1) diblock copolymer (AB) melt, (2) symmetric triblock copolymer (ABA) melt, (3) triblock copolymer (ABA) solution with added electrolyte. For copolymer melts, we varied the composition (that is, the total fraction of A-segments in the system) and the charge density on A blocks and calculated the phase diagram that contains ordered mesophases of lamellar, gyroid, hexagonal, and bcc symmetries, as well as the uniform disordered phase. The phase diagram of charged block copolymer melts in the charge density--system composition coordinates is similar to the classical phase diagram of neutral block copolymer melts, where the composition and the Flory mismatch interaction parameter χ(AB) are used as variables. We found that the transitions between the polyelectrolyte mesophases with the increase of charge density occur in the same sequence, from lamellar to gyroid to hexagonal to bcc to disordered morphologies, as the mesophase transitions for neutral diblocks with the decrease of χ(AB). In a certain range of compositions, the phase diagram for charged triblock copolymers exhibits unexpected features, allowing for transitions from hexagonal to gyroid to lamellar mesophases as the charge density increases. Triblock polyelectrolyte solutions were studied by varying the charge density and solvent concentration at a fixed copolymer composition. Transitions from lamellar to gyroid and gyroid to hexagonal morphologies were observed at lower polymer concentrations than the respective transitions in the similar neutral copolymer, indicating a substantial influence of the charge density on phase behavior.

Interactions of Crosslinked Polyacrylic Acid Polyelectrolyte Gels with Nonionic and Ionic Surfactants.

The morphology and stability of surfactant-loaded polyelectrolyte gels are of great interest for a variety of personal care, cosmetic, and pharmaceutical products. However, the mechanisms of surfactant interactions with gel-forming polymers are poorly understood and experimentally challenging. The aim of this work is to explore in silico the specifics of surfactant absorption within polyelectrolyte gels drawing on the examples of typical non-ionic octaethylene glycol monooctyl ether (C8E8) and anionic sodium dodecyl sulfate (SDS) surfactants and polyacrylic acid modified with hydrophobic sidechains mimicking the practically important Carbopol polymer. Using the systematically parameterized coarse-grained dissipative particle dynamics models, we generate and characterize the equilibrium conformations and swelling of the polymer films in aqueous solutions with the surfactant concentrations varied up to the critical micelle concentration (cmc). We discover the striking difference in interactions of Carbopol-like polymers with nonionic and ionic surfactants under mildly acidic conditions. The sorption of C8E8 within the polymer film is found substantial. As the surfactant concentration increases, the polymer film swells and, close to cmc, becomes unstable due to the formation and growth of water pockets filled with surfactant micelles. Sorption of SDS at the same bulk concentrations is found much lower, with only about 1% of surfactant mass fraction achieved at cmc. As the SDS concentration increases further, a lamellae structure is formed within the film, which remains stable. Reduced swelling and higher stability indicate better prospects of using SDS-type surfactants with Carbopol-based gels in formulations for detergents and personal care products.

Cavitation in metastable liquid nitrogen confined to nanoscale pores.

We studied cavitation in metastable fluids drawing on the example of liquid nitrogen confined to spheroidal pores of specially prepared well-characterized mesoporous silica materials with mean pore diameters ranging from approximately 6 to approximately 35 nm. Cavitation was monitored in the process of evaporation/desorption from fully saturated samples with gradually decreasing vapor pressure at the isothermal conditions. The onset of cavitation was displayed by a sharp step on the desorption isotherm. We found that the vapor pressure at the onset of cavitation depended on the pore size for the samples with pores smaller than approximately 11 nm and remained practically unchanged for the samples with larger pores. We suggest that the observed independence of the cavitation pressure on the size of confinement indicates that the conditions of bubble nucleation in pores larger than approximately 11 nm approach the nucleation conditions in the bulk metastable liquid. To test this hypothesis and to evaluate the nucleation barriers, we performed grand canonical and gauge cell Monte Carlo simulations of nitrogen adsorption and desorption in spherical silica pores ranging from 5.5 to 10 nm in diameter. Simulated and experimental adsorption isotherms were in good agreement. Exploiting the correlation between the experimental cavitation pressure and the simulated nucleation barrier, we found that the nucleation barrier increased almost linearly from approximately 40 to approximately 70 k(B)T in the range of pores from approximately 6 to approximately 11 nm, and varied in diapason of 70-75 k(B)T in larger pores, up to 35 nm. We constructed the dependence of the nucleation barrier on the vapor pressure, which asymptotically approaches the predictions of the classical nucleation theory for the metastable bulk liquid at larger relative pressures (>0.6). Our findings suggest that there is a limit to the influence of the confinement on the onset of cavitation, and thus, cavitation of nanoconfined fluids may be employed to explore cavitation in macroscopic systems.

Adhesion, intake, and release of nanoparticles by lipid bilayers.

Understanding the interactions between nanoparticles (NP) and lipid bilayers (LB), which constitute the foundations of cell membranes, is important for emerging biomedical technologies, as well as for assessing health threats related to nanoparticle commercialization. Applying dissipative particle dynamic simulations, we explore adhesion, intake, and release of hydrophobic nanoparticles by DMPC bilayers. To replicate experimental conditions, we develop a novel simulation setup for modeling membranes at isotension conditions. NP-LB interactions are quantified by the free energy landscape calculated by the ghost tweezers method. NPs are studied z of diameter 2 nm (comparable with the LB hydrophobic core), 4 nm (comparable with the LB thickness) and 8 nm (exceeding the LB thickness). NPs are pre-covered by an adsorbed lipid monolayer. It is shown that NP translocation across LB includes (1) NP intake into the hydrophobic core via merging of the monolayer adsorbed on NP with the outer leaflet of bilayer (2) NP release via formation and rupture of a lipid junction connecting NP and LB. Both stages are associated with free energy barriers. The barrier for the intake stage increases with the NP size and becomes prohibitively high for 8 nm NP. The barriers for the release stage are significantly higher which implies that the release stage controls the translocation rate and dynamics. The release energy barrier of 4 nm NP is found smaller than those for 2 and 8 nm NPs which implies the existence of the optimal NP size for unforced trans-membrane transport. Based on the calculated free energy landscape, the dynamics of unforced transport of NP across LB is evaluated using the Fokker-Planck equation, which mimics NP diffusion along the free energy landscape with multiple attempts to reach the barrier. We found that the number of attempts required for successful translocation scales exponentially with the energy barrier.

Stability of Lipid Coatings on Nanoparticle-Decorated Surfaces.

Lipid membranes supported on solid surfaces and nanoparticles find multiple applications in industrial and biomedical technologies. Here, we explore in silico the mechanisms of the interactions of lipid membranes with nanostructured surfaces with deposited nanoparticles and explain the characteristic particle size dependence of the uniformity and stability of lipid coatings observed in vitro. Simulations are performed to demonstrate the specifics of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) lipid membrane adhesion to hydrophilic and hydrophobic nanoparticles ranging in size from 1.5 to 40 nm using an original coarse-grained molecular dynamics model with implicit solvent and large simulation boxes (scales up to 280 × 154 × 69 nm3). We find that one of the major factors that affects the uniformity and stability of lipid coatings is the disjoining pressure in the water hydration layer formed between the lipid membrane and hydrophilic solid surface. This effect is accounted for by introducing a special long-range lipid-solid interaction potential that mimics the effects of the disjoining pressure in thin water layers. Our simulations reveal the physical mechanisms of interactions of lipid bilayers with solid surfaces that are responsible for the experimentally observed nonmonotonic particle size dependence of the uniformity and stability of lipid coatings: particles smaller than the hydration layer thickness (<2-3 nm) or larger than ∼20 nm are partially or fully enfolded by a lipid bilayer, whereas particles of the intermediate size (5-20 nm) cause membrane perforation and pore formation. In contrast, hydrophobic nanoparticles, which repel the hydration layer, tend to be encapsulated within the hydrophobic interior of the membrane and coated by the lipid monolayer. The proposed model can be further extended and applied to a wide class of systems comprising nanoparticles and nanostructured substrates interacting with lipid and surfactant bilayers and monolayers.

Multicomponent gauge cell method.

The gauge cell Monte Carlo method [Neimark and Vishnyakov, J. Chem. Phys. 122, 234108 (2005)] for calculations of chemical potential in dense and strongly inhomogeneous fluids is extended to multicomponent systems. The system of interest is simulated in a sample cell that is placed in chemical contact with several gauge cells of limited capacity, one gauge cell per component. Thus, each component can be exchanged between the sample cell and the respective gauge cell. The sample and gauge cells are immersed into the thermal bath of a given temperature. The size of the gauge cell controls the level of concentration fluctuations for the respective component in the sample cell. The chemical potentials are rigorously calculated from the equilibrium distribution of particles between the system and the gauges, and the results do not depend on the gauge size. For large systems, the chemical potentials can be accurately estimated from the average densities in the gauge cells. The proposed method was tested against the literature data on the vapor-liquid equilibrium in a binary mixture of subcritical and supercritical fluids and against the grand canonical and Widom insertion Monte Carlo methods for a binary mixture confined to a very narrow spherical pore. The method is specifically suitable for simulations of metastable and labile states in multicomponent confined fluids.

Molecular dynamics simulation of nanoscale distribution and mobility of water and dimethylmethylphosphonate in sulfonated polystyrene.

The interest in a better understanding of the specific interactions of phosphor-organic compounds and water with sulfonated polystyrene (sPS) is motivated by the use of block copolymers as protective membranes against chemical warfare agents. Using classical molecular dynamics simulations, we explored the nanoscale segregation and diffusion of water and nerve gas simulant dimethylmethylphosphonate (DMMP) in sPS neutralized with calcium counterions at different sulfonation and hydration levels. The water content was varied from 15 to 54% of dry polymer weight, and the DMMP content was varied from 0 to 100 wt %. We found that, in the 40% sulfonated polystyrene, water forms well defined aggregates, which grow in size as the hydration increases, reaching approximately 20 A at the maximum water content. In the 100% sulfonated polystyrene, the overall structure of hydrated polymer is more uniform with smaller water aggregates. Diffusion of water at the same number of water molecules per sulfonate group is faster at a lower sulfonation level. The solvation of sPS in water-DMMP binary mixtures was found to differ substantially from Nafion, where DMMP forms a layer between the hydropholic and hydrophobic subphases. In sPS with divalent Ca(2+) counterions, DMMP and water compete for the solvation of the sulfonate group. At high water and DMMP contents, the diffusion of DMMP turned out to be rather fast with a diffusion coefficient of ca. 30% of that of water. At the same time, water diffusion slows down as the DMMP concentration increases. This observation suggests that although sPS is permeable for both solvents, water and DMMP are partially segregated on the scale of 1-2 nm and have different pathways through the system. The nonuniform nanoscale distribution of water and DMMP in sPS is confirmed by analyses of different pair correlation functions. This feature may significantly affect the perm-selective properties of sPS-contained block copolymer membranes.

Interactions of silica nanoparticles in supercritical carbon dioxide.

We report molecular simulation studies on the interaction forces between silica nanoparticles in supercritical carbon dioxide at 318 K. Our goal is to find a better understanding of the interparticle solvation forces during rapid expansion of supercritical solutions. The parameters for interatomic potentials of fluid-fluid and solid-fluid interactions are obtained by fitting our simulations to (i) experimental bulk CO(2) phase diagram at a given temperature and pressure and (ii) CO(2) sorption isotherms on silica at normal boiling and critical temperatures. Our simulations show that the interaction forces between particles and supercritical CO(2) at near-critical pressure of p=69 atm (i.e., slightly below critical condition) reaches a minimum at distances of 0.5-0.8 nm between the outer surfaces of the particles and practically vanishes at distances of approximately 3 nm. The attraction is most prominent for densely hydroxylated particle surfaces that interact strongly with CO(2) via hydrogen bonds. The effective attraction between silica and CO(2) is significantly weaker for dehydroxylated particles. We also compared fluid sorption and interparticle forces between supercritical CO(2) and subcritical nitrogen vapor, and our results showed qualitative similarities, suggesting that the CO(2) configuration between the particles resembles a liquidlike junction.

Nanoparticle-Engendered Rupture of Lipid Membranes.

Tension-induced rupture of 1,2-dimyristoyl- sn-glycero-3-phosphocholine (DMPC) lipid membranes with encapsulated hydrophobic nanoparticles is elucidated using dissipative particle dynamics simulations. The dynamics of hole formation is studied, and a nanoparticle size-dependent relationship is established for the probability of membrane rupture within a given time as a function of the membrane tension. Two mechanisms of hole formation are explored: homogeneous nucleation and heterogeneous nucleation at the nanoparticle surface. While the kinetics of homogeneous nucleation in unloaded membranes complies with the predictions of the classical Deryagin-Gutop (DG) theory, the heterogeneous nucleation causes progressively lower lysis tensions as the particle size increases. The thermodynamics of heterogeneous nucleation is treated by introducing an effective contact angle at the three-phase, solvent-membrane-solid boundary into the DG theory. The proposed approach helps quantitatively interpret the simulation results and predict the membrane stability in real experiments with significantly larger (by many orders of magnitude) observation times and spatial scales.