Source author record

Daniel de las Heras

Daniel de las Heras appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2020arXiv

Flow and structure in nonequilibrium Brownian many-body systems

We present a fundamental classification of forces relevant in nonequilibrium structure formation under collective flow in Brownian many-body systems. The internal one-body force field is systematically split into contributions relevant for the spatial structure and for the coupled motion. We demonstrate that both contributions can be obtained straightforwardly in computer simulations, and present a power functional theory that describes all types of forces quantitatively. Our conclusions and methods are relevant for flow in inertial systems, such as molecular liquids and granular media.

preprint2020arXiv

Non-negative Interfacial Tension in Phase-Separated Active Brownian Particles

We present a microscopic theory for the nonequilibrium interfacial tension $γ_{\rm gl}$ of the free interface between gas and liquid phases of active Brownian particles. The underlying square gradient treatment and the splitting of the force balance in flow and structural contributions is general and applies to inhomogeneous nonequilibrium steady states. We find $γ_{\rm gl}\geq 0$, which opposes claims by Bialké et al. [Phys. Rev. Lett. $\bf{115}$, 098301 (2015)] and delivers the theoretical justification for the widely observed interfacial stability in active Brownian dynamics many-body simulations.

preprint2016arXiv

Particle Conservation in Dynamical Density Functional Theory

We present the exact adiabatic theory for the dynamics of the inhomogeneous density distribution of a classical fluid. Erroneous particle number fluctuations of dynamical density functional theory are absent, both for canonical and grand canonical initial conditions. We obtain the canonical free energy functional, which yields the adiabatic interparticle forces of overdamped Brownian motion. Using an exact and one of the most advanced approximate hard core free energy functionals, we obtain excellent agreement with simulations. The theory applies to finite systems in and out of equilibrium.

preprint2016arXiv

Percolation in binary and ternary mixtures of patchy colloids

We investigate percolation in binary and ternary mixtures of patchy colloidal particles theoretically and using Monte Carlo simulations. Each particle has three identical patches, with distinct species having different types of patch. Theoretically we assume tree-like clusters and calculate the bonding probabilities using Wertheim's first-order perturbation theory for association. For ternary mixtures we find up to eight fundamentally different percolated states. The states differ in terms of the species and pairs of species that have percolated. The strongest gel is a trigel or tricontinuous gel, in which each of the three species has percolated. The weakest gel is a mixed gel in which all of the particles have percolated, but none of the species percolates by itself. The competition between entropy of mixing and internal energy of bonding determines the stability of each state. Theoretical and simulation results are in very good agreement. The only significant difference is the temperature at the percolation threshold, which is overestimated by the theory due to the absence of closed loops in the theoretical description.

preprint2016arXiv

Reentrant Network Formation in Patchy Colloidal Mixtures under Gravity

We study a two-dimensional binary mixture of patchy colloids in sedimentation-diffusion-equilibrium using Monte Carlo simulation and Wertheim's theory. By tuning the buoyant masses of the colloids we can control the gravity-induced sequence of fluid stacks of differing density and percolation properties. We find complex stacking sequences with up to four layers and reentrant network formation, consistently in simulations and theoretically using only the bulk phase diagram as input. Our theory applies to general patchy colloidal mixtures and is relevant to understanding experiments under gravity.

preprint2016arXiv

Temperature (de)activated patchy colloidal particles

We present a new model of patchy particles in which the interaction sites can be activated or deactivated by varying the temperature of the system. We study the thermodynamics of the system by means of Wertheim's first order perturbation theory, and use Flory-Stockmayer theory of polymerization to analyse the percolation threshold. We find a very rich phase behaviour including lower critical points and reentrant percolation.

preprint2016arXiv

The role of sample height in the stacking diagram of colloidal mixtures under gravity

