Free volumes and the anomalous self-diffusivity of attractive colloids

Free volume theories for the dynamics of dense fluids commonly assume (i) that diffusivity increases with average free volume per particle and (ii) that the size distribution of free volumes can be approximated by that of an equivalent hard-sphere reference system. We use molecular simulations to demonstrate that these assumptions break down when one considers concentrated suspensions of particles with short-range attractions. In these systems, self-diffusivity shows nonmonotonic dependencies on both average free volume and the strength of the interparticle attraction. Moreover, when interparticle attractions are strong, the shape of the free volume distribution is qualitatively different than that of the corresponding hard-sphere reference fluid. We propose a conceptual revision to the traditional free volume perspective that takes into account both the size distribution and the persistence time of the free volumes, and we demonstrate that it can qualitatively capture the disparate behaviors of a model fluid with short-range attractions and its hard-sphere reference fluid.     

Probing memory effects in confined fluids via diffusion measurements

Confinement of fluids in porous materials is widely exploited in a variety of technologies, including chemical conversion by heterogeneous catalysis and adsorption separations. Important fundamental phenomena associated with many-molecule interactions occur in such systems, including a remarkably long "memory" of the past when the actual amount of molecules in the pores dramatically depends on the history of how the external conditions have been changed. We demonstrate that the intrinsic diffusivity as measured by NMR serves as an excellent probe of the history-dependent states of the confined fluid. A remarkable feature of our results are differences in diffusivity between out-of-equilibrium states with the same density within the hysteresis loop. This reflects different spatial distributions of the confined fluid that accompany the arrested equilibration of the system in this region.     

Does confining the hard-sphere fluid between hard walls change its average properties?

We use grand canonical transition-matrix Monte Carlo and discontinuous molecular dynamics simulations to generate precise thermodynamic and kinetic data for the equilibrium hard-sphere fluid confined between smooth hard walls. These simulations show that the pronounced inhomogeneous structuring of the fluid normal to the confining walls, often the primary focus of density functional theory studies, has a negligible effect on many of its average properties over a surprisingly broad range of conditions. We present one consequence of this insensitivity to confinement: a simple analytical equation relating the average density of the confined fluid to that of the bulk fluid with equal activity. Nontrivial implications of confinement for average fluid properties do emerge in this system, but only when the fluid is both (i) dense and (ii) confined to a gap smaller than approximately three particle diameters. For this limited set of conditions, we find that "in-phase" oscillatory deviations in excess entropy and self-diffusivity (relative to the behavior of the bulk fluid at the same average density) occur as a function of gap size. These paired thermodynamic/kinetic deviations from bulk behavior appear to reflect the geometric packing frustration that arises when the confined space cannot naturally accommodate an integer number of particle layers.     

Phase behavior of model confined fluids. Influence of substrate-fluid interaction strength

We examine the relationship between the macroscopic phase behavior of nanoconfined fluids and the nature of microscopic interactions between a confining substrate and fluid. Two model slit-pore systems are explored using grand canonical transition-matrix Monte Carlo simulation. One system consists of a square-well fluid interacting with a square-well substrate, and the other contains an embedded point charge model of lysozyme interacting with a mica surface. Fluid phase diagrams are constructed for a broad range of substrate conditions. Our results indicate that one observes a maximum in the critical temperature of the fluid phase envelope upon variation of substrate strength for a given slit width. Both systems studied exhibit such maxima at intermediate wall strength. The physical rationale for this observation suggests that this behavior should be generally expected. We introduce two metrics that enable one to predict conditions that produce maxima in critical temperature. The first is related to the contact angle a fluid develops at a single confining substrate. The second is based upon virial coefficient information and requires knowledge of the substrate-fluid and fluid-fluid interaction potentials only.     

Simulations of the PDF functions for dilute colloidal suspensions of molecular particles flowing in mesopores with rough surface boundaries

Simulations have been carried out to analyze the dynamics of dilute colloidal suspensions of macromolecular particles in solutions flowing in pores, subject to hydrodynamic forces, Brownian motion and stochastic collisions at rough pore boundaries in a two-dimensional spatial frame. A theoretical model is developed and intensively analyzed for the treatment of the mechanical restitution of the particles due to dynamic collisions at these boundaries. In particular we are able to calculate the Probability distribution functions for the spatial positions and the orientations of rod-like particles inside the pores. The results are presented for different widths of pore channels referenced to the size of a rod-like particle. These simulations are general in the sense that they are developed for confining and open pore channels, rough at the nano scale. The simulations also permit calculating the nematic order parameters for colloidal suspensions; the model calculation is applied for dilute colloidal suspensions of carbon nano-tubes in an aqueous single-stranded DNA solution flowing inside pores. Our calculated nematic order results for dilute suspensions of particles of known lengths flowing inside porous systems should indicate, when coupled to birefringence and dichroism experimental results, the possibility to estimate the pore widths for these systems.     

