Classical Liquids in Fractal Dimension

We introduce fractal liquids by generalizing classical liquids of integer dimensions $d = 1, 2, 3$ to a fractal dimension $d_f$. The particles composing the liquid are fractal objects and their configuration space is also fractal, with the same non-integer dimension. Realizations of our generic model system include microphase separated binary liquids in porous media, and highly branched liquid droplets confined to a fractal polymer backbone in a gel. Here we study the thermodynamics and pair correlations of fractal liquids by computer simulation and semi-analytical statistical mechanics. Our results are based on a model where fractal hard spheres move on a near-critical percolating lattice cluster. The predictions of the fractal Percus-Yevick liquid integral equation compare well with our simulation results.

Liquid pair correlations in four spatial dimensions: Theory versus simulation

Using liquid integral equation theory, we calculate the pair correlations of particles that interact via a smooth repulsive pair potential in d = 4 spatial dimensions. We discuss the performance of different closures for the Ornstein-Zernike equation, by comparing the results to computer simulation data. Our results are of relevance to understand crystal and glass formation in high-dimensional systems.

Dynamical density functional theory for the diffusion of injected Brownian particles

While the theory of diffusion of a single Brownian particle in confined geometries is well-established by now, we discuss here the theoretical framework necessary to generalize the theory of diffusion to dense suspensions of strongly interacting Brownian particles. Dynamical density functional theory (DDFT) for classical Brownian particles represents an ideal tool for this purpose. After outlining the basic ingredients to DDFT we show that it can be readily applied to flowing suspensions with time-dependent particle sources. Particle interactions lead to considerable layering in the mean density profiles, a feature that is absent in the trivial case of noninteracting, freely diffusing particles. If the particle injection rate varies periodically in time with a suitable frequency, a resonance in the layering of the mean particle density profile is predicted.

Coupling between bulk- and surface chemistry in suspensions of charged colloids

The ionic composition and pair correlations in fluid phases of realistically salt-free charged colloidal sphere suspensions are calculated in the primitive model. We obtain the number densities of all ionic species in suspension, including low-molecular weight microions, and colloidal macroions with acidic surface groups, from a self-consistent solution of a coupled physicochemical set of nonlinear algebraic equations and non-mean-field liquid integral equations. Here, we study suspensions of colloidal spheres with sulfonate or silanol surface groups, suspended in demineralized water that is saturated with carbon dioxide under standard atmosphere. The only input required for our theoretical scheme are the acidic dissociation constants pKa, and effective sphere diameters of all involved ions. Our method allows for an ab initio calculation of colloidal bare and effective charges, at high numerical efficiency.

Electrokinetic and hydrodynamic properties of charged-particles systems: From small electrolyte ions to large colloids

Dynamic processes in dispersions of charged spherical particles are of importance both in fundamental science, and in technical and bio-medical applications. There exists a large variety of charged-particles systems, ranging from nanometer-sized electrolyte ions to micron-sized charge-stabilized colloids. We review recent advances in theoretical methods for the calculation of linear transport coefficients in concentrated particulate systems, with the focus on hydrodynamic interactions and electrokinetic effects. Considered transport properties are the dispersion viscosity, self- and collective diffusion coefficients, sedimentation coefficients, and electrophoretic mobilities and conductivities of ionic particle species in an external electric field. Advances by our group are also discussed, including a novel mode-coupling-theory method for conduction-diffusion and viscoelastic properties of strong electrolyte solutions. Furthermore, results are presented for dispersions of solvent-permeable particles, and particles with non-zero hydrodynamic surface slip. The concentration-dependent swelling of ionic microgels is discussed, as well as a far-reaching dynamic scaling behavior relating colloidal long- to short-time dynamics.

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, all ions are approximated as non-overlapping hyperspheres with Coulomb interactions. Our spectral HNC solver is based on a Fourier-Bessel transform introduced by Talman [J. Comput. Phys., 29, 35 (1978)], 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 one thousand 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 primitive model results. For moderately large ion size- and charge-asymmetries, we present Molecular Dynamics simulation results that agree well with the approximate HNC pair correlations.

Dynamics of charged gibbsite platelets in the isotropic phase

We report on depolarized and non-depolarized dynamic light scattering, static light scattering, and static viscosity measurements on interacting charged gibbsite platelets suspended in dimethyl sulfoxide (DMSO). The average collective and (long-time) translational self-diffusion coefficients, and the rotational diffusion coefficient, have been measured as functions of the platelet volume fraction phi, up to the isotropic-liquid crystal (I/LC) transition. The non-depolarized intensity autocorrelation function, measured at low scattering wavenumbers, consists of a fast and a slowly decaying mode which we interpret as the orientationally averaged collective and translational self-diffusion coefficients, respectively. Both the rotational and the long-time self-diffusion coefficients decrease very strongly, by more than two orders of magnitude, in going from the very dilute limit to the I/LC transition concentration. A similarly strong decrease, with increasing phi, is observed for the inverse zero-strain limiting static shear viscosity. With increasing phi, increasingly strong shear-thinning is observed, accompanied by a shrinking of the low shear-rate Newtonian plateau. The measured diffusion coefficients are interpreted theoretically in terms of a simple model of effective charged spheres interacting by a screened Coulomb potential, with hydrodynamic interactions included. The disk-like particle shape, and the measured particle radius and thickness polydispersities, enter into the model calculations via the scattering amplitudes. The interaction-induced enhancement of the collective diffusion coefficient by more than a factor of 20 at larger volume fractions is well captured in the effective sphere model, whereas the strong declines both of the experimental translational and rotational self-diffusion coefficients are underestimated.