BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Yan Levin

Yan Levin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softcond-mat.stat-mechastro-ph.GAphysics.chem-phphysics.plasm-phBiological Physicscond-matcond-mat.mes-hallBiomoleculescond-mat.mtrl-scicond-mat.othercond-mat.str-elOther Quantitative Biologyphysics.acc-phphysics.gen-ph

Catalog footprint

What is connected

52works
15topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

52 published item(s)

preprint2022arXiv

Widom insertion method in simulations with Ewald summation

We discuss the application of Widom insertion method for calculation of the chemical potential of individual ions in computer simulations with Ewald summation. Two approaches are considered. In the first approach, an individual ion is inserted into a periodically replicated overall charge neutral system representing an electrolyte solution. In the second approach an inserted ion is also periodically replicated, leading to the violation of the overall charge neutrality. This requires an introduction of an additional neutralizing background. We find that the second approach leads to a much better agreement with the results of grand canonical Monte Carlo simulation for the total chemical potential of a neutral ionic cluster.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Charge regulation of colloidal particles in aqueous solutions

We study charge regulation of colloidal particles inside aqueous electrolyte solutions. To stabilize colloidal suspension against precipitation, colloidal particles are synthesized with either acidic or basic groups on their surface. In contact with water these surface groups undergo proton transfer reaction, resulting in colloidal surface charge. The charge is determined by the condition of local chemical equilibrium between hydronium ions inside the solution and at the colloidal surface. We use a model of Baxter sticky spheres to explicitly calculate the equilibrium dissociation constants and to construct a theory which is able to quantitatively predict the effective charge of colloidal particles with either acidic or basic surface groups. The predictions of the theory for the model are found to be in excellent agreement with the results of Monte Carlo simulations. The theory is further extended to treat colloidal particles with a mixture of both acidic and basic surface groups.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Entropy production and Vlasov equation for self-gravitating systems

The evolution of a self-gravitating system to a non-equilibrium steady state occurs through a process of violent relaxation. In the thermodynamic limit the dynamics of a many body system should be governed by the Vlasov equation. Recently, however, a question was raised regarding the validity of Vlasov equation during the process of violent relaxation. In this paper we will explore the entropy production during the relaxation process using N-body molecular dynamics simulations. We will show that the entropy production time grows as $N^α$, with $α> 0$ and in the limit $N \rightarrow \infty$, entropy will remain constant, consistent with the Vlasov equation. Furthermore, we will show that the mean field dynamics constructed on the basis of the Vlasov equation is in excellent agreement with the full molecular dynamics simulations, justifying the applicability of Vlasov equation during the violent relaxation phase of evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Osmotic stress and pore nucleation in charged biological nanoshells and capsids

A model system is proposed to investigate the chemical equilibrium and mechanical stability of biological spherical-like nanoshells in contact with an aqueous solution with added dissociated electrolyte of a given concentration. The ionic chemical equilibrium across the permeable shell is investigated in the framework of an accurate Density Functional Theory (DFT) that incorporates electrostatic and hardcore correlations beyond the traditional mean-field ({\it e. g.}, Poisson-Boltzmann) limit. The accuracy of the theory is tested by a direct comparison with Monte Carlo (MC) simulations. A simple analytical expression is then deduced which clearly highlights the entropic, electrostatic, and self-energy contributions to the osmotic stress over the shell in terms of the calculated ionic profiles. By invoking a continuum mean-field elastic approach to account for the shell surface stress upon osmotic stretching, the mechanical equilibrium properties of the shell under a wide variety of ionic strengths and surface charges are investigated. The model is further coupled to a continuum mechanical approach similar in structure to a Classical Nucleation Theory (CNT) to address the question of mechanical stability of the shells against a pore nucleation. This allows us to construct a phase diagram which delimits the mechanical stability of capsids for different ionic strengths and shell surface charges.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Thermodynamic collapse in a lattice-gas model for a two-component system of penetrable particles

