Equation of state in the fugacity format for the two-dimensional Coulomb gas

We derive the exact general form of the equation of state, in the fugacity format, for the two-dimensional Coulomb gas. Our results are valid in the conducting phase of the Coulomb gas, for temperatures above the Kosterlitz-Thouless transition. The derivation of the equation of state is based on the knowledge of the general form of the short-distance expansion of the correlation functions of the Coulomb gas. We explicitly compute the expansion up to order $O(\zeta^6)$ in the activity $\zeta$. Our results are in very good agreement with Monte Carlo simulations at very low density

Evolution of quantum systems with a scaling type of time-dependent Hamiltonians

We introduce a new class of quantum models with time-dependent Hamiltonians of a special scaling form. By using a couple of time-dependent unitary transformations, the time evolution of these models is expressed in terms of related systems with time-independent Hamiltonians. The mapping of dynamics can be performed in any dimension, for an arbitrary number of interacting particles and for any type of the scaling interaction potential. The exact solvability of a "dual" time-independent Hamiltonian automatically means the exact solvability of the original problem with model time-dependence.Comment: 9 page

Universal behavior of quantum Green's functions

We consider a general one-particle Hamiltonian H = - \Delta_r + u(r) defined in a d-dimensional domain. The object of interest is the time-independent Green function G_z(r,r') = . Recently, in one dimension (1D), the Green's function problem was solved explicitly in inverse form, with diagonal elements of Green's function as prescribed variables. The first aim of this paper is to extract from the 1D inverse solution such information about Green's function which cannot be deduced directly from its definition. Among others, this information involves universal, i.e. u(r)-independent, behavior of Green'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's function for various shapes of the confining domain boundary.Comment: 46 pages, the shortened version submitted to J. Math. Phy

Exactly solvable model of the 2D electrical double layer

We consider equilibrium statistical mechanics of a simplified model for the ideal conductor electrode in an interface contact with a classical semi-infinite electrolyte, modeled by the two-dimensional Coulomb gas of pointlike $\pm$ unit charges in the stability-against-collapse regime of reduced inverse temperatures $0\le \beta<2$. If there is a potential difference between the bulk interior of the electrolyte and the grounded interface, the electrolyte region close to the interface (known as the electrical double layer) carries some nonzero surface charge density. The model is mappable onto an integrable semi-infinite sine-Gordon theory with Dirichlet boundary conditions. The exact form-factor and boundary state information gained from the mapping provide asymptotic forms of the charge and number density profiles of electrolyte particles at large distances from the interface. The result for the asymptotic behavior of the induced electric potential, related to the charge density via the Poisson equation, confirms the validity of the concept of renormalized charge and the corresponding saturation hypothesis. It is documented on the non-perturbative result for the asymptotic density profile at a strictly nonzero $\beta$ that the Debye-H\"uckel $\beta\to 0$ limit is a delicate issue.Comment: 14 page

Expanded Vandermonde powers and sum rules for the two-dimensional one-component plasma

The two-dimensional one-component plasma (2dOCP) is a system of $N$ mobile particles of the same charge $q$ on a surface with a neutralising background. The Boltzmann factor of the 2dOCP at temperature $T$ can be expressed as a Vandermonde determinant to the power $\Gamma=q^{2}/(k_B T)$. Recent advances in the theory of symmetric and anti-symmetric Jack polymonials provide an efficient way to expand this power of the Vandermonde in their monomial basis, allowing the computation of several thermodynamic and structural properties of the 2dOCP for $N$ values up to 14 and $\Gamma$ equal to 4, 6 and 8. In this work, we explore two applications of this formalism to study the moments of the pair correlation function of the 2dOCP on a sphere, and the distribution of radial linear statistics of the 2dOCP in the plane

Ordering and Demixing Transitions in Multicomponent Widom-Rowlinson Models

We use Monte Carlo techniques and analytical methods to study the phase diagram of multicomponent Widom-Rowlinson models on a square lattice: there are M species all with the same fugacity z and a nearest neighbor hard core exclusion between unlike particles. Simulations show that for M between two and six there is a direct transition from the gas phase at z < z_d (M) to a demixed phase consisting mostly of one species at z > z_d (M) while for M \geq 7 there is an intermediate ``crystal phase'' for z lying between z_c(M) and z_d(M). In this phase, which is driven by entropy, particles, independent of species, preferentially occupy one of the sublattices, i.e. spatial symmetry but not particle symmetry is broken. The transition at z_d(M) appears to be first order for M \geq 5 putting it in the Potts model universality class. For large M the transition between the crystalline and demixed phase at z_d(M) can be proven to be first order with z_d(M) \sim M-2 + 1/M + ..., while z_c(M) is argued to behave as \mu_{cr}/M, with \mu_{cr} the value of the fugacity at which the one component hard square lattice gas has a transition, and to be always of the Ising type. Explicit calculations for the Bethe lattice with the coordination number q=4 give results similar to those for the square lattice except that the transition at z_d(M) becomes first order at M>2. This happens for all q, consistent with the model being in the Potts universality class.Comment: 26 pages, 15 postscript figure

"Screening" of universal van der Waals - Casimir terms by Coulomb gases in a fully-finite two-dimensional geometry

This paper is a continuation of a previous one [Jancovici and Samaj, 2004 J. Stat. Mech. P08006] dealing with classical Casimir phenomena in semi-infinite wall geometries. In that paper, using microscopic Coulomb systems, the long-ranged Casimir force due to thermal fluctuations in conducting walls was shown to be screened by the presence of an electrolyte between the walls into some residual short-ranged force. Here, we aim to extend the study of the screening (cancellation) phenomena to universal Casimir terms appearing in the large-size expansions of the grand potentials for microscopic Coulomb systems confined in fully-finite 2D geometries, in particular the disc geometry. Two cases are solved exactly: the high-temperature (Debye-H\"uckel) limit and the Thirring free-fermion point. Similarities and fundamental differences between fully-finite and semi-infinite geometries are pointed out.Comment: 21 pages, 1 figur

Casimir force between two ideal-conductor walls revisited

The high-temperature aspects of the Casimir force between two neutral conducting walls are studied. The mathematical model of "inert" ideal-conductor walls, considered in the original formulations of the Casimir effect, is based on the universal properties of the electromagnetic radiation in the vacuum between the conductors, with zero boundary conditions for the tangential components of the electric field on the walls. This formulation seems to be in agreement with experiments on metallic conductors at room temperature. At high temperatures or large distances, at least, fluctuations of the electric field are present in the bulk and at the surface of a particle system forming the walls, even in the high-density limit: "living" ideal conductors. This makes the enforcement of the inert boundary conditions inadequate. Within a hierarchy of length scales, the high-temperature Casimir force is shown to be entirely determined by the thermal fluctuations in the conducting walls, modelled microscopically by classical Coulomb fluids in the Debye-H\"{u}ckel regime. The semi-classical regime, in the framework of quantum electrodynamics, is studied in the companion letter by P.R.Buenzli and Ph.A.Martin, cond-mat/0506363, Europhys.Lett.72, 42 (2005).Comment: 7 pages.One reference updated. Domain of validity of eq.(11) correcte

Guest charges in an electrolyte: renormalized charge, long- and short-distance behavior of the electric potential and density profile

We complement a recent exact study by L. Samaj on the properties of a guest charge $Q$ immersed in a two-dimensional electrolyte with charges $+1/-1$. In particular, we are interested in the behavior of the density profiles and electric potential created by the charge and the electrolyte, and in the determination of the renormalized charge which is obtained from the long-distance asymptotics of the electric potential. In Samaj's previous work, exact results for arbitrary coulombic coupling $\beta$ were obtained for a system where all the charges are points, provided $\beta Q<2$ and $\beta < 2$. Here, we first focus on the mean field situation which we believe describes correctly the limit $\beta\to 0$ but $\beta Q$ large. In this limit we can study the case when the guest charge is a hard disk and its charge is above the collapse value $\beta Q>2$. We compare our results for the renormalized charge with the exact predictions and we test on a solid ground some conjectures of the previous study. Our study shows that the exact formulas obtained by Samaj for the renormalized charge are not valid for $\beta Q>2$, contrary to a hypothesis put forward by Samaj. We also determine the short-distance asymptotics of the density profiles of the coions and counterions near the guest charge, for arbitrary coulombic coupling. We show that the coion density profile exhibit a change of behavior if the guest charge becomes large enough ($\beta Q\geq 2-\beta$). This is interpreted as a first step of the counterion condensation (for large coulombic coupling), the second step taking place at the usual Manning--Oosawa threshold $\beta Q=2$