Search

The non-equilibrium steady states of a semi-infinite quasi-one-dimensional univalent binary electrolyte solution, characterised by non-vanishing electric currents, are investigated by means of Poisson-Nernst-Planck (PNP) theory. Exact analytical expressions of the electric field, the charge density and the number density are derived, which depend on the electric current density as a parameter. From a non-equilibrium version of the Grahame equation, which relates the total space charge per cross-sectional area and the corresponding contribution of the electric potential drop, the current-dependent differential capacitance of the diffuse layer is derived. In the limit of vanishing electric current these results reduce to those within Gouy-Chapman theory. It is shown that improperly chosen boundary conditions lead to non-equilibrium steady state solutions of the PNP equations with negative ion number densities. A necessary and sufficient criterion on surface conductivity constitutive relations is formulated which allows one to detect such unphysical solutions., Comment: Corrected d-symbol in Eqs. (31)-(33)

The structures of dilute electrolyte solutions close to non-uniformly charged planar substrates are systematically studied within the entire spectrum of microscopic to macroscopic length scales by means of a unified classical density functional theory (DFT) approach. This is in contrast to previous investigations, which are applicable either to short or to long length scales. It turns out that interactions with microscopic ranges, e.g., due to the hard cores of the fluid molecules and ions, have negligible influence on the formation of non-uniform lateral structures of the electrolyte solutions. This partly justifies the Debye-H\"uckel approximation schemes applied in previous studies of that system. In general, a coupling between the lateral and the normal fluid structures leads to the phenomenology that, upon increasing the distance from the substrate, less details of the lateral non-uniformities contribute to the fluid structure, such that ultimately only large-scale surface features remain relevant. It can be expected that this picture also applies to other fluids characterized by several length scales.

Previous theoretical studies of calamitic (i.e., rod-like) ionic liquid crystals (ILCs) based on an effective one-species model led to indications of a novel smectic-A phase with a layer spacing being much larger than the length of the mesogenic (i.e., liquid-crystal forming) ions. In order to rule out the possibility that this wide smectic-A phase is merely an artifact caused by the one-species approximation, we investigate an extension which accounts explicitly for cations and anions in ILCs. Our present findings, obtained by grand canonical Monte Carlo simulations, show that the phase transitions between the isotropic and the smectic-A phases of the cation-anion system are in qualitative agreement with the effective one-species model used in the preceding studies. In particular, for ILCs with mesogenes (i.e., liquid-crystal forming species) carrying charged sites at their tips, the wide smectic-A phase forms, at low temperatures and within an intermediate density range, in between the isotropic and a hexagonal crystal phase. We find that in the ordinary smectic-A phase the spatial distribution of the counterions of the mesogens is approximately uniform, whereas in the wide smectic-A phase the small counterions accumulate in between the smectic layers. Due to this phenomenology the wide smectic-A phase could be interesting for applications which hinge on the presence of conductivity channels for mobile ions.

Liquid-liquid phase separation has emerged as one of the important paradigms in the chemical physics as well as biophysics of charged macromolecular systems. We elucidate an equilibrium phase separation mechanism based on charge regulation, i.e., protonation-deprotonation equilibria controlled by pH, in an idealized macroion system which can serve as a proxy for simple coacervation. First, a low-density density-functional calculation reveals the dominance of two-particle configurations coupled by ion adsorption on neighboring macroions. Then a binary cell model, solved on the Debye-H\"uckel as well as the full nonlinear Poisson-Boltzmann level, unveils the charge-symmetry breaking as inducing the phase separation between low- and high-density phases as a function of pH. These results can be identified as a charge symmetry broken complex coacervation between chemically identical macroions., Comment: 11 pages, 7 figures

The structure of dilute electrolyte solutions close to a surface carrying a spatially inhomogeneous surface charge distribution is investigated by means of classical density functional theory (DFT) within the approach of fundamental measure theory (FMT). For electrolyte solutions the influence of these inhomogeneities is particularly strong because the corresponding characteristic length scale is the Debye length, which is large compared to molecular sizes. Here a fully three-dimensional investigation is performed, which accounts explicitly for the solvent particles, and thus provides insight into effects caused by ion-solvent coupling. The present study introduces a versatile framework to analyze a broad range of types of surface charge heterogeneities even beyond the linear response regime. This reveals a sensitive dependence of the number density profiles of the fluid components and of the electrostatic potential on the magnitude of the charge as well as on details of the surface charge patterns at small scales.