We study a lattice-gas model of penetrable particles on a square-lattice substrate with same-site and nearest-neighbor interactions. Penetrability implies that the number of particles occupying a single lattice site is unlimited and the model itself is intended as a simple representation of penetrable particles encountered in realistic soft-matter systems. Our specific focus is on a binary mixture, where particles of the same species repel and those of the opposite species attract each other. As a consequence of penetrability and the unlimited occupation of each site, the system exhibits thermodynamic collapse, which in simulations is manifested by an emergence of extremely dense clusters scattered throughout the system with energy of a cluster $E\propto -n^2$ where $n$ is the number of particles in a cluster. After transforming a particle system into a spin system, in the large density limit the Hamiltonian recovers a simple harmonic form, resulting in the discrete Gaussian model used in the past to model the roughening transition of interfaces. For finite densities, due to the presence of a non-harmonic term, the system is approximated using a variational Gaussian model.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Charge neutrality breakdown in confined aqueous electrolytes: theory and simulation

We study, using Density Functional theory and Monte Carlo simulations, aqueous electrolyte solutions between charged infinite planar surfaces, in a contact with a bulk salt reservoir. In agreement with recent experimental observations [Z. Luo et al., Nat. Comm. 6, 6358 (2015)], we find that the confined electrolyte lacks local charge neutrality. We show that a Density Functional Theory (DFT) based on a bulk- HNC expansion properly accounts for strong electrostatic correlations and allows us to accurately calculate the ionic density profiles between the charged surfaces, even for electrolytes containing trivalent counterions. The DFT allows us to explore the degree of local charge neutrality violation, as a function of plate separation and bulk electrolyte concentration, and to accurately calculate the interaction force between the charged surfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Influence of network topology on the swelling of polyelectrolyte nanogels

It is well-known that the swelling behavior of ionic nanogels depends on their cross-link density, however it is unclear how different topologies should affect the response of the polyelectrolyte network. Here we perform Monte Carlo simulations to obtain the equilibrium properties of ionic nanogels as a function of salt concentration $C_s$ and the fraction $f$ of ionizable groups in a polyelectrolyte network formed by cross-links of functionality $z$. Our results indicate that the network with cross-links of low connectivity result in nanogel particles with higher swelling ratios. We also confirm a de-swelling effect of salt on nanogel particles.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Electrolytes between dielectric charged surfaces: Simulations and theory

We present a simulation method to study electrolyte solutions in a dielectric slab geometry using a modified 3D Ewald summation. The method is fast and easy to implement, allowing us to rapidly resum an infinite series of image charges. In the weak coupling limit, we also develop a mean-field theory which allows us to predict the ionic distribution between the dielectric charged plates. The agreement between both approaches, theoretical and simulational, is very good, validating both methods. Examples of ionic density profiles in the strong electrostatic coupling limit are also presented. Finally, we explore the confinement of charge asymmetric electrolytes between neutral surfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Interaction between random heterogeneously charged surfaces in an electrolyte solution

We study, using Monte Carlo simulations, the interaction between infinite heterogeneously charged surfaces inside an electrolyte solution. The surfaces are overall neutral with quenched charged domains. An average over the quenched disorder is performed to obtain the net force. We find that the interaction between the surfaces is repulsive at short distances and is attractive for larger separations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Adsorption of cationic polyions to a hydrophobic surface in the presence of Hofmeister salts

We study, using extensive Monte Carlos simulations, the behavior of cationic polyelectrolytes near hydrophobic surfaces in solutions containing Hofmeister salts. The Hofmeister anions are divided into kosmotropes and chaotropes. Near a hydrophobic surface, the chaotropes lose their solvation sheath and become partially adsorbed to the interface, while the kosmotropes remain strongly hydrated and are repelled from the interface. If the polyelectrolyte solution contains chaotropic anions, a significant adsorption of polyions to the surface is also observed. On the other hand, the kosmotropic anions have only a small influence on the polyion adsorption. These findings can have important implications for exploring the antibacterial properties of cationic polyelectrolytes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Ensemble inequivalence in a mean-field XY model with ferromagnetic and nematic couplings

