Source author record

René van Roij

René van Roij 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

37works
11topics
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

37 published item(s)

preprint2026arXiv

Optical Voltammetry of redox processes inside a nanohole with opto-iontronic microscopy

Cyclic Voltammetry (CV) is the most commonly used method in electrochemistry to characterize electrochemical reactions, usually involving macroscopic electrodes. Here we demonstrate an optical CV technique called Opto-iontronic Microscopy, which is capable of monitoring electrochemical processes at the nanoscale. By integrating optical microscopy with nanohole electrodes, we enhance sensitivity in detecting redox reactions within volumes as small as an attoliter ($(100 \text{~nm})^{3}$). This technique uses total internal reflection illumination, Electric-double-layer modulation, cyclic voltammetry, and lock-in detection, to probe ion dynamics inside nanoholes. We applied this method to study EDL (dis)charging coupled to ferrocenedimethanol (Fc(MeOH)$_2$) redox reactions. Experimental results were validated against a theoretical Poisson-Nernst-Planck-Butler-Volmer model, providing insights into ion concentration changes of reaction species that contribute to the optical contrast. This work opens up opportunities for high-sensitivity, label-free analysis of electrochemical reactions in nanoconfined environments, with potential applications in pure nanocrystal growth and monitoring.

preprint2022arXiv

Absence of anomalous underscreening in highly concentrated aqueous electrolytes confined between smooth silica surfaces

Recent experiments and a series of subsequent theoretical studies suggest the occurrence of universal underscreening in highly concentrated electrolyte solutions. We performed a set of systematic Atomic Force Spectroscopy measurements for aqueous salt solutions in a concentration range from 1 mM to 5 M using chloride salts of various alkali metals as well as mixed concentrated salt solutions (involving both mono- and divalent cations and anions), that mimic concentrated brines typically encountered in geological formations. Experiments were carried out using flat substrates and submicrometer-sized colloidal probes made of smooth oxidized silicon immersed in salt solutions at pH values of 6 and 9 and temperatures of 25 °C and 45 °C. While strong repulsive forces were observed for the smallest tip-sample separations, none of the conditions explored displayed any indication of anomalous long range electrostatic forces as reported for macroscopic mica surfaces. Instead, forces are universally dominated by attractive van der Waals interactions at tip-sample separations of $\approx$2nm and beyond for salt concentrations of 1 M and higher. Complementary calculations based on classical density functional theory for the primitive model support these experimental observations and display a consistent decrease in screening length with increasing ion concentration.

preprint2022arXiv

Direct numerical simulations of the modified Poisson-Nernst-Planck equations for the charging dynamics of cylindrical electrolyte-filled pores

Understanding how electrolyte-filled porous electrodes respond to an applied potential is important to many electrochemical technologies. Here, we consider a model supercapacitor of two blocking cylindrical pores on either side of a cylindrical electrolyte reservoir. A stepwise potential difference $2Φ$ between the pores drives ionic fluxes in the setup, which we study through the modified Poisson-Nernst-Planck equations, solved with finite elements. We focus our discussion on the dominant timescales with which the pores charge and how these timescales depend on three dimensionless numbers. Next to the dimensionless applied potential $Φ$, we consider the ratio $R/R_b$ of the pore's resistance $R$ to the bulk reservoir resistance $R_b$ and the ratio $r_{p}/λ$ of the pore radius $r_p$ to the Debye length $λ$. We compare our data to theoretical predictions by Aslyamov and Janssen ($Φ$), Posey and Morozumi ($R/R_b$), and Henrique, Zuk, and Gupta ($r_{p}/λ$). Through our numerical approach, we delineate the validity of these theories and the assumptions on which they were based.

preprint2022arXiv

Dynamical steady-states of active colloids interacting via chemical fields

