Where the linearized Poisson–Boltzmann cell model fails: Spurious phase separation in charged colloidal suspensions

Effect of prenucleation clusters arising from liquid-liquid phase transition on nucleation in a one-component charged colloidal suspension

It is generally thought that liquid-liquid phase transition (LLPT) in a one-component suspension never or only very rarely happens. If this were true, it would contradict the two nonclassical nucleation models building on either liquid droplets or prenucleation clusters (PNCs). One way out of this paradox is to suppose that LLPT occurs in pathway to nucleation. This study specifies the physical parameters of charged colloids which can bring out LLPT according to the consistent prediction of the DLVO (Derjaguin-Landau-Verwey-Overbeek) potential and the Sogami potential about long-range attraction, and reveals that surface charge is not the only factor to affect attraction, size also plays an essential role. For the first time, we follow exactly the evolution from LLPT to nucleation in which PNCs participate, and characterize pre-ordered liquid-like property of the PNCs and their particle-like and template effect by optical microscopy and light scattering. Furthermore, it is found that when the configuration of the PNCs is changed by a little salt, the pathway to nucleation is altered significantly. Our results demystify LLPT in a one-component suspension and dissolve the paradox, thus extending the range of applicability of the nonclassical nucleation models.

Grand-canonical approach to density functional theory of electrocatalytic systems: Thermodynamics of solid-liquid interfaces at constant ion and electrode potentials

Properties of solid-liquid interfaces are of immense importance for electrocatalytic and electrochemical systems, but modeling such interfaces at the atomic level presents a serious challenge and approaches beyond standard methodologies are needed. An atomistic computational scheme needs to treat at least part of the system quantum mechanically to describe adsorption and reactions, while the entire system is in thermal equilibrium. The experimentally relevant macroscopic control variables are temperature, electrode potential, and the choice of the solvent and ions, and these need to be explicitly included in the computational model as well; this calls for a thermodynamic ensemble with fixed ion and electrode potentials. In this work, a general framework within density functional theory (DFT) with fixed electron and ion chemical potentials in the grand canonical (GC) ensemble is established for modeling electrocatalytic and electrochemical interfaces. Starting from a fully quantum mechanical description of multi-component GC-DFT for nuclei and electrons, a systematic coarse-graining is employed to establish various computational schemes including (i) the combination of classical and electronic DFTs within the GC ensemble and (ii) on the simplest level a chemically and physically sound way to obtain various (modified) Poisson-Boltzmann (mPB) implicit solvent models. The detailed and rigorous derivation clearly establishes which approximations are needed for coarse-graining as well as highlights which details and interactions are omitted in vein of computational feasibility. The transparent approximations also allow removing some of the constraints and coarse-graining if needed. We implement various mPB models within a linear dielectric continuum in the GPAW code and test their capabilities to model capacitance of electrochemical interfaces as well as study different approaches for modeling partly periodic charged systems. Our rigorous and well-defined DFT coarse-graining scheme to continuum electrolytes highlights the inadequacy of current linear dielectric models for treating properties of the electrochemical interface.

The metastable states of foam films containing electrically charged micelles or particles: experiment and quantitative interpretation

The stepwise thinning (stratification) of liquid films containing electrically charged colloidal particles (in our case - surfactant micelles) is investigated. Most of the results are applicable also to films from nanoparticle suspensions. The aim is to achieve agreement between theory and experiment, and to better understand the physical reasons for this phenomenon. To test different theoretical approaches, we obtained experimental data for free foam films from micellar solutions of three ionic surfactants. The theoretical problem is reduced to the interpretation of the experimental concentration dependencies of the step height and of the final film thickness. The surface charges of films and micelles are calculated by means of the charge-regulation model, with a counterion-binding (Stern) constant determined from the fit of surface tension isotherms. The applicability of three models was tested: the Poisson-Boltzmann (PB) model; the jellium-approximation (JA), and the cell model (CM). The best agreement theory/experiment was obtained with the JA model without using any adjustable parameters. Two theoretical approaches are considered. First, in the energy approach the step height is identified with the effective diameter of the charged micelles, which represents an integral of the electrostatic-repulsion energy calculated by the JA model. Second, in the osmotic approach the step height is equal to the inverse cubic root of micelle number density in the bulk of solution. Both approaches are in good agreement with the experiment if the suspension of charged particles (micelles) represents a jellium, i.e. if the particle concentration is uniform despite the field of the electric double layers. The results lead to a convenient method for determining the aggregation number of ionic surfactant micelles from the experimental heights of the steps.