We explore ensemble inequivalence in long-range interacting systems by studying an XY model of classical spins with ferromagnetic and nematic coupling. We demonstrate the inequivalence by mapping the microcanonical phase diagram onto the canonical one, and also by doing the inverse mapping. We show that the equilibrium phase diagrams within the two ensembles strongly disagree within the regions of first-order transitions, exhibiting interesting features like temperature jumps. In particular, we discuss the coexistence and forbidden regions of different macroscopic states in both the phase diagrams.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Equilibrium properties of charged microgels: a Poisson-Boltzmann-Flory approach

The equilibrium properties of ionic microgels are investigated using a combination of the Poisson-Boltzmann and Flory theories. Swelling behavior, density profiles, and effective charges are all calculated in a self-consistent way. Special attention is given to the effects of salinity on these quantities. It is found that the equilibrium microgel size is strongly influenced by the amount of added salt. Increasing the salt concentration leads to a considerable reduction of the microgel volume, which therefore releases its internal material -- solvent molecules and dissociated ions -- into the solution. Finally, the question of charge renormalization of ionic microgels in the context of the cell model is briefly addressed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Ion Specificity and Micellization of Ionic Surfactants: A Monte Carlo Study

We develop a simulation method which allows us to calculate the critical micelle concentrations for ionic surfactants in the presence of different salts. The results are in good agreement with the experimental data. The simulations are performed on a simple cubic lattice. The anionic interactions with the alkyl chains are taken into account based on the previously developed theory of the interfacial tensions of hydrophobic interfaces: the kosmotropic anions do not interact with the hydrocarbon tails of ionic surfactants, while chaotropic anions interact with the alkyl chains through a dispersion potential proportional to the anionic polarizability.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Kosterlitz-Thouless and Potts transitions in a generalized XY model

We present extensive numerical simulations of a generalized XY model with nematic-like terms recently proposed by Poderoso {\it et al} [PRL 106(2011)067202]. Using finite size scaling and focusing on the $q=3$ case, we locate the transitions between the paramagnetic (P), the nematic-like (N) and the ferromagnetic (F) phases. The results are compared with the recently derived lower bounds for the P-N and P-F transitions. While the P-N transition is found to be very close to the lower bound, the P-F transition occurs significantly above the bound. Finally, the transition between the nematic-like and the ferromagnetic phases is found to belong to the 3-states Potts universality class.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nonequilibrium stationary states of 3D self-gravitating systems

Three dimensional self-gravitating systems do not evolve to thermodynamic equilibrium, but become trapped in nonequilibrium quasistationary states. In this Letter we present a theory which allows us to a priori predict the particle distribution in a final quasistationary state to which a self-gravitating system will evolve from an initial condition which is isotropic in particle velocities and satisfies a virial constraint 2K=-U, where K is the total kinetic energy and U is the potential energy of the system.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Yukawa particles in a confining potential

We study the density distribution of repulsive Yukawa particles confined by an external potential. In the weak coupling limit, we show that the mean-field theory is able to accurately account for the particle distribution. In the strong coupling limit, the correlations between the particles become important and the mean-field theory fails. For strongly correlated systems, we construct a density functional theory which provides an excellent description of the particle distribution, without any adjustable parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nonequilibrium phase transitions in systems with long-range interactions

We introduce a generalized Hamiltonian Mean Field Model (gHMF)-XY model with both linear and quadratic coupling between spins and explicit Hamiltonian dynamics. In addition to the usual paramagnetic and ferromagnetic phases, this model also possesses a nematic phase. The gHMF can be solved explicitly using Boltzmann-Gibbs (BG) statistical mechanics, in both canonical and microcanonical ensembles. However, when the resulting microcanonical phase diagram is compared with the one obtained using molecular dynamics simulations, it is found that the two are very different. We will present a dynamical theory which allows us to explicitly calculate the phase diagram obtained using molecular dynamics simulations without any adjustable parameters. The model illustrates the fundamental role played by dynamics as well the inadequacy of BG statistics for systems with long-range forces in the thermodynamic limit.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nonequilibrium Statistical Mechanics of Systems with Long-Range Interactions: Ubiquity of Core-Halo Distributions