Motivated by biological membrane-containing organelles in plants and photosynthetic bacteria, we study charge regulation in a model membrane stack. Considering (de)protonation as the simplest mechanism of charge equilibration between the membranes and with the bathing environment, we uncover a symmetry-broken charge state in the stack with a quasiperiodic effective charge sequence. In the case of a monovalent bathing salt solution our model predicts complex, inhomogeneous charge equilibria depending on the strength of the (de)protonation reaction, salt concentration, intermembrane separation, and their number in the stack. Our results shed light on the basic reorganization mechanism of thylakoid membrane stacks., Comment: 16 pages, 7 figures

We study the transient response of an electrolytic cell subject to a small, suddenly applied temperature increase at one of its two bounding electrode surfaces. An inhomogeneous temperature profile then develops, causing, via the Soret effect, ionic rearrangements towards a state of polarized ionic charge density $q$ and local salt density $c$. For the case of equal cationic and anionic diffusivities, we derive analytical approximations to $q, c$, and the thermovoltage $V_{T}$ for early ($t\ll\tau_{T}$) and late ($t\gg\tau_{T}$) times as compared to the relaxation time $\tau_{T}$ of the temperature. We challenge the conventional wisdom that the typically large Lewis number, the ratio $a/D$ of thermal to ionic diffusivities, of most liquids implies a quickly reached steady-state temperature profile onto which ions relax slowly. Though true for the evolution of $c$, it turns out that $q$ (and $V_{T}$) can respond much faster. Particularly when the cell is much bigger than the Debye length, a significant portion of the transient response of the cell falls in the $t\ll\tau_{T}$ regime, for which our approximated $q$ (corroborated by numerics) exhibits a density wave that has not been discussed before in this context. For electrolytes with unequal ionic diffusivities, $V_{T}$ exhibits a two-step relaxation process, in agreement with experimental data of Bonetti et al. [J. Chem. Phys. 142, 244708 (2015)]., Comment: 13 pages, 5 figures

Ionic liquid crystals (ILCs) are anisotropic mesogenic molecules which additionally carry charges. This combination gives rise to a complex interplay of the underlying (anisotropic) contributions to the pair interactions. It promises interesting and distinctive structural and orientational properties to arise in systems of ILCs, combining properties of liquid crystals and ionic liquids. While previous theoretical studies have focused on the phase behavior of ILCs and the structure of the respective bulk phases, in the present study we provide new results, obtained within density functional theory, concerning (planar) free interfaces between an isotropic liquid $L$ and two types of smectic-A phases ($S_A$ or $S_{AW}$). We discuss the structural and orientational properties of these interfaces in terms of the packing fraction profile $\eta(\vec{r})$ and the orientational order parameter profile $S_2(\vec{r})$ concerning the tilt angle $\alpha$ between the (bulk) smectic layer normal and the interface normal. The asymptotic decay of $\eta(\vec{r})$ and of $S_2(\vec{r})$ towards their values in the isotropic bulk is discussed, too.

The behavior of a uniformly magnetized ferronematic slab is investigated numerically in a situation in which an external magnetic field is applied parallel and antiparallel, respectively, to its initial magnetization direction. The employed numerical method allows one to determine hysteresis curves, from which a critical magnetic field strength (i.e., the one at which the ferronematic sample becomes distorted) as function of the system parameters can be inferred. Two possible mechanisms of switching the magnetization by applying a magnetic field in the antiparallel direction are observed and characterized in terms of the coupling constant between the magnetization and the nematic director as well as in terms of the coupling strength of the nematic liquid crystal and the walls of the slab. Suitably prepared walls allow one to combine both switching mechanisms in one setup, such that one can construct a cell the magnetization of which can be reversibly switched off.

The electrostatic interaction between pairs of spherical or macroscopically long, parallel cylindrical colloids trapped at fluid interfaces is studied theoretically for the case of small inter-particle separations. Starting from the effective interaction between two planar walls and by using the Derjaguin approximation, we address the issue of how the electrostatic interaction between such particles is influenced by their curvatures and by the wetting contact angle at their surfaces. Regarding the influence of curvature, our findings suggest that the discrepancies between linear and nonlinear Poisson-Boltzmann theory, which have been noticed before for planar walls, also occur for spheres and macroscopically long, parallel cylinders, though their magnitude depends on the wetting contact angle. Concerning the influence of the wetting contact angle $\theta$ simple relations are obtained for equally sized particles which indicate that the inter-particle force varies significantly with $\theta$ only within an interval around $90^\circ$. This interval depends on the Debye length of the fluids and on the size of the particles but not on their shape. For unequally sized particles, a more complicated relation is obtained for the variation of the inter-particle force with the wetting contact angle., Comment: 9 pages, 4 figures