Poisson-Boltzmann theory of charged colloids: limits of the cell model for salty suspensions

Thermodynamic properties of charge-stabilized colloidal suspensions and polyelectrolyte solutions are commonly modelled by implementing the mean-field Poisson-Boltzmann (PB) theory within a cell model. This approach models a bulk system by a single macroion, together with counterions and salt ions, confined to a symmetrically shaped, electroneutral cell. While easing numerical solution of the nonlinear PB equation, the cell model neglects microion-induced interactions and correlations between macroions, precluding modelling of macroion ordering phenomena. An alternative approach, which avoids the artificial constraints of cell geometry, exploits the mapping of a macroion-microion mixture onto a one-component model of pseudo-macroions governed by effective interparticle interactions. In practice, effective-interaction models are usually based on linear-screening approximations, which can accurately describe strong nonlinear screening only by incorporating an effective (renormalized) macroion charge. Combining charge renormalization and linearized PB theories, in both the cell model and an effective-interaction (cell-free) model, we compute osmotic pressures of highly charged colloids and monovalent microions, in Donnan equilibrium with a salt reservoir, over a range of concentrations. By comparing predictions with primitive model simulation data for salt-free suspensions, and with predictions from nonlinear PB theory for salty suspensions, we chart the limits of both the cell model and linear-screening approximations in modelling bulk thermodynamic properties. Up to moderately strong electrostatic couplings, the cell model proves accurate for predicting osmotic pressures of deionized (counterion-dominated) suspensions. With increasing salt concentration, however, the relative contribution of macroion interactions to the osmotic pressure grows, leading predictions from the cell and effective-interaction models to deviate. No evidence is found for a liquid-vapour phase instability driven by monovalent microions. These results may guide applications of PB theory to colloidal suspensions and other soft materials.

Colloidal gas-liquid condensation of polystyrene latex particles with intermediate kappa a values (5 to 160, a >> kappa(-1))

Polystyrene latex particles showed gas-liquid condensation under the conditions of large particle radius (a >> kappa(-1)) and intermediate kappa a, where kappa is the Debye-Hückel parameter and a is the particle radius. The particles were dissolved in deionized water containing ethanol from 0 to 77 vol %, settled to the bottom of the glass plate within 1 h, and then laterally moved toward the center of a cell over a 20 h period in reaching a state of equilibrium condensation. All of the suspensions that were 1 and 3 microm in diameter and 0.01-0.20 vol % in concentration realized similar gas-liquid condensation with clear gas-liquid boundaries. In 50 vol % ethanol solvent, additional ethanol was added to enhance the sedimentation force so as to restrict the particles in a monoparticle layer thickness. The coexistence of gas-liquid-solid (crystalline solid) was microscopically recognized from the periphery to the center of the condensates. A phase diagram of the gas-liquid condensation was created as a function of KCl concentration at a particle diameter of 3 microm, 0.10 vol % concentration, and 50:50 water/ethanol solvent at room temperature. The miscibility gap was observed in the concentration range from 1 to 250 microM. There was an upper limit of salt concentration where the phase separation disappeared, showing nearly critical behavior of macroscopic density fluctuation from 250 microM to 1 mM. These results add new experimental evidence to the existence of colloidal gas-liquid condensation and specify conditions of like-charge attraction between particles.

Electroneutrality and phase behavior of colloidal suspensions

Several statistical mechanical theories predict that colloidal suspensions of highly charged macroions and monovalent microions can exhibit unusual thermodynamic phase behavior when strongly deionized. Density-functional, extended Debye-Hückel, and response theories, within mean-field and linearization approximations, predict a spinodal phase instability of charged colloids below a critical salt concentration. Poisson-Boltzmann cell model studies of suspensions in Donnan equilibrium with a salt reservoir demonstrate that effective interactions and osmotic pressures predicted by such theories can be sensitive to the choice of reference system, e.g., whether the microion density profiles are expanded about the average potential of the suspension or about the reservoir potential. By unifying Poisson-Boltzmann and response theories within a common perturbative framework, it is shown here that the choice of reference system is dictated by the constraint of global electroneutrality. On this basis, bulk suspensions are best modeled by density-dependent effective interactions derived from a closed reference system in which the counterions are confined to the same volume as the macroions. Lower-dimensional systems (e.g., monolayers, clusters), depending on the strength of macroion-counterion correlations, may be governed instead by density-independent effective interactions tied to an open reference system with counterions dispersed throughout the reservoir, possibly explaining the observed structural crossover in colloidal monolayers and anomalous metastability of colloidal crystallites.

