Where the linearized Poisson-Boltzmann cell model fails: (I) spurious phase separation in charged colloidal suspensions

We perform a linearization of the Poisson-Boltzmann (PB) density functional for spherical Wigner-Seitz cells that yields Debye-H\"uckel-like equations agreeing asymptotically with the PB results in the weak-coupling (high-temperature) limit. Both the canonical (fixed number of microions) as well as the semi-grand-canonical (in contact with an infinite salt reservoir) cases are considered and discussed in a unified linearized framework. In the canonical case, for sufficiently large colloidal charges the linearized theory predicts the occurrence of a thermodynamical instability with an associated phase separation of the homogeneous suspension into dilute (gas) and dense (liquid) phases. In the semi-grand-canonical case it is predicted that the isothermal compressibility and the osmotic-pressure difference between the colloidal suspension and the salt reservoir become negative in the low-temperature, high-surface charge or infinite-dilution (of polyions) limits. As already pointed out in the literature for the latter case, these features are in disagreement with the exact nonlinear PB solution inside a Wigner-Seitz cell and are thus artifacts of the linearization. By using explicitly gauge-invariant forms of the electrostatic potential we show that these artifacts, although thermodynamically consistent with quadratic expansions of the nonlinear functional and osmotic pressure, may be traced back to the non-fulfillment of the underlying assumptions of the linearization.Comment: 32 pages, 3 PostScript figures, submitted to J. Chem. Phy

First-principles dynamical CPA to finite-temperature magnetism of transition metals

We present here the first-principles dynamical CPA (coherent potential approximation) combined with the tight-binding LMTO LDA+U method towards quantitative calculations of the electronic structure and magnetism at finite temperatures in transition metals and compounds. The theory takes into account the single-site dynamical charge and spin fluctuations using the functional integral technique as well as an effective medium. Numerical results for Fe, Co, and Ni show that the theory explains quantitatively the high-temperature properties such as the effective Bohr magneton numbers and the excitation spectra in the paramagnetic state, and describes the Curie temperatures semiquantitatively.Comment: ICM'09 Proceeding

Polyelectrolyte Bundles

Using extensive Molecular Dynamics simulations we study the behavior of polyelectrolytes with hydrophobic side chains, which are known to form cylindrical micelles in aqueous solution. We investigate the stability of such bundles with respect to hydrophobicity, the strength of the electrostatic interaction, and the bundle size. We show that for the parameter range relevant for sulfonated poly-para-phenylenes (PPP) one finds a stable finite bundle size. In a more generic model we also show the influence of the length of the precursor oligomer on the stability of the bundles. We also point out that our model has close similarities to DNA solutions with added condensing agents, hinting to the possibility that the size of DNA aggregates is under certain circumstances thermodynamically limited.Comment: 10 pages, 8 figure

The osmotic pressure of charged colloidal suspensions: A unified approach to linearized Poisson-Boltzmann theory

We study theoretically the osmotic pressure of a suspension of charged objects (e.g., colloids, polyelectrolytes, clay platelets, etc.) dialyzed against an electrolyte solution using the cell model and linear Poisson-Boltzmann (PB) theory. From the volume derivative of the grand potential functional of linear theory we obtain two novel expressions for the osmotic pressure in terms of the potential- or ion-profiles, neither of which coincides with the expression known from nonlinear PB theory, namely, the density of microions at the cell boundary. We show that the range of validity of linearization depends strongly on the linearization point and proof that expansion about the selfconsistently determined average potential is optimal in several respects. For instance, screening inside the suspension is automatically described by the actual ionic strength, resulting in the correct asymptotics at high colloid concentration. Together with the analytical solution of the linear PB equation for cell models of arbitrary dimension and electrolyte composition explicit and very general formulas for the osmotic pressure ensue. A comparison with nonlinear PB theory is provided. Our analysis also shows that whether or not linear theory predicts a phase separation depends crucially on the precise definition of the pressure, showing that an improper choice could predict an artificial phase separation in systems as important as DNA in physiological salt solution.Comment: 16 pages, 5 figures, REVTeX4 styl

On the fluid-fluid phase separation in charged-stabilized colloidal suspensions

We develop a thermodynamic description of particles held at a fixed surface potential. This system is of particular interest in view of the continuing controversy over the possibility of a fluid-fluid phase separation in aqueous colloidal suspensions with monovalent counterions. The condition of fixed surface potential allows in a natural way to account for the colloidal charge renormalization. In a first approach, we assess the importance of the so called ``volume terms'', and find that in the absence of salt, charge renormalization is sufficient to stabilize suspension against a fluid-fluid phase separation. Presence of salt, on the other hand, is found to lead to an instability. A very strong dependence on the approximations used, however, puts the reality of this phase transition in a serious doubt. To further understand the nature of the instability we next study a Jellium-like approximation, which does not lead to a phase separation and produces a relatively accurate analytical equation of state for a deionized suspensions of highly charged colloidal spheres. A critical analysis of various theories of strongly asymmetric electrolytes is presented to asses their reliability as compared to the Monte Carlo simulations