An effective Landau-like description of ferronematics, i.e., suspensions of magnetic colloidal particles in a nematic liquid crystal (NLC), is developed in terms of the corresponding magnetization and nematic director fields. The study is based on a microscopic model and on classical density functional theory. Ferronematics are susceptible to weak magnetic fields and they can exhibit a ferromagnetic phase, which has been predicted several decades ago and which has recently been found experimentally. Within the proposed effective Landau theory of ferronematics one has quantitative access, e.g., to the coupling between the magnetization of the magnetic colloids and the nematic director of the NLC. On mesoscopic length scales this generates complex response patterns.

The influence of a chemically or electrically heterogeneous distribution of interaction sites at a planar substrate on the number density of an adjacent fluid is studied by means of classical density functional theory (DFT). In the case of electrolyte solutions the effect of this heterogeneity is particularly long ranged, because the corresponding relevant length scale is set by the Debye length which is large compared to molecular sizes. The DFT used here takes the solvent particles explicitly into account and thus captures phenomena, inter alia, due to ion-solvent coupling. The present approach provides closed analytic expressions describing the influence of chemically and electrically nonuniform walls. The analysis of isolated \delta-like interactions, isolated interaction patches, and hexagonal periodic distributions of interaction sites reveals a sensitive dependence of the fluid density profiles on the type of the interaction, as well as on the size and the lateral distribution of the interaction sites.

We revisit a classical problem of theoretical electrochemistry: the response of an electric double layer capacitor (EDLC) subject to a small, suddenly applied external potential. We solve the Debye-Falkenhagen equation to obtain exact expressions for key EDLC quantities: the ionic charge density, the ionic current density, and the electric field. In contrast to earlier works, our results are not restricted to the long-time asymptotics of those quantities. The solutions take the form of infinite sums whose successive terms all decay exponentially with increasingly short relaxation times. Importantly, this set of relaxation times is the same among all aforementioned EDLC quantities; this property is demanded on physical grounds but not generally achieved within approximation schemes. The scaling of the largest relaxation timescale $\tau_{1}$, that determines the long-time decay, is in accordance with earlier results: depending on the Debye length, $\lambda_{D}$, and the electrode separation, $2L$, it amounts to $\tau_{1}\simeq\lambda_{D} L/D$ for $L\gg\lambda_{D}$, and $\tau_{1} \simeq 4 L^2/(\pi^2 D)$ for $L\ll\lambda_{D}$, respectively (with $D$ being the ionic diffusivity)., Comment: 12 pages, 4 figures

The electrostatic interaction between two non-identical, moderately charged colloids situated in close proximity of each other at a fluid interface is studied. By resorting to a well-justified model system, this problem is analytically solved within the framework of linearized Poisson-Boltzmann (PB) density functional theory (DFT). The resulting interaction comprises a surface and a line part, both of which, as functions of the inter-particle separation, show a rich behavior including monotonic as well as non-monotonic variations. In almost all cases, these variations cannot be captured correctly by using the superposition approximation. Moreover, expressions for the surface tensions, the line tensions and the fluid-fluid interfacial tension, which are all independent of the inter-particle separation, are obtained. Our results are expected to be particularly useful for emulsions stabilized by oppositely charged particles., Comment: 12 pages, 5 figures

We derive a general closed expression for the local pressure exerted onto the corrugated walls of a channel confining a fluid medium. When the fluid medium is at equilibrium the local pressure is a functional of the shape of the walls. It is shown that, due to the intrinsic non-local character of the interactions among the particles forming the fluid, the applicability of approximate schemes such as the concept of a surface of tension or morphometric thermodynamics is limited to wall curvatures small compared to the range of particle-particle interactions., Comment: corrections of typos of the previous version and reorganization of part of the material in an additional appendix

The interaction between two chemically identical charge-regulated surfaces is studied using the classical density functional theory. In contrast to common expectations and assumptions, under certain realistic conditions we find a spontaneous emergence of disparate charge densities on the two surfaces. The surface charge densities can differ not only in their magnitude, but quite unexpectedly, even in their sign, implying that the electrostatic interaction between the two chemically identical surfaces can be attractive instead of repulsive. Moreover, an initial symmetry with equal charge densities on both surfaces can also be broken spontaneously upon decreasing the separation between the two surfaces. The origin of this phenomenon is a competition between the adsorption of ions from the solution to the surface and the interaction between the adsorbed ions already on the surface.These findings are fundamental for the understanding of the forces between colloidal objects and, in particular, they are bound to strongly influence the present picture of protein interaction., Comment: 19 pages, 5 figures

