Virial expansion for charged colloids and electrolytes

Using a field-theoretic approach, we derive the first few coefficients of the exact low-density (``virial'') expansion of a binary mixture of positively and negatively charged hard spheres (two-component hard-core plasma, TCPHC). Our calculations are nonperturbative with respect to the diameters $d_+$ and $d_-$ and charge valences $q_+$ and $q_-$ of positive and negative ions. Consequently, our closed-form expressions for the coefficients of the free energy and activity can be used to treat dilute salt solutions, where typically $d_+ \sim d_-$ and $q_+ \sim q_-$, as well as colloidal suspensions, where the difference in size and valence between macroions and counterions can be very large. We show how to map the TCPHC on a one-component hard-core plasma (OCPHC) in the colloidal limit of large size and valence ratio, in which case the counterions effectively form a neutralizing background. A sizable discrepancy with the standard OCPHC with uniform, rigid background is detected, which can be traced back to the fact that the counterions cannot penetrate the colloids. For the case of electrolyte solutions, we show how to obtain the cationic and anionic radii as independent parameters from experimental data for the activity coefficient.Comment: 15 page

Small ball probability and Dvoretzky theorem

Large deviation estimates are by now a standard tool inthe Asymptotic Convex Geometry, contrary to small deviationresults. In this note we present a novel application of a smalldeviations inequality to a problem related to the diameters of random sections of high dimensional convex bodies. Our results imply an unexpected distinction between the lower and the upper inclusions in Dvoretzky Theorem

Surface critical behavior of fluids: Lennard-Jones fluid near weakly attractive substrate

The phase behavior of fluids near weakly attractive substrates is studied by computer simulations of the coexistence curve of a Lennard-Jones (LJ) fluid confined in a slitlike pore. The temperature dependence of the density profiles of the LJ fluid is found to be very similar to the behavior of water near hydrophobic surfaces (Brovchenko et al. J.Phys.: Cond.Matt. v.16, 2004). A universal critical behavior of the local order parameter, defined as the difference between the local densities of the coexisting liquid and vapor phases at some distance z from the pore walls, Deltarho(z) = (rho_l(z) - rho_v(z))/2, is observed in a wide temperature range and found to be consistent with the surface critical behavior of the Ising model. Near the surface the dependence of the order parameter on the reduced temperature tau = (T_c - T)/T_c obeys a scaling law ~ tau^(beta_1) with a critical exponent beta_1 of about 0.8, corresponding to the ordinary surface transition. A crossover from bulk-like to surface-like critical behavior with increasing temperature occurs, when the correlation length is about half the distance to the surface. Relations between the ordinary and normal transitions in Ising systems and the surface critical behavior of fluids are discussed.Comment: 14 pages, 19 figures, submitted to PR

Connectedness percolation of hard convex polygonal rods and platelets

The properties of polymer composites with nanofiller particles change drastically above a critical filler density known as the percolation threshold. Real nanofillers, such as graphene flakes and cellulose nanocrystals, are not idealized disks and rods but are often modeled as such. Here we investigate the effect of the shape of the particle cross section on the geometric percolation threshold. Using connectedness percolation theory and the second-virial approximation, we analytically calculate the percolation threshold of hard convex particles in terms of three single-particle measures. We apply this method to polygonal rods and platelets and find that the universal scaling of the percolation threshold is lowered by decreasing the number of sides of the particle cross section. This is caused by the increase of the surface area to volume ratio with decreasing number of sides.Comment: 7 pages, 3 figures; added references, corrected typo, results unchange

L\'evy-type diffusion on one-dimensional directed Cantor Graphs

L\'evy-type walks with correlated jumps, induced by the topology of the medium, are studied on a class of one-dimensional deterministic graphs built from generalized Cantor and Smith-Volterra-Cantor sets. The particle performs a standard random walk on the sets but is also allowed to move ballistically throughout the empty regions. Using scaling relations and the mapping onto the electric network problem, we obtain the exact values of the scaling exponents for the asymptotic return probability, the resistivity and the mean square displacement as a function of the topological parameters of the sets. Interestingly, the systems undergoes a transition from superdiffusive to diffusive behavior as a function of the filling of the fractal. The deterministic topology also allows us to discuss the importance of the choice of the initial condition. In particular, we demonstrate that local and average measurements can display different asymptotic behavior. The analytic results are compared with the numerical solution of the master equation of the process.Comment: 9 pages, 9 figure

Computer simulation study of the closure relations in hard sphere fluids

We study, using Monte Carlo simulations, the cavity and the bridge functions of various hard sphere fluids: one component system, equimolar additive and non additive binary mixtures. In particular, we numerically check the assumption of local dependency of the bridge functions from the indirect correlation functions, on which most of the existing integral equation theories hinge. We find that this condition can be violated either in the region around the first and second neighbors shell, or inside the hard core, for the systems here considered. The violations manifest themselves clearly in the so called Duh-Haymet plots of the bridge functions versus the indirect correlation functions and become amplified as the coupling of the system increases.Comment: 24 pages, 13 figure

Three-dimensional molecular dynamics simulations of void coalescence during dynamic fracture of ductile metals

Void coalescence and interaction in dynamic fracture of ductile metals have been investigated using three-dimensional strain-controlled multi-million atom molecular dynamics simulations of copper. The correlated growth of two voids during the coalescence process leading to fracture is investigated, both in terms of its onset and the ensuing dynamical interactions. Void interactions are quantified through the rate of reduction of the distance between the voids, through the correlated directional growth of the voids, and through correlated shape evolution of the voids. The critical inter-void ligament distance marking the onset of coalescence is shown to be approximately one void radius based on the quantification measurements used, independent of the initial separation distance between the voids and the strain-rate of the expansion of the system. The interaction of the voids is not reflected in the volumetric asymptotic growth rate of the voids, as demonstrated here. Finally, the practice of using a single void and periodic boundary conditions to study coalescence is examined critically and shown to produce results markedly different than the coalescence of a pair of isolated voids.Comment: Accepted for publication in Physical Review

Osmotic force resisting chain insertion in a colloidal suspension

We consider the problem of inserting a stiff chain into a colloidal suspension of particles that interact with it through excluded volume forces. The free energy of insertion is associated with the work of creating a cavity devoid of colloid and sufficiently large to accomodate the chain. The corresponding work per unit length is the force that resists the entry of the chain into the colloidal suspension. In the case of a hard sphere fluid, this work can be calculated straightforwardly within the scaled particle theory; for solutions of flexible polymers, on the other hand, we employ simple scaling arguments. The forces computed in these ways are shown, for nanometer chain and colloid diameters, to be of the order of tens of pN for solution volume fraction for biophysical processes such as the ejection of DNA from viral capsids into the cell cytoplasm.Comment: 16 pages,3 figures. Accepted for publication in European Physical Journal