Inferring bulk self-assembly properties from simulations of small systems with multiple constituent species and small systems in the grand canonical ensemble

In this paper we generalize a methodology [T. E. Ouldridge, A. A. Louis, and J. P. K. Doye, J. Phys.: Condens. Matter {\bf 22}, 104102 (2010)] for dealing with the inference of bulk properties from small simulations of self-assembling systems of characteristic finite size. In particular, schemes for extrapolating the results of simulations of a single self-assembling object to the bulk limit are established in three cases: for assembly involving multiple particle species, for systems with one species localized in space and for simulations in the grand canonical ensemble. Furthermore, methodologies are introduced for evaluating the accuracy of these extrapolations. Example systems demonstrate that differences in cluster concentrations between simulations of a single self-assembling structure and bulk studies of the same model under identical conditions can be large, and that convergence on bulk results as system size is increased can be slow and non-trivial.Comment: Accepted by J. Chem. Phy

Langlands duality for representations of quantum groups

We establish a correspondence (or duality) between the characters and the crystal bases of finite-dimensional representations of quantum groups associated to Langlands dual semi-simple Lie algebras. This duality may also be stated purely in terms of semi-simple Lie algebras. To explain this duality, we introduce an "interpolating quantum group" depending on two parameters which interpolates between a quantum group and its Langlands dual. We construct examples of its representations, depending on two parameters, which interpolate between representations of two Langlands dual quantum groups.Comment: 37 pages. References added. Accepted for publication in Mathematische Annale

The role of long-range forces in the phase behavior of colloids and proteins

The phase behavior of colloid-polymer mixtures, and of solutions of globular proteins, is often interpreted in terms of a simple model of hard spheres with short-ranged attraction. While such a model yields a qualitative understanding of the generic phase diagrams of both colloids and proteins, it fails to capture one important difference: the model predicts fluid-fluid phase separation in the metastable regime below the freezing curve. Such demixing has been observed for globular proteins, but for colloids it appears to be pre-empted by the appearance of a gel. In this paper, we study the effect of additional long-range attractions on the phase behavior of spheres with short-ranged attraction. We find that such attractions can shift the (metastable) fluid-fluid critical point out of the gel region. As this metastable critical point may be important for crystal nucleation, our results suggest that long-ranged attractive forces may play an important role in the crystallization of globular proteins. However, in colloids, where refractive index matching is often used to switch off long-ranged dispersion forces, gelation is likely to inhibit phase separation.Comment: EURO-LATEX, 6 pages, 2 figure

Core-Softened System With Attraction: Trajectory Dependence of Anomalous Behavior

In the present article we carry out a molecular dynamics study of the core-softened system and show that the existence of the water-like anomalies in this system depends on the trajectory in $P-\rho-T$ space along which the behavior of the system is studied. For example, diffusion and structural anomalies are visible along isotherms as a function of density, but disappears along the isochores and isobars as a function of temperature. On the other hand, the diffusion anomaly may be seen along adiabats as a function of temperature, density and pressure. It should be noted that it may be no signature of a particular anomaly along a particular trajectory, but the anomalous region for that particular anomaly can be defined when all possible trajectories in the same space are examined (for example, signature of diffusion anomaly is evident through the crossing of different isochors. However, there is no signature of diffusion anomaly along a particular isochor). We also analyze the applicability of the Rosenfeld entropy scaling relations to this system in the regions with the water-like anomalies. It is shown that the validity of the Rosenfeld scaling relation for the diffusion coefficient also depends on the trajectory in the $P-\rho-T$ space along which the kinetic coefficients and the excess entropy are calculated.Comment: 16 pages, 21 figures. arXiv admin note: this contains much of the content of arXiv:1010.416

Inversion of Sequence of Diffusion and Density Anomalies in Core-Softened Systems

In this paper we present a simulation study of water-like anomalies in core-softened system introduced in our previous publications. We investigate the anomalous regions for a system with the same functional form of the potential but with different parameters and show that the order of the region of anomalous diffusion and the region of density anomaly is inverted with increasing the width of the repulsive shoulder.Comment: 8 pages, 10 figure

Quasi-binary amorphous phase in a 3D system of particles with repulsive-shoulder interactions

We report a computer-simulation study of the equilibrium phase diagram of a three-dimensional system of particles with a repulsive step potential. Using free-energy calculations, we have determined the equilibrium phase diagram of this system. At low temperatures, we observe a number of distinct crystal phases. However, under certain conditions the system undergoes a glass transition in a regime where the liquid appears thermodynamically stable. We argue that the appearance of this amorphous low-temperature phase can be understood by viewing this one-component system as a pseudo-binary mixture.Comment: 4 pages, 4 figure

Phase diagram of softly repulsive systems: The Gaussian and inverse-power-law potentials

We redraw, using state-of-the-art methods for free-energy calculations, the phase diagrams of two reference models for the liquid state: the Gaussian and inverse-power-law repulsive potentials. Notwithstanding the different behavior of the two potentials for vanishing interparticle distances, their thermodynamic properties are similar in a range of densities and temperatures, being ruled by the competition between the body-centered-cubic (BCC) and face-centered-cubic (FCC) crystalline structures and the fluid phase. We confirm the existence of a reentrant BCC phase in the phase diagram of the Gaussian-core model, just above the triple point. We also trace the BCC-FCC coexistence line of the inverse-power-law model as a function of the power exponent $n$ and relate the common features in the phase diagrams of such systems to the softness degree of the interaction.Comment: 22 pages, 8 figure

The structure of the hard sphere solid

We show that near densest-packing the perturbations of the HCP structure yield higher entropy than perturbations of any other densest packing. The difference between the various structures shows up in the correlations between motions of nearest neighbors. In the HCP structure random motion of each sphere impinges slightly less on the motion of its nearest neighbors than in the other structures.Comment: For related papers see: http://www.ma.utexas.edu/users/radin/papers.htm

A Generalization of Metropolis and Heat-Bath Sampling for Monte Carlo Simulations

For a wide class of applications of the Monte Carlo method, we describe a general sampling methodology that is guaranteed to converge to a specified equilibrium distribution function. The method is distinct from that of Metropolis in that it is sometimes possible to arrange for unconditional acceptance of trial moves. It involves sampling states in a local region of phase space with probability equal to, in the first approximation, the square root of the desired global probability density function. The validity of this choice is derived from the Chapman-Kolmogorov equation, and the utility of the method is illustrated by a prototypical numerical experiment.Comment: RevTeX, 7 pages, 2 table