Ionic liquid crystals (ILCs) are anisotropic mesogenic molecules which carry charges and therefore combine properties of liquid crystals, e.g., the formation of mesophases, and of ionic liquids, such as low melting temperatures and tiny triple-point pressures. Previous density functional calculations have revealed that the phase behavior of ILCs is strongly affected by their molecular properties, i.e., their aspect ratio, the loci of the charges, and their interaction strengths. Here, we report new findings concerning the phase behavior of ILCs as obtained by density functional theory and Monte Carlo simulations. The most important result is the occurrence of a novel, wide smectic-A phase $S_{AW}$, at low temperature, the layer spacing of which is larger than that of the ordinary high-temperature smectic-A phase $S_{A}$. Unlike the ordinary smectic $S_A$ phase, the structure of the $S_{AW}$ phase consists of alternating layers of particles oriented parallel to the layer normal and oriented perpendicular to it.

The wetting of a charged wedge-like wall by an electrolyte solution is investigated by means of classical density functional theory. As in other studies on wedge wetting, this geometry is considered as the most simple deviation from a planar substrate, and it serves as a first step towards more complex confinements of fluids. By focusing on fluids containing ions and surface charges, features of real systems are covered which are not accessible within the vast majority of previous theoretical studies concentrating on simple fluids in contact with uncharged wedges. In particular, the filling transition of charged wedges is necessarily of first order, because wetting transitions of charged substrates are of first order and the barrier in the effective interface potential persists below the wetting transition of a planar wall; hence, critical filling transitions are not expected to occur for ionic systems. The dependence of the critical opening angle on the surface charge, as well as the dependence of the filling height, of the wedge adsorption, and of the line tension on the opening angle and on the surface charge are analyzed in detail.

The salt-induced microheterogeneity (MH) formation in binary liquid mixtures is studied by small-angle X-ray scattering (SAXS) and liquid state theory. Previous experiments have shown that this phenomenon occurs for antagonistic salts, whose cations and anions prefer different components of the solvent mixture. However, so far the precise mechanism leading to the characteristic length scale of MHs remained unclear. Here, it is shown that MHs can be generated by the competition of short-ranged interactions and long-ranged monopole-dipole interactions. The experimental SAXS patterns can be quantitatively reproduced by fitting to the derived correlation functions without assuming any specific model. The dependency of the MH structure with respect to ionic strength and temperature is analyzed. Close to the demixing phase transition, critical-like behavior occurs with respect to the spinodal line in the phase diagram.

Density functional theory is used to describe electrolyte solutions in contact with electrodes of planar or spherical shape. For the electrolyte solutions we consider the so-called civilized model, in which all species present are treated on equal footing. This allows us to discuss the features of the electric double layer in terms of the differential capacitance. The model provides insight into the microscopic structure of the electric double layer, which goes beyond the mesoscopic approach studied in the accompanying paper. This enables us to judge the relevance of microscopic details, such as the radii of the particles forming the electrolyte solutions or the dipolar character of the solvent particles, and to compare the predictions of various models. Similar to the preceding paper, a general behavior is observed for small radii of the electrode in that in this limit the results become independent of the surface charge density and of the particle radii. However, for large electrode radii non-trivial behaviors are observed. Especially the particle radii and the surface charge density strongly influence the capacitance. From the comparison with the Poisson-Boltzmann approach it becomes apparent that the shape of the electrode determines whether the microscopic details of the full civilized model have to be taken into account or whether already simpler models yield acceptable predictions., Comment: accepted for publication in the Journal of Chemical Physics

The effective electrostatic interaction between a pair of colloids, both of them located close to each other at an electrolyte interface, is studied by employing the full, nonlinear Poisson-Boltzmann (PB) theory within classical density functional theory. Using a simplified yet appropriate model, all contributions to the effective interaction are obtained exactly, albeit numerically. The comparison between our results and those obtained within linearized PB theory reveals that the latter overestimates these contributions significantly at short inter-particle separations. Whereas the surface contributions to the linear and the nonlinear PB results differ only quantitatively, the line contributions show qualitative differences at short separations. Moreover, a dependence of the line contribution on the solvation properties of the two adjacent fluids is found, which is absent within the linear theory. Our results are expected to enrich the understanding of effective interfacial interactions between colloids.

The concept of a non-equilibrium interfacial tension, defined via the work required to deform the system such that the interfacial area is changed while the volume is conserved, is investigated theoretically in the context of the relaxation of an initial perturbation of a colloidal fluid towards the equilibrium state. The corresponding general formalism is derived for systems with planar symmetry and applied to fluid models of colloidal suspensions and polymer solutions. It is shown that the non-equilibrium interfacial tension is not necessarily positive, that negative non-equilibrium interfacial tensions are consistent with strictly positive equilibrium interfacial tensions, and that the sign of the interfacial tension can influence the morphology of density perturbations during relaxation.

The three-phase contact line formed by the intersection of a liquid-vapor interface of an electrolyte solution with a charged planar substrate is studied in terms of classical density functional theory applied to a lattice model. The influence of the substrate charge density and of the ionic strength of the solution on the intrinsic structure of the three-phase contact line and on the corresponding line tension is analyzed. We find a negative line tension for all values of the surface charge density and of the ionic strength considered. The strength of the line tension decreases upon decreasing the contact angle via varying either the temperature or the substrate charge density.