Phase separation of charge-stabilized colloids: a Gibbs ensemble Monte Carlo simulation study

The fluid phase behavior of charge-stabilized colloidal suspensions is explored by applying a variant of the Gibbs ensemble Monte Carlo simulation method to a coarse-grained one-component model with implicit microions and solvent. The simulations take as input linear-response approximations for the effective electrostatic interactions--a hard-sphere-Yukawa pair potential and a one-body volume energy. The conventional Gibbs ensemble trial moves are supplemented by exchange of (implicit) salt between coexisting phases, with acceptance probabilities influenced by the state dependence of the effective interactions. Compared with large-scale simulations of the primitive model, with explicit microions, our computationally practical simulations of the one-component model closely match the pressures and pair distribution functions at moderate electrostatic couplings. For macroion valences and couplings within the linear-response regime, deionized aqueous suspensions with monovalent microions exhibit separation into macroion-rich and macroion-poor fluid phases below a critical salt concentration. The resulting pressures and phase diagrams are in excellent agreement with predictions of a variational free energy theory based on the same model.

Electro-osmotic equilibria for a semipermeable shell filled with a solution of polyions

The authors study theoretically the electrostatic equilibria for a shell filled with a suspension of polyions (e.g., colloids, polyelectrolytes, etc.) and immersed in an infinite salt-free reservoir. The shell is treated as impermeable for polyions, but allowing free diffusion of counterions. From the solution of the linearized Poisson-Boltzmann equation we obtain the distribution of the potential and concentration profiles for polyions. The authors then derive explicit formulas for the excess electro-osmotic pressure of a polyion solution exerted by the shell. This is shown to be due to a concentration of polyions at the inner shell boundary and can be very different from the pressure of a corresponding bulk polyion solution.

Rayleigh instability of charged aggregates: Role of the dimensionality, ionic strength, and dielectric contrast

We extended a previous analysis of the classical Rayleigh instability of spherical charged droplets in the presence of neutralizing monovalent counterions [M. Deserno, Eur. Phys. J. E 6, 163 (2001)], by generalizing the problem for suspensions of aggregates with D-dimensional symmetry, corresponding for D = 2 to infinite (rodlike) cylindrical charged bundles and for D = 3 to spherical charged droplets. In addition, we include the effects of added monovalent salt and of dielectric contrast between the charged aggregate and the surrounding solvent. The electrostatic energy taking the microion screening into account is estimated using uniform profiles within the framework of the cell model. We verify the robustness of these results by also considering Debye-Hückel-type microion profiles that are obtained by the minimization of a linearized Poisson-Boltzmann free-energy functional. In the case when the microions can enter inside the charged aggregates, we confirm the occurrence of a discontinuous phase change between aggregates of finite size and an infinite aggregate, which takes place at a collapse temperature that depends on their volume fraction phi and on the salt content. Decrease of phi shifts the phase-change temperature toward higher values, while salt addition has an opposite effect. We obtain analytical expressions for the phase-separation line in the asymptotic limit of infinite dilution (phi-->0), showing that the collapse temperature depends logarithmically on phi . As an application for D = 3 we discuss the stability of the pearl-necklace structures of flexible polyelectrolytes in poor solvents. The case D = 2 is applied to the problem of finite bundle sizes of stiff polyelectrolytes that attract each other-via, e.g., multivalent counterions-leading to an effective surface tension.

Shear viscosity of claylike colloids in computer simulations and experiments

Dense suspensions of small strongly interacting particles are complex systems that are rarely understood on the microscopic level. We investigate properties of dense suspensions and sediments of small spherical Al2O3 particles in a shear cell by means of a combined molecular-dynamics and stochastic rotation dynamics 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 with experimental data, the relation between the bulk pH value and surface charge of spherical colloidal particles is modeled by Debye-Hückel theory in conjunction with a 2 pK charge regulation model.

Phase separation in charge-stabilized colloidal suspensions: influence of nonlinear screening

