Simulating critical dynamics in liquid mixtures: Short-range and long-range contributions

Hydrodynamic effects in kinetics of phase separation in binary fluids: Critical versus off-critical compositions

Via hydrodynamics-preserving molecular dynamics simulations we study growth phenomena in a phase-separating symmetric binary mixture model. We quench high-temperature homogeneous configurations to state points inside the miscibility gap, for various mixture compositions. For compositions at the symmetric or critical value we capture the rapid linear viscous hydrodynamic growth due to advective transport of material through tubelike interconnected domains. For state points very close to any of the branches of the coexistence curve, the growth in the system, following nucleation of disconnected droplets of the minority species, occurs via a coalescence mechanism. Using state-of-the-art techniques, we have identified that these droplets, between collisions, exhibit diffusive motion. The value of the exponent for the power-law growth, related to this diffusive coalescence mechanism, has been estimated. While the exponent nicely agrees with that for the growth via the well-known Lifshitz-Slyozov particle diffusion mechanism, the amplitude is stronger. For the intermediate compositions we observe initial rapid growth that matches the expectations for viscous or inertial hydrodynamic pictures. However, at later times these types of growth cross over to the exponent that is decided by the diffusive coalescence mechanism.

Active particles in explicit solvent: Dynamics of clustering for alignment interaction

We study the dynamics of clustering in systems containing active particles that are immersed in an explicit solvent. For this, we have adopted a hybrid simulation method, consisting of molecular dynamics and multiparticle collision dynamics. In our model, the overlap-avoiding passive interaction of an active particle with another active particle or a solvent particle has been taken care of via variants of the Lennard-Jones potential. Dynamic interactions among the active particles have been incorporated via a Vicsek-like alignment rule in self-propulsion that facilitates clustering. We quantify the effects of activity and importance of hydrodynamics on the dynamics of clustering via variations of relevant system parameters. Results are obtained for low overall density of active particles, for which the state point is close to the vapor branch of the coexistence curve, and thus the morphology consists of disconnected clusters. In such a situation, the mechanism of growth switches among particle diffusion, diffusive coalescence, and ballistic aggregation, depending upon the presence or absence of active and hydrodynamic interactions providing different kinds of mobilities to the clusters. Corresponding growth laws have been quantified and discussed in the context of appropriate theoretical pictures. Our results suggest that multiparticle collision dynamics is an effective method for the investigation of hydrodynamic phenomena in phase-separating active matter systems.

How do clusters in phase-separating active matter systems grow? A study for Vicsek activity in systems undergoing vapor-solid transition

Via molecular dynamics simulations, we have studied the kinetics of vapor-"solid" phase transition in an active matter model in which self-propulsion is introduced via the well-known Vicsek rule. The overall density of the particles is chosen in such a way that the evolution morphology consists of disconnected clusters that are defined as regions of high density of particles. Our focus has been on understanding the influence of the above-mentioned self-propulsion on structure and growth of these clusters by comparing the results with those for the passive limit of the model that also exhibits vapor-"solid" transition. While in the passive case growth occurs due to a standard diffusive mechanism, the Vicsek activity leads to very rapid growth, via a process that is practically equivalent to the ballistic aggregation mechanism. The emerging growth law in the latter case has been accurately estimated and explained by invoking information on velocity and structural aspects of the clusters into a relevant theory. Some of these results are also discussed with reference to a model for active Brownian particles.

Critical Casimir interactions between colloids around the critical point of binary solvents

Critical Casimir interactions between colloidal particles arise from the confinement of fluctuations of a near-critical solvent in the liquid gap between closely-spaced particles. So far, the comparison of theoretical predictions and experimental measurements of critical Casimir forces (CCFs) has focused on the critical solvent composition, while it has been lacking for off-critical compositions. We address this issue by investigating CCFs between spherical colloidal particles around the critical point of a binary solvent through a combination of experiments, previous Ising Monte Carlo simulation results and field-theoretical methods. By measuring the correlation length of the near-critical solvent and the pair potentials of the particles in terms of radial distribution functions and by determining the second virial coefficient, we test in detail theoretical predictions. Our results indicate that the critical Casimir theory gives quantitative correct predictions for the interaction potential between particles in a near critical binary mixture if weak preferential adsorption of the particle surface is taken into account.

The Effect of Alkali Halides on the Critical Exponents of the 2,6-Dimethylpyridine-Water System

The diffusion coefficients and shear viscosities of 2,6-dimethylpyridine-water mixtures of critical composition have been measured without and with small amounts of alkali halides added. The data have been analyzed in terms of power law behavior. Deviations from power law behavior indicate a coupling between the critical fluctuations in the local concentrations and the formation of mesoscopic molecular aggregates. The critical exponent of the fluctuation correlation length, the shear viscosity exponent, and the critical exponent of the relaxation rate of fluctuations have been evaluated to show noticeable influences from the salts. The correlation length exponent indicates a suppression of the critical fluctuations, whereas the viscosity exponent rather points at the activation of some extra fluctuations. No clear evidence has been obtained that the added salts affect the critical behavior, and thus cause the opposed effects in the exponents, directly by the long-range ionic fields. Alternatively, the ions may have an influence on the aggregate formation which in turn could modify the critical fluctuations either by reducing the region of true critical exponents or by affecting the critical behavior due to the presence as well as the formation and disintegration kinetics of the multimolecular structures.

Finite-size scaling study of shear viscosity anomaly at liquid-liquid criticality

We study the equilibrium dynamics of a symmetrical binary Lennard-Jones fluid mixture near its consolute criticality. Molecular dynamics simulation results for the shear viscosity, η, from a microcanonical ensemble are compared with those from a canonical ensemble with various thermostats. It is observed that the Nosé-Hoover thermostat is a good candidate for this purpose, and is therefore adopted for the quantification of the critical singularity of η, to avoid the temperature fluctuations (or even drifts) that are often encountered in microcanonical simulations. Via a finite-size scaling analysis of our simulation data we have been able to confirm that the shear viscosity exhibits a weak critical singularity in agreement with the theoretical predictions.