Dynamical arrest transition of a bidisperse two-patchy colloidal dispersion: A dynamic Monte Carlo study.

During the last few years, patchy colloidal dispersions have emerged as ideal candidates of glass-formers of systems composed of particles that interact with non-isotropic potentials. However, from the computational point of view, the characterization of their dynamical properties close to the glass transition via any kind of molecular dynamics simulation technique can be very difficult due to the slowing down of both the rotational and translational dynamics. Although a plethora of dynamical techniques have been developed to account for the dynamics of patchy colloids, new and complementary simulation techniques are required to explore, much faster and more efficiently, the dynamical arrest transition of patchy colloidal dispersions when computer simulation consists of a large number of particles and, due to the slow particle dynamics at the glass transition, an extended time window is explicitly required.

Modeling the structure and thermodynamics of multicomponent and polydisperse hard-sphere dispersions with continuous potentials.

A model system of identical particles interacting via a hard-sphere potential is essential in condensed matter physics; it helps to understand in and out of equilibrium phenomena in complex fluids, such as colloidal dispersions. Yet, most of the fixed time-step algorithms to study the transport properties of those systems have drawbacks due to the mathematical nature of the interparticle potential. Because of this, mapping a hard-sphere potential onto a soft potential has been recently proposed [Báez et al.

Patchy colloidal gels under the influence of gravity.

In this contribution, gravitational effects in gel-forming patchy colloidal systems are studied. We focus on how the gel structure is modified by gravity. Through Monte Carlo computer simulations of gel-like states recently identified by the rigidity percolation criterion [J.

Phase separation and dynamical arrest of protein solutions dominated by short-range attractions.

The interplay of liquid-liquid phase separation (LLPS) and dynamical arrest can lead to the formation of gels and glasses, which is relevant for such diverse fields as condensed matter physics, materials science, food engineering, and the pharmaceutical industry. In this context, protein solutions exhibit remarkable equilibrium and non-equilibrium behaviors. In the regime where attractive and repulsive forces compete, it has been demonstrated, for example, that the location of the dynamical arrest line seems to be independent of ionic strength, so that the arrest lines at different ionic screening lengths overlap, in contrast to the LLPS coexistence curves, which strongly depend on the salt concentration.

Determining depletion interactions by contracting forces.

Depletion forces are fundamental for determining the phase behavior of a vast number of materials and colloidal dispersions and have been used for the manipulation of in- and out-of-equilibrium thermodynamic states. The entropic nature of depletion forces is well understood; however, most theoretical approaches, and also molecular simulations, work quantitatively at moderate size ratios in much diluted systems since large size asymmetries and high particle concentrations are difficult to deal with. The existing approaches for integrating out the degrees of freedom of the depletant species may fail under these extreme physical conditions.

Clusters in colloidal dispersions with a short-range depletion attraction: Thermodynamic identification and morphology.

Hypothesis: Particle aggregation is ubiquitous for many colloidal systems, and drives the phase separation or the formation of materials with a highly heterogeneous large-scale structure, such as gels, porous media and attractive glasses. While the macroscopic properties of such materials strongly depend on the shape and size of these particle aggregates, the morphology and underlining aggregation physical mechanisms are far from being fully understood. Recently, it has been proposed that for reversible colloidal aggregation, the cluster morphology in the case of colloids interacting with short-range attractive forces is determined by a single variable, namely, the reduced second virial coefficient, B.

Soft representation of the square-well and square-shoulder potentials to be used in Brownian and molecular dynamics simulations.

The discrete hard-sphere (HS), square-well (SW), and square-shoulder (SS) potentials have become the battle horse of molecular and complex fluids because they contain the basic elements to describe the thermodynamic, structural, and transport properties of both types of fluids. The mathematical simplicity of these discrete potentials allows us to obtain some analytical results despite the nature and complexity of the modeled systems. However, the divergent forces arising at the potential discontinuities may lead to severe issues when discrete potentials are used in computer simulations with uniform time steps.

Location of the gel-like boundary in patchy colloidal dispersions: Rigidity percolation, structure, and particle dynamics.

During the past decade, there has been a hot debate about the physical mechanisms that determine when a colloidal dispersion approaches the gel transition. However, there is still no consensus on a possible unique route that leads to the conditions for the formation of a gel-like state. Based on gel states identified in experiments, Valadez-Pérez et al.

Thermodynamics, static properties and transport behaviour of fluids with competing interactions.

Competing interaction fluids have become ideal model systems to study a large number of phenomena, for example, the formation of intermediate range order structures, condensed phases not seen in fluids driven by purely attractive or repulsive forces, the onset of particle aggregation under in- and out-of-equilibrium conditions, which results in the birth of reversible and irreversible aggregates or clusters whose topology and morphology depend additionally on the thermodynamic constrictions, and a particle dynamics that has a strong influence on the transport behaviour and rheological properties of the fluid. In this contribution, we study a system of particles interacting through a potential composed by a continuous succession of a short-ranged square-well (SW), an intermediate-ranged square-shoulder and a long-ranged SW. This potential model is chosen to systematically analyse the contribution of every component of the interaction potential on the phase behaviour, the microstructure, the morphology of the resulting aggregates and the transport phenomena of fluids described by competing interactions.

Fluorescent organic particle doped polymer-based gel dosimeter for neutron detection.

The purpose of this work is to develop a material capable of detecting neutrons produced by photodisintegration in a linear accelerator for its medical use. In this study, we have developed a gel-like material doped with fluorescent organic particles. PPO at 1 wt% is used as primary dopant and POPOP as secondary one at 0.

Potential-invariant network structures in Asakura-Oosawa mixtures with very short attraction range.

We systematically investigated the structure and aggregate morphology of gel networks formed by colloid-polymer mixtures with a moderate colloid volume fraction and different values of the polymer-colloid size ratio, always in the limit of short-range attraction. Using the coordinates obtained from confocal microscopy experiments, we determined the radial, angular, and nearest-neighbor distribution functions together with the cluster radius of gyration as a function of size ratio and polymer concentration. The analysis of the structural correlations reveals that the network structure becomes increasingly less sensitive to the potential strength with the decreasing polymer-colloid size ratio.

Competing interactions in the depletion forces of ternary colloidal mixtures.

Depletion interactions between colloidal particles surrounded by smaller depletants are typically characterized by a strong attraction at contact and a moderately repulsive barrier in front of it that extends at distances similar to the size of the depletants; the appearance and height of the barrier basically depend on the concentration and, therefore, the correlation between depletants. From a thermodynamic point of view, the former can drive the system to phase separation or toward non-equilibrium states, such as gel-like states, but its effects on both local and global properties may be controlled by the latter, which acts as a kind of entropic gate. However, the latter has not been entirely analyzed and understood within the context of colloidal mixtures mainly driven by entropy.

Generalized equation of state for fluids: From molecular liquids to colloidal dispersions.

In this work, a new parameterization for the Statistical Association Fluid Theory for potentials of Variable Range (SAFT-VR) is coupled to the discrete potential theory to represent the thermodynamic properties of several fluids, ranging from molecular liquids to colloidal-like dispersions. In this way, this version of the SAFT-VR approach can be straightforwardly applied to any kind of either simple or complex fluid. In particular, two interaction potentials, namely, the Lennard-Jones and the hard-core attractive Yukawa potentials, are discretized to study the vapor-liquid equilibrium properties of both molecular and complex liquids, respectively.

Concentration and size effects on the size-selective particle purification method using the critical Casimir force.

Critical Casimir force (CCF) is a solvent fluctuation introduced interaction between particles dispersed in a binary solvent. Recently, it has been demonstrated that the CCF induced attraction between particles can trigger particle size-sensitive aggregation, and has thus been used as an efficient way to purify nanoparticles by size. Here, combining small angle neutron scattering and dynamic light scattering, we investigate the effects of size and concentration on this particle size separation method.

Assessment of the Wolf method using the Stillinger-Lovett sum rules: From strong electrolytes to weakly charged colloidal dispersions.

The Ewald method has been the cornerstone in molecular simulations for modeling electrostatic interactions of charge-stabilized many-body systems. In the late 1990s, Wolf and collaborators developed an alternative route to describe the long-range nature of electrostatic interactions; from a computational perspective, this method provides a more efficient and straightforward way to implement long-range electrostatic interactions than the Ewald method. Despite these advantages, the validity of the Wolf potential to account for the electrostatic contribution in charged fluids remains controversial.

Machine learning for condensed matter physics.

Condensed matter physics (CMP) seeks to understand the microscopic interactions of matter at the quantum and atomistic levels, and describes how these interactions result in both mesoscopic and macroscopic properties. CMP overlaps with many other important branches of science, such as chemistry, materials science, statistical physics, and high-performance computing. With the advancements in modern machine learning (ML) technology, a keen interest in applying these algorithms to further CMP research has created a compelling new area of research at the intersection of both fields.

Spherical harmonic projections of the static structure factor of the dipolar hard sphere model: Theory vs simulations.

We investigate the static correlations of a dipolar fluid in terms of the irreducible coefficients of the spherical harmonic expansion of the static structure factor. To this end, we develop a theoretical framework based on a soft-core version of Wertheim's solution of the mean spherical approximation (MSA), which renders the analytical determination of such coefficients possible. The accuracy of this approximation is tested by a comparison against the results obtained with the assistance of extensive molecular dynamics simulations at different regimes of concentration and temperature.

Arrested dynamics of the dipolar hard sphere model.

We report the combined results of molecular dynamics simulations and theoretical calculations concerning various dynamical arrest transitions in a model system representing a dipolar fluid, namely, N (soft core) rigid spheres interacting through a truncated dipole-dipole potential. By exploring different regimes of concentration and temperature, we find three distinct scenarios for the slowing down of the dynamics of the translational and orientational degrees of freedom: at low (η = 0.2) and intermediate (η = 0.

Glassy dynamics in asymmetric binary mixtures of hard spheres.

We perform a systematic and detailed study of the glass transition in highly asymmetric binary mixtures of colloidal hard spheres, combining differential dynamic microscopy experiments, event-driven molecular dynamics simulations, and theoretical calculations, exploring the whole state diagram and determining the self-dynamics and collective dynamics of both species. Two distinct glassy states involving different dynamical arrest transitions are consistently described, namely, a double glass with the simultaneous arrest of the self-dynamics and collective dynamics of both species, and a single glass of large particles in which the self-dynamics of the small species remains ergodic. In the single-glass scenario, spatial modulations in the collective dynamics of both species occur due to the structure of the large spheres, a feature not observed in the double-glass domain.

Using the second virial coefficient as physical criterion to map the hard-sphere potential onto a continuous potential.

In the same sense as in the extended law of corresponding states [M. Noro and D. Frenkel, J.

Reversible Aggregation and Colloidal Cluster Morphology: The Importance of the Extended Law of Corresponding States.

Cluster morphology of spherical particles interacting with a short-range attraction has been extensively studied due to its relevance to many applications, such as the large-scale structure in amorphous materials, phase separation, protein aggregation, and organelle formation in cells. Although it was widely accepted that the range of the attraction solely controls the fractal dimension of clusters, recent experimental results challenged this concept by also showing the importance of the strength of attraction. Using Monte Carlo simulations, we conclusively demonstrate that it is possible to reduce the dependence of the cluster morphology to a single variable, namely, the reduced second virial coefficient, B_{2}^{*}, linking the local properties of colloidal systems to the extended law of corresponding states.

Friction and diffusion of a nano-colloidal disk in a two-dimensional solvent with a liquid-liquid transition.

We report on the friction and diffusion of a single mobile nano-colloidal disk, whose size and mass are one and two orders of magnitude, respectively, greater than the molecules of the host solvent; all particles are restricted to move in a two-dimensional space. Using molecular dynamics simulations, the variation of the transport coefficients as a function of the thermodynamic state of the supporting fluid, in particular, around those states in the neighbourhood of the liquid-liquid phase coexistence, is investigated. The diffusion coefficient is determined through the fit of the mean-square displacement at long times and with the Green-Kubo relationship for the velocity autocorrelation function, whereas the friction coefficient is computed from the correlation of the fluctuating force.

Counterion accumulation effects on a suspension of DNA molecules: Equation of state and pressure-driven denaturation.

The study of the effects associated with the electrostatic properties of DNA is of fundamental importance to understand both its molecular properties at the single molecule level, like the rigidity of the chain, and its interaction with other charged bio-molecules, including other DNA molecules; such interactions are crucial to maintain the thermodynamic stability of the intra-cellular medium. In the present work, we combine the Poisson-Boltzmann mean-field theory with an irreversible thermodynamic approximation to analyze the effects of counterion accumulation inside DNA on both the denaturation profile of the chain and the equation of state of the suspension. To this end, we model the DNA molecule as a porous charged cylinder immersed in an aqueous solution.

Single file dynamics in soft materials.

The term single file (SF) dynamics refers to the motion of an assembly of particles through a channel with cross-sections comparable to the particles' diameter. Single file diffusion (SFD) is then the diffusion of a tagged particle in a single file, i.e.

Structure of colloidal gels at intermediate concentrations: the role of competing interactions.

Colloidal gels formed by colloid-polymer mixtures with an intermediate volume fraction (ϕ ≈ 0.4) are investigated by confocal microscopy. In addition, we have performed Monte Carlo simulations based on a simple effective pair potential that includes a short-range attractive contribution representing depletion interactions, and a longer-range repulsive contribution describing the electrostatic interactions due to the presence of residual charges.