The phase behavior of charge-stabilized colloidal suspensions is modeled by a combination of response theory for electrostatic interparticle interactions and variational theory for free energies. Integrating out degrees of freedom of the microions (counterions, salt ions), the macroion-microion mixture is mapped onto a one-component system governed by effective macroion interactions. Linear response of microions to the electrostatic potential of the macroions results in a screened-Coulomb (Yukawa) effective pair potential and a one-body volume energy, while nonlinear response modifies the effective interactions [A. R. Denton, Phys. Rev. E 70, 031404 (2004)]. The volume energy and effective pair potential are taken as input to a variational free energy, based on thermodynamic perturbation theory. For both linear and first-order nonlinear effective interactions, a coexistence analysis applied to aqueous suspensions of highly charged macroions and monovalent microions yields bulk separation of macroion-rich and macroion-poor phases below a critical salt concentration, in qualitative agreement with predictions of related linearized theories [R. van Roij, M. Dijkstra, and J.-P. Hansen, Phys. Rev. E 59, 2010 (1999); P. B. Warren, J. Chem. Phys. 112, 4683 (2000)]. It is concluded that nonlinear screening can modify phase behavior but does not necessarily suppress bulk phase separation of de-ionized suspensions.

Volume terms for charged colloids: a grand-canonical treatment

We present a study of thermodynamic properties of suspensions of charged colloids on the basis of linear Poisson-Boltzmann theory. We calculate the effective Hamiltonian of the colloids by integrating out the ionic degrees of freedom grand canonically. This procedure not only yields the well-known pairwise screened-Coulomb interaction between the colloids, but also additional volume terms that affect the phase behavior and the thermodynamic properties, such as the osmotic pressure. These calculations are greatly facilitated by the grand-canonical character of our treatment of the ions and allow for relatively fast computations compared to earlier studies in the canonical ensemble. Moreover, the present derivation of the volume terms are relatively simple, make a direct connection with Donnan equilibrium, yield an explicit expression for the effective screening constant, and allow for extensions to include, for instance, nonlinear effects.

Long-ranged electrostatic repulsion and crystallization of emulsion droplets in an ultralow dielectric medium supercritical carbon dioxide

Electrostatic repulsion stabilizes micrometer-sized water droplets with spacings greater than 10 microm in an ultralow dielectric medium, CO2 (epsilon = 1.5), at elevated pressures. The morphology of the water/CO2 emulsion is characterized by optical microscopy and laser diffraction as a function of height. The counterions, stabilized with a nonionic, highly branched, stubby hydrocarbon surfactant, form an extremely thick double layer with a Debye screening length of 8.9 microm. As a result of the balance between electrostatic repulsion and the downward force due to gravity, the droplets formed a hexagonal crystalline lattice at the bottom of the high-pressure cell with spacings of over 10 microm. The osmotic pressure, calculated by solving the Poisson-Boltzmann equation in the framework of the Wigner-Seitz cell model, is in good agreement with that determined from the sedimentation profile measured by laser diffraction. Thus, the long-ranged stabilization of the emulsion may be attributed to electrostatic stabilization. The ability to form new types of colloids in CO2 with electrostatic stabilization is beneficial because steric stabilization is often unsatisfactory because of poor solvation of the stabilizers.

Phase separation, interface properties, and charge density waves in a simplified model for a macroion suspension

A simplified density functional theory for a macroion suspension is examined, where the correlation free energy corresponds to the macroion self-energy, treated within a linearized or Debye-Hückel approximation. The model possesses a miscibility gap (liquid-liquid phase separation) at low ionic strength. Within the gap, density profiles, electrical structure, and surface tension are calculated for the interface between coexisting phases, using a variational approximation. Additionally, structure factors are calculated for the homogeneous system. As one approaches the critical points, the structure factors can diverge at a nonzero wave vector, signaling the onset of charge density wave phases. Although the quantitative results should be treated with care, the results may be indicative of the rich phenomenology that can arise in asymmetric charged systems.

Gas-solid coexistence in highly charged colloidal suspensions

Aqueous suspensions of highly charged polystyrene particles with different volume fractions have been investigated for structural ordering and phase behavior using static light scattering (SLS) and confocal laser scanning microscope (CLSM). Under deionized conditions, suspensions of high-charge-density colloidal particles remained disordered whereas suspensions of relatively low charge density showed crystallization by exhibiting iridescence for the visible light. Though for the unaided eye crystallized suspensions appeared homogeneous, SLS measurements and CLSM observations have revealed their inhomogeneous nature in the form of the coexistence of voids with dense ordered regions. CLSM investigations on disordered suspensions showed their inhomogeneous nature in the form coexistence of voids with dense disordered (amorphous) regions. Our studies on highly charged colloids confirm the occurrence of gas-solid transition and are in accordance with predictions of Monte Carlo simulations using a pair-potential having a long-range attractive term [Mohanty, P. S.; Tata, B. V. R. J. Colloid Interface Sci. 2003, 264, 101]. On the basis of our experimental and simulation results, we argue that the reported reentrant disordered state [Yamanaka et al. Phys. Rev. Lett. 1998, 80, 5806 and Toyotama et al. Langmuir 2003, 19, 3236] in charged colloids observed at high charge densities is a gas-solid coexistence state.