Systems with long-range (LR) forces, for which the interaction potential decays with the interparticle distance with an exponent smaller than the dimensionality of the embedding space, remain an outstanding challenge to statistical physics. The internal energy of such systems lacks extensivity and additivity. Although the extensivity can be restored by scaling the interaction potential with the number of particles, the non-additivity still remains. Lack of additivity leads to inequivalence of statistical ensembles. Before relaxing to thermodynamic equilibrium, isolated systems with LR forces become trapped in out-of-equilibrium quasi-stationary state (qSS), the lifetime of which diverges with the number of particles. Therefore, in thermodynamic limit LR systems will not relax to equilibrium. The qSSs are attained through the process of collisionless relaxation. Density oscillations lead to particle-wave interactions and excitation of parametric resonances. The resonant particles escape from the main cluster to form a tenuous halo. Simultaneously, this cools down the core of the distribution and dampens out the oscillations. When all the oscillations die out the ergodicity is broken and a qSS is born. In this report, we will review a theory which allows us to quantitatively predict the particle distribution in the qSS. The theory is applied to various LR interacting systems, ranging from plasmas to self-gravitating clusters and kinetic spin models.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Simulation and Theory of Ions at Atmospherically Relevant Aqueous Liquid-Air Interfaces

Chemistry occurring at or near the surfaces of aqueous droplets and thin films in the atmosphere influences air quality and climate. Molecular dynamics simulations are becoming increasingly useful for gaining atomic-scale insight into the structure and reactivity of aqueous interfaces in the atmosphere. Here we review simulation studies of atmospherically relevant aqueous liquid-air interfaces, with an emphasis on ions that play important roles in the chemistry of atmospheric aerosols. In addition to surveying results from simulation studies, we discuss challenges to the refinement and experimental validation of the methodology for simulating ion adsorption to the air-water interface, and recent advances in elucidating the driving forces for adsorption. We also review the recent development of a dielectric continuum theory that is capable of reproducing simulation and experimental data on ion behavior at aqueous interfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Symmetry Breaking in d-Dimensional Self-gravitating Systems

Systems with long-range interactions, such as self-gravitating clusters and magnetically confined plasmas, do not relax to the usual Boltzmann-Gibbs thermodynamic equilibrium, but become trapped in quasi-stationary states (QSS) the life time of which diverges with the number of particles. The QSS are characterized by the lack of ergodicity which can result in a symmetry broken QSS starting from a spherically symmetric particle distribution. We will present a theory which allows us to quantitatively predict the instability threshold for spontaneous symmetry breaking for a class of d-dimensional self-gravitating systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The double-layer of penetrable ions: an alternative route to charge reversal

We investigate a double-layer of penetrable ions near a charged wall. We find a new mechanism for charge reversal that occurs in the weak-coupling regime and, accordingly, the system is suitable for the mean-field analysis. The penetrability is achieved by smearing-out the ionic charge inside a sphere, so there is no need to introduce non-electrostatic forces and the system in the low coupling limit can be described by a modified version of the Poisson-Boltzmann equation. The predictions of the theory are compared with the Monte Carlo simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Electrochemical Surface Potential Due to Classical Point Charge Models Drives Anion Adsorption to the Air-Water Interface

We demonstrate that the driving forces for ion adsorption to the air-water interface for point charge models results from both cavitation and a term that is of the form of a negative electrochemical surface potential. We carefully characterize the role of the free energy due to the electrochemical surface potential computed from simple empirical models and its role in ionic adsorption within the context of dielectric continuum theory. Our research suggests that the electrochemical surface potential due to point charge models provides anions with a significant driving force to the air-water interface. This is contrary to the results of ab initio simulations that indicate that the average electrostatic surface potential should favor the desorption of anions at the air-water interface. The results have profound implications for the studies of ionic distributions in the vicinity of hydrophobic surfaces and proteins.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Topology of Collisionless Relaxation