The static dielectric properties of dense ionic fluids, e.g., room temperature ionic liquids (RTILs) and inorganic fused salts, are investigated on different length scales by means of grandcanonical Monte Carlo simulations. A generally applicable scheme is developed which allows one to approximately decompose the electric susceptibility of dense ionic fluids into the orientation and the distortion polarization contribution. It is shown that at long range the well-known plasma-like perfect screening behavior occurs, which corresponds to a diverging distortion susceptibility, whereas at short range orientation polarization dominates, which coincides with that of a dipolar fluid of attached cation-anion pairs. This observation suggests that the recently debated interpretation of RTILs as dilute electrolyte solutions might not be simply a yes-no-question but it might depend on the considered length scale.

The electrostatic interaction in ionic fluids is well-known to give rise to a characteristic phase behavior and structure. Sometimes its long range is proposed to single out the electrostatic potential over other interactions with shorter ranges. Here the importance of the range for the phase behavior and the structure of ionic fluids is investigated by means of grandcanonical Monte Carlo simulations of the lattice restricted primitive model (LRPM). The long-ranged electrostatic interaction is compared to various types of short-ranged potentials obtained by sharp and/or smooth cut-off schemes. Sharply cut off electrostatic potentials are found to lead to a strong dependence of the phase behavior and the structure on the cut-off radius. However, when combined with a suitable additional smooth cut-off, the short-ranged LRPM is found to exhibit quantitatively the same phase behavior and structure as the conventional long-ranged LRPM. Moreover, the Stillinger-Lovett perfect screening property, which is well-known to be generated by the long-ranged electrostatic potential, is also fulfilled by short-ranged LRPMs with smooth cut-offs. By showing that the characteristic phase behavior and structure of ionic fluids can also be found in systems with short-ranged potentials, one can conclude that the decisive property of the electrostatic potential in ionic fluids is not the long range but rather the valency dependence.

Several model fluids in contact with planar, spherical, and cylindrical walls are investigated for small number densities within density functional theory. The dependence of the solid-fluid interfacial tension on the curvature of spherical and cylindrical walls is examined and compared with the corresponding expression derived within the framework of morphometric thermodynamics. Particular attention is paid to the implications of the choice of the interface location, which underlies the definition of the interfacial tension. We find that morphometric thermodynamics is never exact for the considered systems and that its quality as an approximation depends sensitively on the choice of the interface location.

The electrostatic interaction between colloidal particles trapped at the interface between two immiscible electrolyte solutions is studied in the limit of small inter-particle distances. Within an appropriate model exact analytic expressions for the electrostatic potential as well as for the surface and line interaction energies are obtained. They demonstrate that the widely used superposition approximation, which is commonly applied to large distances between the colloidal particles, fails qualitatively at small distances and is quantitatively unreliable even at large distances. Our results contribute to an improved description of the interaction between colloidal particles trapped at fluid interfaces., Comment: Submitted

For wetting films in dilute electrolyte solutions close to charged walls we present analytic expressions for their effective interface potentials. The analysis of these expressions renders the conditions under which corresponding wetting transitions can be first- or second-order. Within mean field theory we consider two models, one with short- and one with long-ranged solvent-solvent and solvent-wall interactions. The analytic results reveal in a transparent way that wetting transitions in electrolyte solutions, which occur far away from their critical point (i.e., the bulk correlation length is less than half of the Debye length) are always first-order if the solvent-solvent and solvent-wall interactions are short-ranged. In contrast, wetting transitions close to the bulk critical point of the solvent (i.e., the bulk correlation length is larger than the Debye length) exhibit the same wetting behavior as the pure, i.e., salt-free, solvent. If the salt-free solvent is governed by long-ranged solvent-solvent as well as long-ranged solvent-wall interactions and exhibits critical wetting, adding salt can cause the occurrence of an ion-induced first-order thin-thick transition which precedes the subsequent continuous wetting as for the salt-free solvent., Comment: Submitted