Confinement, entropy, and single-particle dynamics of equilibrium hard-sphere mixtures

We use discontinuous molecular dynamics and grand-canonical transition-matrix Monte Carlo simulations to explore how confinement between parallel hard walls modifies the relationships between packing fraction, self-diffusivity, partial molar excess entropy, and total excess entropy for binary hard-sphere mixtures. To accomplish this, we introduce an efficient algorithm to calculate partial molar excess entropies from the transition-matrix Monte Carlo simulation data. We find that the species-dependent self-diffusivities of confined fluids are very similar to those of the bulk mixture if compared at the same, appropriately defined, packing fraction up to intermediate values, but then deviate negatively from the bulk behavior at higher packing fractions. On the other hand, the relationships between self-diffusivity and partial molar excess entropy (or total excess entropy) observed in the bulk fluid are preserved under confinement even at relatively high packing fractions and for different mixture compositions. This suggests that the excess entropy, calculable from classical density functional theories of inhomogeneous fluids, can be used to predict some of the nontrivial dynamical behaviors of fluid mixtures in confined environments.     

Extending the simple weighted density approximation for a hard-sphere fluid to a Lennard-Jones fluid II. Application

A simple weighted density approximation (SWDA) was extended to nonuniform Lennard-Jones fluids by following the spirit of a partitioned density function theory [S. Zhou, Phys. Rev. E 68 (2003) 061201] and mapping the hard-core part onto an effective hard-sphere fluid whose higher order terms beyond the second order of the functional perturbation expansion are treated by the SWDA. The resultant DFT formalism performs well for Lennard-Jones fluids under the influence of diverse external fields. With the present DFT formalism, we investigate in detail the structure and adsorption properties of a low-density LJ gas in a spherical cavity with a wall consisting of hard-sphere or LJ particles. It was found that when the cavity wall exerts an attractive external potential on the LJ particles in the cavity, the excess adsorption decreases as the temperature increases, while when the cavity wall exerts a hard repulsive external potential on the LJ particles in the cavity, the excess adsorption increases as the temperature increases.     

高分子链在球内表面上的吸附/排空转变

应用自洽场理论(SCFT)研究了受限于球内的高分子溶液的结构,重点关注高分子链在受限壁附近的行为.根据自洽场理论数值计算结果,讨论了球半径、高分子与球限制壁的相互作用、高分子平均浓度等因素对球内高分子浓度分布的影响.从高分子浓度分布和吸附/排空层厚度可以发现,在一定的条件下,受限的高分子在受限壁上会发生吸附/排空转变.吸附/排空转变与受限球大小、高分子链长和平均浓度,以及高分子链与受限壁之间相互作用都有关系.理论预测发生吸附/排空转变时的高分子与球限制壁的临界相互作用参数与链长的倒数成线性关系,且斜率与球半径有关.限制球越小,要发生吸附/排空转变,需要高分子与球之间有更大的临界吸引能.     

A simulation study of electrostatic effects on mixed ionic micelles confined between two parallel charged plates

Confined colloidal systems have been the subject of extensive theoretical and experimental research, and the recent observation of long-range like-charge attraction in such systems has only highlighted their peculiar behavior. On the other hand, surfactant solutions are often used in small confined space, yet their behavior in confinement has received relatively little attention. A distinct feature of confined self-assembling systems is that the aggregates are capable of adjusting their composition, size, and shape in response to their external environment, which may lead to very different phase characteristics compared to bulk solutions. The primary objective of this study is to explore the effects of varying micelle composition on the structural behavior of a confined mixed ionic micellar solution. Mesoscale canonical Monte Carlo simulations were used to probe the structure of the confined solution, while a molecular-thermodynamic model was used to systematically account for the change in micelle size as we varied its composition. Significant micelle ordering was found under certain conditions, which implies that large deviations from the minimum-energy micelle configuration may not be entropically favorable. Accumulation of micelles along the midplane was observed when the confining walls are weakly charged, suggesting that micelle shape transformation should be considered in more detail. On the other hand, with high wall charge density, apparent attraction was found between like-charged micelles and wall. These findings point to the need for a more quantitative theoretical treatment in describing surfactant self-assembly in confined geometries.     

