Researcher profile

L. Samaj

L. Samaj contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
8topics
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

17 published item(s)

preprint2019arXiv

Electric double layers with surface charge modulations: Novel exact Poisson-Boltzmann solutions

Poisson-Boltzmann theory is the cornerstone for soft matter electrostatics. We provide novel exact analytical solutions to this non-linear mean-field approach, for the diffuse layer of ions in the vicinity of a planar or a cylindrical macroion. While previously known solution are for homogeneously charged objects, the cases worked out exhibit a modulated surface charge --or equivalently surface potential-- on the macroion (wall) surface. In addition to asymptotic features at large distances from the wall, attention is paid to the fate of the contact theorem, relating the contact density of ions to the local wall charge density. For salt-free systems (counterions only), we make use of results pertaining to the two-dimensional Liouville equation, supplemented by an inverse approach. When salt is present, we invoke the exact two-soliton solution to the 2D sinh-Gordon equation. This leads to inhomogeneous charge patterns, that are either localized or periodic in space. Without salt, the electrostatic signature of a charge pattern on the macroion fades exponentially with distance for a planar macroion, while it decays as an inverse power-law for a cylindrical macroion. With salt, our study is limited to the planar geometry, and reveals that pattern screening is exponential.

preprint2016arXiv

Mean-field beyond mean-field: the single particle view for moderately to strongly coupled charged fluids

In a counter-ion only charged fluid, Coulomb coupling is quantified by a single dimensionless parameter. Yet, the theoretical treatment of moderately to strongly coupled charged fluids is a difficult task, central to the understanding of a wealth of soft matter problems, including biological systems. We show that the corresponding coupling regime can be remarkably well described by a single particle treatment, which, at variance with previous works, takes due account of inter-ionic interactions. To this end, the prototypical problem of a planar charged dielectric interface is worked out. Testing our predictions against Monte Carlo simulation data reveals an excellent agreement.

preprint2016arXiv

The rich polymorphic behavior of Wigner bilayers

Self-assembly into target structures is an efficient material design strategy. Combining analytical calculations and computational techniques of evolutionary and Monte Carlo types, we report about a remarkable structural variability of Wigner bilayer ground states, when charges are confined between parallel charged plates. Changing the inter-layer separation, or the plate charge asymmetry, a cascade of ordered patterns emerges. At variance with the symmetric case phenomenology, the competition between commensurability features and charge neutralization leads to long range attraction, appearance of macroscopic charges, exotic phases, and non conventional phase transitions with distinct critical indices, offering the possibility of a subtle, but precise and convenient control over patterns.

preprint2015arXiv

Effective charge of cylindrical and spherical colloids immersed in an electrolyte: the quasi-planar limit

We consider the non-linear Poisson-Boltzmann theory for a single cylindrical or spherical macro-ion in symmetric 1:1, together with asymmetric 1:2 and 2:1 electrolytes. We focus on the regime where $κa $, the ratio of the macro-ion radius $a$ over the inverse Debye length in the bulk electrolyte, is large. Analyzing the structure of the analytical expansion emerging from a multiple scale analysis, we uncover a hidden structure for the electrostatic potential. This structure, which appears after a heuristic resummation, suggests a new and convenient expansion scheme that we present and work out in detail. We show that novel exact results can thereby be obtained, in particular pertaining to effective charge properties, in complete agreement with the direct numerical solution to the problem.

preprint2014arXiv

Counter-ions between or at asymmetrically charged walls: 2D free-fermion point

This work contributes to the problem of determining effective interaction between asymmetrically (likely or oppositely) charged objects whose total charge is neutralized by mobile pointlike counter-ions of the same charge, the whole system being in thermal equilibrium. The problem is formulated in two spatial dimensions with logarithmic Coulomb interactions. The charged objects correspond to two parallel lines at distance d, with fixed line charge densities. Two versions of the model are considered: the standard "unconstrained" one with particles moving freely between the lines and the "constrained" one with particles confined to the lines. We solve exactly both systems at the free-fermion coupling and compare the results for the pressure (i.e. the force between the lines per unit length of one of the lines) with the mean-field Poisson-Boltzmann solution. For the unconstrained model, the large-d asymptotic behaviour of the free-fermion pressure differs from that predicted by the mean-field theory. For the constrained model, the asymptotic pressure coincides with the attractive van der Waals-Casimir fluctuational force. For both models, there are fundamental differences between the cases of likely-charged and oppositely-charged lines, the latter case corresponding at large distances d to a capacitor.

preprint2013arXiv

Ground state structure of a bilayer Wigner crystal with repulsive dielectric images

