DNA-condensation, redissolution and mesocrystals induced by tetravalent counterions

The distance-resolved effective interaction potential between two parallel DNA molecules is calculated by computer simulations with explicit tetravalent counterions and monovalent salt. Adding counterions first yields an attractive minimum in the potential at short distances which then disappears in favor of a shallower minimum at larger separations. The resulting phase diagram includes a DNA-condensation and redissolution transition and a stable mesocrystal with an intermediate lattice constant for high counterion concentration.Comment: 4 pages, 4 figure

DNA condensation and redissolution: Interaction between overcharged DNA molecules

The effective DNA-DNA interaction force is calculated by computer simulations with explicit tetravalent counterions and monovalent salt. For overcharged DNA molecules, the interaction force shows a double-minimum structure. The positions and depths of these minima are regulated by the counterion density in the bulk. Using two-dimensional lattice sum and free energy perturbation theories, the coexisting phases for DNA bundles are calculated. A DNA-condensation and redissolution transition and a stable mesocrystal with an intermediate lattice constant for high counterion concentration are obtained.Comment: 26 pages, 10 figure

Phase diagram of aggregation of oppositely charged colloids in salty water

Aggregation of two oppositely charged colloids in salty water is studied. We focus on the role of Coulomb interaction in strongly asymmetric systems in which the charge and size of one colloid is much larger than the other one. In the solution, each large colloid (macroion) attracts certain number of oppositely charged small colloids ($Z$-ion) to form a complex. If the concentration ratio of the two colloids is such that complexes are not strongly charged, they condense in a macroscopic aggregate. As a result, the phase diagram in a plane of concentrations of two colloids consists of an aggregation domain sandwiched between two domains of stable solutions of complexes. The aggregation domain has a central part of total aggregation and two wings corresponding to partial aggregation. A quantitative theory of the phase diagram in the presence of monovalent salt is developed. It is shown that as the Debye-H\"{u}ckel screening radius $r_s$ decreases, the aggregation domain grows, but the relative size of the partial aggregation domains becomes much smaller. As an important application of the theory, we consider solutions of long double-helix DNA with strongly charged positive spheres (artificial chromatin). We also consider implications of our theory for in vitro experiments with the natural chromatin. Finally, the effect of different shapes of macroions on the phase diagram is discussed.Comment: 10 pages, 9 figures. The text is rewritten, but results are not change

Temperature Effects on Threshold Counterion Concentration to Induce Aggregation of fd Virus

We seek to determine the mechanism of like-charge attraction by measuring the temperature dependence of critical divalent counterion concentration ($\rm{C_{c}}$) for the aggregation of fd viruses. We find that an increase in temperature causes $\rm{C_c}$ to decrease, primarily due to a decrease in the dielectric constant ($\epsilon$) of the solvent. At a constant $\epsilon$, $\rm{C_c}$ is found to increase as the temperature increases. The effects of $T$ and $\epsilon$ on $\rm {C_{c}}$ can be combined to that of one parameter: Bjerrum length ($l_{B}$). $\rm{C_{c}}$ decreases exponentially as $l_{B}$ increases, suggesting that entropic effect of counterions plays an important role at the onset of bundle formation.Comment: 12 pages, 3 figure

Recursive graphs with small-world scale-free properties

We discuss a category of graphs, recursive clique trees, which have small-world and scale-free properties and allow a fine tuning of the clustering and the power-law exponent of their discrete degree distribution. We determine relevant characteristics of those graphs: the diameter, degree distribution, and clustering parameter. The graphs have also an interesting recursive property, and generalize recent constructions with fixed degree distributions.Comment: 4 pages, 2 figure

Symmetries of Electrostatic Interaction between DNA Molecules

We study a model for pair interaction $U$ of DNA molecules generated by the discrete dipole moments of base-pairs and the charges of phosphate groups, and find noncommutative group of eighth order ${\cal S}$ of symmetries that leave $U$ invariant. We classify the minima using group ${\cal S}$ and employ numerical methods for finding them. The minima may correspond to several cholesteric phases, as well as phases formed by cross-like conformations of molecules at an angle close to $\rm{90}^{o}$, "snowflake phase". The results depend on the effective charge $Q$ of the phosphate group which can be modified by the polycations or the ions of metals. The snowflake phase could exist for $Q$ above the threshold $Q_C$. Below $Q_C$ there could be several cholesteric phases. Close to $Q_C$ the snowflake phase could change into the cholesteric one at constant distance between adjacent molecules.Comment: 13 pages, 4 figure

Molecular dynamics simulation of multivalent ion mediated DNA attraction

All atom molecular dynamics simulations with explicit water were done to study the interaction between two parallel double-stranded DNA molecules in the presence of the multivalent counterions putrescine (2+), spermidine (3+), spermine (4+) and cobalt hexamine (3+). The inter-DNA interaction potential is obtained with the umbrella sampling technique. The attractive force is rationalized in terms of the formation of ion bridges, i.e. multivalent ions which are simultaneously bound to the two opposing DNA molecules. The lifetime of the ion bridges is short on the order of a few nanoseconds.Comment: 4 pages, 5 figures, to be published in Physical Review Letter

Attractive Interactions Between Rod-like Polyelectrolytes: Polarization, Crystallization, and Packing

We study the attractive interactions between rod-like charged polymers in solution that appear in the presence of multi-valence counterions. The counterions condensed to the rods exhibit both a strong transversal polarization and a longitudinal crystalline arrangement. At short distances between the rods, the fraction of condensed counterions increases, and the majority of these occupy the region between the rods, where they minimize their repulsive interactions by arranging themselves into packing structures. The attractive interaction is strongest for multivalent counterions. Our model takes into account the hard-core volume of the condensed counterions and their angular distribution around the rods. The hard core constraint strongly suppresses longitudinal charge fluctuations.Comment: 4 figures, uses revtex, psfig and epsf. The new version contains a different introduction, and the bibliography has been expande

Stabilization and Controlled Association of Inorganic Nanoparticles using Block Copolymers

We report on the structural properties of mixed aggregates made from rare-earth inorganic nanoparticles (radius 20 Angstroms) and polyelectrolyte-neutral block copolymers in aqueous solutions. Using scattering experiments and Monte Carlo simulations, we show that these mixed aggregates have a hierarchical core-shell microstructure. The core is made of densely packed nanoparticles and it is surrounded by a corona of neutral chains. This microstructure results from a process of controlled association and confers to the hybrid aggregates a remarkable colloidal stability.Comment: 14 pages, 5 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