Using extensive molecular dynamics simulations we explore the fine-grained phase space structure of systems with long-range interactions. We find that if the initial phase space particle distribution has no holes, the final stationary distribution will also contain a compact simply connected region. The microscopic holes created by the filamentation of the initial distribution function are always restricted to the outer regions of the phase space. In general, for complex multilevel distributions it is very difficult to a priori predict the final stationary state without solving the full dynamical evolution. However, we show that for multilevel initial distributions satisfying a generalized virial condition, it is possible to predict the particle distribution in the final stationary state using Casimir invariants of the Vlasov dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Vortex distribution in a confining potential

We study a model of interacting vortices in a type II superconductor. In the weak coupling limit, we constructed a mean-field theory which allows us to accurately calculate the vortex density distribution inside a confining potential. In the strong coupling limit, the correlations between the particles become important and the mean-field theory fails. Contrary to recent suggestions, this does not imply failure of the Boltzmann-Gibbs statistical mechanics, as we clearly demonstrate by comparing the results of Molecular Dynamics and Monte Carlo simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A close look into the excluded volume effects within a double layer

We explore the effect of steric interaction on the ionic density distribution near a charged hard wall. For weakly charged walls, small particles, and monovalent ions the mean-field Poisson-Boltzmann equation provides an excellent description of the density profiles. For large ions and large surface charges, however, deviations appear. To explore these, we use the density functional theory. We find that local density functionals are not able to account for steric interactions near a wall. Based on the weighted density approximation we derive a simple analytical expression for the contact electrostatic potential which allows us to analytically calculate the differential capacitance of the double layer.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Equation of state of charged colloidal suspensions and its dependence on the thermodynamic route

The thermodynamic properties of highly charged colloidal suspensions in contact with a salt reservoir are investigated in the framework of the Renormalized Jellium Model (RJM). It is found that the equation of state is very sensitive to the particular thermodynamic route used to obtain it. Specifically, the osmotic pressure calculated within the RJM using the contact value theorem can be very different from the pressure calculated using the Kirkwood-Buff fluctuation relations. On the other hand, Monte Carlo (MC) simulations show that both the effective pair potentials and the correlation functions are accurately predicted by the RJM. It is suggested that the lack of self-consistency in the thermodynamics of the RJM is a result of neglected electrostatic correlations between the counterions and coions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ergodicity Breaking and Parametric Resonances in Systems with Long-Range Interactions

We explore the mechanism responsible for the ergodicity breaking in systems with long-range forces. In thermodynamic limit such systems do not evolve to the Boltzmann-Gibbs equilibrium, but become trapped in an out-of-equilibrium quasi-stationary-state. Nevertheless, we show that if the initial distribution satisfies a specific constraint - a generalized virial condition - the quasi- stationary-state is very close to ergodic and can be described by Lynden-Bell statistics. On the other hand if the generalized virial condition is violated, parametric resonances are excited, leading to chaos and ergodicity breaking.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ion fluxes through nano-pores and transmembrane channels

We introduce an implicit solvent Molecular Dynamics approach for calculating ionic fluxes through narrow nano-pores and transmembrane channels. The method relies on a dual-control- volume grand-canonical molecular dynamics (DCV-GCMD) simulation and the analytical solution for the electrostatic potential inside a cylindrical nano-pore recently obtained by Levin [Europhys. Lett., 76, 163 (2006)]. The theory is used to calculate the ionic fluxes through an artificial trans-membrane c hannel which mimics the antibacterial gramicidin A channel. Both current-voltage and current-concentration relations are calculated under various experimental conditions. We show that our results are comparable to the characteristics associated to the gramicidin A pore, specially the existence of two binding sites inside the pore and the observed saturation in the current-concentration profiles.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Polarizable Surfaces: Weak and Strong Coupling Regimes

We study the ionic distribution near a charged surface. A new method for performing Monte Carlo simulations in this geometry is discussed. A theory is then presented that allows us to accurately reproduce the density profiles obtained in the simulations. In the weak-coupling regime, a theory accounts for the ion-image interactions, leading to a modified Poisson-Boltzmann equation. When the correlations between the ions are significant, a strong-coupling theory is used to calculate the density profiles near the surface and a Poisson-Boltzmann equation with a renormalized boundary condition to account for the counterion distribution in the far-field.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Comment on: Thermostatistics of Overdamped Motion of Interacting Particles [arXiv:1008.1421]