We study the dynamical steady-states of a monolayer of chemically active self-phoretic colloids as a function of packing fraction and self-propulsion speed by means of Brownian dynamics simulations. We focus on the case that a chemical field induces competing attractive positional and repulsive orientational interactions. Analyzing the distribution of cluster size and local density as well as the hexatic order parameter, we distinguish four distinct dynamical states which include collapsed, active gas, dynamical clustering, and motility-induced phase-separated states. The long-range chemical field-induced interactions shift the onset of motility-induced phase separation (MIPS) to very low packing fractions at intermediate self-propulsion speeds. We also find that the fraction of particles in the largest clusters is a suitable order parameter characterizing the dynamical phase transitions from an active gas or dynamical clustering steady-state to a phase-separated state upon increase of the packing fraction. The order parameter changes discontinuously when going from an active gas to a MIPS-like state at intermediate self-propulsion speeds, whereas it changes continuously at larger activities where the system undergoes a transition from a dynamical clustering state to MIPS-like state.

preprint2021arXiv

Statistical-physics-inspired model for intrinsic fluctuations driving supply and demand in markets

We propose a simple statistical-physics-inspired model for the effect of intrinsic fluctuations on supply and demand in markets. The model consists of agents that trade in two types of goods of which the total number is separately conserved. The relative preference of an individual agent for the two types of goods is determined by a utility that is identical for all agents. Market supply and demand curves are computed and compared for various motivated choices of the distribution of goods over the agents. In particular, we compare the "mean-field" case, in which all agents have the same number of goods and that is akin to the economics textbook case, to the case of Boltzmann-Gibbs distributed goods, in which agents have a fluctuating number of goods. We find that the resulting equilibrium prices are not equal for these two approaches, especially when a large fraction of the agents can neither buy nor sell.

preprint2020arXiv

A Blessing and a Curse: How a Supercapacitor's Large Capacitance Causes its Slow Charging

The development of novel electrolytes and electrodes for supercapacitors is hindered by a gap of several orders of magnitude between experimentally measured and theoretically predicted charging timescales. Here, we propose an electrode model, containing many parallel stacked electrodes, that explains the slow charging dynamics of supercapacitors. At low applied potentials, the charging behavior of this model is described well by an equivalent circuit model. Conversely, at high potentials, charging dynamics slow down and evolve on two relaxation time scales: a generalized $RC$ time and a diffusion time, which, interestingly, become similar for porous electrodes. The charging behavior of the stack-electrode model presented here helps to understand the charging dynamics of porous electrodes and qualitatively agrees with experimental time scales measured with porous electrodes.

preprint2016arXiv

A Landau-de Gennes theory for hard colloidal rods: defects and tactoids

We construct a phenomenological Landau-de Gennes theory for hard colloidal rods by performing an order parameter expansion of the chemical-potential dependent grand potential. By fitting the coefficients to known results of Onsager theory, we are not only able to describe the isotropic-nematic phase transition as function of density, including the well-known density jump, but also the isotropic-nematic planar interface. The resulting theory is applied in calculations of the isotropic core size in a radial hedgehog defect, the density dependence of linear defects of hard rods in square confinement, and the formation of a nematic droplet in an isotropic background.

preprint2016arXiv

Alternating strings and clusters in suspensions of charged colloids

We report the formation of alternating strings and clusters in a binary suspension of repulsive charged colloids with double layers larger than the particle size. Within a binary cell model we include many-body and charge-regulation effects under the assumption of a constant surface potential, and consider their repercussions on the two-particle interaction potential. We find that the formation of induced dipoles close to a charge-reversed state may explain the formation of these structures. Finally, we will touch upon the formation of dumbbells and small clusters in a one-component system, where the effective electrostatic interaction is always repulsive.

preprint2016arXiv

Density-induced reentrant melting of colloidal Wigner crystals

Electrostatic repulsions can drive crystallization in many-particle systems. For charged colloidal systems, the phase boundaries as well as crystal structure are highly tunable by experimental parameters such as salt concentration and pH. By using projections of the colloid-ion mixture to a system of (soft) repulsive spheres and the one-component plasma (OCP), we study the hitherto unexplained experimentally observed reentrant melting of electrostatically repelling colloids upon increasing the colloid density. Our study shows that the surface chemistry should involve a competition between adsorption of cations and anions to explain the observed density-induced reentrant melting.

