BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Marjolein Dijkstra

Marjolein Dijkstra appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softcond-mat.mtrl-scicond-mat.stat-mechcond-mat.mes-hallphysics.comp-phBiological Physicsphysics.chem-phmath-phmath.MPnlin.AOphysics.soc-ph

Catalog footprint

What is connected

51works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

51 published item(s)

preprint2022arXiv

Exploiting anisotropic particle shape to electrostatically assemble colloidal molecules with high yield and purity

Hypothesis: Colloidal molecules with anisotropic shapes and interactions are powerful model systems for deciphering the behavior of real molecules and building units for creating materials with designed properties. While many strategies for their assembly have been developed, they typically yield a broad distribution or are limited to a specific type. We hypothesize that the shape and relative sizes of colloidal particles can be exploited to efficiently direct their assembly into colloidal molecules of desired valence. Experiments: We exploit electrostatic self-assembly of negatively charged spheres made from either polystyrene or silica onto positively charged hematite cubes. We thoroughly analyze the role of the shape and size ratio of particles on the cluster size and yield of colloidal molecules. Findings: Using a combination of experiments and simulations, we demonstrate that cubic particle shape is crucial to generate high yields of distinct colloidal molecules over a wide variety of size ratios. We find that electrostatic repulsion between the satellite spheres is important to leverage the templating effect of the cubes, leading the spheres to preferentially assemble on the facets rather than the edges and corners of the cube. Furthermore, we reveal that our protocol is not affected by the specific choice of the material of the colloidal particles. Finally, we show that the permanent magnetic dipole moment of the hematite cubes can be utilized to separate colloidal molecules from non-assembled satellite particles. Our simple and effective strategy might be extended to other templating particle shapes, thereby greatly expanding the library of colloidal molecules that can be achieved with high yield and purity.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Kinetic phase diagram for nucleation and growth of competing crystal polymorphs in charged colloidal suspensions

We determine the full kinetic phase diagram for nucleation and growth of crystal phases in a suspension of charged colloids. We calculate nucleation barrier heights for face-centered cubic (fcc) and body-centered cubic (bcc) crystal phases for varying screening lengths and supersaturations using the seeding approach in extensive simulations. Using classical nucleation theory, we determine for the entire metastable fluid region the crystal polymorph with the lowest nucleation barrier. Surprisingly, we find a regime close to the triple point where metastable bcc can form due to a lower nucleation barrier, even though fcc is the stable phase. For higher supersaturation, we find that the difference in barrier heights decreases and we observe a mix of hexagonal close-packed (hcp), fcc and bcc structures in the growth of crystalline seeds as well as in spontaneously formed crystals.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Machine-learning effective many-body potentials for anisotropic particles using orientation-dependent symmetry functions

Spherically-symmetric atom-centered descriptors of atomic environments have been widely used for constructing potential or free energy surfaces of atomistic and colloidal systems and to characterize local structures using machine learning techniques. However, when particle shapes are non-spherical, as in the case of rods and ellipsoids, standard spherically-symmetric structure functions alone produce imprecise descriptions of local environments. In order to account for the effects of orientation, we introduce two- and three-body orientation-dependent particle-centered descriptors for systems composed of rod-like particles. To demonstrate the suitability of the proposed functions, we use an efficient feature selection scheme and simple linear regression to construct coarse-grained many-body interaction potentials for computationally-efficient simulations of model systems consisting of colloidal particles with anisotropic shape: mixtures of colloidal rods and nonadsorbing polymer, hard rods enclosed by an elastic microgel shell, and ligand-stabilized nanorods. We validate the machine-learning (ML) effective many-body potentials based on orientation-dependent symmetry functions by using them in direct coexistence simulations to map out the phase behavior of colloidal rods and non-adsorbing polymer. We find good agreement with results obtained from simulations of the true binary mixture, demonstrating that the effective interactions are well-described by the orientation-dependent ML potentials.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Steering self-organisation through confinement

Self-organisation is the spontaneous emergence of spatio-temporal structures and patterns from the interaction of smaller individual units. Examples are found across many scales in very different systems and scientific disciplines, from physics, materials science and robotics to biology, geophysics and astronomy. Recent research has highlighted how self-organisation can be both mediated and controlled by confinement. Confinement occurs through interactions with boundaries, and can function as either a catalyst or inhibitor of self-organisation. It can then become a means to actively steer the emergence or suppression of collective phenomena in space and time. Here, to provide a common framework for future research, we examine the role of confinement in self-organisation and identify overarching scientific challenges across disciplines that need to be addressed to harness its full scientific and technological potential. This framework will not only accelerate the generation of a common deeper understanding of self-organisation but also trigger the development of innovative strategies to steer it through confinement, with impact, e.g., on the design of smarter materials, tissue engineering for biomedicine and crowd management.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Predicting the phase behavior of mixtures of active spherical particles

An important question in the field of active matter is whether or not it is possible to predict the phase behavior of these systems. Here, we study the phase coexistence of binary mixtures of torque-free active Brownian particles, for both systems with purely repulsive interactions and systems with attractions. Using Brownian dynamics simulations, we show that phase coexistences can be predicted quantitatively for these systems by measuring the pressure and "reservoir densities". Specifically, in agreement with previous literature, we find that the coexisting phases are in mechanical equilibrium, i.e. the two phases have the same pressure. Importantly, we also demonstrate that the coexisting phases are in chemical equilibrium by bringing each phase into contact with particle reservoirs, and showing that for each species these reservoirs are characterized by the same density for both phases. Using this requirement of mechanical and chemical equilibrium we accurately construct the phase boundaries from properties which can be measured purely from the individual coexisting phases. This result highlights that torque-free active Brownian systems follow simple coexistence rules, thus shedding new light on their thermodynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Speeding up dynamics by tuning the non-commensurate size of rod-like particles in a smectic phase

