Connectivity, dynamics, and structure in a tetrahedral network liquid

Structure and elasticity of model disordered, polydisperse, and defect-free polymer networks

The elasticity of disordered and polydisperse polymer networks is a fundamental problem of soft matter physics that is still open. Here, we self-assemble polymer networks via simulations of a mixture of bivalent and tri- or tetravalent patchy particles, which result in an exponential strand length distribution analogous to that of experimental randomly cross-linked systems. After assembly, the network connectivity and topology are frozen and the resulting system is characterized. We find that the fractal structure of the network depends on the number density at which the assembly has been carried out, but that systems with the same mean valence and same assembly density have the same structural properties. Moreover, we compute the long-time limit of the mean-squared displacement, also known as the (squared) localization length, of the cross-links and of the middle monomers of the strands, showing that the dynamics of long strands is well described by the tube model. Finally, we find a relation connecting these two localization lengths at high density and connect the cross-link localization length to the shear modulus of the system.

Extended Wertheim theory predicts the anomalous chain length distributions of divalent patchy particles under extreme confinement

Colloidal patchy particles with divalent attractive interaction can self-assemble into linear polymer chains. Their equilibrium properties in 2D and 3D are well described by Wertheim's thermodynamic perturbation theory which predicts a well-defined exponentially decaying equilibrium chain length distribution. In experi- mental realizations, due to gravity, particles sediment to the bottom of the suspension forming a monolayer of particles with a gravitational height smaller than the particle diameter. In accordance with experiments, an anomalously high monomer concentration is observed in simulations which is not well understood. To account for this observation, we interpret the polymerization as taking place in a highly confined quasi-2D plane and extend the Wertheim thermodynamic perturbation theory by defining addition reactions constants as functions of the chain length. We derive the theory, test it on simple square well potentials, and apply it to the experimental case of synthetic colloidal patchy particles immersed in a binary liquid mixture that are described by an accurate effective critical Casimir patchy particle potential. The important interaction parameters entering the theory are explicitly computed using the integral method in combination with Monte Carlo sampling. Without any adjustable parameter, the predictions of the chain length distribution are in excellent agreement with explicit simulations of self-assembling particles. We discuss generality of the approach, and its application range.

The physics of empty liquids: from patchy particles to water

Empty liquids represent a wide class of materials whose constituents arrange in a random network through reversible bonds. Many key insights on the physical properties of empty liquids have originated almost independently from the study of colloidal patchy particles on one side, and a large body of theoretical and experimental research on water on the other side. Patchy particles represent a family of coarse-grained potentials that allows for a precise control of both the geometric and the energetic aspects of bonding, while water has arguably the most complex phase diagram of any pure substance, and a puzzling amorphous phase behavior. It was only recently that the exchange of ideas from both fields has made it possible to solve long-standing problems and shed new light on the behavior of empty liquids. Here we highlight the connections between patchy particles and water, focusing on the modelling principles that make an empty liquid behave like water, including the factors that control the appearance of thermodynamic and dynamic anomalies, the possibility of liquid-liquid phase transitions, and the crystallization of open crystalline structures.

A temperature-dependent critical Casimir patchy particle model benchmarked onto experiment

Synthetic colloidal patchy particles immersed in a binary liquid mixture can self-assemble via critical Casimir interactions into various superstructures, such as chains and networks. Up to now, there are no quantitatively accurate potential models that can simulate and predict this experimentally observed behavior precisely. Here, we develop a protocol to establish such a model based on a combination of theoretical Casimir potentials and angular switching functions. Using Monte Carlo simulations, we optimize several material-specific parameters in the model to match the experimental chain length distribution and persistence length. Our approach gives a systematic way to obtain accurate potentials for critical Casimir induced patchy particle interactions and can be used in large-scale simulations.

From colloidal particles to photonic crystals: advances in self-assembly and their emerging applications