Bulk phase separation is responsible for the occurrence of stacks of different layers in sedimentation of colloidal mixtures. A recently proposed theory (de las Heras and Schmidt 2013 Soft Matter 9 8636) establishes a unique connection between the bulk phase behaviour and sedimentation-diffusion-equilibrium. The theory constructs a stacking diagram of all possible sequences of stacks under gravity in the limit of very high (infinite) sample heights. Here, we study the stacking diagrams of colloidal mixtures at finite sample height, h. We demonstrate that h plays a vital role in sedimentation-diffusion-equilibrium of colloidal mixtures. The region of the stacking diagram occupied by a given sequence of stacks depends on h. Hence, two samples with different heights but identical colloidal concentrations can develop different stacking sequences. In addition, the stacking diagrams for different heights can be qualitatively different since some stacking sequences occur only in a given interval of sample heights. We use the theory to investigate the stacking diagrams of both model bulk systems and mixtures of patchy particles that differ either by the number or by the types of patches.

preprint2016arXiv

Topological protection of multiparticle dissipative transport

Topological protection allows robust transport of localized phenomena such as quantum information, solitons, and dislocations. The transport can be either dissipative or non-dissipative. Here, we experimentally demonstrate and theoretically explain the topologically protected dissipative motion of colloidal particles above a periodic hexagonal magnetic pattern. By driving the system with periodic modulation loops of an external and spatially homogeneous magnetic field, we achieve total control over the motion of diamagnetic and paramagnetic colloids. We can transport simultaneously and independently each type of colloid along any of the six crystallographic directions of the pattern via adiabatic or deterministic ratchet motion. Both types of motion are topologically protected. As an application, we implement an automatic topologically protected quality control of a chemical reaction between functionalized colloids. Our results are relevant to other systems with the same symmetry.

preprint2016arXiv

Topologically protected colloidal transport above a square magnetic lattice

We theoretically study the motion of magnetic colloidal particles above a magnetic pattern and compare the predictions with Brownian dynamics simulations. The pattern consists of alternating square domains of positive and negative magnetization. The colloidal motion is driven by periodic modulation loops of an external magnetic field. There exist loops that induce topologically protected colloidal transport between two different unit cells of the pattern. The transport is very robust against internal and external perturbations. Theory and simulations are in perfect agreement.

preprint2015arXiv

Capillary and winding transitions in a confined cholesteric liquid crystal

We consider a Lebwohl-Lasher model of chiral particles confined in a planar cell (slit pore) with different boundary conditions, and solve it using mean-field theory. The phase behaviour of the system with respect to temperature and pore width is studied. Two phenomena are observed: (i) an isotropic-cholesteric transition which exhibits an oscillatory structure with respect to pore width, and (ii) an infinite set of winding transitions caused by commensuration effects between cholesteric pitch and pore width. The latter transitions have been predicted and analysed by other authors for cholesterics confined in a fixed pore and subject to an external field promoting the uniaxial nematic phase; here we induce winding transitions solely from geometry by changing the pore width at zero external field (a setup recently explored in Atomic-Force Microscopy experiments). In contrast with previous studies, we obtain the phase diagrams in the temperature vs pore width plane, including the isotropic-cholesteric transition, the winding transitions and their complex relationship. In particular, the structure of winding transitions terminates at the capillary isotropic-cholesteric transition via triple points where the confined isotropic phase coexists with two cholesterics with different helix indices. For symmetric and asymmetric monostable plate anchorings the phase diagram are qualitatively similar.

preprint2015arXiv

Ordering in granular rod monolayers driven far from thermodynamic equilibrium

The orientational order in vertically agitated granular rod monolayers is investigated experimentally and compared quantitatively with equilibrium Monte Carlo simulations and density functional theory. At sufficiently high number density, short rods form a tetratic state and long rods form a uniaxial nematic state. The length-to-width ratio at which the order changes from tetratic to uniaxial is around $7.3$ in both experiments and simulations. This agreement illustrates the universal aspects of the ordering of rod-shaped particles across equilibrium and nonequilibrium systems. Moreover, the assembly of granular rods into ordered states is found to be independent of the agitation frequency and strength, suggesting that the detailed nature of energy injection into such a nonequilibrium system does not play a crucial role.