Using simulations, we study the diffusion of rod-like guest particles in a smectic environment of rod-like host particles. We find that the dynamics of guest rods across smectic layers changes from a fast nematic-like diffusion to a slow hopping-type dynamics via an intermediate switching regime by varying the length of the guest rods with respect to the smectic layer spacing. We determine the optimal rod length that yields the fastest and the slowest diffusion in a lamellar environment. We show that this behavior can be rationalized by a complex 1D effective periodic potential exhibiting two energy barriers, resulting in a varying preferred mean position of the guest particle in the smectic layer. The interplay of these two barriers controls the dynamics of the guest particles yielding a slow, an intermediate and a fast diffusion regime depending on the particle length.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Binary icosahedral quasicrystals of hard spheres in spherical confinement

The influence of geometry on the local and global packing of particles is important to many fundamental and applied research themes such as the structure and stability of liquids, crystals and glasses. Here, we show by experiments and simulations that a binary mixture of hard-sphere-like particles crystallizing into the MgZn2 Laves phase in bulk, spontaneously forms 3D icosahedral quasicrystals in slowly drying droplets. Moreover, the local symmetry of 70-80% of the particles changes to that of the MgCu2 Laves phase. Both of these findings are significant for photonic applications. If the stoichiometry deviates from that of the Laves phase, our experiments show that the crystallization of MgZn2 is hardly affected by the spherical confinement. Our simulations show that the quasicrystals nucleate away from the spherical boundary and grow along five-fold symmetric structures. Our findings not only open the way for particle-level studies of nucleation and growth of 3D quasicrystals, but also of binary crystallization.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The role of topological defects in the two-stage melting and elastic behavior of active Brownian particles

We find that crystalline states of repulsive active Brownian particles at high activity melt into a hexatic state but this transition is not driven by an unbinding of bound dislocation pairs as suggested by the Kosterlitz-Thouless-Halperin-Nelson-Young (KTHNY) theory. Upon reducing the density, the crystalline state melts into a high-density hexatic state devoid of any defects. Decreasing the density further, the dislocations proliferate and introduce plasticity in the system, nevertheless maintaining the hexatic state, but eventually melting into a fluid state. Remarkably, the elastic constants of active solids are equal to those of their passive counterparts, as the swim contribution to the stress tensor is negligible in the solid state. The sole effect of activity is that the stable solid regime shifts to higher densities. Furthermore, discontinuities in the elastic constants as a function of density correspond to changes in the defect concentrations rather than to the solid-hexatic transition.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Unsupervised learning for local structure detection in colloidal systems

We introduce a simple, fast, and easy to implement unsupervised learning algorithm for detecting different local environments on a single-particle level in colloidal systems. In this algorithm, we use a vector of standard bond-orientational order parameters to describe the local environment of each particle. We then use a neural-network-based autoencoder combined with Gaussian mixture models in order to autonomously group together similar environments. We test the performance of the method on snapshots of a wide variety of colloidal systems obtained via computer simulations, ranging from simple isotropically interacting systems, to binary mixtures, and even anisotropic hard cubes. Additionally, we look at a variety of common self-assembled situations such as fluid-crystal and crystal-crystal coexistences, grain boundaries, and nucleation. In all cases, we are able to identify the relevant local environments to a similar precision as "standard", manually-tuned and system-specific, order parameters. In addition to classifying such environments, we also use the trained autoencoder in order to determine the most relevant bond orientational order parameters in the systems analyzed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Phase diagrams of charged colloidal rods: can a uniaxial charge distribution break chiral symmetry?

We construct phase diagrams for charged rodlike colloids within the second-virial approximation as a function of rod concentration, salt concentration, and colloidal charge. Besides the expected isotropic-nematic transition, we also find parameter regimes with a coexistence between a nematic and a second, more highly aligned nematic phase including an isotropic-nematic-nematic triple point and a nematic-nematic critical point, which can all be explained in terms of the twisting effect. We compute the Frank elastic constants to see if the twist elastic constant can become negative, which would indicate the possibility of a cholesteric phase spontaneously forming. Although the twisting effect reduces the twist elastic constant, we find that it always remains positive. In addition, we find that for finite aspect-ratio rods the twist elastic constant is also always positive, such that there is no evidence of chiral symmetry breaking due to a uniaxial charge distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Removing grain boundaries from three-dimensional colloidal crystals using active dopants

Using computer simulations we explore how grain boundaries can be removed from three-dimensional colloidal crystals by doping with a small fraction of active colloids. We show that for sufficient self-propulsion, the system is driven into a crystal-fluid coexistence. In this phase separated regime, the active dopants become mobile and spontaneously gather at the grain boundaries. The resulting surface melting and recrystallization of domains result in the motion of the grain boundaries over time and lead to the formation of a large single crystal. However, when the self-propulsion is too low to cause a phase separation, we observe no significant enhancement of grain growth.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stabilizing the Hexagonal Close Packed Structure of Hard Spheres with Polymers : Phase diagram, Structure, and Dynamics

We study the phase behaviour of a binary mixture of colloidal hard spheres and freely-jointed chains of beads using Monte Carlo simulations. Recently Panagiotopoulos and coworkers predicted [Nat. Commun. 5, 4472 (2014)] that the hexagonal close packed (HCP) structure of hard spheres can be stabilized in such a mixture due to the interplay between polymer and the void structure in the crystal phase. Their predictions were based on estimates of the free-energy penalty for adding a single hard polymer chain in the HCP and the competing face centered cubic (FCC) phase. Here we calculate the phase diagram using free-energy calculations of the full binary mixture and find a broad fluid-solid coexistence region and a metastable gas-liquid coexistence region. For the colloid-monomer size ratio considered in this work, we find that the HCP phase is only stable in a small window at relatively high polymer reservoir packing fractions, where the coexisting HCP phase is nearly close packed. Additionally we investigate the structure and dynamic behaviour of these mixtures.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Vapour-Liquid Coexistence of an Active Lennard-Jones fluid

