Heterogeneous crystallization in colloids and complex plasmas: the role of binary mobilities

Both charged colloidal suspensions and complex (dusty) plasmas represent classical many-body strongly coupled Coulomb systems. Here we discuss their basic properties and focus on their heterogeneous crystallization from an undercooled melt. In particular, a model with different mobilities is proposed which is realizable in binary mixtures of charged particles. Within this binary-mobility model, the crystallization behaviour near a structured wall is explored by Brownian dynamics computer simulations. As a result, the propagation velocity of the crystal-fluid interface is a nonmonotonic function of the mobility ratio (if expressed in terms of an averaged mobility).     

Effect of geometrical confinement on the interaction between charged colloidal suspensions

The effective interaction between charged colloidal particles confined between two planar like-charged walls is investigated using computer simulations of the primitive model describing asymmetric electrolytes. In detail, we calculate the effective force acting onto a single macroion and onto a macroion pair in the presence of slitlike confinement. For moderate Coulomb coupling, we find that this force is repulsive. Under strong-coupling conditions, however, the sign of the force depends on the distance to the plates and on the interparticle distance. In particular, the particle-plate interaction becomes strongly attractive for small distances which may explain the occurrence of colloidal crystalline layers near the plates observed in recent experiments.     

Highly asymmetric electrolytes in the primitive model: Hypernetted chain solution in arbitrary spatial dimensions

The pair‐correlation functions for fluid ionic mixtures in arbitrary spatial dimensions are computed in hypernetted chain (HNC) approximation. In the primitive model (PM), all ions are approximated as nonoverlapping hyperspheres with Coulomb interactions. Our spectral HNC solver is based on a Fourier‐Bessel transform introduced by Talman (J. Comput. Phys. 1978, 29, 35), with logarithmically spaced computational grids. Numeric efficiency for arbitrary spatial dimensions is a commonly exploited virtue of this transform method. Here, we highlight another advantage of logarithmic grids, consisting in efficient sampling of pair‐correlation functions for highly asymmetric ionic mixtures. For three‐dimensional fluids, ion size and charge‐ratios larger than 1000 can be treated, corresponding to hitherto computationally not accessed micrometer‐sized colloidal spheres in 1‐1 electrolyte. Effective colloidal charge numbers are extracted from our PM results. For moderately large ion size and charge‐asymmetries, we present molecular dynamics simulation results that agree well with the approximate HNC pair correlations. © 2013 Wiley Periodicals, Inc.     

Doubled heterogeneous crystal nucleation in sediments of hard sphere binary-mass mixtures

Crystallization during the sedimentation process of a binary colloidal hard spheres mixture is explored by Brownian dynamics computer simulations. The two species are different in buoyant mass but have the same interaction diameter. Starting from a completely mixed system in a finite container, gravity is suddenly turned on, and the crystallization process in the sample is monitored. If the Peclet numbers of the two species are both not too large, crystalline layers are formed at the bottom of the cell. The composition of lighter particles in the sedimented crystal is non-monotonic in the altitude: it is first increasing, then decreasing, and then increasing again. If one Peclet number is large and the other is small, we observe the occurrence of a doubled heterogeneous crystal nucleation process. First, crystalline layers are formed at the bottom container wall which are separated from an amorphous sediment. At the amorphous-fluid interface, a secondary crystal nucleation of layers is identified. This doubled heterogeneous nucleation can be verified in real-space experiments on colloidal mixtures.     

Simulation study of the equilibrium morphology in ionomers with different architectures

The influence of the side chain and the backbone segment length on the shape of the hydrophilic clusters in an ionomer with sulfonate‐terminated side chains is investigated by molecular dynamics computer simulations. Different ionomer architectures, ranging from Nafion‐like molecules to ionomer structures with extremely long side chains and backbone segments, are analyzed. It is shown that the size and the sulfonate population of hydrophilic clusters are differently affected by the variation of the backbone length in homologous membranes. Although the clusters swell when the backbone length is increased, they contain fewer head groups, and the average sulfonate–sulfonate separation in these swollen clusters becomes larger. An increase in the equivalent weight (EW) of homologous membranes thus hinders proton transport through hydrophilic channels. We also have shown that larger sulfonate clusters are formed in modified membranes with longer side chains. Compared to the case of membranes with shorter side chains but the same EW, those with longer side chains have a higher proton diffusivity. © 2010 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys, 2010     