We study the ground-state structures of identical classical point charges with Coulomb interactions, confined between two symmetric parallel charged walls. For the well understood homogeneous dielectric case with no electrostatic images, the charges evenly condense on the opposite walls, thereby forming a bilayer Wigner crystal; five structures compete upon changing the inter-wall separation. Here, we consider a dielectric jump between the walls and a solvent in which charges are immersed, implying repulsive images. Using recently developed series representations of lattice sums for Coulomb law, we derive the complete phase diagram. In contrast to the homogeneous dielectric case, the particles remain in a hexagonal Wigner monolayer up to a certain distance between the walls. Beyond this distance, a bifurcation occurs to a sequence of Wigner bilayers, each layer having a nonzero spacing from the nearest wall. Another new phenomenon is that the ground-state energy as a function of the wall separation exhibits a global minimum.

preprint2012arXiv

Critical phenomena and phase sequence in classical bilayer Wigner crystal at zero temperature

We study the ground-state properties of a system of identical classical Coulombic point particles, evenly distributed between two equivalently charged parallel plates at distance $d$; the system as a whole is electroneutral. It was previously shown that upon increasing d from 0 to infinity, five different structures of the bilayer Wigner crystal become energetically favored, starting from a hexagonal lattice (phase I, d=0) and ending at a staggered hexagonal lattice (phase V, d -> infinity). In this paper, we derive new series representations of the ground-state energy for all five bilayer structures. The derivation is based on a sequence of transformations for lattice sums of Coulomb two-particle potentials plus the neutralizing background, having their origin in the general theory of Jacobi theta functions. The new series provide convenient starting points for both analytical and numerical progress. Its convergence properties are indeed excellent: Truncation at the fourth term determines in general the energy correctly up to 17 decimal digits. The accurate series representations are used to improve the specification of transition points between the phases and to solve a controversy in previous studies. In particular, it is shown both analytically and numerically that the hexagonal phase I is stable only at d=0, and not in a finite interval of small distances between the plates as was anticipated before. The expansions of the structure energies around second-order transition points can be done analytically, which enables us to show that the critical behavior is of the Ginzburg-Landau type, with a mean-field critical index beta=1/2 for the growth of the order parameters.

preprint2012arXiv

Ground state of classical bilayer Wigner crystals

We study the ground state structure of electronic-like bilayers, where different phases compete upon changing the inter-layer separation or particle density. New series representations with exceptional convergence properties are derived for the exact Coulombic energies under scrutiny. The complete phase transition scenario --including critical phenomena-- can subsequently be worked out in detail, thereby unifying a rather scattered or contradictory body of literature, hitherto plagued by the inaccuracies inherent to long range interaction potentials.

preprint2012arXiv

Like-charge colloidal attraction: a simple argument

By a length scale analysis, we study the equilibrium interactions between two like-charge planes confining neutralising counter-ions. At large Coulombic couplings, approaching the two charged bodies leads to an unbinding of counter-ions, a situation that is amenable to an exact treatment. This phenomenon is the key to attractive effective interactions. A particular effort is made for pedagogy, keeping equations and formalism to a minimum.

preprint2011arXiv

Strong-coupling theory for a polarizable planar colloid

We propose a strong-coupling analysis of a polarizable planar interface, in the spirit of a recently introduced Wigner-Crystal formulation. The system is made up of two moieties: a semi-infinite medium (z<0) with permittivity epsilon&#39; while the other half space in z>0 is occupied by a solution with permittivity epsilon, and mobile counter-ions (no added electrolyte). The interface at z=0 bears a uniform surface charge. The counter-ion density profile is worked out explicitly for both repulsive and attractive dielectric image cases.

preprint2011arXiv

Survival probability (heat content) and the lowest eigenvalue of Dirichlet Laplacian

We study the survival probability of a particle diffusing in a two-dimensional domain, bounded by a smooth absorbing boundary. The short-time expansion of this quantity depends on the geometric characteristics of the boundary, whilst its long-time asymptotics is governed by the lowest eigenvalue of the Dirichlet Laplacian defined on the domain. We present a simple algorithm for calculation of the short-time expansion for an arbitrary &#34;star-shaped&#34; domain. The coefficients are expressed in terms of powers of boundary curvature, integrated around the circumference of the domain. Based on this expansion, we look for a Padé interpolation between the short-time and the long-time behavior of the survival probability, i.e. between geometric characteristics of the boundary and the lowest eigenvalue of the Dirichlet Laplacian.

preprint2010arXiv

Counter-ions at Charged Walls: Two Dimensional Systems