preprint2014arXiv

Domain walls in two-dimensional nematics confined in a small circular cavity

Using Monte Carlo simulation, we study a fluid of two-dimensional hard rods inside a small circular cavity bounded by a hard wall, from the dilute regime to the high-density, layering regime. Both planar and homeotropic anchoring of the nematic director can be induced at the walls through a free-energy penalty. The circular geometry creates frustration in the nematic phase and a polar-symmetry configuration with a distorted director field plus two $+1/2$ disclinations is created. At higher densities, a quasi-uniform structure is observed with a (minimal) director distortion which is relaxed via the formation of orientational domain walls. This novel structure is not predicted by elasticity theory and is similar to the step-like structures observed in three-dimensional hybrid slit pores. We speculate that the formation of domain walls is a general mechanism to relax elastic stresses in conditions of strong surface anchoring and severe spatial confinement.

preprint2014arXiv

Sedimentation stacking diagram of binary colloidal mixtures and bulk phases in the plane of chemical potentials

We give a full account of a recently proposed theory that explicitly relates the bulk phase diagram of a binary colloidal mixture to its phase stacking phenomenology under gravity [Soft Matter 9, 8636 (2013)]. As we demonstrate, the full set of possible phase stacking sequences in sedimentation-diffusion equilibrium originates from straight lines (sedimentation paths) in the chemical potential representation of the bulk phase diagram. From the analysis of various standard topologies of bulk phase diagrams, we conclude that the corresponding sedimentation stacking diagrams can be very rich, even more so when finite sample height is taken into account. We apply the theory to obtain the stacking diagram of a mixture of nonadsorbing polymers and colloids. We also present a catalog of generic phase diagrams in the plane of chemical potentials in order to facilitate the practical application of our concept, which also generalizes to multi component mixtures.

preprint2014arXiv

Superadiabatic forces in Brownian many-body dynamics

Theoretical approaches to nonequilibrium many-body dynamics generally rest upon an adiabatic assumption, whereby the true dynamics is represented as a sequence of equilibrium states. Going beyond this simple approximation is a notoriously difficult problem. For the case of classical Brownian many-body dynamics we present a simulation method that allows to isolate and precisely evaluate superadiabatic correlations and the resulting forces. Application of the method to a system of one-dimensional hard particles reveals the importance for the dynamics, as well as the complexity, of these nontrivial out-of-equilibrium contributions. Our findings help clarify the status of dynamical density functional theory and provide a rational basis for the development of improved theories.

preprint2013arXiv

Computing the phase diagram of binary mixtures: a patchy particle case study

We investigate the phase behaviour of 2D mixtures of bi-functional and three-functional patchy particles and 3D mixtures of bi-functional and tetra-functional patchy particles by means of Monte Carlo simulations and Wertheim theory. We start by computing the critical points of the pure systems and then we investigate how the critical parameters change upon lowering the temperature. We extend the Successive Umbrella Sampling method to mixtures to make it possible to extract information about the phase behaviour of the system at a fixed temperature for the whole range of densities and compositions of interest.

preprint2013arXiv

Phase stacking diagram of colloidal mixtures under gravity

The observation of stacks of distinct layers in a colloidal or liquid mixture in sedimentation-diffusion equilibrium is a striking consequence of bulk phase separation. Drawing quantitative conclusions about the phase diagram is, however, very delicate. Here we introduce the Legendre transform of the chemical potential representation of the bulk phase diagram to obtain a unique stacking diagram of all possible stacks under gravity. Simple bulk phase diagrams generically lead to complex stacking diagrams. We apply the theory to a binary hard core platelet mixture with only two-phase bulk coexistence, and find that the stacking diagram contains six types of stacks with up to four distinct layers. These results can be tested experimentally in colloidal platelet mixtures. In general, an extended Gibbs phase rule determines the maximum number of sedimented layers to be $3+2(n_b-1)+n_i$, where $n_b$ is the number of binodals and $n_i$ is the number of their inflection points.