Monte Carlo Studies on the Deformation Behavior of Entangled Polymer Chains

In this work we present a continuous three‐dimensional bond‐fluctuation model (CBFM) in combination with the method of confined self‐avoiding walks (CSAW). This method enables us to analyse both the compression and the stretching regime of the deformation process of two entangled polymer chains. We studied the deformation behavior in respect to the number of entanglements, the distance of the end monomers on the confining surfaces and the orientation of the entanglements to the deformation axis. Our analysis of the behavior of the force and the structural properties of the systems during the deformation process leads us to the assumption that the entanglements act as one topological crosslink with variable strength.     

Wall Effects on Spheres Settling Through Non-Newtonian Fluid Media in Cylindrical Tubes

Experimental study is performed to understand and quantify the wall and eccentric retardation effects on spheres settling in shear thinning and shear-thinning viscoelastic fluids over a wide range of diameter ratios (0.02 < λ < 0.9). The four-parameter Carreau viscosity equation has been chosen to represent the apparent viscosity-shear rate of polyacrylamide solutions. Two new wall factor corrections are presented with excellent agreement compared to experimental data.

The terminal settling velocity of a sphere in bounded fluid is significantly reduced by the presence of confining boundaries, named wall retardation effect that decreases due to the shear-thinning behavior of power law fluids, which is weaken further by the elastic effect of viscoelastic fluids. The wall factors of spheres settling in viscoelastic fluids increase at low ξ up to 50, followed by a horizontal confidence region (0.7 ≤ f ≤ 1) at high ξ. In this region, the wall factor is mainly dominated by fluids’ elasticity, which is more distinguished for small spheres. As the settling spheres approach to the wall (b/R → 1), the neighboring wall exert more intensive retardation that reduce the terminal settling velocity greatly when b/R > 0.6 in pure shear-thinning fluids, and the extra retardation effect of nearby wall increases at high concentration due to the enhanced non-Newtonian property. In contrast, the eccentric effect on settling velocity in viscoelastic fluids is cut down greatly by the fluid's elasticity, which is negligible.     

Structural transition of nematic liquid crystal in cylindrical capillary as a result of the annihilation of two point defects

We study the annihilation of hedgehog-antihedgehog defects in confined nematic liquid crystals using Brownian molecular dynamics simulations. After the collision, merging of defects, and building a loop disclination structure, system can experience a structural transition into another nematic structure, triggered by a nucleation of loop disclination structure. In our rough theoretical approach we calculate the size of the emerged loop structure as the function of the typical size of the confining cavity. Attention is paid also to the dynamics of the loop structure after collision.     

Local structure and density fluctuations in confined fluids

Spatial confinement modifies the microscopic structure of dense fluids, thereby inducing for example structural forces between the confining walls. However, confinement also modifies the fluids' density fluctuations, resulting in more elusive but equally important effects. In this brief review it is shown that both of these phenomena are naturally analyzed using the confined fluid's pair densities, which have recently become also experimentally accessible. Two particular topics are discussed, namely, the mechanisms of oscillatory density profiles and ensuing solvation forces in dense confined fluids as well as the behavior of liquids in solvophobic confinement.     

Confinement induced critical micelle concentration shift

In this work, extensive lattice Monte Carlo simulations were performed to investigate the influence of confinement on critical micelle concentration (CMC). It is found that the CMC of surfactants in a confined space is shifted from its bulk value, and the shift is affected by the presence of the confining boundaries, which induces both the finite size effect and the wall-surfactant interaction. In general, for strongly confined system (the system with narrow pore size), the finite size effect dominates the CMC shift because the confined space cannot accommodate fully developed micelles, and the rapid increase of the entropic loss due to the decrease of the pore size results in the rapid increase of CMC. In contrast, for a weakly confined space, the CMC shift depends on the interaction between the walls and surfactants. For the systems with two weakly hydrophilic surfaces, the local density depletion of the surfactants near the walls results in lower CMCs than the bulk value, and the CMC shifts to a higher value as the pore size increases. For the systems with moderately hydrophilic surfaces, the shifts of CMCs show a similar behavior as those for weakly hydrophilic surfaces, but the CMCs are near their bulk values in the range of weak confinement. For the systems with strongly attractive wall-surfactant interactions, the strong adsorption also results in lower CMCs than their bulk value, but the CMCs decrease with the increase of pore size.     