We study a three-dimensional system of self-propelled Lennard-Jones particles using Brownian Dynamics simulations. Using recent theoretical results for active matter, we calculate the pressure and report equations of state for the system. Additionally, we chart the vapour-liquid coexistence and show that the coexistence densities can be well described using simple power laws. Lastly, we demonstrate that our out-of-equilibrium system shows deviations from both the law of rectilinear diameters and the law of corresponding states.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Novel Chiral Phase of Achiral Hard Triangles and an Entropy-Driven Demixing of Enantiomers

We investigate the phase behavior of a system of hard equilateral and right-angled triangles in two dimensions using Monte Carlo simulations. Hard equilateral triangles undergo a continuous isotropic-triatic liquid crystal phase transition at packing fraction $ϕ=0.7$. Similarly, hard right-angled isosceles triangles exhibit a first-order phase transition from an isotropic fluid phase to a rhombic liquid crystal phase with a coexistence region $ϕ\in \left[0.733,0.782\right]$. Both these liquid crystal phases undergo a continuous phase transition to their respective close-packed crystal structures at high pressures. Although the particles and their close-packed crystals are both achiral, the solid phases of equilateral and right-angled triangles exhibit spontaneous chiral symmetry breaking at sufficiently high packing fractions. The colloidal triangles rotate either in clockwise or anti-clockwise direction with respect to one of the lattice vectors for packing fractions higher than $ϕ_χ$. As a consequence, these triangles spontaneously form a regular lattice of left- or right-handed chiral holes which are surrounded by six triangles in the case of equilateral triangles and four or eight triangles for right-angled triangles. Moreover, our simulations show a spontaneous entropy-driven demixing transition of the right- and left-handed "enantiomers".

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Entropy-driven formation of prolate and oblate cholesteric phases by computer simulations

Predicting the macroscopic chiral behaviour of cholesteric liquid crystals from the microscopic chirality of the particles is highly non-trivial, even when the chiral interactions are purely entropic in nature. Here we introduce a novel chiral hard-particle model, namely particles with a twisted polyhedral shape and obtain, for the first time, a stable fully-entropy-driven cholesteric phase by computer simulations. By slightly modifying the triangular base of the particle, we are able to switch from a left-handed prolate to a right-handed oblate cholesteric using the same right-handed twisted particle model. Furthermore, we find qualitative agreement with an Onsager-like theory, suggesting that the latter can be used as a quick tool to scan the huge parameter space associated to the microscopic chirality. Our results unveil how the competition between particle biaxiality and chirality is reflected into the nematic self-organization and new theoretical challenges on the self-assembly of chiral particles can be undertaken.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Phase Separation and Self-Assembly in a Fluid of Mickey Mouse Particles

Recent developments in the synthesis of colloidal particles allow for control over shape and inter-particle interaction. One example, among others, is the so-called "Mickey Mouse" (MM) particle for which the self-assembly properties have been previously studied yielding a stable cluster phase together with elongated, tube-like structures. Here, we investigate under which conditions a fluid of Mickey Mouse particles can yield phase separation and how the self-assembly behaviour affects the gas-liquid coexistence. We vary the distance between the repulsive and the attractive lobes (bond length), and the interaction range, and follow the evolution of the gas-liquid (GL) coexistence curve. We find that upon increasing the bond length distance the binodal line shifts to lower temperatures, and that the interaction range controls the transition between phase separation and self-assembly of clusters. Upon further reduction of the interaction range and temperature, the clusters assume an increasingly ordered tube-like shape, ultimately matching the one previously reported in literature. These results are of interest when designing particle shape and particle-particle interaction for self-assembly processes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cholesterics of colloidal helices: Predicting the macroscopic pitch from the particle shape and thermodynamic state

Building a general theoretical framework to describe the microscopic origin of macroscopic chirality in (colloidal) liquid crystals is a long-standing challenge. Here, we combine classical density functional theory with Monte Carlo calculations of virial-type coefficients, to obtain the equilibrium cholesteric pitch as a function of thermodynamic state and microscopic details. Applying the theory to hard helices, we observe both right- and left-handed cholesteric phases that depend on a subtle combination of particle geometry and system density. In particular, we find that entropy alone can even lead to a (double) inversion in the cholesteric sense of twist upon changing the packing fraction. We show how the competition between single-particle properties (shape) and thermodynamics (local alignment) dictates the macroscopic chiral behavior. Moreover, by expanding our free-energy functional we are able to assess, quantitatively, Straley's theory of weak chirality, used in several earlier studies. Furthermore, by extending our theory to different lyotropic and thermotropic liquid-crystal models, we analyze the effect of an additional soft interaction on the chiral behavior of the helices. Finally, we provide some guidelines for the description of more complex chiral phases, like twist-bend nematics. Our results provide new insights on the role of entropy in the microscopic origin of this state of matter.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Critical Casimir Forces and Colloidal Phase Transitions in a Near-Critical Solvent : A Simple Model Reveals a Rich Phase Diagram

From experimental studies it is well-known that colloidal particles suspended in a near-critical binary solvent exhibit interesting aggregation phenomena, often associated with colloidal phase transitions, and assumed to be driven by long-ranged solvent mediated (SM) interactions (critical Casimir forces), set by the (diverging) correlation length of the solvent. We present the first simulation and theoretical study of an explicit model of a ternary mixture that mimics this situation. Both the effective SM pair interactions and the full ternary phase diagram are determined for Brownian discs suspended in an explicit two-dimensional supercritical binary liquid mixture. Gas-liquid and fluid-solid transitions are observed in a region that extends well-away from criticality of the solvent reservoir. We discuss to what extent an effective pair-potential description can account for the phase behavior we observe. Our study provides a fresh perspective on how critical fluctuations of the solvent might influence colloidal self-assembly.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Phase Behavior of a Family of Truncated Hard Cubes

