Source author record

Ladislav Šamaj

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

8works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

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

8 published item(s)

preprint2022arXiv

Electric double layers with modulated surface charge density: Exact 2D results

Electric double layers (EDL) with counterions only, say electrons with the elementary charge $-e$, in thermal equilibrium at the inverse temperature $β$ are considered. In particular, we study the effect of the surface charge modulation on the particle number density profile and the pressure. The mobile particles are constrained to the surface of a 2D cylinder and immersed in vacuum (no dielectric image charges). An EDL corresponds to the end-circle of the cylinder which carries a (fixed) position-dependent line charge density. The geometries of one single EDL and two EDLs at distance $d$ are considered; the particle density profile is studied for both geometries, the effective interaction of two EDLs is given by the particle pressure on either of the line walls. For any coupling constant $Γ\equiv βe^2 = 2\times {\rm integer}$, there exists a mapping of the 2D one-component Coulomb system onto the 1D interacting anticommuting-field theory defined on a chain of sites. Using specific transformation symmetries of anticommuting variables, the contact value theorem is generalized to the EDL with modulated line charge density. For the free-fermion coupling $Γ=2$ it is shown that, under certain conditions, the matrix of interaction strengths between anticommuting variables diagonalizes itself which permits one to obtain exact formulas for the particle density profile as well as the pressure. The obtained results confirm the previous indications of weak-coupling and Monte Carlo observations that the surface charge inhomogeneity implies an enhancement of the counterion density at the contact with the charged line and a diminution of the pressure between two parallel lines in comparison with the uniformly charged ones (with the same mean charge densities).

preprint2020arXiv

Full nonuniversality of the symmetric 16-vertex model on the square lattice

We consider the symmetric two-state 16-vertex model on the square lattice whose vertex weights are invariant under any permutation of adjacent edge states. The vertex-weight parameters are restricted to a critical manifold which is self-dual under the gauge transformation. The critical properties of the model are studied numerically by using the Corner Transfer Matrix Renormalization Group method. Accuracy of the method is tested on two exactly solvable cases: the Ising model and a specific version of the Baxter 8-vertex model in a zero field that belong to different universality classes. Numerical results show that the two exactly solvable cases are connected by a line of critical points with the polarization as the order parameter. There are numerical indications that critical exponents vary continuously along this line in such a way that the weak universality hypothesis is violated.

preprint2019arXiv

Short-distance symmetry of pair correlations in two-dimensional jellium

We consider the two-dimensional one-component plasma (jellium) of mobile pointlike particles with the same charge $e$, interacting pairwisely by the logarithmic Coulomb potential and immersed in a fixed neutralizing background charge density. Particles are in thermal equilibrium at the inverse temperature $β$, the only relevant dimensionless parameter is the coupling constant $Γ\equivβe^2$. In the bulk fluid regime and for any value of the coupling constant $Γ=2\times{\rm integer}$, Šamaj and Percus [J. Stat. Phys. {\bf 80}, 811--824 (1995)] have derived an infinite sequence of sum rules for the coefficients of the short-distance expansion of particle pair correlation function. In the context of the equivalent fractional quantum Hall effect, by using specific methods of quantum geometry Haldane [PRL {\bf 107}, 116801 (2011) and arXiv:1112.0990v2] derived a self-dual relation for the Landau-level guiding-center structure factor. In this paper, we establish the relation between the guiding-center structure factor and the pair correlation function of jellium particles. It is shown that the self-dual formula, which provides an exact relation between the pair correlation function and its Fourier component, comes directly from the short-distance symmetry of the bulk jellium. The short-distant symmetry of pair correlations is extended to the semi-infinite geometry of a rectilinear plain hard wall with a fixed surface charge density, constraining particles to a half-space. The symmetry is derived for the original jellium model as well as its simplified version with no background charge (charged wall surface with ``counter-ions only''). The obtained results are checked at the exactly solvable free-fermion coupling $Γ=2$.

preprint2016arXiv

Critical properties of the eight-vertex model in a field

The general eight-vertex model on a square lattice is studied numerically by using the Corner Transfer Matrix Renormalization Group method. The method is tested on the symmetric (zero-field) version of the model, the obtained dependence of critical exponents on model's parameters is in agreement with Baxter's exact solution and weak universality is verified with a high accuracy. It was suggested longtime ago that the symmetric eight-vertex model is a special exceptional case and in the presence of external fields the eight-vertex model falls into the Ising universality class. We confirm numerically this conjecture in a subspace of vertex weights, except for two specific combinations of vertical and horizontal fields for which the system still exhibits weak universality.

preprint2016arXiv

Reentrant disorder-disorder transitions in generalized multicomponent Widom-Rowlinson models