Density functional approach to the adsorption of spherical molecules on a surface modified with attached short chains

A density functional and Monte Carlo simulation study of end-grafted polymers immersed by simple fluids is presented. The polymer molecules are modeled as freely jointed tangent hard spheres with the end segments linked to the surface. The authors analyze an influence of the chain length, the grafting density, and a nature of solvent on the brush structure. Adsorption of hard-sphere mixtures on the modified surface is also discussed. The theory precisely approximates simulation data.     

Model for the free-volume distributions of equilibrium fluids

We introduce and test via molecular simulation a simple model for predicting the manner in which interparticle interactions and thermodynamic conditions impact the single-particle free-volume distributions of equilibrium fluids. The model suggests a scaling relationship for the density-dependent behavior of the hard-sphere system. It also predicts how the second virial coefficients of fluids with short-range attractions affect their free-volume distributions.     

Confined polyelectrolytes: The complexity of a simple system

The interaction between polyelectrolytes and counterions in confined situations and the mutual relationship between chain conformation and ion condensation is an important issue in several areas. In the biological field, it assumes particular relevance in the understanding of the packaging of nucleic acids, which is crucial in the design of gene delivery systems. In this work, a simple coarse‐grained model is used to assess the cooperativity between conformational change and ion condensation in spherically confined backbones, with capsides permeable to the counterions. It is seen that the variation on the degree of condensation depends on counterion valence. For monovalent counterions, the degree of condensation passes through a minimum before increasing as the confining space diminishes. In contrast, for trivalent ions, the overall tendency is to decrease the degree of condensation as the confinement space also decreases. Most of the particles reside close to the spherical wall, even for systems in which the density is higher closer to the cavity center. This effect is more pronounced, when monovalent counterions are present. Additionally, there are clear variations in the charge along the concentric layers that cannot be totally ascribed to polyelectrolyte behavior, as shown by decoupling the chain into monomers. If both chain and counterions are confined, the formation of a counterion rich region immediately before the wall is observed. Spool and doughnut‐like structures are formed for stiff chains, within a nontrivial evolution with increasing confinement. © 2015 Wiley Periodicals, Inc.     

Close-packed structures and phase diagram of soft spheres in cylindrical pores

It is shown for a model system consisting of spherical particles confined in cylindrical pores that the first ten close-packed phases are in one-to-one correspondence with the first ten ways of folding a triangular lattice, each being characterized by a roll-up vector like the single-walled carbon nanotube. Phase diagrams in pressure-diameter and temperature-diameter planes are obtained by inherent-structure calculation and molecular dynamics simulation. The phase boundaries dividing two adjacent phases are infinitely sharp in the low-temperature limit but are blurred as temperature is increased. Existence of such phase boundaries explains rich, diameter-sensitive phase behavior unique for cylindrically confined systems.     

Coarse-grained simulations of ABA amphiphilic triblock copolymer solutions in thin films

We study a coarse-grained model of A(10)-B(20)-A(10) amphiphilic triblock copolymers in aqueous solution under confinement. We focus on the influence of the wall interaction on the morphology of the ensuing self-assembled structures. We also study the dynamics of the polymers. All our simulations are confined between two walls. We study three different combinations of walls: hydrophobic and hydrophobic, hydrophobic and hydrophilic, hydrophilic and hydrophilic. We moreover elucidate the concentration influence. The conformation and behavior of the copolymer in strongly confined systems depend on the type of wall interaction and concentration.     

A cubic perturbed hard-chain equation of state

A cubic equation of state is developed on the basis of perturbation theory. The equation is an association of three segments: the hard-sphere, the hard-chain, and the attraction. The expression for each segment was invoked from approximations of computer simulations of rigorous molecular theories of fluids, but compromised to some extent accuracy and theory for simplicity. This model equation is shown to be potentially capable of describing the PVT behavior of real fluids. As limiting cases, the new equation is reduced to expressions for the hard-sphere and the hard-body fluids. It also represents square-well fluids when the hard-chain contribution is eliminated. The square-well equation was found satisfactory in conforming with the molecular simulation results for square-well fluids and their mixtures.