No Arabic abstract
We compute the fourth virial coefficient of a binary nonadditive hard-sphere mixture over a wide range of deviations from diameter additivity and size ratios. Hinging on this knowledge, we build up a $y$ expansion [B. Barboy and W. N. Gelbart, J. Chem. Phys. {bf 71}, 3053 (1979)] in order to trace the fluid-fluid coexistence lines which we then compare with the available Gibbs-ensemble Monte Carlo data and with the estimates obtained through two refined integral-equation theories of the fluid state. We find that in a regime of moderately negative nonadditivity and largely asymmetric diameters, relevant to the modelling of sterically and electrostatically stabilized colloidal mixtures, the fluid-fluid critical point is unstable with respect to crystallization.
Using scaled-particle theory for binary mixtures of two-dimensional hard particles with rotational freedom, we analyse the stability of nematic phases and the demixing phase behaviour of a variety of mixtures, focussing on cases where at least one of the components consists of hard rectangles or hard squares. A pure fluid of hard rectangles may exhibit, aside from the usual uniaxial nematic phase, an additional (tetratic) oriented phase, possessing two directors, which is the analogue of the biaxial or cubatic phases in three- dimensional fluids. There is computer simulation evidence that the tetratic phase might be stable with respect to phases with spatial order for rectangles with low aspect ratios. As hard rectangles are mixed with other particles not possessing stable tetratic order by themselves, the tetratic phase is destabilised, via a first- or second-order phase transition, to uniaxial nematic or isotropic phases; for hard rectangles of low aspect ratio tetratic order persists in a relatively large range of volume fractions. The order of these transitions depends on the particle geometry, dimensions and thermodynamic conditions of the mixture. The second component of the mixture has been chosen to be hard discs or disco-rectangles, the geometry of which is different from that of rectangles, leading to packing frustration and demixing behaviour, or simply rectangles of different aspect ratio. These mixtures may be good candidates for observing thermodynamically stable tetratic phases in monolayers of hard particles. Finally, demixing between fluid (isotropic--tetratic or tetratic--tetratic) phases is seen to occur in mixtures of hard squares of different sizes when the size ratio is sufficiently large.
We study a binary mixture of polar chiral (counterclockwise or clockwise) active particles in a two-dimensional box with periodic boundary conditions. Beside the excluded volume interactions between particles, particles are also subject to the polar velocity alignment. From the extensive Brownian dynamics simulations, it is found that the particle configuration (mixing or demixing) is determined by the competition between the chirality difference and the polar velocity alignment. When the chirality difference competes with the polar velocity alignment, the clockwise particles aggregate in one cluster and the counterclockwise particles aggregate in the other cluster, thus particles are demixed and can be separated. However, when the chirality difference or the polar velocity alignment is dominated, particles are mixed. Our findings could be used for the experimental pursuit of the separation of binary mixtures of chiral active particles.
Obtaining a rigorous and reliable method for linking computer simulations of polymer blends and composites at different length scales of interest is a highly desirable goal in soft matter physics. In this paper a multiscale modeling procedure is presented for the efficient calculation of the static structural properties of binary homopolymer blends. The procedure combines computer simulations of polymer chains on two different length scales, using a united atom representation for the finer structure and a highly coarse-grained approach on the meso-scale, where chains are represented as soft colloidal particles interacting through an effective potential. A method for combining the structural information by inverse mapping is discussed, allowing for the efficient calculation of partial correlation functions, which are compared with results from full united atom simulations. The structure of several polymer mixtures is obtained in an efficient manner for several mixtures in the homogeneous region of the phase diagram. The method is then extended to incorporate thermal fluctuations through an effective chi parameter. Since the approach is analytical, it is fully transferable to numerous systems.
Using density functionals from fundamental measure theory, phase diagrams and crystal-fluid surface tensions in additive and nonadditive (Asakura-Oosawa model) two-dimensional hard disk mixtures are determined for the whole range of size ratios $q$ between disks, assuming random disorder in the crystal phase. The fluid-crystal transitions are first-order due to the assumption of a periodic unit cell in the density functional calculations. Qualitatively, the shape of the phase diagrams is similar to the case of three-dimensional hard sphere mixtures. For the nonadditive case, a broadening of the fluid-crystal coexistence region is found for small $q$ whereas for higher $q$ a vapor--fluid transition intervenes. In the additive case, we find a sequence of spindle type, azeotropic and eutectic phase diagrams upon lowering $q$ from 1 to 0.6. The transition from azeotropic to eutectic is different from the three-dimensional case. Surface tensions in general become smaller (up to a factor 2) upon addition of a second species and they are rather small. The minimization of the functionals proceeds without restrictions and optimized graphics card routines are used.
We study the thermodynamic stability of fluid-fluid phase separation in binary nonadditive mixtures of hard-spheres for moderate size ratios. We are interested in elucidating the role played by small amounts of nonadditivity in determining the stability of fluid-fluid phase separation with respect to the fluid-solid phase transition. The demixing curves are built in the framework of the modified-hypernetted chain and of the Rogers-Young integral equation theories through the calculation of the Gibbs free energy. We also evaluate fluid-fluid phase equilibria within a first-order thermodynamic perturbation theory applied to an effective one-component potential obtained by integrating out the degrees of freedom of the small spheres. A qualitative agreement emerges between the two different approaches. We also address the determination of the freezing line by applying the first-order thermodynamic perturbation theory to the effective interaction between large spheres. Our results suggest that for intermediate size ratios a modest amount of nonadditivity, smaller than earlier thought, can be sufficient to drive the fluid-fluid critical point into the thermodinamically stable region of the phase diagram. These findings could be significant for rare-gas mixtures in extreme pressure and temperature conditions, where nonadditivity is expected to be rather small.