Aqueous suspensions of charged spherical colloids: dependence of the surface charge on ionic strength, acidity, and colloid concentration

We theoretically investigate the dependence of the surface charge developed on charged spherical colloids upon several environmental parameters: the ionic strength of the monovalent added electrolyte, acidity (stabilized by a pH buffer solution), and colloid concentration. In the framework of the mean-field Poisson-Boltzmann spherical cell model, we include the charged colloid-microion correlations into the buffer equation, and we allow for the specific binding of ions to the ionizable groups on the colloid surface. Theoretical predictions are compared to the results obtained under the planar-symmetry Gouy-Chapman approximation and analyzed for the experimental conditions of an aqueous dispersion of the phospholipid dimyristoyl phosphatidylglycerol (DMPG). Experimental measurements of the partition ratio of an aqueous soluble cationic spin label on buffered dispersions of polyanionic unilamellar vesicles of DMPG in the presence of added monovalent salt are theoretically interpreted in terms of ion partition due to electrostatic interactions. We show that the specific binding of the probe must be admitted to explain the experimental results.

Osmotic pressure of aqueous chondroitin sulfate solution: a molecular modeling investigation

The osmotic pressure of chondroitin sulfate (CS) solution in contact with an aqueous 1:1 salt reservoir of fixed ionic strength is studied using a recently developed coarse-grained molecular model. The effects of sulfation type (4- vs. 6-sulfation), sulfation pattern (statistical distribution of sulfate groups along a chain), ionic strength, CS intrinsic stiffness, and steric interactions on CS osmotic pressure are investigated. At physiological ionic strength (0.15 M NaCl), the sulfation type and pattern, as measured by a standard statistical description of copolymerization, are found to have a negligible influence on CS osmotic pressure, which depends principally on the mean volumetric fixed charge density. The intrinsic backbone stiffness characteristic of polysaccharides such as CS, however, is demonstrated to contribute significantly to its osmotic pressure behavior, which is similar to that of a solution of charged rods for the 20-disaccharide chains considered. Steric excluded volume is found to play a negligible role in determining CS osmotic pressure at physiological ionic strength due to the dominance of repulsive intermolecular electrostatic interactions that maintain chains maximally spaced in that regime, whereas at high ionic-strength steric interactions become dominant due to electrostatic screening. Osmotic pressure predictions are compared to experimental data and to well-established theoretical models including the Donnan theory and the Poisson-Boltzmann cylindrical cell model.

Equilibrium bundle size of rodlike polyelectrolytes with counterion-induced attractive interactions

Multivalent counterions can induce an effective attraction between like-charged rodlike polyelectrolytes, leading to the formation of polyelectrolyte bundles. In this paper, we calculate the equilibrium bundle size using a simple model in which the attraction between polyelectrolytes (assumed to be pairwise additive) is treated phenomenologically. If the counterions are pointlike, they almost completely neutralize the charge of the bundle, and the equilibrium bundle size diverges. When the counterions are large, however, steric and short-range electrostatic interactions prevent charge neutralization of the bundle, thus forcing the equilibrium bundle size to be finite. We also show that if the attractive interactions between the rods become frustrated as the bundle grows, finite-size bundles can be obtained with pointlike counterions.

Contribution of the Nernst potential to stiffness constants: the asymmetrical case

Inside biological membranes, one of the fundamental functions of active proteins such as pumps is to generate some electrochemical gradient across the membrane and then, to establish a new stationary state. The membrane electric potential generated by activity modifies the stiffness constants of the membrane. A spontaneous curvature appears if the inner and outer Debye lengths are different. The corresponding characteristic radius falls in the range from 0.08 microm to 50 microm. The bending elastic modulus is always increased. This effect is only noticeable in the limit of large Debye length from 0.5 microm to 0.09 microm. For a Nernst potential of 100 mV and a Debye length of 0.2 microm, the bending modulus can reach 40k(B)T.

Effective charge of colloidal particles