In continuation of our work in [A.P. Gantapara et al., Phys. Rev. Lett. 111, 015501 (2013)], we investigate here the thermodynamic phase behavior of a family of truncated hard cubes, for which the shape evolves smoothly from a cube via a cuboctahedron to an octahedron. We used Monte Carlo simulations and free-energy calculations to establish the full phase diagram. This phase diagram exhibits a remarkable richness in crystal and mesophase structures, depending sensitively on the precise particle shape. In addition, we examined in detail the nature of the plastic crystal (rotator) phases that appear for intermediate densities and levels of truncation. Our results allow us to probe the relation between phase behavior and building-block shape and to further the understanding of rotator phases. Furthermore, the phase diagram presented here should prove instrumental for guiding future experimental studies on similarly-haped nanoparticles and the creation of new materials.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Self-assembly of "Mickey Mouse" shaped colloids into tube-like structures: experiments and simulations

The self-assembly of anisotropic patchy particles with triangular shape was studied by experiments and computer simulations. The colloidal particles were synthesized in a two-step seeded emulsion polymerization process, and consist of a central smooth lobe connected to two rough lobes at an angle of $\sim$90$^{\circ}$, resembling the shape of a "Mickey Mouse" head. Due to the difference in overlap volume, adding an appropriate depletant induces an attractive interaction between the smooth lobes of the colloids only, while the two rough lobes act as steric constraints. The essentially planar geometry of the "Mickey Mouse" particles is a first geometric deviation of dumbbell shaped patchy particles. This new geometry is expected to form one-dimensional tube-like structures rather than spherical, essentially zero-dimensional micelles. At sufficiently strong attractions, we indeed find tube-like structures with the sticky lobes at the core and the non-sticky lobes pointing out as steric constraints that limit the growth to one direction, providing the tubes with a well-defined diameter but variable length both in experiments and simulations. In the simulations, we found that the internal structure of the tubular fragments could either be straight or twisted into so-called Bernal spirals.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Self-Assembly of Patchy Colloidal Dumbbells

We employ Monte Carlo simulations to investigate the self-assembly of patchy colloidal dumbbells interacting via a modified Kern-Frenkel potential by probing the system concentration and dumbbell shape. We consider dumbbells consisting of one attractive sphere with diameter $σ_1$ and one repulsive sphere with diameter $σ_2$ and center-to-center distance $d$ between the spheres. For three different size ratios, we study the self-assembled structures for different separations $l = 2d/(σ_1+σ_2)$ between the two spheres. In particular, we focus on structures that can be assembled from the homogeneous fluid, as these might be of interest in experiments. We use cluster order parameters to classify the shape of the formed structures. When the size of the spheres is almost equal, $q=σ_2/σ_1=1.035$, we find that, upon increasing $l$, spherical micelles are transformed to elongated micelles and finally to vesicles and bilayers. For size ratio $q=1.25$ we observe a continuously tunable transition from spherical to elongated micelles upon increasing the sphere separation. For size ratio $q=0.95$ we find bilayers and vesicles, plus faceted polyhedra and liquid droplets. Our results identify key parameters to create colloidal vesicles with attractive dumbbells in experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Two-stage melting induced by dislocations and grain boundaries in Monolayers of Hard Spheres

Melting in two-dimensional systems has remained controversial as theory, simulations, and experiments show contrasting results. One issue that obscures this discussion is whether or not theoretical predictions on strictly 2D systems describe those of quasi-2D experimental systems, where out-of-plane fluctuations may alter the melting mechanism. Using event-driven Molecular Dynamics simulations, we find that the peculiar two-stage melting scenario of a continuous solid-hexatic and a first-order hexatic-liquid transition as observed for a truly 2D system of hard disks [Bernard and Krauth, Phys. Rev. Lett. 107, 155704 (2011)] persists for a quasi-2D system of hard spheres with out-of-plane particle motions as high as half the particle diameter. By calculating the renormalized Young's modulus, we show that the solid-hexatic transition is of the Kosterlitz-Thouless type, and occurs via dissociation of bound dislocation pairs. In addition, we find a first-order hexatic-liquid transition that seems to be driven by a spontaneous proliferation of grain boundaries.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Orientational order of carbon nanotube guests in a nematic host suspension of colloidal viral rods

In order to investigate the coupling between the degrees of alignment of elongated particles in binary nematic dispersions, surfactant stabilized single-wall carbon nanotubes (CNTs) have been added to nematic suspensions of colloidal rodlike viruses in aqueous solution.We have independently measured the orientational order parameter of both components of the guest-host system by means of polarized Raman spectroscopy and by optical birefringence, respectively. Our system allows us therefore to probe the regime where the guest particles (CNTs) are shorter and thinner than the fd virus host particles. We show that the degree of order of the CNTs is systematically smaller than that of the fd virus particles for the whole nematic range. These measurements are in good agreement with predictions of an Onsager-type second-viral theory, which explicitly includes the flexibility of the virus particles, and the polydispersity of the CNTs.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Pushing the glass transition towards random close packing using self-propelled hard spheres

Although the concept of random close packing with an almost universal packing fraction of ~ 0.64 for hard spheres was introduced more than half a century ago, there are still ongoing debates. The main difficulty in searching the densest packing is that states with packing fractions beyond the glass transition at ~ 0.58 are inherently non-equilibrium systems, where the dynamics slows down with a structural relaxation time diverging with density; hence, the random close packing is inaccessible. Here we perform simulations of self-propelled hard spheres, and we find that with increasing activity the relaxation dynamics can be sped up by orders of magnitude. The glass transition shifts to higher packing fractions upon increasing the activity, allowing the study of sphere packings with fluid-like dynamics at packing fractions close to random close packing. Our study opens new possibilities of investigating dense packings and the glass transition in systems of hard particles.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Depletion-induced biaxial nematic states of boardlike particles