The phenomenon of electrowetting, i.e., the dependence of the macroscopic contact angle of a fluid on the electrostatic potential of the substrate, is analyzed in terms of the density functional theory of wetting. It is shown that electrowetting is not an electrocapillarity effect, i.e., it cannot be consistently understood in terms of the variation of the substrate-fluid interfacial tension with the electrostatic substrate potential, but it is related to the depth of the effective interface potential. The key feature, which has been overlooked so far and which occurs naturally in the density functional approach is the structural change of a fluid if it is brought into contact with another fluid. These structural changes occur in the present context as the formation of finite films of one fluid phase in between the substrate and the bulk of the other fluid phase. The non-vanishing Donnan potentials (Galvani potential differences) across such film-bulk fluid interfaces, which generically occur due to an unequal partitioning of ions as a result of differences of solubility contrasts, lead to correction terms in the electrowetting equation, which become relevant for sufficiently small substrate potentials. Whereas the present density functional approach confirms the commonly used electrocapillarity-based electrowetting equation as a good approximation for the cases of metallic electrodes or electrodes coated with a hydrophobic dielectric in contact with an electrolyte solution and an ion-free oil, a significantly reduced tendency for electrowetting is predicted for electrodes coated with a dielectric which is hydrophilic or which is in contact with two immiscible electrolyte solutions., Comment: Submitted

The possibility to extract properties of an interface between two immiscible liquids, e.g., electrolyte solutions or polyelectrolyte multilayers, by means of impedance spectroscopy is investigated theoretically within a dynamic density functional theory which is equivalent to the Nernst-Planck-Poisson theory. A novel approach based on a two-step fitting procedure of an equivalent circuit to impedance spectra is proposed which allows to uniquely separate bulk and interfacial elements. Moreover, the proposed method avoids overfitting of the bulk properties of the two liquids in contact and underfitting of the interfacial properties, as they might occur for standard one-step procedures. The key idea is to determine the bulk elements of the equivalent circuit in a first step by fitting corresponding sub-circuits to the spectra of uniform electrolyte solutions, and afterwards fitting the full equivalent circuit with fixed bulk elements to the impedance spectrum containing the interface. This approach is exemplified for an equivalent circuit which leads to a physically intuitive qualitative behavior as well as to quantitively realistic values of the interfacial elements. The proposed method is robust such that it can be expected to be applicable to a wide class of systems with liquid-liquid interfaces.

The time-dependent structure, interfacial tension, and evaporation of an oversaturated colloid-rich (liquid) phase in contact with an undersaturated colloid-poor (vapor) phase of a colloidal dispersion is investigated theoretically during the early-stage relaxation, where the interface is relaxing towards a local equilibrium state while the bulk phases are still out of equilibrium. Since systems of this type exhibit a clear separation of colloidal and solvent relaxation time scales with typical times of interfacial tension measurements in between, they can be expected to be suitable for analogous experimental studies, too. The major finding is that, irrespective of how much the bulk phases differ from two-phase coexistence, the interfacial structure and the interfacial tension approach those at two-phase coexistence during the early-stage relaxation process. This is a surprising observation since it implies that the relaxation towards global equilibrium of the interface is not following but preceding that of the bulk phases. Scaling forms for the local chemical potential, the flux, and the dissipation rate exhibit qualitatively different leading order contributions depending on whether an equilibrium or a non-equilibrium system is considered. The degree of non-equilibrium between the bulk phases is found to not influence the qualitative relaxation behavior (i.e., the values of power-law exponents), but to determine the quantitative deviation of the observed quantities from their values at two-phase coexistence. Whereas the underlying dynamics differs between colloidal and molecular fluids, the behavior of quantities such as the interfacial tension approaching the equilibrium values during the early-stage relaxation process, during which non-equilibrium conditions of the bulk phases are not changed, can be expected to occur for both types of systems., Comment: Submitted

Wetting of a charged substrate by an electrolyte solution is investigated by means of classical density functional theory applied to a lattice model. Within the present model the pure, i.e., salt-free solvent, for which all interactions are of the nearest-neighbor type only, exhibits a second-order wetting transition for all strengths of the substrate-particle and the particle-particle interactions for which the wetting transition temperature is nonzero. The influences of the substrate charge density and of the ionic strength on the wetting transition temperature and on the order of the wetting transition are studied. If the substrate is neutral, the addition of salt to the solvent changes neither the order nor the transition temperature of the wetting transition of the system. If the surface charge is nonzero, upon adding salt this continuous wetting transition changes to first-order within the wide range of substrate surface charge densities and ionic strengths studied here. As the substrate surface charge density is increased, at fixed ionic strength, the wetting transition temperature decreases and the prewetting line associated with the first-order wetting transition becomes longer. This decrease of the wetting transition temperature upon increasing the surface charge density becomes more pronounced by decreasing the ionic strength., Comment: submitted

The influence of ions on the bulk phase behavior of binary liquid mixtures acting as their solvents and on the corresponding interfacial structures close to a planar wall is investigated by means of density functional theory based on local descriptions of the effective interactions between ions and their solvents. The bilinear coupling approximation (BCA), which has been used in numerous previous related investigations, is compared with a novel local density approximation (LDA) for the ion-solvent interactions. It turns out that within BCA the bulk phase diagrams, the two-point correlation functions, and critical adsorption exhibit qualitative features which are not compatible with the available experimental data. These discrepancies do not occur within the proposed LDA. Further experimental investigations are suggested which assess the reliability of the proposed LDA. This approach allows one to obtain a consistent and rather general understanding of the effects of ions on solvent properties. From our analysis we infer in particular that there can be an experimentally detectable influence of ions on binary liquid mixtures due to steric effects but not due to charge effects.

