Structure and equation of state of interaction site models for disc-shaped lamellar colloids

Early stages of aggregation in fluid mixtures of dimers and spheres: a theoretical and simulation study

We use Monte Carlo simulation and the Reference Interaction Site Model (RISM) theory of molecular fluids to investigate a simple model of colloidal mixture consisting of dimers, made up of two tangent hard monomers of different size, and hard spheres. In addition to steric repulsion, the two species interact via a square-well attraction only between small monomers and spheres. Recently, we have characterized the low-temperature regime of this mixture by Monte Carlo, reporting on the spontaneous formation of a wide spectrum of supramolecular aggregates [Prestipino et al., J. Phys. Chem. B, 2019, 123, 9272]. Here we focus on a regime of temperatures where, on cooling, the appearance of local inhomogeneties first, and the early stages of aggregation thereafter, are observed. In particular, we find signatures of aggregation in the onset of a low-wavevector peak in the structure factors of the mixture, as computed by both theory and simulation. Then, we link the structural information to the microscopic arrangement through a detailed cluster analysis of Monte Carlo configurations. In this regard, we devise a novel method to compute the maximum distance for which two spheres can be regarded as bonded together, a crucial issue in the proper identification of fluid aggregates. The RISM theory provides relatively accurate structural and thermodynamic predictions in comparison with Monte Carlo, but with slightly degrading performances as the fluid progresses inside the locally inhomogeneous phase. Our study certifies the efficacy of the RISM approach as a useful complement to numerical simulation for a reasoned analysis of aggregation properties in colloidal mixtures.

Isotropic-nematic phase transition of uniaxial variable softness prolate and oblate ellipsoids

Onsager's theory of the isotropic-nematic phase separation of rod shaped particles is generalized to include particle softness and attractions in the anisotropic interparticle force field. The procedure separates a scaled radial component from the angular integral part, the latter being treated in essentially the same way as in the original Onsager formulation. Building on previous treatments of more idealised hard-core particle models, this is a step toward representing more realistic rod-like systems and also allowing temperature (and in principle specific chemical factors) to be included at a coarse grained level in the theory. The focus of the study is on the coexisting concentrations and associated coexistence properties. Prolate and oblate ellipsoids are considered in both the small and very large aspect ratio limits. Approximations to the terms in the angular integrals derived assuming the very large (prolate) and very small (oblate) aspect ratios limits are compared with the formally exact treatment. The approximation for the second virial coefficient matches the exact solution for aspect ratios above about 20 for the prolate ellipsoids and less than ca. 0.05 for the oblate ellipsoids from the numerical evaluation of the angular integrals. The temperature dependence of the coexistence density could be used to help determine the interaction potential of two molecules. The method works at temperatures above a certain threshold temperature where the second virial coefficient is positive.

Development of molecular closures for the reference interaction site model theory with application to square-well and Lennard-Jones homonuclear diatomics

Inspired by significant improvements obtained for the performances of the polymer reference interaction site model (PRISM) theory of the fluid phase when coupled with 'molecular closures' (Schweizer and Yethiraj 1993 J. Chem. Phys. 98 9053), we exploit a matrix generalization of this concept, suitable for the more general RISM framework. We report a preliminary test of the formalism, as applied to prototype square-well homonuclear diatomics. As for the structure, comparison with Monte Carlo shows that molecular closures are slightly more predictive than their 'atomic' counterparts, and thermodynamic properties are equally accurate. We also devise an application of molecular closures to models interacting via continuous, soft-core potentials, by using well established prescriptions in liquid state perturbation theories. In the case of Lennard-Jones dimers, our scheme definitely improves over the atomic one, providing semi-quantitative structural results, and quite good estimates of internal energy, pressure and phase coexistence. Our finding paves the way to a systematic employment of molecular closures within the RISM framework to be applied to more complex systems, such as molecules constituted by several non-equivalent interaction sites.

Generalized Onsager theory for strongly anisometric patchy colloids

The implications of soft "patchy" interactions on the orientational disorder-order transition of strongly elongated colloidal rods and flat disks is studied within a simple Onsager-van der Waals density functional theory. The theory provides a generic framework for studying the liquid crystal phase behaviour of highly anisometric cylindrical colloids which carry a distinct geometrical pattern of repulsive or attractive soft interactions localized on the particle surface. In this paper, we apply our theory to the case of charged rods and disks for which the local electrostatic interactions can be described by a screened-Coulomb potential. We consider infinitely thin rod like cylinders with a uniform line charge and infinitely thin discotic cylinders with several distinctly different surface charge patterns. Irrespective of the backbone shape the isotropic-nematic phase diagrams of charged colloids feature a generic destabilization of nematic order at low ionic strength, a dramatic narrowing of the biphasic density region, and a reentrant phenomenon upon reducing the electrostatic screening. The low screening regime is characterized by a complete suppression of nematic order in favor of positionally ordered liquid crystal phases.