preprint2016arXiv

Harvesting vibrational energy with liquid-bridged electrodes: thermodynamics in mechanically and electrically driven RC-circuits

We theoretically study a vibrating pair of parallel electrodes bridged by a (deformed) liquid droplet, which is a recently developed microfluidic device to harvest vibrational energy. The device can operate with various liquids, including liquid metals, electrolytes, as well as ionic liquids. We numerically solve the Young-Laplace equation for all droplet shapes during a vibration period, from which the time-dependent capacitance follows that serves as input for an equivalent circuit model. We first investigate two existing energy harvesters (with a constant and a vanishing bias potential), for which we explain an open issue related to their optimal electrode separations, which is as small as possible or as large as possible in the two cases, respectively. Then we propose a new engine with a time-dependent bias voltage, with which the harvested work and the power can be increased by orders of magnitude at low vibration frequencies and by factors 2-5 at high frequencies, where frequencies are to be compared to the inverse RC-time of the circuit.

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.

preprint2016arXiv

Sedimentation stacking diagrams of binary mixtures of thick and thin hard rods

We use Onsager theory and the local density approximation to study sedimentation-diffusion equilibrium density profiles of binary mixtures of thick and thin hard rods. We construct stacking diagrams for three diameter ratios, and find that even a simple spindle-shaped phase diagram with only isotropic-nematic demixing can lead to counter-intuitive stacking sequences such as an isotropic phase sandwiched between two nematic phases. For the most complex phase diagram considered here, we find sixteen distinct stacking sequences, including several with five sedimented layers. By adding sedimentation paths to composition-pressure and density-density phase diagrams and calculating density and composition profiles, we show that conclusions about bulk phase diagrams of binary mixtures on the basis of sedimentation-diffusion equilibria should be drawn warily.

preprint2016arXiv

Tuning colloid-interface interactions by salt partitioning

We show that the interaction of an oil-dispersed colloidal particle with an oil-water interface is highly tunable from attractive to repulsive, either by varying the sign of the colloidal charge via charge regulation, or by varying the difference in hydrophilicity between the dissolved cations and anions. In addition, we investigate the yet unexplored interplay between the self-regulated colloidal surface charge distribution with the planar double layer across the oil-water interface and the spherical one around the colloid. Our findings explain recent experiments and have direct relevance for tunable Pickering emulsions.

preprint2015arXiv

Boosting capacitive blue-energy and desalination devices with waste heat

We show that sustainably harvesting 'blue' energy from the spontaneous mixing process of fresh and salty water can be boosted by varying the water temperature during a capacitive mixing process. Our modified Poisson-Boltzmann calculations predict a strong temperature dependence of the electrostatic potential of a charged electrode in contact with an adjacent aqueous 1:1 electrolyte. We propose to exploit this dependence to boost the efficiency of capacitive blue engines, which are based on cyclically charging and discharging nanoporous supercapacitors immersed in salty and fresh water, respectively [D. Brogioli, Phys. Rev. Lett. 103, 058501 (2009)]. We show that the energy output of blue engines can be increased by a factor of order two if warm (waste-heated) fresh water is mixed with cold sea water. Moreover, the underlying physics can also be used to optimize the reverse process of capacitive desalination of water.

preprint2015arXiv

Fundamental measure theory for the electric double layer: implications for blue-energy harvesting and water desalination