With the aim of investigating the stability conditions of biaxial nematic liquid crystals, we study the effect of adding a non-adsorbing ideal depletant on the phase behavior of colloidal hard boardlike particles. We take into account the presence of the depletant by introducing an effective depletion attraction between a pair of boardlike particles. At fixed depletant fugacity, the stable liquid crystal phase is determined through a mean-field theory with restricted orientations. Interestingly, we predict that for slightly elongated boardlike particles a critical depletant density exists, where the system undergoes a direct transition from an isotropic liquid to a biaxial nematic phase. As a consequence, by tuning the depletant density, an easy experimental control parameter, one can stabilize states of high biaxial nematic order even when these states are unstable for pure systems of boardlike particles.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Effect of bond length fluctuations on crystal nucleation of hard bead chains

We study the spontaneous nucleation and crystallization of linear and cyclic chains of flexibly connected hard spheres using extensive molecular dynamics simulations. To this end, we present a novel event-driven molecular dynamics simulation method, which is easy to implement and very efficient. We find that the nucleation rates are predominately determined by the number of bonds per sphere in the system, rather than the precise details of the chain topology, chain length, and polymer composition. O ur results thus show that the crystal nucleation rate of bead chains can be enhanced by adding monomers to the system. In addition, we find that the resulting crystal n uclei contain significantly more face-centered-cubic than hexagonal-close-packed ordered particles. More surprisingly, the resulting crystal nuclei possess a range of crystal morphologies including structures with a five-fold symmetry.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Electrostatic Interactions between Janus Particles

