Tuning the liquid-liquid transition by modulating the hydrogen bond angular flexibility in a model for water

We propose a simple extension of the well known ST2 model for water [F.H. Stillinger and A. Rahman, J. Chem. Phys. {\bf 60}, 1545 (1974)] that allows for a continuous modification of the hydrogen bond angular flexibility. We show that the bond flexibility affects the relative thermodynamic stability of the liquid and of the hexagonal (or cubic) ice. On increasing flexibility, the liquid-liquid critical point, which in the original ST2 model is located in the no-man's land (i. e. the region where ice is the thermodynamically stable phase) progressively moves to a temperature where the liquid is more stable than ice. Our study definitively proves that the liquid-liquid transition in ST2 is a genuine phenomenon, of high relevance in all tetrahedral network-forming liquids, including water.Comment: Accepted in Phys. Rev. Let

Phase diagram of the ST2 model of water

We evaluate the free energy of the fluid and crystal phases for the ST2 potential [F.H. Stillinger and A. Rahman, J. Chem. Phys. 60, 1545 (1974)] with reaction field corrections for the long-range interactions. We estimate the phase coexistence boundaries in the temperature-pressure plane, as well as the gas-liquid critical point and gas-liquid coexistence conditions. Our study frames the location of the previously identified liquid-liquid critical point relative to the crystalline phase boundaries, and opens the way for exploring crystal nucleation in a model where the metastable liquid-liquid critical point is computationally accessible

Understanding tetrahedral liquids through patchy colloids

We investigate the structural properties of a simple model for tetrahedral patchy colloids in which the patch width and the patch range can be tuned independently. For wide bond angles, a fully bonded network can be generated by standard Monte Carlo or molecular dynamics simulations of the model, providing a neat method for generating defect-free random tetrahedral networks. This offers the possibility of focusing on the role of the patch angular width on the structure of the fully bonded network. The analysis of the fully bonded configurations as a function of the bonding angle shows how the bonding angle controls the system compressibility, the strength of the pre-peak in the structure factor and ring size distribution. Comparison with models of liquid water and silica allows us to find the best mapping between these continuous potentials and the colloidal one. Building on previous studies focused on the connection between angular range and crystallization, the mapping makes it possible to shed new light on the glass-forming ability of network-forming tetrahedral liquids.Comment: 10 pages, 6 figure

Gelling by Heating

We introduce a simple model, a binary mixture of patchy particles, which has been designed to form a gel upon heating. Due to the specific nature of the particle interactions, notably the number and geometry of the patches as well as their interaction energies, the system is a fluid both at high and at low temperatures, whereas at intermediate temperatures the system forms a solid-like disordered open network structure, i.e. a gel. Using molecular dynamics we investigate the static and dynamic properties of this system

Modeling of many-body interactions between elastic spheres through symmetry functions

Simple models for spherical particles with a soft shell have been shown to self-assemble into numerous crystal phases and even quasicrystals. However, most of these models rely on a simple pairwise interaction, which is usually a valid approximation only in the limit of small deformations, i.e. low densities. In this work, we consider a many-body yet simple model for the evaluation of the elastic energy associated with the deformation of a spherical shell. The resulting energy evaluation, however, is relatively expensive for direct use in simulations. We significantly reduce the associated numerical cost by fitting the potential using a set of symmetry functions. We propose a method for selecting a suitable set of symmetry functions that capture the most relevant features of the particle environment in a systematic manner. The fitted interaction potential is then used in Monte Carlo simulations to draw the phase diagram of the system in two dimensions. The system is found to form both a fluid and a hexagonal crystal phase.Comment: 10 pages, 9 figure