Heterogeneous crystallization of hard and soft spheres near flat and curved walls

Crystallization represents a long-standing problem in statistical physics and is of great relevance for many practical and industrial applications. It often occurs in the presence of container walls or impurities, which are usually unavoidable or might even be desirable to facilitate crystallization by exploiting heterogeneous nucleation. Heterogeneous nucleation relies on a seed. Here we discuss the role of the seed and concentrate on a very generic situation, namely crystallization of hard and soft colloidal spheres in the presence of flat or curved hard walls. Curvature serves as a simple means to introduce a tunable mismatch between the seed-induced crystal lattice and the thermodynamically-favoured lattice. The mismatch induces distortions and elastic stress, which accumulate while the crystallite grows. This has an important consequence: once the crystallite reaches a critical size, it detaches from the seed allowing it to relax. The relaxed crystal continues to grow in the bulk, but crystallization ceases before reaching the seed, which now represents an impurity. Therefore, while seeds favour nucleation, any mismatch, like the seed curvature or an incommensurate structure, induces unfavourable distortions and can lead to the detachment of the crystallite. An additional mechanism to relax distortions is available to soft spheres, which can exploit their interaction potential and possibly deform. The different multi-step processes have been investigated by confocal microscopy, which provides particle-level information, and compared to computer simulations and theoretical results.     

Nonlinear optics of nontoxic nanomaterials

We studied the linear and nonlinear optical properties of halloysite nanotubes using Z-scan technique. Halloysite is alumina silicate clay rolled into 50 nm diameter hollow cylinders, where the silica layer is at the outer surface of the tube and alumina layer is in the inner surface. Optical absorption spectra show an absorption peak around 600 nm. Open aperture Z-scan measurements using 3 ns laser pulses at 532 nm reveal two-photon induced absorption. The closed aperture Z-scan indicates a positive nonlinear refractive index. As these nanotubes are nontoxic and biocompatible, they have advantages over conventional carbon nanotubes for biomedical applications.     

Overscreening in colloidal systems: the attractive branch of pair macroion interaction

The effective pair macroion interaction for colloidal systems, obtained in general form by Trigger, Podloubny and Schram (TPS potential) is considered for the case of a strong Coulomb interaction between the macroion and counterion subsystems of a colloid. For this potential a good agreement with the experiments on the shear modules as function of the volume fraction is found.     

Role of electrostatic forces in cluster formation in a dry ionomer

This simulation study investigates the dependence of the structure of dry Nafion-like ionomers on the electrostatic interactions between the components of the molecules. In order to speed equilibration, a procedure was adopted which involved detaching the side chains from the backbone and cutting the backbone into segments, and then reassembling the macromolecule by means of a strong imposed attractive force between the cut ends of the backbone, and between the nonionic ends of the side chains and the midpoints of the backbone segments. Parameters varied in this study include the dielectric constant, the free volume, side chain length, and strength of head group interactions. A series of coarse-grained mesoscale simulations shows the morphology to depend sensitively on the ratio of the strength of the dipole-dipole interactions between the side-chain acidic end groups to the strength of the other electrostatic components of the Hamiltonian. Examples of the two differing morphologies proposed by Gierke and co-workers [J. Polym. Sci., Polym. Phys. Ed. 19, 1687 (1981); Macromolecules 15, 101 (1982); J. Membr. Sci. 13, 307 (1982)] and by Gebel [Fuel Cells 5, 261 (2005); Macromolecules 37, 7772 (2004)] emerge from our simulations.     

Simulation study of poled low-water ionomers with different architectures

The role of the ionomer architecture in the formation of ordered structures in poled membranes is investigated by molecular dynamics computer simulations. It is shown that the length of the sidechain L(s) controls both the areal density of cylindrical aggregates N(c) and the diameter of these cylinders in the poled membrane. The backbone segment length L(b) tunes the average diameter D(s) of cylindrical clusters and the average number of sulfonates N(s) in each cluster. A simple empirical formula is noted for the dependence of the number density of induced rod-like aggregates on the sidechain length L(s) within the parameter range considered in this study.