BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Yan Levin

Yan Levin contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.softcond-mat.stat-mechphysics.chem-phastro-ph.GABiological PhysicsBiomoleculescond-mat.mtrl-sciphysics.acc-phphysics.gen-phphysics.plasm-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
23works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

23 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
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

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

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
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

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

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