The subtle interplay between critical phenomena and electrostatics is investigated by considering the effective force acting on two parallel walls confining a near-critical binary liquid mixture with added salt. The ion-solvent coupling can turn a non-critical repulsive electrostatic force into an attractive one upon approaching the critical point. However, the effective force is eventually dominated by the critical Casimir effect, the universal properties of which are not altered by the presence of salt. This observation allows a consistent interpretation of recent experimental data., Comment: Submitted

We argue that the extremely low vapor pressures of room temperature ionic liquids near their triple points are due to the combination of strong ionic characters and of low melting temperatures., Comment: Initially submitted manuscript of article M. Bier and S. Dietrich, Mol. Phys. 108, 211 (2010) [Corrigendum: Mol. Phys. 108, 1413 (2010)]

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., Comment: Submitted

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., Comment: 9 pages, 4 figures, accepted by JPCM for proceedings of LMC7

The interfacial tension of a liquid droplet surrounded by another liquid in the presence of microscopic ions is studied as a function of the droplet radius. An analytical expression for the interfacial tension is obtained within a linear Poisson-Boltzmann theory and compared with numerical results from non-linear Poisson-Boltzmann theory. The excess liquid-liquid interfacial tension with respect to the pure, salt-free liquid-liquid interfacial tension is found to decompose into a curvature-independent part due to short-ranged interfacial effects and a curvature-dependent electrostatic contribution. Several curvature-dependent regimes of different scaling of the electrostatic excess interfacial tension are identified. Symmetry relations of the interfacial tension upon swapping droplet and bulk liquid are found to hold in the low-curvature limit, which, e.g., lead to a sign change of the excess Tolman length. For some systems a low-curvature expansion up to second order turns out to be applicable if and only if the droplet size exceeds the Debye screening length in the droplet, independent of the Debye length in the bulk., Comment: submitted

We study the spontaneous charging and the crystallization of spherical micron-sized water-droplets dispersed in oil by numerically solving, within a Poisson-Boltzmann theory in the geometry of a spherical cell, for the density profiles of the cations and anions in the system. We take into account screening, ionic Born self-energy differences between oil and water, and partitioning of ions over the two media. We find that the surface charge density of the droplet as induced by the ion partitioning is significantly affected by the droplet curvature and by the finite density of the droplets. We also find that the salt concentration and the dielectric constant regime in which crystallization of the water droplets is predicted is enhanced substantially compared to results based on the planar oil-water interface, thereby improving quantitative agreement with recent experiments., Comment: 10 pages, 7 figures, submitted for publication

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., Comment: submitted

It is theoretically shown that the excess liquid-liquid interfacial tension between two electrolyte solutions as a function of the ionic strength I behaves asymptotically as O(- I^0.5) for small I and as O(+- I) for large I. The former regime is dominated by the electrostatic potential due to an unequal partitioning of ions between the two liquids whereas the latter regime is related to a finite interfacial thickness. The crossover between the two asymptotic regimes depends sensitively on material parameters suggesting that, depending on the actual system under investigation, the experimentally accessible range of ionic strengths can correspond to either the small or the large ionic strength regime. In the limiting case of a liquid-gas surface where ion partitioning is absent, the image chage interaction can dominate the surface tension for small ionic strength I such that an Onsager-Samaras limiting law O(- I ln(I)) is expected. The proposed picture is consistent with more elaborate models and published measurements., Comment: Accepted for publication in Physical Review Letters

Nonequilibrium steady states in an open system connecting two reservoirs of platelike colloidal particles are investigated by means of a recently proposed phenomenological dynamic density functional theory [M. Bier and R. van Roij, Phys. Rev. E 76, 021405 (2007)]. The platelike colloidal particles are approximated within the Zwanzig model of restricted orientations, which exhibits an isotropic-nematic bulk phase transition. Inhomogeneities of the local chemical potential generate a diffusion current which relaxes to a nonvanishing value if the two reservoirs coupled to the system sustain different chemical potentials. The relaxation process of initial states towards the steady state turns out to comprise two regimes: a smoothening of initial steplike structures followed by an ultimate relaxation of the slowest diffusive mode. The position of a nonequilibrium interface and the particle current of steady states depend nontrivially on the structure of the reservoirs due to the coupling between translational and orientational degrees of freedom of the fluid.

