Interaction regimes for oppositely charged plates with multivalent counterions

Within a mean-field treatment of the interaction between two oppositely charged plates in a salt-free solution, the distance at which a transition from an attractive to a repulsive regime appears can be computed analytically. The mean-field description, however, breaks down under strong Coulombic couplings, which can be achieved at room temperature with multivalent counterions and highly charged surfaces. Making use of the contact theorem and simple physical arguments, we propose explicit expressions for the equation of state in several situations at short distances. The possibility of Bjerrum pair formation is addressed and is shown to have profound consequences on the interactions. To complete the picture, we consider the large-distance limit, from which schematic phase diagram discriminating attractive from repulsive regions can be proposed

Charge ordering in quarter-filled ladder systems coupled to the lattice

We investigate charge ordering in the presence of electron-phonon coupling for quarter-filled ladder systems by using Exact Diagonalization. As an example we consider NaV2O5 using model parameters obtained from first-principles band-structure calculations. The relevant Holstein coupling to the lattice considerably reduces the critical value of the nearest-neighbor Coulomb repulsion at which formation of the zig-zag charge-ordered state occurs, which is then accompanied by a static lattice distortion. Energy and length of a kink-like excitation on the background of the distorted lattice are calculated. Spin and charge spectra on ladders with and without static distortion are obtained, and the charge gap and the effective spin-spin exchange parameter J are extracted. J agrees well with experimental results. Analysis of the dynamical Holstein model, restricted to a small number of phonons, shows that low frequency lattice vibrations increase the charge order, accompanied by dynamically produced zig-zag lattice distortions.Comment: 11 pages, 17 figures, revised version as to appear in Phys. Rev.

A numerical study of a binary Yukawa model in regimes characteristic of globular proteins in solutions

The main goal of this paper is to assess the limits of validity, in the regime of low concentration and strong Coulomb coupling (high molecular charges), for a simple perturbative approximation to the radial distribution functions (RDF), based upon a low-density expansion of the potential of mean force and proposed to describe protein-protein interactions in a recent Small-Angle-Scattering (SAS) experimental study. A highly simplified Yukawa (screened Coulomb) model of monomers and dimers of a charged globular protein ($\beta$-lactoglobulin) in solution is considered. We test the accuracy of the RDF approximation, as a necessary complementary part of the previous experimental investigation, by comparison with the fluid structure predicted by approximate integral equations and exact Monte Carlo (MC) simulations. In the MC calculations, an Ewald construction for Yukawa potentials has been used to take into account the long-range part of the interactions in the weakly screened cases. Our results confirm that the perturbative first-order approximation is valid for this system even at strong Coulomb coupling, provided that the screening is not too weak (i.e., for Debye length smaller than monomer radius). A comparison of the MC results with integral equation calculations shows that both the hypernetted-chain (HNC) and the Percus-Yevick (PY) closures have a satisfactory behavior under these regimes, with the HNC being superior throughout. The relevance of our findings for interpreting SAS results is also discussed.Comment: Physical Review E, in press (2005

Spin models for orientational ordering of colloidal molecular crystals

Two-dimensional colloidal suspensions exposed to periodic external fields exhibit a variety of molecular crystalline phases. There two or more colloids assemble at lattice sites of potential minima to build new structural entities, referred to as molecules. Using the strength of the potential and the filling fraction as control parameter, phase transition to unconventional orientationally ordered states can be induced. We introduce an approach that focuses at the discrete set of orientational states relevant for the phase ordering. The orientationally ordered states are mapped to classical spin systems. We construct effective hamiltonians for dimeric and trimeric molecules on triangular lattices suitable for a statistical mechanics discussion. A mean-field analysis produces a rich phase behavior which is substantiated by Monte Carlo simulations.Comment: 19 pages, 21 figures; misplacement of Fig.3 fixe

Shear Viscosity of Clay-like Colloids in Computer Simulations and Experiments

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 figure

Effect of Polydispersity and Anisotropy in Colloidal and Protein Solutions: an Integral Equation Approach

Application of integral equation theory to complex fluids is reviewed, with particular emphasis to the effects of polydispersity and anisotropy on their structural and thermodynamic properties. Both analytical and numerical solutions of integral equations are discussed within the context of a set of minimal potential models that have been widely used in the literature. While other popular theoretical tools, such as numerical simulations and density functional theory, are superior for quantitative and accurate predictions, we argue that integral equation theory still provides, as in simple fluids, an invaluable technique that is able to capture the main essential features of a complex system, at a much lower computational cost. In addition, it can provide a detailed description of the angular dependence in arbitrary frame, unlike numerical simulations where this information is frequently hampered by insufficient statistics. Applications to colloidal mixtures, globular proteins and patchy colloids are discussed, within a unified framework.Comment: 17 pages, 7 figures, to appear in Interdiscip. Sci. Comput. Life Sci. (2011), special issue dedicated to Prof. Lesser Blu

Attraction between DNA molecules mediated by multivalent ions

The effective force between two parallel DNA molecules is calculated as a function of their mutual separation for different valencies of counter- and salt ions and different salt concentrations. Computer simulations of the primitive model are used and the shape of the DNA molecules is accurately modelled using different geometrical shapes. We find that multivalent ions induce a significant attraction between the DNA molecules whose strength can be tuned by the averaged valency of the ions. The physical origin of the attraction is traced back either to electrostatics or to entropic contributions. For multivalent counter- and monovalent salt ions, we find a salt-induced stabilization effect: the force is first attractive but gets repulsive for increasing salt concentration. Furthermore, we show that the multivalent-ion-induced attraction does not necessarily correlate with DNA overcharging.Comment: 51 pages and 13 figure

Effective interaction between helical bio-molecules

The effective interaction between two parallel strands of helical bio-molecules, such as deoxyribose nucleic acids (DNA), is calculated using computer simulations of the "primitive" model of electrolytes. In particular we study a simple model for B-DNA incorporating explicitly its charge pattern as a double-helix structure. The effective force and the effective torque exerted onto the molecules depend on the central distance and on the relative orientation. The contributions of nonlinear screening by monovalent counterions to these forces and torques are analyzed and calculated for different salt concentrations. As a result, we find that the sign of the force depends sensitively on the relative orientation. For intermolecular distances smaller than $6\AA$ it can be both attractive and repulsive. Furthermore we report a nonmonotonic behaviour of the effective force for increasing salt concentration. Both features cannot be described within linear screening theories. For large distances, on the other hand, the results agree with linear screening theories provided the charge of the bio-molecules is suitably renormalized.Comment: 18 pages, 18 figures included in text, 100 bibliog