Over the last three decades, photonic crystals (PhCs) have attracted intense interests thanks to their broad potential applications in optics and photonics. Generally, these structures can be fabricated via either "top-down" lithographic or "bottom-up" self-assembly approaches. The self-assembly approaches have attracted particular attention due to their low cost, simple fabrication processes, relative convenience of scaling up, and the ease of creating complex structures with nanometer precision. The self-assembled colloidal crystals (CCs), which are good candidates for PhCs, have offered unprecedented opportunities for photonics, optics, optoelectronics, sensing, energy harvesting, environmental remediation, pigments, and many other applications. The creation of high-quality CCs and their mass fabrication over large areas are the critical limiting factors for real-world applications. This paper reviews the state-of-the-art techniques in the self-assembly of colloidal particles for the fabrication of large-area high-quality CCs and CCs with unique symmetries. The first part of this review summarizes the types of defects commonly encountered in the fabrication process and their effects on the optical properties of the resultant CCs. Next, the mechanisms of the formation of cracks/defects are discussed, and a range of versatile fabrication methods to create large-area crack/defect-free two-dimensional and three-dimensional CCs are described. Meanwhile, we also shed light on both the advantages and limitations of these advanced approaches developed to fabricate high-quality CCs. The self-assembly routes and achievements in the fabrication of CCs with the ability to open a complete photonic bandgap, such as cubic diamond and pyrochlore structure CCs, are discussed as well. Then emerging applications of large-area high-quality CCs and unique photonic structures enabled by the advanced self-assembly methods are illustrated. At the end of this review, we outlook the future approaches in the fabrication of perfect CCs and highlight their novel real-world applications.

Can oppositely charged polyelectrolyte stars form a gel? A simulational study

We present a Langevin molecular dynamics study of an equimolar mixture of monodispersed oppositely charged di-block four-armed polyelectrolyte stars. We used an implicit solvent coarse-grained representation of the polyelectrolyte stars, and varied the length of the terminal charged blocks that reside on each arm. By varying the polymer concentration, we computed PV diagrams and determined the free-swelling equilibrium concentration with respect to a pure water reservoir as a function of the charged block length. We investigated various structural properties of the resulting equilibrium structures, like the number of ionic bonds, dangling arms, isolated stars, and cluster sizes. The ionic bonds featured a broad distribution of the number of arms involved and also displayed a distribution of net charges peaked around the neutral ionic bond. Our main result is that for charged block length equal to 4 and 5 ionized beads the resulting macro-aggregate spans the box and forms a network phase. Furthermore, we investigated the restructuring dynamics of ionic bonds; the results suggested both short bond lifetimes and a high frequency of ballistic association/dissociation events. Bonds result strong enough to yield a stable gel phase, but they are still weak enough to allow network restructuring under thermal fluctuations.

Non-Gaussian subdiffusion of single-molecule tracers in a hydrated polymer network

Single molecule tracking experiments inside a hydrated polymer network have shown that the tracer motion is subdiffusive due to the viscoelastic environment inside the gel-like network. This property can be related to the negative autocorrelation of the instantaneous displacements at short times. Although the displacements of the individual tracers exhibit Gaussian statistics, the displacement distribution of all the trajectories combined from different spatial locations of the polymer network exhibits a non-Gaussian distribution. Here, we analyze many individual tracer trajectories to show that the central portion of the non-Gaussian distribution can be well approximated by an exponential distribution that spreads sublinearly with time. We explain all these features seen in the experiment by a generalized Langevin model for an overdamped particle with algebraically decaying correlations. We show that the degree of non-Gaussianity can change with the extent of heterogeneity, which is controlled in our model by the experimentally observed distributions of the motion parameters.

Understanding, predicting, and tuning the fragility of vitrimeric polymers

Fragility is an empirical property that describes how abruptly a glass-forming material solidifies upon supercooling. The degree of fragility carries important implications for the functionality and processability of a material, as well as for our fundamental understanding of the glass transition. However, the microstructural properties underlying fragility still remain poorly understood. Here, we explain the microstructure-fragility link in vitrimeric networks, a novel type of high-performance polymers with unique bond-swapping functionality and unusual glass-forming behavior. Our results are gained from coarse-grained computer simulations and first-principles mode-coupling theory (MCT) of star-polymer vitrimers. We first demonstrate that the vitrimer fragility can be tuned over an unprecedentedly broad range, from fragile to strong and even superstrong behavior, by decreasing the bulk density. Remarkably, this entire phenomenology can be reproduced by microscopic MCT, thus challenging the conventional belief that existing first-principles theories cannot account for nonfragile behaviors. Our MCT analysis allows us to rationally identify the microstructural origin of the fragile-to-superstrong crossover, which is rooted in the sensitivity of the static structure factor to temperature variations. On the molecular scale, this behavior stems from a change in dominant length scales, switching from repulsive excluded-volume interactions to intrachain attractions as the vitrimer density decreases. Finally, we develop a simplified schematic MCT model which corroborates our microscopically founded conclusions and which unites our findings with earlier MCT studies. Our work sheds additional light on the elusive structure-fragility link in glass-forming matter and provides a first-principles-based platform for designing amorphous materials with an on-demand dynamic response.

