Static and Dynamic Critical Behavior of a Symmetrical Binary Fluid: A Computer Simulation

A symmetrical binary, A+B Lennard-Jones mixture is studied by a combination of semi-grandcanonical Monte Carlo (SGMC) and Molecular Dynamics (MD) methods near a liquid-liquid critical temperature $T_c$. Choosing equal chemical potentials for the two species, the SGMC switches identities (${\rm A} \to {\rm B} \to {\rm A}$) to generate well-equilibrated configurations of the system on the coexistence curve for $T<T_c$ and at the critical concentration, $x_c=1/2$, for $T>T_c$. A finite-size scaling analysis of the concentration susceptibility above $T_c$ and of the order parameter below $T_c$ is performed, varying the number of particles from N=400 to 12800. The data are fully compatible with the expected critical exponents of the three-dimensional Ising universality class. The equilibrium configurations from the SGMC runs are used as initial states for microcanonical MD runs, from which transport coefficients are extracted. Self-diffusion coefficients are obtained from the Einstein relation, while the interdiffusion coefficient and the shear viscosity are estimated from Green-Kubo expressions. As expected, the self-diffusion constant does not display a detectable critical anomaly. With appropriate finite-size scaling analysis, we show that the simulation data for the shear viscosity and the mutual diffusion constant are quite consistent both with the theoretically predicted behavior, including the critical exponents and amplitudes, and with the most accurate experimental evidence.Comment: 35 pages, 13 figure

Double Binds and Double Blinds: Evaluation Tactics in Critically Oriented HCI

Critically oriented researchers within Human-Computer Interaction (HCI) have fruitfully intersected design and critical analysis to engage users and designers in reflection on underlying values, assumptions and dominant practices in technology. To successfully integrate this work within the HCI community, critically oriented researchers have tactically engaged with dominant practices within HCI in the design and evaluation of their work. This paper draws attention to the ways that tactical engagement with aspects of HCI evaluation methodology shapes and bears consequences for critically oriented research. We reflect on three of our own experiences evaluating critically oriented designs and trace challenges that we faced to the ways that sensibilities about generalizable knowledge are manifested in HCI evaluation methodology. Drawing from our own experiences, as well as other influential critically oriented design projects in HCI, we articulate some of the trade-offs involved in consciously adopting or not adopting certain normative aspects of HCI evaluation. We argue that some forms of this engagement can hamstring researchers from pursuing their intended research goals and have consequences beyond specific research projects to affect the normative discourse in the field as a whole

Phase diagrams of Janus fluids with up-down constrained orientations

A class of binary mixtures of Janus fluids formed by colloidal spheres with the hydrophobic hemispheres constrained to point either up or down are studied by means of Gibbs ensemble Monte Carlo simulations and simple analytical approximations. These fluids can be experimentally realized by the application of an external static electrical field. The gas-liquid and demixing phase transitions in five specific models with different patch-patch affinities are analyzed. It is found that a gas-liquid transition is present in all the models, even if only one of the four possible patch-patch interactions is attractive. Moreover, provided the attraction between like particles is stronger than between unlike particles, the system demixes into two subsystems with different composition at sufficiently low temperatures and high densities.Comment: 10 pages, 6 figure

Molecular Dynamics Simulation of Heat-Conducting Near-Critical Fluids

Using molecular dynamics simulations, we study supercritical fluids near the gas-liquid critical point under heat flow in two dimensions. We calculate the steady-state temperature and density profiles. The resultant thermal conductivity exhibits critical singularity in agreement with the mode-coupling theory in two dimensions. We also calculate distributions of the momentum and heat fluxes at fixed density. They indicate that liquid-like (entropy-poor) clusters move toward the warmer boundary and gas-like (entropy-rich) regions move toward the cooler boundary in a temperature gradient. This counterflow results in critical enhancement of the thermal conductivity

Critical Dynamics in a Binary Fluid: Simulations and Finite-size Scaling

We report comprehensive simulations of the critical dynamics of a symmetric binary Lennard-Jones mixture near its consolute point. The self-diffusion coefficient exhibits no detectable anomaly. The data for the shear viscosity and the mutual-diffusion coefficient are fully consistent with the asymptotic power laws and amplitudes predicted by renormalization-group and mode-coupling theories {\it provided} finite-size effects and the background contribution to the relevant Onsager coefficient are suitably accounted for. This resolves a controversy raised by recent molecular simulations.Comment: 4 pages, 4 figure

Histogram Reweighting Method for Dynamic Properties

The histogram reweighting technique, widely used to analyze Monte Carlo data, is shown to be applicable to dynamic properties obtained from Molecular Dynamics simulations. The theory presented here is based on the fact that the correlation functions in systems in thermodynamic equilibrium are averages over initial conditions of functions of the trajectory of the system in phase-space, the latter depending on the volume, the total number of particles and the classical Hamiltonian. Thus, the well-known histogram reweighting method can almost straightforwardly be applied to reconstruct the probability distribution of initial states at different thermodynamic conditions, without extra computational effort. Correlation functions and transport coefficients are obtained with this method from few simulation data sets.Comment: 4 pages, 3 figure

Effects of patch size and number within a simple model of patchy colloids

We report on a computer simulation and integral equation study of a simple model of patchy spheres, each of whose surfaces is decorated with two opposite attractive caps, as a function of the fraction $\chi$ of covered attractive surface. The simple model explored --- the two-patch Kern-Frenkel model --- interpolates between a square-well and a hard-sphere potential on changing the coverage $\chi$. We show that integral equation theory provides quantitative predictions in the entire explored region of temperatures and densities from the square-well limit $\chi = 1.0$ down to $\chi \approx 0.6$. For smaller $\chi$, good numerical convergence of the equations is achieved only at temperatures larger than the gas-liquid critical point, where however integral equation theory provides a complete description of the angular dependence. These results are contrasted with those for the one-patch case. We investigate the remaining region of coverage via numerical simulation and show how the gas-liquid critical point moves to smaller densities and temperatures on decreasing $\chi$. Below $\chi \approx 0.3$, crystallization prevents the possibility of observing the evolution of the line of critical points, providing the angular analog of the disappearance of the liquid as an equilibrium phase on decreasing the range for spherical potentials. Finally, we show that the stable ordered phase evolves on decreasing $\chi$ from a three-dimensional crystal of interconnected planes to a two-dimensional independent-planes structure to a one-dimensional fluid of chains when the one-bond-per-patch limit is eventually reached.Comment: 26 pages, 11 figures, J. Chem. Phys. in pres

Critical dynamics of an isothermal compressible non-ideal fluid

A pure fluid at its critical point shows a dramatic slow-down in its dynamics, due to a divergence of the order-parameter susceptibility and the coefficient of heat transport. Under isothermal conditions, however, sound waves provide the only possible relaxation mechanism for order-parameter fluctuations. Here we study the critical dynamics of an isothermal, compressible non-ideal fluid via scaling arguments and computer simulations of the corresponding fluctuating hydrodynamics equations. We show that, below a critical dimension of 4, the order-parameter dynamics of an isothermal fluid effectively reduces to "model A," characterized by overdamped sound waves and a divergent bulk viscosity. In contrast, the shear viscosity remains finite above two dimensions. Possible applications of the model are discussed.Comment: 19 pages, 7 figures; v3: minor corrections and clarifications; as published in Phys. Rev.