preprint2012arXiv

Bulk fluid phase behaviour of colloidal platelet-sphere and platelet-polymer mixtures

Using a geometry-based fundamental measure density functional theory, we calculate bulk fluid phase diagrams of colloidal mixtures of vanishingly thin hard circular platelets and hard spheres. We find isotropic-nematic phase separation with strong broadening of the biphasic region upon increasing the pressure. In mixtures with large size ratio of platelet and sphere diameter, there is also demixing between two nematic phases with differing platelet concentrations. We formulate a fundamental-measure density functional for mixtures of colloidal platelets and freely overlapping spheres, which represent ideal polymers, and use it to obtain phase diagrams. We find that for low platelet-polymer size ratio in addition to isotropic-nematic and nematic-nematic phase coexistence, platelet-polymer mixtures also display isotropic-isotropic demixing. In contrast, we do not find isotropic-isotropic demixing in hard core platelet-sphere mixtures for the size ratios considered.

preprint2012arXiv

Floating nematic phase in colloidal platelet-sphere mixtures

The phase behaviour of colloidal dispersions is interesting for fundamental reasons and for technological applications such as photonic crystals and electronic paper. Sedimentation, which in everyday life is relevant from blood analysis to the shelf life of paint, is a means to determine phase boundaries by observing distinct layers in samples that are in sedimentation-diffusion equilibrium. However, disentangling the effects due to interparticle interactions, which generate the bulk phase diagram, from those due to gravity is a complex task. Here we show that a line in the space of chemical potentials $μ_i$, where $i$ labels the species, represents a sedimented sample and that each crossing of this sedimentation path with a binodal generates an interface under gravity. Complex phase stacks can result, such as the sandwich of a floating nematic layer between top and bottom isotropic phases that we observed in a mixture of silica spheres and gibbsite platelets.

preprint2012arXiv

Properties of patchy colloidal particles close to a surface: a Monte Carlo and density functional study

We investigate the behavior of a patchy particle model close to a hard-wall via Monte Carlo simulation and density functional theory (DFT). Two DFT approaches, based on the homogeneous and inhomogeneous versions of Wertheim's first order perturbation theory for the association free energy are used. We evaluate, by simulation and theory, the equilibrium bulk phase diagram of the fluid and analyze the surface properties for two isochores, one of which is close to the liquid side of the gas-liquid coexistence curve. We find that the density profile near the wall crosses over from a typical high-temperature adsorption profile to a low-temperature desorption one, for the isochore close to coexistence. We relate this behavior to the properties of the bulk network liquid and find that the theoretical descriptions are reasonably accurate in this regime. At very low temperatures, however, an almost fully bonded network is formed, and the simulations reveal a second adsorption regime which is not captured by DFT. We trace this failure to the neglect of orientational correlations of the particles, which are found to exhibit surface induced orientational order in this regime.

preprint2012arXiv

Two-dimensional nematics in bulk and confined geometries

Two-dimensional nematics possess peculiar properties that have been studied recently using computer simulation and various theoretical models. Here we review our own contribution to the field using density-functional theory, and present some preliminary simulation results on confined two-dimensional nematics. First we discuss the possible stable bulk phases and phase diagrams and the relation between phases and particle geometry. We then explore the adsorption properties on a single substrate and the confinement effects that arise when the fluid is confined between parallel walls. Next, confinement in circular cavities is presented; this geometry allows us to measure some properties of the simplest defects that arise in two-dimensional nematics. Finally, preliminary Monte Carlo simulation results of confined nematics in circular geometry are shown.

preprint2011arXiv

Bicontinuous and mixed gels in binary mixtures of patchy colloidal particles