In a recent paper, Phys. Rev. Lett. 105 260601 (2010) [arXiv:1008.1421], Andrade et al., argued that classical particles confined in a parabolic trap at T=0 distribute themselves in accordance with the Tsallis statistics. To prove their point the authors performed molecular dynamics simulations. Here we show that the model of Andrade et al. can be solved exactly. The distribution of particles at T=0 has nothing to do with the Tsallis entropy and is determined simply by the force balance.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Core-halo distribution in the Hamiltonian Mean-Field Model

We study a paradigmatic system with long-range interactions: the Hamiltonian Mean-Field Model (HMF). It is shown that in the thermodynamic limit this model does not relax to the usual equilibrium Maxwell-Boltzmann distribution. Instead, the final stationary state has a peculiar core-halo structure. In the thermodynamic limit, HMF is neither ergodic nor mixing. Nevertheless, we find that using dynamical properties of Hamiltonian systems, it is possible to quantitatively predict both the spin distribution and the velocity distribution functions in the final stationary state, without any adjustable parameters. We also show that HMF undergoes a non-equilibrium first-order phase transition between paramagnetic and ferromagnetic states.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Effects of the dielectric discontinuity on the counterion distribution in a colloidal suspension

We introduce a new method for simulating colloidal suspensions with spherical colloidal particles of dielectric constant different from the surrounding medium. The method uses exact calculation of the Green function to obtain the ion-ion interaction potential in the presence of a dielectric discontinuity at the surface of the colloidal particle. The new method is orders of magnitude faster than the traditional approaches based on series expansions of the interaction potential.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Ion specificity and the theory of stability of colloidal suspensions

A theory is presented which allow us to accurately calculate the critical coagulation concentration (CCC) of hydrophobic colloidal suspensions. For positively charged particles the CCC's follow the Hofmeister (lyotropic) series. For negatively charged particles the series is reversed. We find that strongly polarizable chaotropic anions are driven towards the colloidal surface by electrostatic and hydrophobic forces. Within approximately one ionic radius from the surface, the chaotropic anions loose part of their hydration sheath and become strongly adsorbed. The kosmotropic anions, on the other hand, are repelled from the hydrophobic surface. The theory is quantitatively accurate without any adjustable parameters. We speculate that the same mechanism is responsible for the Hofmeister series that governs stability of protein solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Rejoinder on: Thermostatistics of Overdamped Motion of Interacting Particles

In their Reply [arXiv:1104.5036] to our Comment [arXiv:1104.0697] Andrade et al. state that we have "chosen to categorically dismiss their elaborate and solid conceptual approach without employing any concepts or tools from Statistical Mechanics". Because of the one page restriction imposed by PRL, we could not address all of the flaws of the original paper. Therefore, we are grateful to Andrade et al., for giving us an opportunity to further elaborate on our Comment. In this Rejoinder we show that the model studied by Andrade et al. has nothing to do with the Tsallis entropy either at T=0 or at any other temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Statistical Mechanics of 1d Self-Gravitating Systems: The Core-Halo Distribution

We study, using both theory and simulations, a system of self-gravitating sheets. A new statistical mechanics theory - free of any adjustable parameters - is derived to quantitatively predict the final stationary state achieved by this system after the process of collisionless relaxation is completed. The theory shows a very good agreement with the numerical simulations. The model sheds new light on the general mechanism of relaxation of self-gravitating systems and may help us to understand cold matter distribution in the Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Weak and Strong Coupling Theories for Polarizable Colloids and Nano-Particles

A theory is presented which allows us to accurately calculate the density profile of monovalent and multivalent counterions in suspensions of polarizable colloids or nano-particles. In the case of monovalent ions, we derive a weak-coupling theory that explicitly accounts for the ion-image interaction, leading to a modified Poisson-Boltzmann equation. For suspensions with multivalent counterions, a strong-coupling theory is used to calculate the density profile near the colloidal surface and a Poisson-Boltzmann equation with a renormalized boundary condition to account for the counterion distribution in the far-field. All the results are compared with the Monte Carlo simulations, showing an excellent agreement between the theory and the simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Statistical Mechanics of Unbound Two Dimensional Self-Gravitating Systems