Liquid crystal phase behaviour of attractive disc-like particles

We employ a generalized van der Waals-Onsager perturbation theory to construct a free energy functional capable of describing the thermodynamic properties and orientational order of the isotropic and nematic phases of attractive disc particles. The model mesogen is a hard (purely repulsive) cylindrical disc particle decorated with an anisotropic square-well attractive potential placed at the centre of mass. Even for isotropic attractive interactions, the resulting overall inter-particle potential is anisotropic, due to the orientation-dependent excluded volume of the underlying hard core. An algebraic equation of state for attractive disc particles is developed by adopting the Onsager trial function to characterize the orientational order in the nematic phase. The theory is then used to represent the fluid-phase behaviour (vapour-liquid, isotropic-nematic, and nematic-nematic) of the oblate attractive particles for varying values of the molecular aspect ratio and parameters of the attractive potential. When compared to the phase diagram of their athermal analogues, it is seen that the addition of an attractive interaction facilitates the formation of orientationally-ordered phases. Most interestingly, for certain aspect ratios, a coexistence between two anisotropic nematic phases is exhibited by the attractive disc-like fluids.

Hindrance function for sedimentation and creaming of colloidal disks

We report the first measurement of the hindrance function for sedimentation or creaming of disk-shaped colloids via the analytical centrifugation. Disks align with the external flow right above a volume fraction of a few percent, and this effect is extremely sensitive to the aspect ratio of disks. Due to this alignment effect, disk sedimentation or creaming demonstrates distinct trends in dilute and semidilute regions.

Interaction of cylindrical polymer brushes in dilute and semi-dilute solution

We present a systematic study of flexible cylindrical brush-shaped macromolecules in a good solvent by small-angle neutron scattering (SANS), static light scattering (SLS), and by dynamic light scattering (DLS) in dilute and semi-dilute solution. The SLS and SANS data extrapolated to infinite dilution lead to the shape of the polymer that can be modeled in terms of a worm-like chain with a contour length of 380 nm and a persistence length of 17.5 nm. SANS data taken at higher polymer concentration were evaluated by using the polymer reference interaction site model (PRISM). We find that the persistence length reduce from 17.5 nm at infinite dilution to 5.3 nm at the highest concentration (volume fraction 0.038). This is comparable with the decrease of the persistence length in semi-dilute concentration predicted theoretically for polyelectrolytes. This finding reveals a softening of stiffness of the polymer brushes caused by their mutual interaction.

Nonequilibrium steady states in fluids of platelike colloidal particles

Nonequilibrium steady states in an open system connecting two reservoirs of platelike colloidal particles are investigated by means of a recently proposed phenomenological dynamic density functional theory [M. Bier and R. van Roij, Phys. Rev. E 76, 021405 (2007)]. The platelike colloidal particles are approximated within the Zwanzig model of restricted orientations, which exhibits an isotropic-nematic bulk phase transition. Inhomogeneities of the local chemical potential generate a diffusion current which relaxes to a nonvanishing value if the two reservoirs coupled to the system sustain different chemical potentials. The relaxation process of initial states towards the steady state turns out to comprise two regimes: a smoothening of initial steplike structures followed by an ultimate relaxation of the slowest diffusive mode. The position of a nonequilibrium interface and the particle current of steady states depend nontrivially on the structure of the reservoirs due to the coupling between translational and orientational degrees of freedom of the fluid.

Relaxation dynamics in fluids of platelike colloidal particles

The relaxation dynamics of a model fluid of platelike colloidal particles is investigated by means of a phenomenological dynamic density functional theory. The model fluid approximates the particles within the Zwanzig model of restricted orientations. The driving force for time dependence is expressed completely by gradients of the local chemical potential, which in turn is derived from a density functional-hydrodynamic interactions are not taken into account. These approximations are expected to lead to qualitatively reliable results for low densities like those within the isotropic-nematic two-phase region. The formalism is applied to model an initially spatially homogeneous stable or metastable isotropic fluid which is perturbed by switching a two-dimensional array of Gaussian laser beams. Switching on the laser beams leads to an accumulation of colloidal particles in the beam centers. If the initial chemical potential and the laser power are large enough, a preferred orientation of particles occurs, breaking the symmetry of the laser potential. After switching off the laser beams again, the system can follow different relaxation paths: It either relaxes back to the homogeneous isotropic state or it forms an approximately elliptical high-density core which is elongated perpendicular to the dominating orientation in order to minimize the surface free energy. For large supersaturations of the initial isotropic fluid, the high-density cores of neighboring laser beams of the two-dimensional array merge into complex superstructures.

