Mode-coupling theory predictions for a limited valency attractive square-well model

Recently we have studied, using numerical simulations, a limited valency model, i.e. an attractive square well model with a constraint on the maximum number of bonded neighbors. Studying a large region of temperatures $T$ and packing fractions $\phi$, we have estimated the location of the liquid-gas phase separation spinodal and the loci of dynamic arrest, where the system is trapped in a disordered non-ergodic state. Two distinct arrest lines for the system are present in the system: a {\it (repulsive) glass} line at high packing fraction, and a {\it gel} line at low $\phi$ and $T$. The former is essentially vertical ($\phi$-controlled), while the latter is rather horizontal ($T$-controlled) in the $(\phi-T)$ plane. We here complement the molecular dynamics results with mode coupling theory calculations, using the numerical structure factors as input. We find that the theory predicts a repulsive glass line -- in satisfactory agreement with the simulation results -- and an attractive glass line which appears to be unrelated to the gel line.Comment: 12 pages, 6 figures. To appear in J. Phys. Condens. Matter, special issue: "Topics in Application of Scattering Methods for Investigation of Structure and Dynamics of Soft Condensed Matter", Fiesole, November 200

A spherical model with directional interactions: I. Static properties

We introduce a simple spherical model whose structural properties are similar to the ones generated by models with directional interactions, by employing a binary mixture of large and small hard spheres, with a square-well attraction acting only between particles of different size. The small particles provide the bonds between the large ones. With a proper choice of the interaction parameters, as well as of the relative concentration of the two species, it is possible to control the effective valence. Here we focus on a specific choice of the parameters which favors tetrahedral ordering and study the equilibrium static properties of the system in a large window of densities and temperatures. Upon lowering the temperature we observe a progressive increase in local order, accompanied by the formation of a four-coordinated network of bonds. Three different density regions are observed: at low density the system phase separates into a gas and a liquid phase; at intermediate densities a network of fully bonded particles develops; at high densities -- due to the competition between excluded volume and attractive interactions -- the system forms a defective network. The very same behavior has been previously observed in numerical studies of non-spherical models for molecular liquids, such as water, and in models of patchy colloidal particles. Differently from these models, theoretical treatments devised for spherical potentials, e.g. integral equations and ideal mode coupling theory for the glass transition can be applied in the present case, opening the way for a deeper understanding of the thermodynamic and dynamic behavior of low valence molecules and particles.Comment: 11 pages, 11 figure

Energy landscape of a simple model for strong liquids

We calculate the statistical properties of the energy landscape of a minimal model for strong network-forming liquids. Dynamics and thermodynamic properties of this model can be computed with arbitrary precision even at low temperatures. A degenerate disordered ground state and logarithmic statistics for the energy distribution are the landscape signatures of strong liquid behavior. Differences from fragile liquid properties are attributed to the presence of a discrete energy scale, provided by the particle bonds, and to the intrinsic degeneracy of topologically disordered networks.Comment: Revised versio

Non-Gaussian energy landscape of a simple model for strong network-forming liquids: accurate evaluation of the configurational entropy