Capacitive mixing (CAPMIX) and capacitive deionization (CDI) are promising candidates for harvesting clean, renewable energy and for the energy efficient production of potable water, respectively. Both CAPMIX and CDI involve water-immersed porous carbon (supercapacitors) electrodes at voltages of the order of hundreds of millivolts, such that counter-ionic packing is important for the electric double layer (EDL) which forms near the surface of these porous materials. Thus, we propose a density functional theory (DFT) to model the EDL, where the White-Bear mark II fundamental measure theory functional is combined with a mean-field Coulombic and a mean spherical approximation-type correction to describe the interplay between dense packing and electrostatics, in good agreement with molecular dynamics simulations. We discuss the concentration-dependent potential rise due to changes in the chemical potential in capacitors in the context of an over-ideal theoretical description and its impact on energy harvesting and water desalination. Compared to less elaborate mean-field models our DFT calculations reveal a higher work output for blue-energy cycles and a higher energy demand for desalination cycles.

preprint2015arXiv

Self-propulsion mechanism of active Janus particles in near-critical binary mixtures

Gold-capped Janus particles immersed in a near-critical binary mixture can be propelled using illumination. We employ a non-isothermal diffuse interface approach to investigate the self-propulsion mechanism of a single colloid. We attribute the motion to body forces at the edges of a micronsized droplet that nucleates around the particle. Thus, the often-used concept of a surface velocity cannot account for the self-propulsion. The particle's swimming velocity is related to the droplet shape and size, which is determined by a so-called critical isotherm. Two distinct swimming regimes exist, depending on whether the droplet partially or completely covers the particle. Interestingly, the dependence of the swimming velocity on temperature is non-monotonic in both regimes.

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.

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.

preprint2014arXiv

The electric double layer has a life of its own

Using molecular dynamics simulations with recently developed importance sampling methods, we show that the differential capacitance of a model ionic liquid based double-layer capacitor exhibits an anomalous dependence on the applied electrical potential. Such behavior is qualitatively incompatible with standard mean-field theories of the electrical double layer, but is consistent with observations made in experiment. The anomalous response results from structural changes induced in the interfacial region of the ionic liquid as it develops a charge density to screen the charge induced on the electrode surface. These structural changes are strongly influenced by the out-of-plane layering of the electrolyte and are multifaceted, including an abrupt local ordering of the ions adsorbed in the plane of the electrode surface, reorientation of molecular ions, and the spontaneous exchange of ions between different layers of the electrolyte close to the electrode surface. The local ordering exhibits signatures of a first-order phase transition, which would indicate a singular charge-density transition in a macroscopic limit.

preprint2013arXiv

Charge fluctuations in nano-scale capacitors

The fluctuations of the charge on an electrode contain information on the microscopic correlations within the adjacent fluid and their effect on the electronic properties of the interface. We investigate these fluctuations using molecular dynamics simulations in a constant-potential ensemble with histogram reweighting techniques. This approach offers in particular an efficient, accurate and physically insightful route to the differential capacitance that is broadly applicable. We demonstrate these methods with three different capacitors: pure water between platinum electrodes, and a pure as well as a solvent-based organic electrolyte each between graphite electrodes. The total charge distributions with the pure solvent and solvent-based electrolytes are remarkably Gaussian, while in the pure ionic liquid the total charge distribution displays distinct non-Gaussian features, suggesting significant potential-driven changes in the organization of the interfacial fluid.

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.

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.

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.

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.

preprint2012arXiv

Statistical thermodynamics of supercapacitors and blue engines

We study the thermodynamics of electrode-electrolyte systems, for instance supercapacitors filled with an ionic liquid or blue-energy devices filled with river- or sea water. By a suitable mapping of thermodynamic variables, we identify a strong analogy with classical heat engines. We introduce several Legendre transformations and Maxwell relations. We argue that one should distinguish between the differential capacity at constant ion number and at constant ion chemical potential, and derive a relation between them that resembles the standard relation between heat capacity at constant volume and constant pressure. Finally, we consider the probability distribution of the electrode charge at a given electrode potential, the standard deviation of which is given by the differential capacity.

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.

preprint2011arXiv

DNA translocation through nanopores with salt gradients: The role of osmotic flow

Recent experiments of translocation of double stranded DNA through nanopores [M. Wanunu \textit{et al.} Nature Nanotech. {\bf 5}, 160 (2010)] reveal that the DNA capture rate can be significantly influenced by a salt gradient across the pore. We show that osmotic flow combined with electrophoretic effects can quantitatively explain the experimental data on the salt-gradient dependence of the capture rate.

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.