Structure factor and thermodynamics of rigid dendrimers in solution

The "polymer reference interaction site model" (PRISM) integral equation theory is used to determine the structure factor of rigid dendrimers in solution. The theory is quite successful in reproducing experimental structure factors for various dendrimer concentrations. In addition, the structure factor at vanishing scattering vector is calculated via the compressibility equation using scaled particle theory and fundamental measure theory. The results as predicted by both theories are systematically smaller than the experimental and PRISM data for platelike dendrimers.

Surface properties of fluids of charged platelike colloids

Surface properties of mixtures of charged platelike colloids and salt in contact with a charged planar wall are studied within density functional theory. The particles are modeled by hard cuboids with their edges constrained to be parallel to the Cartesian axes corresponding to the Zwanzig model [J. Chem. Phys. 39, 1714 (1963)] and the charges of the particles are concentrated at their centers. The density functional applied is an extension of a recently introduced functional for charged platelike colloids. It provides a qualitative approach because it does not determine the relation between the actual and the effective charges entering into the model. Technically motivated approximations, such as using the Zwanzig model, are expected not to influence the results qualitatively. Analytically and numerically calculated bulk and surface phase diagrams exhibit first-order wetting for sufficiently small macroion charges and isotropic bulk order as well as first-order drying for sufficiently large macroion charges and nematic bulk order. The asymptotic wetting and drying behaviors are investigated by means of effective interface potentials which turn out to be asymptotically the same as for a suitable neutral system governed by isotropic nonretarded dispersion forces. Wetting and drying points as well as predrying lines and the corresponding critical points have been located numerically. A crossover from monotonic to nonmonotonic electrostatic potential profiles upon varying the surface charge density has been observed. Nonmonotonic electrostatic potential profiles are equivalent to the occurrence of charge inversion. Due to the presence of both the Coulomb interactions and the hard-core repulsions, the surface potential and the surface charge do not vanish simultaneously, i.e., the point of zero charge and the isoelectric point of the surface do not coincide.

The Solution Structure of Stilbenoid Dendrimers: A Small-Angle Scattering Study

The spatial structure of a stilbenoid dendrimer is investigated by small-angle X-ray scattering (SAXS) and small-angle neutron scattering (SANS) in dilute solution. All measurements are performed in toluene. The dendrimer consists of a stilbenoid scaffold with appended hexyloxy chains. SAXS is mainly sensitive to the dendrimer scaffold whereas SANS intensity, measured in fully deuterated toluene, derives from the solute molecules. The resulting SAXS and SANS intensities are analyzed by comparison with various models. It is found that the model of a circular disk gives the best description of the data. SAXS data demonstrate that the stilbenoid scaffold is flat as expected for benzene rings conjugated through vinylene units. Thus, it can be described by a circular disk with a radius of 1.6 nm and a thickness of 0.7 nm. SANS, on the other hand, shows that the hexyloxy chains are not confined in the plane defined by the core. This is derived from modeling the SANS data with a much thicker circular disk (radius: 2.4 nm, thickness: 1.8 nm). The structure factor S(q), describing the interaction of the dendrimers at higher concentrations, is modeled quantitatively with the "polymer reference interaction site model" (PRISM) integral equation formalism for hard plates such as particles. Here the structural data obtained from the analysis of the SANS data are used so that no new adjustable parameter is necessary for this description.

Effective interaction of charged platelets in aqueous solution: investigations of colloid laponite suspensions by static light scattering and small-angle x-ray scattering

We study dilute aqueous solutions of charged disklike mineral particles (laponite) by a combination of static light scattering (SLS) and small-angle x-ray scattering (SAXS). Laponite solutions are known to form gels above a certain critical concentration that must be described as nonequilibrium states. Here we focus on the investigation by SLS and SAXS at concentrations below gelation (c<0.016 g/L) and at low concentrations of added salt (0.001M and 0.005M). Thus, we have obtained the scattering function of single Laponite platelets as well as the structure factor describing their interaction at finite concentration. A detailed analysis of the combined sets of data proves that the solutions are in a well-defined equilibrium state. Moreover, this analysis demonstrates the internal consistency and accuracy of the scattering functions obtained at finite concentrations. We find that laponite particles interact through an effective pair potential that is attractive on short range but repulsive on longer range. This finding demonstrates that Laponite solutions exhibit only a limited stability at the concentration of added salt used herein. Raising the ionic strength to 0.005M already leads to slow flocculation as is evidenced from the enhanced scattering intensity at smallest scattering angles. All data strongly suggest that the gelation occurring at higher concentration is related to aggregation.