We present a numerical study of the statistical properties of the potential energy landscape of a simple model for strong network-forming liquids. The model is a system of spherical particles interacting through a square well potential, with an additional constraint that limits the maximum number of bonds, $N_{\rm max}$, per particle. Extensive simulations have been carried out as a function of temperature, packing fraction, and $N_{\rm max}$. The dynamics of this model are characterized by Arrhenius temperature dependence of the transport coefficients and by nearly exponential relaxation of dynamic correlators, i.e. features defining strong glass-forming liquids. This model has two important features: (i) landscape basins can be associated with bonding patterns; (ii) the configurational volume of the basin can be evaluated in a formally exact way, and numerically with arbitrary precision. These features allow us to evaluate the number of different topologies the bonding pattern can adopt. We find that the number of fully bonded configurations, i.e. configurations in which all particles are bonded to $N_{\rm max}$ neighbors, is extensive, suggesting that the configurational entropy of the low temperature fluid is finite. We also evaluate the energy dependence of the configurational entropy close to the fully bonded state, and show that it follows a logarithmic functional form, differently from the quadratic dependence characterizing fragile liquids. We suggest that the presence of a discrete energy scale, provided by the particle bonds, and the intrinsic degeneracy of fully bonded disordered networks differentiates strong from fragile behavior.Comment: Final version. Journal of Chemical Physics 124, 204509 (2006

Effect of bond lifetime on the dynamics of a short-range attractive colloidal system

We perform molecular dynamics simulations of short-range attractive colloid particles modeled by a narrow (3% of the hard sphere diameter) square well potential of unit depth. We compare the dynamics of systems with the same thermodynamics but different bond lifetimes, by adding to the square well potential a thin barrier at the edge of the attractive well. For permanent bonds, the relaxation time $\tau$ diverges as the packing fraction $\phi$ approaches a threshold related to percolation, while for short-lived bonds, the $\phi$-dependence of $\tau$ is more typical of a glassy system. At intermediate bond lifetimes, the $\phi$-dependence of $\tau$ is driven by percolation at low $\phi$, but then crosses over to glassy behavior at higher $\phi$. We also study the wavevector dependence of the percolation dynamics.Comment: Revised; 9 pages, 9 figure

Confirmation of Anomalous Dynamical Arrest in attractive colloids: a molecular dynamics study

Previous theoretical, along with early simulation and experimental, studies have indicated that particles with a short-ranged attraction exhibit a range of new dynamical arrest phenomena. These include very pronounced reentrance in the dynamical arrest curve, a logarithmic singularity in the density correlation functions, and the existence of `attractive' and `repulsive' glasses. Here we carry out extensive molecular dynamics calculations on dense systems interacting via a square-well potential. This is one of the simplest systems with the required properties, and may be regarded as canonical for interpreting the phase diagram, and now also the dynamical arrest. We confirm the theoretical predictions for re-entrance, logarithmic singularity, and give the first direct evidence of the coexistence, independent of theory, of the two coexisting glasses. We now regard the previous predictions of these phenomena as having been established.Comment: 15 pages,15 figures; submitted to Phys. Rev.

Effect of Chain Polydispersity on the Elasticity of Disordered Polymer Networks

Due to their unique structural and mechanical properties, randomly cross-linked polymer networks play an important role in many different fields, ranging from cellular biology to industrial processes. In order to elucidate how these properties are controlled by the physical details of the network (e.g., chain-length and end-to-end distributions), we generate disordered phantom networks with different cross-linker concentrations C and initial densities ρinit and evaluate their elastic properties. We find that the shear modulus computed at the same strand concentration for networks with the same C, which determines the number of chains and the chain-length distribution, depends strongly on the preparation protocol of the network, here controlled by ρinit. We rationalize this dependence by employing a generic stress-strain relation for polymer networks that does not rely on the specific form of the polymer end-to-end distance distribution. We find that the shear modulus of the networks is a nonmonotonic function of the density of elastically active strands, and that this behavior has a purely entropic origin. Our results show that if short chains are abundant, as it is always the case for randomly cross-linked polymer networks, the knowledge of the exact chain conformation distribution is essential for correctly predicting the elastic properties. Finally, we apply our theoretical approach to literature experimental data, qualitatively confirming our interpretations

Multiple glass transitions in star polymer mixtures: Insights from theory and simulations

The glass transition in binary mixtures of star polymers is studied by mode coupling theory and extensive molecular dynamics computer simulations. In particular, we have explored vitrification in the parameter space of size asymmetry $\delta$ and concentration $\rho_2$ of the small star polymers at fixed concentration of the large ones. Depending on the choice of parameters, three different glassy states are identified: a single glass of big polymers at low $\delta$ and low $\rho_2$, a double glass at high $\delta$ and low $\rho_2$, and a novel double glass at high $\rho_2$ and high $\delta$ which is characterized by a strong localization of the small particles. At low $\delta$ and high $\rho_2$ there is a competition between vitrification and phase separation. Centered in the $(\delta, \rho_2)$-plane, a liquid lake shows up revealing reentrant glass formation. We compare the behavior of the dynamical density correlators with the predictions of the theory and find remarkable agreement between the two.Comment: 15 figures, to be published in Macromolecule

Static and Dynamic Anomalies in a Repulsive Spherical Ramp Liquid: Theory and Simulation

We compare theoretical and simulation results for static and dynamic properties for a model of particles interacting via a spherically symmetric repulsive ramp potential. The model displays anomalies similar to those found in liquid water, namely, expansion upon cooling and an increase of diffusivity upon compression. In particular, we calculate the phase diagram from the simulation and successfully compare it with the phase diagram obtained using the Rogers-Young (RY) closure for the Ornstein-Zernike equation. Both the theoretical and the numerical calculations confirm the presence of a line of isobaric density maxima, and lines of compressibility minima and maxima. Indirect evidence of a liquid-liquid critical point is found. Dynamic properties also show anomalies. Along constant temperature paths, as the density increases, the dynamics alternates several times between slowing down and speeding up, and we associate this behavior with the progressive structuring and de-structuring of the liquid. Finally we confirm that mode coupling theory successfully predicts the non-monotonic behavior of dynamics and the presence of multiple glass phases, providing strong evidence that structure (the only input of mode coupling theory) controls dynamics.Comment: Static and Dynamic Anomalies in a Repulsive Spherical Ramp Liquid: Theory and Simulatio

Glasslike Arrest in Spinodal Decomposition as a Route to Colloidal Gelation

Colloid-polymer mixtures can undergo spinodal decomposition into colloid-rich and colloid-poor regions. Gelation results when interconnected colloid-rich regions solidify. We show that this occurs when these regions undergo a glass transition, leading to dynamic arrest of the spinodal decomposition. The characteristic length scale of the gel decreases with increasing quench depth, and the nonergodicity parameter exhibits a pronounced dependence on scattering vector. Mode coupling theory gives a good description of the dynamics, provided we use the full static structure as input.Comment: 14 pages, 4 figures; replaced with published versio