We study equilibrium statistical mechanics of classical point counter-ions, formulated on 2D Euclidean space with logarithmic Coulomb interactions (infinite number of particles) or on the cylinder surface (finite particle numbers), in the vicinity of a single uniformly charged line (one single double-layer), or between two such lines (interacting double-layers). The weak-coupling Poisson-Boltzmann theory, which applies when the coupling constant Gamma is small, is briefly recapitulated (the coupling constant is defined as Gamma = beta e^2 where beta is the inverse temperature, and e the counter-ion charge). The opposite strong-coupling limit (Gamma -> infinity) is treated by using a recent method based on an exact expansion around the ground-state Wigner crystal of counter-ions. The weak- and strong-coupling theories are compared at intermediary values of the coupling constant Gamma=2 gamma (gamma=1,2,3), to exact results derived within a 1D lattice representation of 2D Coulomb systems in terms of anti-commuting field variables. The models (density profile, pressure) are solved exactly for any particles numbers N at Gamma=2 and up to relatively large finite N at Gamma=4 and 6. For the one-line geometry, the decay of the density profile at asymptotic distance from the line undergoes a fundamental change with respect to the mean-field behavior at Gamma=6. The like-charge attraction regime, possible in the strong coupling limit but precluded at mean-field level, survives for Gamma=4 and 6, but disappears at Gamma=2.

preprint2010arXiv

Strong-Coupling Theory of Counter-ions at Charged Plates

We present an analytical approach to the strong coupling regime of similarly and highly charged plates in the presence of counter-ions. The procedure is physically transparent and based on an exact expansion around the ground state formed by the two-dimensional Wigner crystal of counter-ions. The one plate problem is worked out, together with the two plates situation. Unlike previous approaches, the expansion is free of divergences, and is shown to be in excellent agreement with available data of Monte-Carlo simulations under strong Coulombic couplings. The present results shed light on the like-charge attraction regime.

preprint2005arXiv

Saturation of Electrostatic Potential: Exactly Solvable 2D Coulomb Models

We test the concepts of renormalized charge and potential saturation, introduced within the framework of highly asymmetric Coulomb mixtures, on exactly solvable Coulomb models. The object of study is the average electrostatic potential induced by a unique &#34;guest&#34; charge immersed in a classical electrolyte, the whole system being in thermal equilibrium at some inverse temperature $β$. The guest charge is considered to be either an infinite hard wall carrying a uniform surface charge or a charged colloidal particle. The systems are treated as two-dimensional; the electrolyte is modelled by a symmetric two-component plasma of point-like $\pm e$ charges with logarithmic Coulomb interactions. Two cases are solved exactly: the Debye-Hückel limit $βe^2\to 0$ and the Thirring free-fermion point $βe^2 = 2$. The results at the free-fermion point can be summarized as follows: (i) The induced electrostatic potential exhibits the asymptotic behavior, at large distances from the guest charge, whose form is different from that obtained in the Debye-Hückel (linear Poisson-Boltzmann) theory. This means that the concept of renormalized charge, developed within the nonlinear Poisson-Boltzmann theory to describe the screening effect of the electrolyte cloud, fails at the free-fermion point. (ii) In the limit of an infinite bare charge, the induced electrostatic potential saturates at a finite value in every point of the electrolyte region. This fact confirms the previously proposed hypothesis of potential saturation.

preprint2004arXiv

Translation Symmetry Breaking in the One-Component Plasma on the Cylinder

The two-dimensional one-component plasma, i.e. the system of point-like charged particles embedded in a homogeneous neutralizing background, is studied on the surface of a cylinder of finite circumference, or equivalently in a semiperiodic strip of finite width. The model has been solved exactly by Choquard et al. at the free-fermion coupling $Γ=2$: in the thermodynamic limit of an infinitely long strip, the particle density turns out to be a nonconstant periodic function in space and the system exhibits long-range order of the Wigner-crystal type. The aim of this paper is to describe, qualitatively as well as quantitatively, the crystalline state for a larger set of couplings $Γ=2γ$ ($γ=1,2,...$ a positive integer) when the plasma is mappable onto a one-dimensional fermionic theory. The fermionic formalism, supplemented by some periodicity assumptions, reveals that the density profile results from a hierarchy of Gaussians with a uniform variance but with different amplitudes. The number and spatial positions of these Gaussians within an elementary cell depend on the particular value of $γ$. Analytic results are supported by the exact solution at $γ=1$ ($Γ=2$) and by exact finite-size calculations at $γ=2,3$.

preprint2003arXiv

Universal behavior of quantum Green&#39;s functions

We consider a general one-particle Hamiltonian H = - Δ_r + u(r) defined in a d-dimensional domain. The object of interest is the time-independent Green function G_z(r,r&#39;) = < r | (z-H)^{-1} | r&#39; >. Recently, in one dimension (1D), the Green&#39;s function problem was solved explicitly in inverse form, with diagonal elements of Green&#39;s function as prescribed variables. The first aim of this paper is to extract from the 1D inverse solution such information about Green&#39;s function which cannot be deduced directly from its definition. Among others, this information involves universal, i.e. u(r)-independent, behavior of Green&#39;s function close to the domain boundary. The second aim is to extend the inverse formalism to higher dimensions, especially to 3D, and to derive the universal form of Green&#39;s function for various shapes of the confining domain boundary.