We investigate the thermodynamics and percolation regimes of model binary mixtures of patchy colloidal particles. The particles of each species have three sites of two types, one of which promotes bonding of particles of the same species while the other promotes bonding of different species. We find up to four percolated structures at low temperatures and densities: two gels where only one species percolates, a mixed gel where particles of both species percolate but neither species percolates separately, and a bicontinuous gel where particles of both species percolate separately forming two interconnected networks. The competition between the entropy and the energy of bonding drives the stability of the different percolating structures. Appropriate mixtures exhibit one or more connectivity transitions between the mixed and bicontinuous gels, as the temperature and/or the composition changes.

preprint2011arXiv

Phase diagrams of binary mixtures of patchy colloids with distinct numbers and types of patches: The empty fluid regime

We investigate the effect of distinct bonding energies on the onset of criticality of low functionality fluid mixtures. We focus on mixtures of particles with two and three patches as this includes the mixture where 'empty' fluids were originally reported. In addition to the number of patches, the species differ in the type of patches or bonding sites. For simplicity, we consider that the patches on each species are identical: one species has 3 patches of type $A$ and the other 2 patches of type $B$. We have found a rich phase behaviour with closed miscibility gaps, liquid-liquid demixing and negative azeotropes. Liquid-liquid demixing was found to preempt the 'empty' fluid regime, of these mixtures, when the $AB$ bonds are weaker than the $AA$ or $BB$ bonds. By contrast, mixtures in this class exhibit 'empty' fluid behaviour when the $AB$ bonds are stronger than at least one of the other two. Mixtures with bonding energies $ε_{BB}=ε_{AB}$ and $ε_{AA}<ε_{BB}$, were found to exhibit an unusual negative azeotrope.

preprint2011arXiv

Phase diagrams of binary mixtures of patchy colloids with distinct numbers of patches: The network fluid regime

We calculate the network fluid regime and phase diagrams of binary mixtures of patchy colloids, using Wertheim's first order perturbation theory and a generalization of Flory-Stockmayer's theory of polymerization. The colloids are modelled as hard spheres with the same diameter and surface patches of the same type, $A$. The only difference between species is the number of their patches -or functionality-, $f_A^{(1)}$ and $f_A^{(2)}$ (with $f_A^{(2)}>f_A^{(1)}$). We have found that the difference in functionality is the key factor controlling the behaviour of the mixture in the network (percolated) fluid regime. In particular, when $f_A^{(2)}\ge2f_A^{(1)}$ the entropy of bonding drives the phase separation of two network fluids which is absent in other mixtures. This changes drastically the critical properties of the system and drives a change in the topology of the phase diagram (from type I to type V) when $f_A^{(1)}>2$. The difference in functionality also determines the miscibility at high (osmotic) pressures. If $f_A^{(2)}-f_A^{(1)}=1$ the mixture is completely miscible at high pressures, while closed miscibility gaps at pressures above the highest critical pressure of the pure fluids are present if $f_A^{(2)}-f_A^{(1)}>1$. We argue that this phase behaviour is driven by a competition between the entropy of mixing and the entropy of bonding, as the latter dominates in the network fluid regime.

preprint2010arXiv

Surface and smectic layering transitions in binary mixtures of parallel hard rods

The surface phase behavior of binary mixtures of colloidal hard rods in contact with a solid substrate (hard wall) is studied, with special emphasis on the region of the phase diagram that includes the smectic A phase. The colloidal rods are modelled as hard cylinders of the same diameter and different lengths, in the approximation of perfect alignment. A fundamental--measure density functional is used to obtain equilibrium density profiles and thermodynamic properties such as surface tensions and adsorption coefficients. The bulk phase diagram exhibits nematic-smectic and smectic-smectic demixing, with smectic phases having different compositions; in some cases they are microfractionated. The calculated surface phase diagram of the wall-nematic interface shows a very rich phase behavior, including layering transitions and complete wetting at high pressures, whereby an infinitely thick smectic film grows at the wall via an infinite sequence of stepwise first--order layering transitions. For lower pressures complete wetting also obtains, but here the smectic film grows in a continuous fashion. Finally, at very low pressures, the wall-nematic interface exhibits critical adsorption by the smectic phase, due to the second-order character of the bulk nematic-smectic transition.