We study, using both theory and molecular dynamics simulations, the relaxation dynamics of a microcanonical two dimensional self-gravitating system. After a sufficiently large time, a gravitational cluster of N particles relaxes to the Maxwell-Boltzmann distribution. The time to reach the thermodynamic equilibrium, however, scales with the number of particles. In the thermodynamic limit, $N\to\infty$ at fixed total mass, equilibrium state is never reached and the system becomes trapped in a non-ergodic stationary state. An analytical theory is presented which allows us to quantitatively described this final stationary state, without any adjustable parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Surface tensions and surface potentials of acid solutions

A theory is presented which allows us to quantitatively calculate the excess surface tension of acid solutions. The H^+, in the form of hydronium ion, is found to be strongly adsorbed to the solution-air interface. To account for the electrostatic potential difference measured experimentally, it is necessary to assume that the hydronium ion is oriented with its hydrogens pointing into the bulk water. The theory is quantitatively accurate for surface tensions and is qualitative for electrostatic potential difference across the air-water interface.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Surface tensions, surface potentials and the Hofmeister series of electrolyte solutions

A theory is presented which allows us to accurately calculate the surface tensions and the surface potentials of electrolyte solutions. Both the ionic hydration and the polarizability are taken into account. We find a good correlation between the Jones-Dole viscosity $B$-coefficient and the ionic hydration near the air-water interface. The kosmotropic anions such as fluoride, iodate, sulfate and carbonate, are found to be strongly hydrated and are repelled from the interface. The chaotropic anions such as perchlorate, iodide, chlorate and bromide are found to be significantly adsorbed to the interface. Chloride and bromate anions become weakly hydrated in the interfacial region. The sequence of surface tensions and surface potentials is found to follow the Hofmeister ordering. The theory, with only one adjustable parameter, quantitatively accounts for the surface tensions of 10 sodium salts for which there is experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The renormalized Jellium model of colloidal suspensions with multivalent counterions

An extension of the renormalized Jellium model which allows to study colloidal suspensions containing trivalent counterions is proposed. The theory is based on a modified Poisson-Boltzmann equation which incorporates the effects of counterion correlations near the colloidal surfaces using a new boundary condition. The renormalized charges, the counterion density profiles, and osmotic pressures can be easily calculated using the modified renormalized Jellium model. The results are compared with the ones obtained using the traditional Wigner-Seitz (WS) cell approximation also with a new boundary condition. We find that while the thermodynamic functions obtained within the renormalized Jellium model are in a good agreement with their WS counterpart, the effective charges predicted by the two theories can be significantly different.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Polarizable ions at interfaces

A non-perturbative theory is presented which allows to calculate the solvation free energy of polarizable ions near a water-vapor and water-oil interfaces. The theory predicts that larger halogen anions are adsorbed at the interface, while the alkali metal cations are repelled from it. The density profiles calculated theoretically are similar to those obtained using the molecular dynamics simulations with polarizable force fields.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Relaxation and Emittance Growth of a Thermal Charged-Particle Beam

We present a theory that allows us to accurately calculate the distribution functions and the emittance growth of a thermal charged-particle beam after it relaxes to equilibrium. The theory can be used to obtain the fraction of particles, which will evaporate from the beam to form a halo. The calculated emittance growth is found to be in excellent agreement with the simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Electromagnetic instability of the Thomson Problem

The classical Thomson problem of $n$ charged particles confined to the surface of a sphere of radius $a$ is analyzed within the Darwin approximation of electrodynamics. For $n<n_c(a)$ the ground state corresponds to a hexagonal Wigner crystal with a number of topological defects. However, if $n>n_c(a)$ the Wigner lattice is unstable with respect to small perturbations and the ground state becomes spontaneously magnetized for finite $n$.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Pore dynamics of osmotically stressed vesicles

