Researcher profile

Markus Bier

Markus Bier contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
28works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

28 published item(s)

preprint2023arXiv

Non-equilibrium steady states of electrolyte interfaces

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.

preprint2020arXiv

Heterogeneous surface charge confining an electrolyte solution

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.

preprint2020arXiv

The role of counterions in ionic liquid crystals

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.

preprint2016arXiv

Poisson-Boltzmann study of the effective electrostatic interaction between colloids at an electrolyte interface

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.

preprint2015arXiv

Implications of interface conventions for morphometric thermodynamics

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.

preprint2015arXiv

Non-equilibrium interfacial tension during relaxation

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.

preprint2015arXiv

Static dielectric properties of dense ionic fluids

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.

preprint2015arXiv

Thermal and structural properties of ionic fluids

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.

preprint2015arXiv

Three-phase contact line and line tension of electrolyte solutions in contact with charged substrates

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.

preprint2014arXiv

Density functional theory of electrowetting

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.

preprint2014arXiv

Electrostatic interaction between colloidal particles trapped at an electrolyte interface

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.

preprint2014arXiv

Order of wetting transitions in electrolyte solutions

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.

preprint2013arXiv

Impedance spectroscopy of ions at liquid-liquid interfaces

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.

preprint2013arXiv

Non-equilibrium interfaces in colloidal fluids

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.

preprint2013arXiv

Wetting in electrolyte solutions

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.

preprint2012arXiv

Local theory for ions in binary liquid mixtures

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.

preprint2011arXiv

Electrostatic interactions in critical solvents

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.

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.

preprint2008arXiv

Curvature dependence of the electrolytic liquid-liquid interfacial tension

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.

preprint2008arXiv

Liquid-liquid interfacial tension of electrolyte solutions

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.

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.

preprint2008arXiv

Spontaneous Charging and Crystallization of Water Droplets in Oil

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.

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

Nonequilibrium steady states in fluids of platelike colloidal particles

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.

preprint2007arXiv

Relaxation dynamics in fluids of platelike colloidal particles

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.

preprint2006arXiv

Surface properties of fluids of charged platelike colloids

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.

preprint2005arXiv

Free Isotropic-Nematic Interfaces in Fluids of Charged Platelike Colloids

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ückel 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.