In the lattice version of the multicomponent Widom-Rowlinson (WR) model, each site can be either empty or singly occupied by one of $M$ different particles, all species having the same fugacity $z$. The only nonzero interaction potential is a nearest-neighbor hard-core exclusion between unlike particles. For $M<M_0$ with some minimum $M_0$ dependent on the lattice structure, as $z$ increases from 0 to $\infty$ there is a direct transition from the disordered (gas) phase to a demixed (liquid) phase with one majority component at $z>z_d(M)$. If $M\ge M_0$, there is an intermediate ordered "crystal phase" (composed of two nonequivalent even and odd sublattices) for $z$ lying between $z_c(M)$ and $z_d(M)$ which is driven by entropy. We generalize the multicomponent WR model by replacing the hard-core exclusion between unlike particles by more realistic large (but finite) repulsion. The model is solved exactly on the Bethe lattice with an arbitrary coordination number. The numerical calculations, based on the corner transfer matrix renormalization group, are performed for the two-dimensional square lattice. The results for $M=4$ indicate that the second-order phase transitions from the disordered gas to the demixed phase become of first order, for an arbitrarily large finite repulsion. The results for $M\ge M_0$ show that, as the repulsion weakens, the region of crystal phase diminishes itself. For weak enough repulsions, the direct transition between the crystal and demixed phases changes into a separate pair of crystal-gas and gas-demixed transitions; this is an example of a disorder-disorder reentrant transition via an ordered crystal phase. If the repulsion between unlike species is too weak, the crystal phase disappears from the phase diagram. It is shown that the generalized WR model belongs to the Ising universality class.

preprint2016arXiv

Separation of Wigner structures for 2D equimolar binary mixtures of Coulomb particles

We study the lowest energy configurations of an equimolar binary mixture of classical pointlike particles with charges $Q_1$ and $Q_2$, such that $q=Q_2/Q_1\in [0,1]$. The particles interact pairwisely via 3D Coulomb potential and are confined to a 2D plane with a homogeneous neutralizing background charge density. In a recent paper by M. Antlanger and G. Kahl [Cond. Mat. Phys. {\bf 16}, 43501 (2013)], using numerical computations based on evolutionary algorithm, six fully mixed structures were identified for $0\le q\lesssim 0.59$, while the separation of $Q_1$ and $Q_2$ pure hexagonal phases minimizes the energy for $0.59\lesssim q<1$. Here, we introduce a novel structure which consists in the separation of two phases, the pure hexagonal one formed by a fraction of particles with the larger charge $Q_1$ and the other fixed one containing different numbers of $Q_1$ and $Q_2$ charges. Using an analytic method based on an expansion of the interaction energy in Misra functions we show that this novel structure provides the lowest energy in two intervals of $q$ values, $0<q\lesssim 0.04707$ and $0.58895\lesssim q\lesssim 0.61367$. This fact might inspire numerical methods, for both Coulomb and Yukawa interactions, to test more general separations which go beyond the separation of two pure phases.

preprint2015arXiv

Counter-Ions Near a Charged Wall: Exact Results for Disc and Planar Geometries

Macromolecules, when immersed in a polar solvent like water, become charged by a fixed surface charge density which is compensated by ``counter-ions'' moving out of the surface. Such classical particle systems exhibit poor screening properties at any temperature and the trivial bulk regime (far away from the charged surface) with no particles, so the validity of standard Coulomb sum rules is questionable. In the present paper, we concentrate on the two-dimensional version of the model with the logarithmic interaction potential. We go from the finite disc to the semi-infinite planar geometry. The system is exactly solvable for two values of the coupling constant $Γ$: in the Poisson-Boltzmann mean-field limit $Γ\to 0$ and at the free-fermion point $Γ=2$. We show that the finite-size expansion of the free energy does not contain universal term as is usual for Coulomb fluids. For the coupling constant being an arbitrary positive even integer, using an anticommuting representation of the partition function and many-body densities we derive a sequence of sum rules. As a result, the contact density of counter-ions at the wall is available for the disc. The amplitude function, which characterizes the asymptotic inverse-power law behavior of the two-body density along the wall, is found to be related to the particle density profile. The dielectric susceptibility tensor, calculated exactly for an arbitrary coupling and the particle number, exhibits the anticipated disc value in the thermodynamic limit, in spite of zero contribution from the bulk region. Some of the results obtained in the Poisson-Boltzmann limit are generalized to an arbitrary Euclidean dimension.

preprint2014arXiv

Taking one charge off a two-dimensional Wigner crystal

A planar array of identical charges at vanishing temperature forms a Wigner crystal with hexagonal symmetry. We take off one (reference) charge in a perpendicular direction, hold it fixed, and search for the ground state of the whole system. The planar projection of the reference charge should then evolve from a six-fold coordination (center of a hexagon) for small distances to a three-fold arrangement (center of a triangle), at large distances $d$ from the plane. The aim of this paper is to describe the corresponding non-trivial lattice transformation. For that purpose, two numerical methods (direct energy minimization and Monte Carlo simulations), together with an analytical treatment, are presented. Our results indicate that the $d=0$ and $d\to\infty$ limiting cases extend for finite values of $d$ from the respective starting points into two sequences of stable states, with intersecting energies at some value $d_t$; beyond this value the branches continue as metastable states.