Electrolytic depletion interactions

We consider the interactions between two uncharged planar macroscopic surfaces immersed in an electrolyte solution which are induced by interfacial selectivity. These forces are taken into account by introducing a depletion free-energy density functional, in addition to the usual mean-field Poisson-Boltzmann functional. The minimization of the total free-energy functional yields the density profiles of the microions and the electrostatic potential. The disjoining pressure is obtained by differentiation of the total free energy with respect to the separation of the surfaces, holding the range and strength of the depletion forces constant. We find that the induced interaction between the two surfaces is always repulsive for sufficiently large separations, and becomes attractive at shorter separations. The nature of the induced interactions changes from attractive to repulsive at a distance corresponding to the range of the depletion forces.Comment: 17 pages, 4 Postscript figures, submitted to Physical Review

Effective Interactions and Volume Energies in Charged Colloids: Linear Response Theory

Interparticle interactions in charge-stabilized colloidal suspensions, of arbitrary salt concentration, are described at the level of effective interactions in an equivalent one-component system. Integrating out from the partition function the degrees of freedom of all microions, and assuming linear response to the macroion charges, general expressions are obtained for both an effective electrostatic pair interaction and an associated microion volume energy. For macroions with hard-sphere cores, the effective interaction is of the DLVO screened-Coulomb form, but with a modified screening constant that incorporates excluded volume effects. The volume energy -- a natural consequence of the one-component reduction -- contributes to the total free energy and can significantly influence thermodynamic properties in the limit of low-salt concentration. As illustrations, the osmotic pressure and bulk modulus are computed and compared with recent experimental measurements for deionized suspensions. For macroions of sufficient charge and concentration, it is shown that the counterions can act to soften or destabilize colloidal crystals.Comment: 14 pages, including 3 figure

Hydration interactions: aqueous solvent effects in electric double layers

A model for ionic solutions with an attractive short-range pair interaction between the ions is presented. The short-range interaction is accounted for by adding a quadratic non-local term to the Poisson-Boltzmann free energy. The model is used to study solvent effects in a planar electric double layer. The counter-ion density is found to increase near the charged surface, as compared with the Poisson-Boltzmann theory, and to decrease at larger distances. The ion density profile is studied analytically in the case where the ion distribution near the plate is dominated only by counter-ions. Further away from the plate the density distribution can be described using a Poisson-Boltzmann theory with an effective surface charge that is smaller than the actual one.Comment: 11 Figures in 13 files + LaTex file. 20 pages. Accepted to Phys. Rev. E. Corrected typos and reference

Influence of solvent granularity on the effective interaction between charged colloidal suspensions

We study the effect of solvent granularity on the effective force between two charged colloidal particles by computer simulations of the primitive model of strongly asymmetric electrolytes with an explicitly added hard sphere solvent. Apart from molecular oscillating forces for nearly touching colloids which arise from solvent and counterion layering, the counterions are attracted towards the colloidal surfaces by solvent depletion providing a simple statistical description of hydration. This, in turn, has an important influence on the effective forces for larger distances which are considerably reduced as compared to the prediction based on the primitive model. When these forces are repulsive, the long-distance behaviour can be described by an effective Yukawa pair potential with a solvent-renormalized charge. As a function of colloidal volume fraction and added salt concentration, this solvent-renormalized charge behaves qualitatively similar to that obtained via the Poisson-Boltzmann cell model but there are quantitative differences. For divalent counterions and nano-sized colloids, on the other hand, the hydration may lead to overscreened colloids with mutual attraction while the primitive model yields repulsive forces. All these new effects can be accounted for through a solvent-averaged primitive model (SPM) which is obtained from the full model by integrating out the solvent degrees of freedom. The SPM was used to access larger colloidal particles without simulating the solvent explicitly.Comment: 14 pages, 16 craphic

The Critical Behaviour of the Spin-3/2 Blume-Capel Model in Two Dimensions

The phase diagram of the spin-3/2 Blume-Capel model in two dimensions is explored by conventional finite-size scaling, conformal invariance and Monte Carlo simulations. The model in its $\tau$-continuum Hamiltonian version is also considered and compared with others spin-3/2 quantum chains. Our results indicate that differently from the standard spin-1 Blume-Capel model there is no multicritical point along the order-disorder transition line. This is in qualitative agreement with mean field prediction but in disagreement with previous approximate renormalization group calculations. We also presented new results for the spin-1 Blume-Capel model.Comment: latex 18 pages, 4 figure