We present a theory for pore dynamics of osmotically stressed vesicles. When a liposome with an internal concentration of solute is placed inside a solute-depleted medium, an osmotic flow of solvent through the lipid bilayer leads to swelling of vesicle and to increase in membrane surface tension. This can result in membrane rupture and opening of thermal pores. Depending on the internal concentration of solute and the size of the vesicle, pores can close rapidly or be long-lived. We find that the life span of the long-lived pores scales non-trivially with the size of the liposome. Closure of the long-lived pore is followed by a rapid flicker-like opening and closing of short-lived pores. Our model is consistent with the observation of long-lived pores in red blood cell ghosts.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Density-functional theory for attraction between like-charged plates

We study the interactions between two negatively charged macroscopic surfaces confining positive counterions. A density-functional approach is introduced which, besides the usual mean-field interactions, takes into account the correlations in the positions of counterions. The excess free energy is derived in the framework of the Debye-H{ü}ckel theory of the one-component plasma, with the homogeneous density replaced by a weighted density. The minimization of the total free energy yields the density profile of the microions. The pressure is calculated and compared with the simulations and the results derived from integral equations theories. We find that the interaction between the two plates becomes attractive when their separation distance is sufficiently small and the surface charge density is larger than a threshold value.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

The one-component plasma: a conceptual approach

The one-component plasma (OCP) represents the simplest statistical mechanical model of a Coulomb system. For this reason, it has been extensively studied over the last forty years. The advent of the integral equations has resulted in a dramatic improvement in our ability to carry out numerical calculations, but came at the expense of a physical insight gained in a simpler analytic theory. In this paper we present an extension of the Debye-Hueckel (DH) theory to the OCP. The theory allows for analytic calculations of all the thermodynamic functions, as well as the structure factor. The theory explicitly satisfies the Stillinger-Lovett and, for small couplings, the compressibility sum rules, implying its internal self consistency.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Debye-Hueckel-Bjerrum theory for charged colloids

We formulate an extension of the Debye-Hueckel-Bjerrum theory [M. E. Fisher and Y. Levin, Phys. Rev. Lett. 71, 3826 (1993)] to the fluid state of a highly asymmetric charged colloid. Allowing for the formation of clusters consisting of one polyion and n condensed counterions, the total Helmholtz free energy of the colloidal suspension is constructed. The thermodynamic properties, such as the cluster-density distribution and the pressure, are obtained by the minimization of the free energy under the constraints of fixed number of polyions and counterions. In agreement with the current experimental and Monte Carlo results, no evidence of any phase transition is encountered.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Polyelectrolyte Solutions with Multivalent Salts

We investigate the thermodynamic properties of a polyelectrolyte solution in a presence of {\it multivalent} salts. The polyions are modeled as rigid cylinders with the charge distributed uniformly along the major axis. The solution, besides the polyions, contain monovalent and divalent counterions as well as monovalent coions. The strong electrostatic attraction existing between the polyions and the counterions results in formation of clusters consisting of one polyion and a number of associated monovalent and divalent counterions. The theory presented in the paper allows us to explicitly construct the Helmholtz free energy of a polyelectrolyte solution. The characteristic cluster size, as well as any other thermodynamic property can then be determined by an appropriate operation on the free energy.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Sine-Gordon mean field theory of a Coulomb Gas

Sine-Gordon field theory is used to investigate the phase diagram of a neutral Coulomb gas. A variational mean field free energy is constructed and the corresponding phase diagrams in two (2d) and three dimensions (3d) are obtained. When analyzed in terms of chemical potential, the Sine-Gordon theory predicts the phase diagram topologically identical with the Monte Carlo simulations and a recently developed Debye-Hückel-Bjerrum (DHBj) theory. In 2d we find that the infinite order Kosterlitz-Thouless line terminates in a tricritical point, after which the metal-insulator transition becomes first order. However, when the transformation from chemical potential to the density is made the whole of the insulating phase is mapped onto zero density.

0 citations0 reviewsNo code linkedOpen access