The relaxation dynamics of a model fluid of platelike colloidal particles is investigated by means of a phenomenological dynamic density functional theory. The model fluid approximates the particles within the Zwanzig model of restricted orientations. The driving force for time-dependence is expressed completely by gradients of the local chemical potential which in turn is derived from a density functional -- hydrodynamic interactions are not taken into account. These approximations are expected to lead to qualitatively reliable results for low densities as those within the isotropic-nematic two-phase region. The formalism is applied to model an initially spatially homogeneous stable or metastable isotropic fluid which is perturbed by switching a two-dimensional array of Gaussian laser beams. Switching on the laser beams leads to an accumulation of colloidal particles in the beam centers. If the initial chemical potential and the laser power are large enough a preferred orientation of particles occurs breaking the symmetry of the laser potential. After switching off the laser beams again the system can follow different relaxation paths: It either relaxes back to the homogeneous isotropic state or it forms an approximately elliptical high-density core which is elongated perpendicular to the dominating orientation in order to minimize the surface free energy. For large supersaturations of the initial isotropic fluid the high-density cores of neighboring laser beams of the two-dimensional array merge into complex superstructures., Comment: low-resolution figures due to file size restrictions, revised version

Surface properties of mixtures of charged platelike colloids and salt in contact with a charged planar wall are studied within density functional theory. The particles are modeled by hard cuboids with their edges constrained to be parallel to the Cartesian axes corresponding to the Zwanzig model and the charges of the particles are concentrated in their centers. The density functional applied is an extension of a recently introduced functional for charged platelike colloids. Analytically and numerically calculated bulk and surface phase diagrams exhibit first-order wetting for sufficiently small macroion charges and isotropic bulk order as well as first-order drying for sufficiently large macroion charges and nematic bulk order. The asymptotic wetting and drying behavior is investigated by means of effective interface potentials which turn out to be asymptotically the same as for a suitable neutral system governed by isotropic nonretarded dispersion forces. Wetting and drying points as well as predrying lines and the corresponding critical points have been located numerically. A crossover from monotonic to non-monotonic electrostatic potential profiles upon varying the surface charge density has been observed. Due to the presence of both the Coulomb interactions and the hard-core repulsions, the surface potential and the surface charge do not vanish simultaneously, i.e., the point of zero charge and the isoelectric point of the surface do not coincide., Comment: 14 pages, submitted

Dense suspensions of small strongly interacting particles are complex systems, which are rarely understood on the microscopic level. We investigate properties of dense suspensions and sediments of small spherical Al_2O_3 particles in a shear cell by means of a combined Molecular Dynamics (MD) and Stochastic Rotation Dynamics (SRD) simulation. We study structuring effects and the dependence of the suspension's viscosity on the shear rate and shear thinning for systems of varying salt concentration and pH value. To show the agreement of our results to experimental data, the relation between bulk pH value and surface charge of spherical colloidal particles is modeled by Debye-Hueckel theory in conjunction with a 2pK charge regulation model., Comment: 15 pages, 8 figures

In a recent article (cond-mat/0510119) it has been argued that the Jarzynski equation is violated for adiabatic stretching processes of a three dimensional rotor system. Here we want to show that the reasoning is not correct. Rather, the Jarzynski equation is fulfilled for this adiabatically stretched rotor., Comment: 2 pages, no figures

It is shown that ''Theorem 1'' of the article ''Copenhagen Interpretation of Quantum Mechanics Is Incorrect'' by G.-L. Li and V.O.K. Li (see quant-ph/0509089) is false. Therefore the assertion expressed in the title of that article is untenable., Comment: 1/2 page, no figures

Bulk properties and free interfaces of mixtures of charged platelike colloids and salt are studied within density-functional theory. The particles are modeled by hard cuboids with their edges constrained to be parallel to the artesian axes corresponding to the Zwanzig model. The charges of the particles are concentrated in their center. The density functional is derived by functional integration of an extension of the Debye-H\"uckel pair distribution function with respect to the interaction potential. For sufficiently small macroion charges, the bulk phase diagrams exhibit one isotropic and one nematic phase separated by a first-order phase transition. With increasing platelet charge, the isotropic and nematic binodals are shifted to higher densities. The Donnan potential between the coexisting isotropic and nematic phases is inferred from bulk structure calculations. Non-monotonic density and nematic order parameter profiles are found at a free interface interpolating between the coexisting isotropic and nematic bulk phases. Moreover, electrically charged layers form at the free interface leading to monotonically varying electrostatic potential profiles. Both the widths of the free interfaces and the bulk correlation lengths are approximately given by the Debye length. For fixed salt density, the interfacial tension decreases upon increasing the macroion charge., Comment: 11 pages, submitted to J. Chem. Phys