No Arabic abstract
The structure of polydisperse hard sphere fluids, in the presence of a wall, is studied by the Rosenfeld density functional theory. Within this approach, the local excess free energy depends on only four combinations of the full set of density fields. The case of continuous polydispersity thereby becomes tractable. We predict, generically, an oscillatory size segregation close to the wall, and connect this, by a perturbation theory for narrow distributions, with the reversible work for changing the size of one particle in a monodisperse reference fluid.
The solid-solid coexistence of a polydisperse hard sphere system is studied by using the Monte Carlo simulation. The results show that for large enough polydispersity the solid-solid coexistence state is more stable than the single-phase solid. The two coexisting solids have different composition distributions but the same crystal structure. Moreover, there is evidence that the solid-solid transition terminates in a critical point as in the case of the fluid-fluid transition.
The coupling-parameter method, whereby an extra particle is progressively coupled to the rest of the particles, is applied to the sticky-hard-sphere fluid to obtain its equation of state in the so-called chemical-potential route ($mu$ route). As a consistency test, the results for one-dimensional sticky particles are shown to be exact. Results corresponding to the three-dimensional case (Baxters model) are derived within the Percus-Yevick approximation by using different prescriptions for the dependence of the interaction potential of the extra particle on the coupling parameter. The critical point and the coexistence curve of the gas-liquid phase transition are obtained in the $mu$ route and compared with predictions from other thermodynamics routes and from computer simulations. The results show that the $mu$ route yields a general better description than the virial, energy, compressibility, and zero-separation routes.
We compute the rheological properties of inelastic hard spheres in steady shear flow for general shear rates and densities. Starting from the microscopic dynamics we generalise the Integration Through Transients (textsc{itt}) formalism to a fluid of dissipative, randomly driven granular particles. The stress relaxation function is computed approximately within a mode-coupling theory---based on the physical picture, that relaxation of shear is dominated by slow structural relaxation, as the glass transition is approached. The transient build-up of stress in steady shear is thus traced back to transient density correlations which are computed self-consistently within mode-coupling theory. The glass transition is signalled by the appearance of a yield stress and a divergence of the Newtonian viscosity, characterizing linear response. For shear rates comparable to the structural relaxation time, the stress becomes independent of shear rate and we observe shear thinning, while for the largest shear rates Bagnold scaling, i.e., a quadratic increase of shear stress with shear rate, is recovered. The rheological properties are qualitatively similar for all values of $varepsilon$, the coefficient of restitution; however, the magnitude of the stress as well as the range of shear thinning and thickening show significant dependence on the inelasticity.
We demonstrate that the time evolution of the van Hove dynamical pair correlation function is governed by adiabatic forces that arise from the free energy and by superadiabatic forces that are induced by the flow of the van Hove function. The superadiabatic forces consist of drag, viscous, and structural contributions, as occur in active Brownian particles, in liquids under shear and in lane forming mixtures. For hard sphere liquids we present a power functional theory that predicts these universal force fields in quantitative agreement with our Brownian dynamics simulation results.
We recently found that crystallization of monodisperse hard spheres from the bulk fluid faces a much higher free energy barrier in four than in three dimensions at equivalent supersaturation, due to the increased geometrical frustration between the simplex-based fluid order and the crystal [J.A. van Meel, D. Frenkel, and P. Charbonneau, Phys. Rev. E 79, 030201(R) (2009)]. Here, we analyze the microscopic contributions to the fluid-crystal interfacial free energy to understand how the barrier to crystallization changes with dimension. We find the barrier to grow with dimension and we identify the role of polydispersity in preventing crystal formation. The increased fluid stability allows us to study the jamming behavior in four, five, and six dimensions and compare our observations with two recent theories [C. Song, P. Wang, and H. A. Makse, Nature 453, 629 (2008); G. Parisi and F. Zamponi, Rev. Mod. Phys, in press (2009)].