No Arabic abstract
We report on a large scale computer simulation study of crystal nucleation in hard spheres. Through a combined analysis of real and reciprocal space data, a picture of a two-step crystallization process is supported: First dense, amorphous clusters form which then act as precursors for the nucleation of well-ordered crystallites. This kind of crystallization process has been previously observed in systems that interact via potentials that have an attractive as well as a repulsive part, most prominently in protein solutions. In this context the effect has been attributed to the presence of metastable fluid-fluid demixing. Our simulations, however, show that a purely repulsive system (that has no metastable fluid-fluid coexistence) crystallizes via the same mechanism.
We use numerical simulations to study the crystallization of monodisperse systems of hard aspherical particles. We find that particle shape and crystallizability can be easily related to each other when particles are characterized in terms of two simple and experimentally accessible order parameters: one based on the particle surface-to-volume ratio, and the other on the angular distribution of the perturbations away from the ideal spherical shape. We present a phase diagram obtained by exploring the crystallizability of 487 different particle shapes across the two-order-parameter spectrum. Finally, we consider the physical properties of the crystalline structures accessible to aspherical particles, and discuss limits and relevance of our 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)].
The smallest maximum kissing-number Voronoi polyhedron of 3d spheres is the icosahedron and the tetrahedron is the smallest volume that can show up in Delaunay tessalation. No periodic lattice is consistent with either and hence these dense packings are geometrically frustrated. Because icosahedra can be assembled from almost perfect tetrahedra, the terms icosahedral and polytetrahedral packing are often used interchangeably, which leaves the true origin of geometric frustration unclear. Here we report a computational study of freezing of 4d hard spheres, where the densest Voronoi cluster is compatible with the symmetry of the densest crystal, while polytetrahedral order is not. We observe that, under otherwise comparable conditions, crystal nucleation in 4d is less facile than in 3d. This suggest that it is the geometrical frustration of polytetrahedral structures that inhibits crystallization.
We report the emergence of a giant Mpemba effect in the uniformly heated gas of inelastic rough hard spheres: The initially hotter sample may cool sooner than the colder one, even when the initial temperatures differ by more than one order of magnitude. In order to understand this behavior, it suffices to consider the simplest Maxwellian approximation for the velocity distribution in a kinetic approach. The largeness of the effect stems from the fact that the rotational and translational temperatures, which obey two coupled evolution equations, are comparable. Our theoretical predictions agree very well with molecular dynamics and direct simulation Monte Carlo data.
Using super-heterodyne Doppler velocimetry with multiple scattering correction, we extend the opti-cally accessible range of concentrations in experiments on colloidal electro-kinetics. We here meas-ured the electro-phoretic mobility and the DC conductivity of aqueous charged sphere suspensions covering about three orders of magnitude in particle concentrations and transmissions as low as 40%. The extended concentration range for the first time allows the demonstration of a non-monotonic con-centration dependence of the mobility for a single particle species. Our observations reconcile previ-ous experimental observations made on other species over restricted concentration ranges. We com-pare our results to state of the art theoretical calculations using a constant particle charge and the carefully determined experimental boundary conditions as input. In particular, we consider so-called realistic salt free conditions, i.e. we respect the release of counter-ions by the particles, the solvent hydrolysis and the formation of carbonic acid from dissolved neutral CO2. We also compare to previ-ous results obtained under similarly well-defined conditions. This allows identification of three dis-tinct regions of differing density dependence. An ascent during the built up of double layer overlap which is not expected by theory, an extended plateau region in quantitative agreement with theoretical expectation based on a constant effective charge and a sudden decrease which occurs way before the expected gradual decrease. Our observations suggest a relation of the non-monotonic behavior to a decrease of particle charge, and we tentatively discuss possibly underlying mechanisms.