In this paper we study the electrostatic properties of `Janus' spheres with unequal charge densities on both hemispheres. We introduce a method to compare primitive-model Monte Carlo simulations of the ionic double layer with predictions of (mean-field) nonlinear Poisson-Boltzmann theory. We also derive practical DLVO-like expressions that describe the Janus-particle pair interactions by mean-field theory. Using a large set of parameters, we are able to probe the range of validity of the Poisson-Boltzmann approximation, and thus of DLVO-like theories, for such particles. For homogeneously charged spheres this range corresponds well to the range that was predicted by field-theoretical studies of homogeneously charged flat surfaces. Moreover, we find similar ranges for colloids with a Janus-type charge distribution. The techniques and parameters we introduce show promise for future studies of an even wider class of charged-patterned particles.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Phase diagram of colloidal hard superballs: from cubes via spheres to octahedra

The phase diagram of colloidal hard superballs, of which the shape interpolates between cubes and octahedra via spheres, is determined by free-energy calculations in Monte Carlo simulations. We discover not only a stable face-centered cubic (fcc) plastic crystal phase for near-spherical particles, but also a stable body-centered cubic (bcc) plastic crystal close to the octahedron shape, and in fact even coexistence of these two plastic crystals with a substantial density gap. The plastic fcc and bcc crystals are, however, unstable in the cube and octahedron limit, suggesting that the rounded corners of superballs play an important role in stablizing the rotator phases. In addition, we observe a two-step melting phenomenon for hard octahedra, in which the Minkowski crystal melts into a metastable bcc plastic crystal before melting into the fluid phase.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Predicting Patchy Particle Crystals: Variable Box Shape Simulations and Evolutionary Algorithms

We consider several patchy particle models that have been proposed in literature and we investigate their candidate crystal structures in a systematic way. We compare two different algorithms for predicting crystal structures: (i) an approach based on Monte Carlo simulations in the isobaric-isothermal ensemble and (ii) an optimization technique based on ideas of evolutionary algorithms. We show that the two methods are equally successful and provide consistent results on crystalline phases of patchy particle systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Vacancy-stabilized crystalline order in hard cubes

We examine the effect of vacancies on the phase behavior and structure of systems consisting of hard cubes using event-driven molecular dynamics and Monte Carlo simulations. We find a first-order phase transition between a fluid and a simple cubic crystal phase that is stabilized by a surprisingly large number of vacancies, reaching a net vacancy concentration of ~6.4% near bulk coexistence. Remarkably, we find that vacancies increase the positional order in the system. Finally, we show that the vacancies are delocalized and therefore hard to detect.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dense Regular Packings of Irregular Non-Convex Particles

We present a new numerical scheme to study systems of non-convex, irregular, and punctured particles in an efficient manner. We employ this method to analyze regular packings of odd-shaped bodies, not only from a nanoparticle but also both from a computational geometry perspective. Besides determining close-packed structures for many shapes, we also discover a new denser configuration for Truncated Tetrahedra. Moreover, we consider recently synthesized nanoparticles and colloids, where we focus on the excluded volume interactions, to show the applicability of our method in the investigation of their crystal structures and phase behavior. Extensions to the presented scheme include the incorporation of soft particle-particle interactions, the study of quasicrystalline systems, and random packings.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Frustration of the isotropic-columnar phase transition of colloidal hard platelets by a transient cubatic phase

Using simulations and theory, we show that the cubatic phase is metastable for three model hard platelets. The locally favored structures of perpendicular particle stacks in the fluid prevent the formation of the columnar phase through geometric frustration resulting in vitrification. Also, we find a direct link between structure and dynamic heterogeneities in the cooperative rotation of particle stacks, which is crucial for the devitrification process. Finally, we show that the life time of the glassy cubatic phase can be tuned by surprisingly small differences in particle shape.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Ion association in low-polarity solvents: comparisons between theory, simulation, and experiment

The association of ions in electrolyte solutions at very low concentration and low temperature is studied using computer simulations and quasi-chemical ion-pairing theory. The specific case of the restricted primitive model (charged hard spheres) is considered. Specialised simulation techniques are employed that lead to efficient sampling of the arrangements and distributions of clusters and free ions, even at conditions corresponding to nanomolar solutions of simple salts in solvents with dielectric constants in the range 5-10, as used in recent experimental work on charged-colloid sus- pensions. A direct comparison is effected between theory and simulation using a variety of clustering criteria and theoretical approximations. It is shown that conventional distance-based cluster criteria can give erroneous results. A reliable set of theoretical and simulation estimators for the degree of association is proposed. The ion-pairing theory is then compared to experimental results for salt solutions in low-polarity solvents. The agreement is excellent, and on this basis some calculations are made for the screening lengths which will figure in the treatment of colloid-colloid interactions in such solutions. The accord with available experimental results is complete.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Measuring Colloidal Forces from Particle Position Deviations inside an Optical Trap

We measure interaction forces between pairs of charged PMMA colloidal particles suspended in a relatively low-polar medium (5 $\lesssim ε\lesssim$ 8) directly from the deviations of particle positions inside two time-shared optical traps. The particles are confined to optical point traps; one is held in a stationary trap and the other particle is brought closer in small steps while tracking the particle positions using confocal microscopy. From the observed particle positions inside the traps we calculate the interparticle forces using an ensemble-averaged particle displacement-force relationship. The force measurements are confirmed by independent measurements of the different parameters using electrophoresis and a scaling law for the liquid-solid phase transition. When increasing the salt concentration by exposing the sample to UV light, the force measurements agree well with the classical DLVO theory assuming a constant surface potential. On the other hand, when adding tetrabutylammonium chloride (TBAC) to vary the salt concentration, surface charge regulation seems to play an important role.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Phase diagram and effective shape of semi-flexible colloidal rods and biopolymers

We study suspensions of semi-flexible colloidal rods and biopolymers using an Onsager-type second-virial functional for a segmented-chain model. For suspensions of thin and thick fd virus particles we calculate phase diagrams in quantitative agreement with experimental observations, and we find their effective state-point dependent shape to be much shorter and thicker than the actual shape. We also calculate the stretching of worm-like micelles in a host fd virus solution, again finding agreement with experiments. For both systems, our results show that the fd virus stiffness can play a key role in system behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Collective Diffusion of Colloidal Hard Rods in Smectic Liquid Crystals: Effect of Particle Anisotropy

We study the layer-to-layer diffusion in smectic-A liquid crystals of colloidal hard rods with different length-to-diameter ratios using computer simulations. The layered arrangement of the smectic phase yields a hopping-type diffusion due to the presence of permanent barriers and transient cages. Remarkably, we detect stringlike clusters composed of inter-layer rods moving cooperatively along the nematic director. Furthermore, we find that the structural relaxation in equilibrium smectic phases shows interesting similarities with that of out-of-equilibrium supercooled liquids, although there the particles are kinetically trapped in transient rather than permanent cages. Additionally, at fixed packing fraction we find that the barrier height increases with increasing particle anisotropy, and hence the dynamics is more heterogeneous and non-Gaussian for longer rods, yielding a lower diffusion coefficient along the nematic director and smaller clusters of inter-layer particles that move less cooperatively. At fixed barrier height, the dynamics becomes more non-Gaussian and heterogeneous for longer rods that move more collectively giving rise to a higher diffusion coefficient along the nematic director.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Computer simulations of the restricted primitive model at very low temperature and density

The problem of successfully simulating ionic fluids at low temperature and low density states is well known in the simulation literature: using conventional methods, the system is not able to equilibrate rapidly due to the presence of strongly associated cation-anion pairs. In this manuscript we present a numerical method for speeding up computer simulations of the restricted primitive model (RPM) at low temperatures (around the critical temperature) and at very low densities (down to $10^{-10}σ^{-3}$, where $σ$ is the ion diameter). Experimentally, this regime corresponds to typical concentrations of electrolytes in nonaqueous solvents. As far as we are aware, this is the first time that the RPM has been equilibrated at such extremely low concentrations. More generally, this method could be used to equilibrate other systems that form aggregates at low concentrations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Crystal nucleation of colloidal hard dumbbells

Using computer simulations we investigate the homogeneous crystal nucleation in suspensions of colloidal hard dumbbells. The free energy barriers are determined by Monte Carlo simulations using the umbrella sampling technique. We calculate the nucleation rates for the plastic crystal and the aperiodic crystal phase using the kinetic prefactor as determined from event driven molecular dynamics simulations. We find good agreement with the nucleation rates determined from spontaneous nucleation events observed in event driven molecular dynamics simulations within error bars of one order of magnitude. We study the effect of aspect ratio of the dumbbells on the nucleation of plastic and aperiodic crystal phases and we also determine the structure of the critical nuclei. Moreover, we find that the nucleation of the aligned CP1 crystal phase is strongly suppressed by a high free energy barrier at low supersaturations and slow dynamics at high supersaturations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Crystal nucleation of hard spheres using molecular dynamics, umbrella sampling and forward flux sampling: A comparison of simulation techniques

Over the last number of years several simulation methods have been introduced to study rare events such as nucleation. In this paper we examine the crystal nucleation rate of hard spheres using three such numerical techniques: molecular dynamics, forward flux sampling and a Bennett-Chandler type theory where the nucleation barrier is determined using umbrella sampling simulations. The resulting nucleation rates are compared with the experimental rates of Harland and Van Megen [J. L. Harland and W. van Megen, Phys. Rev. E 55, 3054 (1997)], Sinn et al. [C. Sinn et al., Prog. Colloid Polym. Sci. 118, 266 (2001)] and Schatzel and Ackerson [K. Schatzel and B.J. Ackerson, Phys. Rev. E, 48, 3766 (1993)] and the predicted rates for monodisperse and 5% polydisperse hard spheres of Auer and Frenkel [S. Auer and D. Frenkel, Nature 409, 1020 (2001)]. When the rates are examined in long-time diffusion units, we find agreement between all the theoretically predicted nucleation rates, however, the experimental results display a markedly different behaviour for low supersaturation. Additionally, we examined the pre-critical nuclei arising in the molecular dynamics, forward flux sampling, and umbrella sampling simulations. The structure of the nuclei appear independent of the simulation method, and in all cases, the nuclei contain on average significantly more face-centered-cubic ordered particles than hexagonal-close-packed ordered particles.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Glassy dynamics, spinodal fluctuations, and the kinetic limit of hard-rod nucleation

Using simulations we identify three dynamic regimes in supersaturated isotropic fluid states of short hard rods: (i) for moderate supersaturations we observe nucleation of multi-layered crystalline clusters; (ii) at higher supersaturation, we find nucleation of small crystallites which arrange into long-lived locally favored structures that get kinetically arrested, while (iii) at even higher supersaturation the dynamic arrest is due to the conventional cage-trapping glass transition. For longer rods we find that the formation of the (stable) smectic phase out of a supersaturated isotropic state is strongly suppressed by an isotropic-nematic spinodal instability that causes huge spinodal-like orientation fluctuations with nematic clusters diverging in size. Our results show that glassy dynamics and spinodal instabilities set kinetic limits to nucleation in a highly supersaturated hard-rod fluids.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Heterogeneous Dynamics in Columnar Liquid Crystals of Parallel Hard Rods

In the wake of previous studies on the rattling-and-jumping diffusion in smectic liquid crystal phases of colloidal rods, we analyze here for the first time the heterogeneous dynamics in columnar phases. More specifically, we perform computer simulations to investigate the relaxation dynamics of a binary mixture of perfectly aligned hard spherocylinders. We detect that the columnar arrangement of the system produces free-energy barriers the particles should overcome to jump from one column to another, thus determining a hopping-type diffusion. This phenomenon accounts for the non-Gaussian inter-column diffusion and shows a two-step structural relaxation which is remarkably analogous to that of out-of-equilibrium glass-forming systems and gels. Surprisingly enough, slight deviations from the behavior of simple liquids due to transient cages is also observed in the direction perpendicular to this plane, where the system is usually referred to as liquid-like.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Phase behavior and structure of colloidal bowl-shaped particles: simulations

We study the phase behavior of bowl-shaped particles using computer simulations. These particles were found experimentally to form a meta-stable worm-like fluid phase in which the bowl-shaped particles have a strong tendency to stack on top of each other [M.Marechal et al, Nano Letters 10, 1907 (2010)]. In this work, we show that the transition from the low-density fluid to the worm-like phase has an interesting effect on the equation of state. The simulation results also show that the worm-like fluid phase transforms spontaneously into a columnar phase for bowls that are sufficiently deep. Furthermore, we describe the phase behavior as obtained from free energy calculations employing Monte Carlo simulations. The columnar phase is stable for bowl shapes ranging from infinitely thin bowls to surprisingly shallow bowls. Aside from a large region of stability for the columnar phase, the phase diagram features four novel crystal phases and a region where the stable fluid contains worm-like stacks.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Phase Diagrams of Binary Mixtures of Oppositely Charged Colloids

Phase diagrams of binary mixtures of oppositely charged colloids are calculated theoretically. The proposed mean-field-like formalism interpolates between the limits of a hard-sphere system at high temperatures and the colloidal crystals which minimize Madelung-like energy sums at low temperatures. Comparison with computer simulations of an equimolar mixture of oppositely charged, equally sized spheres indicate semi-quantitative accuracy of the proposed formalism. We calculate global phase diagrams of binary mixtures of equally sized spheres with opposite charges and equal charge magnitude in terms of temperature, pressure, and composition. The influence of the screening of the Coulomb interaction upon the topology of the phase diagram is discussed. Insight into the topology of the global phase diagram as a function of the system parameters leads to predictions on the preparation conditions for specific binary colloidal crystals.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Phase diagrams of colloidal spheres with a constant zeta-potential

We study suspensions of colloidal spheres with a constant zeta-potential within Poisson-Boltzmann theory, quantifying the discharging of the spheres with increasing colloid density and decreasing salt concentration. We use the calculated renormalized charge of the colloids to determine their pairwise effective screened-Coulomb repulsions. Bulk phase diagrams in the colloid concentration-salt concentration representation follow, for various zeta-potentials, by a mapping onto published fits of phase boundaries of point-Yukawa systems. Although the resulting phase diagrams do feature face-centered cubic (fcc) and body-centered cubic (bcc) phases, they are dominated by the (re-entrant) fluid phase due to the colloidal discharging with increasing colloid concentration and decreasing salt concentration.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A Triangular Tessellation Scheme for the Adsorption Free Energy at the Liquid-Liquid Interface: Towards Non-Convex Patterned Colloids

We introduce a new numerical technique, namely triangular tessellation, to calculate the free energy associated with the adsorption of a colloidal particle at a flat interface. The theory and numerical scheme presented here are sufficiently general to handle non-convex patchy colloids with arbitrary surface patterns characterized by a wetting angle, e.g., amphiphilicity. We ignore interfacial deformation due to capillary, electrostatic, or gravitational forces, but the method can be extended to take such effects into account. It is verified that the numerical method presented is accurate and sufficiently stable to be applied to more general situations than presented in this paper. The merits of the tessellation method prove to outweigh those of traditionally used semi-analytic approaches, especially when it comes to generality and applicability.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Non-Gaussian Dynamics in Smectic Liquid Crystals of Parallel Hard Rods

Using computer simulations, we studied the diffusion and structural relaxation in equilibrium smectic liquid crystal bulk phases of parallel hard spherocylinders. These systems exhibit a non-Gaussian layer-to-layer diffusion due to the presence of periodic barriers and transient cages, and show remarkable similarities with the behavior of out-of-equilibrium supercooled liquids. We detect a very slow inter-layer relaxation dynamics over the whole density range of the stable smectic phase which spans a time interval of four time decades. The intrinsic nature of the layered structure yields a hopping-type diffusion which becomes more heterogeneous for higher packing fractions. In contrast, the in-layer dynamics is typical of a dense fluid with a relatively fast decay. Our results on the dynamic behavior agree well with that observed in systems of freely rotating hard rods, but differ quantitavely, as the height of the periodic barriers reduces to zero at the nematic-smectic transition for aligned rods, while it remains finite for freely rotating rods.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Random close packing of polydisperse hard spheres

We study jammed configurations of hard spheres as a function of compression speed using an event-driven molecular dynamics algorithm. We find that during the compression, the pressure follows closely the metastable liquid branch until the system gets arrested into a glass state as the relaxation time exceeds the compression speed. Further compression yields a jammed configuration that can be regarded as the infinite pressure configuration of that glass state. Consequently, we find that the density of jammed packings varies from 0.638 to 0.658 for polydisperse hard spheres and from 0.635 to 0.645 for pure hard spheres upon decreasing the compression rate. This demonstrates that the density at which the systems falls out of equilibrium determines the density at which the system jams at infinite pressure. In addition, we give accurate data for the jamming density as a function of compression rate and size polydispersity.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Self diffusion of particles in complex fluids: temporary cages and permanent barriers

We study the self diffusion of individual particles in dense (non-)uniform complex fluids within dynamic density functional theory and explicitly account for their coupling to the temporally fluctuating background particles. Applying the formalism to rod-like particles in uniaxial nematic and smectic liquid crystals, we find correlated diffusion in different directions: The temporary cage formed by the neighboring particles competes with permanent barriers in periodic inhomogeneous systems such as the lamellar smectic state and delays self diffusion of particles even in uniform systems. We compare our theory with recent experimental data on the self diffusion of fluorescently labelled filamentous virus particles in aqueous dispersions in the smectic phase and find qualitative agreement. This demonstrates the importance of explicitly dealing with the time-dependent self-consistent molecular field that every particle experiences.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

A Case Study of Sedimentation of Charged Colloids: The Primitive Model and the Effective One-Component Approach

Sedimentation-diffusion equilibrium density profiles of suspensions of charge-stabilized colloids are calculated theoretically and by Monte Carlo simulation, both for a one-component model of colloidal particles interacting through pairwise screened-Coulomb repulsions and for a three-component model of colloids, cations, and anions with unscreened-Coulomb interactions. We focus on a state point for which experimental measurements are available [C.P. Royall et al., J. Phys.: Cond. Matt. {\bf 17}, 2315 (2005)]. Despite the apparently different picture that emerges from the one- and three-component model (repelling colloids pushing each other to high altitude in the former, versus a self-generated electric field that pushes the colloids up in the latter), we find similar colloidal density profiles for both models from theory as well as simulation, thereby suggesting that these pictures represent different view points of the same phenomenon. The sedimentation profiles obtained from an effective one-component model by MC simulations and theory, together with MC simulations of the multi-component primitive model are consistent among themselves, but differ quantitatively from the results of a theoretical multi-component description at the Poisson-Boltzmann level. We find that for small and moderate colloid charge the Poisson-Boltzmann theory gives profiles in excellent agreement with the effective one-component theory if a smaller effective charge is used. We attribute this discrepancy to the poor treatment of correlations in the Poisson-Boltzmann theory.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Evidence for out-of-equilibrium crystal nucleation in suspensions of oppositely charged colloids

We report a numerical study of the rate of crystal nucleation in a binary suspension of oppositely charged colloids. Two different crystal structures compete in the thermodynamic conditions under study. We find that the crystal phase that nucleates is metastable and, more surprisingly, its nucleation free energy barrier is not the lowest one. This implies that, during nucleation, there is insufficient time for sub-critical nuclei to relax to their lowest free-energy structure. Such behavior is in direct contradiction with the common assumption that the phase that crystallizes most readily is the one with the lowest free-energy barrier for nucleation. The phenomenon that we describe should be relevant for crystallization experiments where competing solid structures are not connected by an easy transformation.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Phase behaviour of charged colloidal sphere dispersions with added polymer chains

We study the stability of mixtures of highly screened repulsive charged spheres and non-adsorbing ideal polymer chains in a common solvent using free volume theory. The effective interaction between charged colloids in an aqueous salt solution is described by a screened-Coulomb pair potential, which supplements the pure hard-sphere interaction. The ideal polymer chains are treated as spheres that are excluded from the colloids by a hard-core interaction, whereas the interaction between two ideal chains is set to zero. In addition, we investigate the phase behaviour of charged colloid-polymer mixtures in computer simulations, using the two-body (Asakura-Oosawa pair potential) approximation to the effective one-component Hamiltonian of the charged colloids. Both our results obtained from simulations and from free volume theory show similar trends. We find that the screened-Coulomb repulsion counteracts the effect of the effective polymer-mediated attraction. For mixtures of small polymers and relatively large charged colloidal spheres, the fluid-crystal transition shifts to significantly larger polymer concentrations with increasing range of the screened-Coulomb repulsion. For relatively large polymers, the effect of the screened-Coulomb repulsion is weaker. The resulting fluid-fluid binodal is only slightly shifted towards larger polymer concentrations upon increasing the range of the screened-Coulomb repulsion. In conclusion, our results show that the miscibility of dispersions containing charged colloids and neutral non-adsorbing polymers increases, upon increasing the range of the screened-Coulomb repulsion, or upon lowering the salt concentration, especially when the polymers are small compared to the colloids.

0 citations0 reviewsNo code linkedOpen access