A new dynamical definition of the effective colloidal charge, especially suitable for the Monte Carlo and Molecular-dynamics simulations, is proposed. It is shown that for aqueous colloidal suspensions containing monovalent counterions the "dynamical" effective charge agrees perfectly with the "statistical" effective charge calculated using the Alexander prescription. In the case of multivalent ions, the effective charge behaves in a qualitatively different way from the predictions of the Poisson-Boltzmann theory.

Nonlinear screening and effective electrostatic interactions in charge-stabilized colloidal suspensions

A nonlinear response theory is developed and applied to electrostatic interactions between spherical macroions, screened by surrounding microions, in charge-stabilized colloidal suspensions. The theory describes leading-order nonlinear response of the microions (counterions, salt ions) to the electrostatic potential of the macroions and predicts microion-induced effective many-body interactions between macroions. A linear response approximation [A.R. Denton, Phys. Rev. E 62, 3855 (2000)] yields an effective pair potential of screened-Coulomb (Yukawa) form, as well as a one-body volume energy, which contributes to the free energy. Nonlinear response generates effective many-body interactions and essential corrections to both the effective pair potential and the volume energy. By adopting a random-phase approximation (RPA) for the response functions, and thus neglecting microion correlations, practical expressions are derived for the effective pair and triplet potentials and for the volume energy. Nonlinear screening is found to weaken repulsive pair interactions, induce attractive triplet interactions, and modify the volume energy. Numerical results for monovalent microions are in good agreement with available ab initio simulation data and demonstrate that nonlinear effects grow with increasing macroion charge and concentration and with decreasing salt concentration. In the dilute limit of zero macroion concentration, leading-order nonlinear corrections vanish. Finally, it is shown that nonlinear response theory, when combined with the RPA, is formally equivalent to the mean-field Poisson-Boltzmann theory and that the linear response approximation corresponds, within integral-equation theory, to a linearized hypernetted-chain closure.

Screening of spherical colloids beyond mean field: a local density functional approach

We study the counterion distribution around a spherical macroion and its osmotic pressure in the framework of the recently developed Debye-Hückel-hole-cavity (DHHC) theory. This is a local density functional approach which incorporates correlations into Poisson-Boltzmann theory by adding a free energy correction based on the one-component plasma. We compare the predictions for ion distribution and osmotic pressure obtained by the full theory and by its zero temperature limit with Monte Carlo simulations. They agree excellently for weakly developed correlations and give the correct trend for stronger ones. In all investigated cases the DHHC theory and its computationally simpler zero temperature limit yield better results than the Poisson-Boltzmann theory.

Renormalized jellium model for charge-stabilized colloidal suspensions

We introduce a renormalized jellium model to calculate the equation of state for charged colloidal suspensions. An almost perfect agreement with Monte Carlo simulations is found. Our self-consistent approach naturally allows to define the effective charge of particles at finite colloidal density. Although this quantity may differ significantly from its counterpart obtained from the standard Poisson-Boltzmann cell approach, the osmotic pressures for both models are in good agreement. We argue that by construction, the effective charge obtained using the jellium approximation is more appropriate to the study of colloidal interactions. We also discuss a possibility of a fluid-fluid critical point and show how our equation of state can be used to shed light on the surprising results found in recent sedimentation experiments.

Where the linearized Poisson-Boltzmann cell model fails: the planar case as a prototype study

The linearized Poisson-Boltzmann (PB) approximation is investigated for the classical problem of two infinite, uniformly charged planes in electrochemical equilibrium with an infinite monovalent salt reservoir. At the nonlinear level, we obtain an explicit expression of the associated electrostatic contribution to the semi-grand-canonical potential. The linearized osmotic-pressure difference between the interplane region and the salt reservoir becomes negative in the low-temperature, large-separation, or high-surface charge limits, in disagreement with the exact (at mean-field level) nonlinear PB solution. We show that these artifacts--although thermodynamically consistent with quadratic expansions of the nonlinear functional--can be traced back to the nonfulfillment of the underlying assumptions of the linearization. Explicit comparison between the analytical expressions of the exact nonlinear solution and the corresponding linearized equations allows us to show that the linearized results are asymptotically exact in the weak-coupling and counterionic ideal-gas limits, but always fail otherwise, predicting negative osmotic-pressure differences. By taking appropriate limits of the full nonlinear PB solution, we provide asymptotic expressions for the semi-grand-canonical potential and the osmotic-pressure difference that involve only elementary functions, which cover the complementary region where the linearized theory breaks down.