preprint2010arXiv

Charge regulation and ionic screening of patchy surfaces

The properties of surfaces with charge-regulated patches are studied using non-linear Poisson-Boltzmann theory. Using a mode expansion to solve the non-linear problem efficiently, we reveal the charging behaviour of Debye-length sized patches. We find that patches charge up to higher charge densities if their size is relatively small and if the patches are well separated. The numerical results are used to construct a basic analytical model which predicts the average surface charge density on surfaces with patchy chargeable groups.

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.

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.

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.

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.

preprint2009arXiv

Inflation of the screening length induced by Bjerrum pairs

Within a modified Poisson-Boltzmann theory we study the effect of Bjerrum pairs on the typical length scale $1/\barκ$ over which electric fields are screened in electrolyte solutions, taking into account a simple association-dissociation equilibrium between free ions and Bjerrum pairs. At low densities of Bjerrum pairs, this length scale is well approximated by the Debye length $1/κ\propto 1/\sqrt{ρ_\mathrm{s}}$, with $ρ_\mathrm{s}$ the free ion density. At high densities of Bjerrum pairs, however, we find $1/\barκ\propto \sqrt{ρ_\mathrm{s}}$ which is significantly larger than $1/κ$ due to the enhanced effective permittivity of the electrolyte, caused by the polarization of Bjerrum pairs. We argue that this mechanism may explain the recently observed anomalously large colloid-free zones between an oil-dispersed colloidal crystal and a colloidal monolayer at the oil-water interface.

preprint2009arXiv

Poisson-Boltzmann cell model for heterogeneously charged colloids

We introduce the Poisson-Boltzmann cell model for spherical colloidal particles with a heterogeneous surface charge distribution. This model is obtained by generalizing existing cell models for mixtures of homogeneously charged colloidal spheres. Our new model has similar features as Onsager's second-virial theory for liquid crystals, but it predicts no orientational ordering if there is no positional ordering. This implies that all phases of heterogeneously charged colloids that are liquid-like with respect to translational degrees of freedom are also isotropic with respect to particle orientation.

preprint2008arXiv

Stability of additive-free water-in-oil emulsions

We calculate ion distributions near a planar oil-water interface within non-linear Poisson-Boltzmann theory, taking into account the Born self-energy of the ions in the two media. For unequal self-energies of cations and anions, a spontaneous charge separation is found such that the water and oil phase become oppositely charged, in slabs with a typical thickness of the Debye screening length in the two media. From the analytical solutions, the corresponding interfacial charge density and the contribution to the interfacial tension is derived, together with an estimate for the Yukawa-potential between two spherical water droplets in oil. The parameter regime is explored where the plasma coupling parameter exceeds the crystallization threshold, i.e. where the droplets are expected to form crystalline structures due to a strong Yukawa repulsion, as recently observed experimentally. Extensions of the theory that we discuss briefly include numerical calculations on spherical water droplets in oil, and analytical calculations of the linear PB-equation for a finite oil-water interfacial width.

preprint2007arXiv

Destructed double-layer and ionic charge separation near the oil-water interface

We study suspensions of hydrophobic charged colloidal spheres dispersed in a demixed oil-water mixture by means of a modified Poisson-Boltzmann theory, taking into account image charge effects and partitioning of the monovalent ions. We find that the ion's aversion for oil can destroy the double layers of the oil-dispersed colloids. This affects the salt-concentration dependence of the colloidal adsorption to the oil-water interface qualitatively. The theory also predicts a narrow range of the oil-dielectric constant in which micron-sized water-in-oil droplets acquire enough charge to crystallize at volume fractions as small as $\sim 10^{-3}$ in the absence of colloids. These findings explain recent observations [M.E. Leunissen {\em et al.}, Proc. Nat. Ac. Sci {\bf 104}, 2585 (2007)].