q-Independent Slow Dynamics in Atomic and Molecular Systems

Investigating million-atom systems for very long simulation times, we demonstrate that the collective density-density correlation time (τ_{α}) in simulated supercooled water and silica becomes wave-vector independent (q^{0}) when the probing wavelength is several times larger than the interparticle distance. The q independence of the collective density-density correlation functions, a feature clearly observed in light-scattering studies of some soft-matter systems, is thus a genuine feature of many (but not all) slow-dynamics systems, either atomic, molecular, or colloidal. Indeed, we show that when the dynamics of the density fluctuations includes particle-type diffusion, as in the case of the Lennard-Jones binary-mixture model, the q^{0} regime does not set in and the relaxation time continues to scale as τ_{α}∼q^{-2} even at small q.

Two time scales for self and collective diffusion near the critical point in a simple patchy model for proteins with floating bonds

Using dynamic Monte Carlo and Brownian dynamics, we investigate a floating bond model in which particles can bind through mobile bonds. The maximum number of bonds (here fixed to 4) can be tuned by appropriately choosing the repulsive, nonadditive interactions among bonds and particles. We compute the static and dynamic structure factor (intermediate scattering function) in the vicinity of the gas-liquid critical point. The static structure exhibits a weak tetrahedral network character. The intermediate scattering function shows a temporal decay deviating from a single exponential, which can be described by a double exponential decay where the two time scales differ approximately by one order of magnitude. This time scale separation is robust over a range of wave numbers. The analysis of clusters in real space indicates the formation of noncompact clusters and shows a considerable stretch in the instantaneous size distribution when approaching the critical point. The average time evolution of the largest subcluster of given initial clusters with 10 or more particles also shows a double exponential decay. The observation of two time scales in the intermediate scattering function at low packing fractions is consistent with similar findings in globular protein solutions with trivalent metal ions that act as bonds between proteins.

How to simulate patchy particles⋆

Patchy particles is the name given to a large class of systems of mesoscopic particles characterized by a repulsive core and a discrete number of short-range and highly directional interaction sites. Numerical simulations have contributed significantly to our understanding of the behaviour of patchy particles, but, although simple in principle, advanced simulation techniques are often required to sample the low temperatures and long time-scales associated with their self-assembly behaviour. In this work we review the most popular simulation techniques that have been used to study patchy particles, with a special focus on Monte Carlo methods. We cover many of the tools required to simulate patchy systems, from interaction potentials to biased moves, cluster moves, and free-energy methods. The review is complemented by an educationally oriented Monte Carlo computer code that implements all the techniques described in the text to simulate a well-known tetrahedral patchy particle model.

Dynamics of network fluids

Network fluids are structured fluids consisting of chains and branches. They are characterized by unusual physical properties, such as, exotic bulk phase diagrams, interfacial roughening and wetting transitions, and equilibrium and nonequilibrium gels. Here, we provide an overview of a selection of their equilibrium and dynamical properties. Recent research efforts towards bridging equilibrium and non-equilibrium studies are discussed, as well as several open questions.

Path-integral formalism for stochastic resetting: Exactly solved examples and shortcuts to confinement

We study the dynamics of overdamped Brownian particles diffusing in conservative force fields and undergoing stochastic resetting to a given location at a generic space-dependent rate of resetting. We present a systematic approach involving path integrals and elements of renewal theory that allows us to derive analytical expressions for a variety of statistics of the dynamics such as (i) the propagator prior to first reset, (ii) the distribution of the first-reset time, and (iii) the spatial distribution of the particle at long times. We apply our approach to several representative and hitherto unexplored examples of resetting dynamics. A particularly interesting example for which we find analytical expressions for the statistics of resetting is that of a Brownian particle trapped in a harmonic potential with a rate of resetting that depends on the instantaneous energy of the particle. We find that using energy-dependent resetting processes is more effective in achieving spatial confinement of Brownian particles on a faster time scale than performing quenches of parameters of the harmonic potential.