Finite-Size Effects in Non-Neutral Two-Dimensional Coulomb Fluids
Thermodynamic potential of a neutral two-dimensional (2D) Coulomb fluid, confined to a large domain with a smooth boundary, exhibits at any (inverse) temperature a logarithmic finite-size correction term whose universal prefactor depends only on the Euler number of the domain and the conformal anomaly number . A minimal free boson conformal field theory, which is equivalent to the 2D symmetric two-component plasma of elementary charges at coupling constant , was studied in the past. It was shown that creating a non-neutrality by spreading out a charge at infinity modifies the anomaly number to . Here, we study the effect of non-neutrality on the finite-size expansion of the free energy for another Coulomb fluid, namely the 2D one-component plasma (jellium) composed of identical pointlike -charges in a homogeneous background surface charge density. For the disk geometry of the confining domain we find that the non-neutrality induces the same change of the anomaly number in the finite-size expansion. We derive this result first at the free-fermion coupling and then, by using a mapping of the 2D one-component plasma onto an anticommuting field theory formulated on a chain, for an arbitrary coupling constant.
Keywords:Coulomb fluidFinite-size correctionCentral chargeConformal field theory
According to the principle of conformal invariance Affleck86 ; Blote86 ; Cardy88a , two-dimensional (2D) systems of constituents with short-range interactions, confined to a large domain with a smooth boundary, exhibit at the critical point universal finite-size properties. In particular, for a disk domain of radius , with the Euler number , the (dimensionless) free energy has the large- expansion
where is the inverse temperature and is a length scale. The bulk specific free energy and the surface tension are non-universal, i.e., dependent on microscopic details of the given model and on the temperature. The logarithmic term is universal, depending only on the conformal anomaly number (central charge) of the critical theory. Like for instance, the massless Gaussian theory with a field , defined by the Hamiltonian
has Cardy88b .
In this paper, we study 2D classical (i.e., non-quantum) systems of charged pointlike particles interacting pairwisely by the long-range Coulomb potential , where the dielectric constants of the medium in which particles move and of the walls are equal to (vacuum in Gauss units), for simplicity. The logarithmic potential is the solution of the 2D Poisson equation
and in three dimensions corresponds to the effective interaction of infinitely long parallel charged lines, perpendicular to the plane. Two types of Coulomb systems are of special interest. The one-component plasma (OCP), or jellium, is a system of identical mobile (say elementary) charges immersed in a homogeneous neutralizing background charge density. The symmetric two-component plasma (TCP), or the Coulomb gas, is a system of oppositely charged mobile particles with no background. While the OCP is usually studied in the canonical ensemble, the TCP is treated within the grand-canonical ensemble, with the enforced neutrality condition for each microscopic configuration. In both cases, the thermodynamics and the particle correlation functions are determined by the only dimensionless parameter, the coupling constant . The 2D OCP is exactly solvable at by mapping onto free fermions Jancovici81 . The 2D TCP undergoes the collapse of positive-negative pairs of pointlike charges at coupling , which corresponds to the exactly solvable free-fermion point of the equivalent Thirring model Cornu87 ; Gaudin85 . For reviews about exactly solvable 2D Coulomb systems with various geometries of confining domains, see Refs. Forrester98 ; Jancovici92 .
At arbitrary temperature of the conducting regime, the long-range tail of the Coulomb potential induces screening. Due to the screening effect, the bulk particle correlations exhibit short-range (exponential or even Gaussian) decay which indicates a non-criticality. On the other hand, the same screening phenomenon causes that the induced electrical-field correlations are long-ranged Jancovici95 ; Lebowitz84 . As a result, the free energy (or the grand potential) of any Coulomb system exhibits a universal finite-size correction of type (1.1). Since the logarithmic Coulomb potential is the inverse of the Laplacian operator , in the functional sense it is associated with the Gaussian Hamiltonian (1.2) rewritten as . Although the conformal anomaly number of Coulomb systems is expected to be related to the Gaussian one, there is a change of sign, namely
The explicit checks of the universal finite-size behavior with were done, for both OCP and TCP, at the exactly solvable coupling . The studied cases involve Coulomb gases with periodic boundary conditions Forrester91 , confined to a domain by plain hard walls Jancovici94 , by ideal-conductor boundaries Jancovici96 and by ideal-dielectric walls Jancovici01 ; Tellez01 . For Coulomb systems living on a surface of the sphere, a direct derivation of the universal finite-size correction was performed for any coupling by combining stereographic projection of the sphere onto an infinite plane with linear response theory (TCP, Ref. Jancovici00a ) or with density functional approach (OCP, Ref. Jancovici00c ). The prefactor to the universal correction term was related to the second moment of the short-range part of the planar direct correlation function. Using a renormalized Mayer expansion Deutsch74 ; Friedman62 and observing a cancellation property of specific families of renormalized diagrams, this second moment was evaluated for both TCP Jancovici00b and OCP Kalinay00 . Another explicit derivation of the universal logarithmic term at an arbitrary coupling constant was made for the TCP in a disk geometry Samaj02 . All obtained results confirm the conformal prediction (1.1) for a critical system as we had .
All that has been said holds for neutral Coulomb systems. Recently, Ferrero and Téllez Ferrero14 studied the 2D TCP of charges at the exactly solvable collapse point, confined to a disk of radius . They fixed a “guest” hard-core impurity of charge at the disk origin and found that the large- expansion for the grand potential of type (1.1) applies if the conformal anomaly number is modified to
presence of the impurity has no effect on bulk and boundary terms. It has been noted in the same reference that the same modification of the anomaly number occurs when the minimal free boson conformal field theory Dotsenko04 ; Ginsparg89 , equivalent formally to the 2D TCP, is deformed by spreading out a charge at infinity, creating in this way a non-neutral system. The equivalence of the two systems, one with a charge fixed at the disk origin and one with the same charge spread out at infinity, can be explained by noting that a charge fixed at the disk origin is screened at microscopic distances by counterions from the TCP and thus the universal large- correction looks like the one for a non-neutral system with the net charge . Notice that the definition of in Refs. Dotsenko04 ; Ferrero14 has opposite sign with respect to the standard notation used in this paper. At an arbitrary , conformal field theory Dotsenko04 ; Ginsparg89 yields
With this conformal number, the logarithmic correction in (1.1) is no longer universal, its prefactor depends on both the coupling constant and the excess charge.
In this paper, we concentrate on the non-neutral 2D OCP confined to the disk. We perform a microscopic derivation of the finite-size expansion of the free energy and find that it is of type (1.1) with the same anomaly number (1.6) as found for the 2D TCP in the context of conformal field theories. We would like to note that in contrast to the previous studies Jancovici03 ; Levesque00 ; Martin80 of charge fluctuations in finite Coulomb systems in a much larger overall neutral system or with particle reservoir put at infinity, here we consider non-neutral systems within the constrained canonical ensemble with the fixed net charge and look for the dependence of finite-size corrections of the free energy on this excess charge.
The paper is organized as follows. In Sect. 2, we derive the total energy of mobile charged particles and the fixed background in presence of non-neutrality. Section 3 summarizas the general mapping of the 2D OCP onto the theory of anticommuting variables on a one-dimensional chain. The exactly solvable free-fermion case is the subject of Sect. 4. The free-fermion result is generalized to an arbitrary coupling in Sect. 5. In Sect. 6, we give a short recapitulation and some concluding remarks.
2 Derivation of Boltzmann factor
Let us consider the 2D OCP inside the disk of radius . There are pointlike particles of charge immersed in a uniform background charge density , where the background density is given by
The net charge of the system is thus and the neutral case corresponds to . The particle density is given by
the particle and background densities coincide in the thermodynamic limit .
The potential induced by the background charge density reads as
This potential satisfies the Poisson equation with circular symmetry
The potential at the disk boundary is given by the total background charge placed at the origin, i.e., . The solution of Eq. (2.5), supplemented by this boundary condition, reads
To derive the total Boltzmann factor, we proceed in close analogy with neutral systems, see e.g. Refs. DiFrancesco94 ; Forrester98 ; Sari76 . The total Coulomb interaction energy of the particle-background system consists of three parts:
the particle-particle interaction
the particle-background interaction
the background-background interaction
With the definition of the coupling constant , the corresponding Boltzmann factor reads as
From the expression for the canonical partition function
the free energy is given by , i.e.,
is the configuration integral with the circularly symmetric one-body Boltzmann factor
3 Mapping onto 1D fermions
There exist two basic approaches how to treat the configuration integral (2.13). The first one is based on the expansion of the power of van der monde determinants in the basis of Jack polynomials Salazar16 ; Tellez99 ; Tellez12 . The second one is based on the mapping onto a 1D lattice fermion system; this formalism has been introduced in Ref. Samaj95 and developed further in Refs. Samaj00 ; Samaj04 ; Samaj11 ; Samaj14 ; Samaj15 ; Samaj16 . The relation between the two methods was established in Ref. Grimaldo15 . Here, we apply the mapping onto 1D fermions.
For a positive integer, the configuration integral (2.13) can be expressed in terms of two sets of anticommuting variables each with components , defined on a discrete chain of sites , as follows
Here, and the action involves pair interactions of composite operators
i.e. the products of anticommuting variables from a given set with the fixed sum of site indices. The elements of the diagonal interaction matrix are given by
For the present one-body Boltzmann factor (2.14), the interaction strengths can be rewritten as follows
The one-body density of particles at point is defined by
where is the microscopic density of particles at point and denotes the statistical average over canonical ensemble. The one-body density can be obtained from the configuration integral (2.13) in the standard way as the functional derivative:
Within the formalism of anticommuting variables, the circularly symmetric one-body density is expressible explicitly as
denotes averaging over anticommuting variables.
4 The free-fermion point
Since the correlators
the particle density (3.7) takes the form
According to (3.4), the interaction strengths are expressible as
the expression (2.12) for the free energy can be rewritten as
Now we aim at estimating the crucial sum . Since it holds
we get the important relation
where is the particle number density at , i.e., at the wall contact. We can write
This formula completes the determination of the free energy (4.6).
Let us now consider the large- limit when the background density is related to the particle density via . We start with the sum representation (4.10). With the aid of formulas presented in Appendix A of Ref. DiFrancesco94 , we find that
Using the asymptotic formula Erdelyi
being the error function, we obtain Jancovici94
Finally, the integral
depends on the statistical average which is finite for any , and therefore in the limit contributes to . It is natural to suppose that in the limit the statistical averages do not depend on the finite charge imbalance within an infinite domain and so the particle density at the disk border tends to its half-space “wall contact” value Jancovici82 , independent of . Thus we have
which is the term of order . To summarize,
Inserting this large- formula into (4.6) and taking into account that
we arrive at the finite-size expansion of type (1.1) with
We see that for the coupling the central charge depends on in the same way as for the 2D TCP at the same coupling constant , see Eq. (1.5).
5 Arbitrary with integer
For an arbitrary integer , there exist certain analogies with the exactly solvable case.
In the configuration integral (3.1), we perform the unitary transformation of anticommuting variables, say ’s,
which keeps the composite form of -operators. Here, and are free as-yet unspecified parameters. Taking into account the corresponding Jacobian, this transformation modifies to the form
we eliminate the -dependent prefactor in (3.4) and in this way pass from to . We obtain that
To obtain a convenient representation of , we differentiate this quantity with respect to :
As concerns the correlators , we return back to the original anticommuting ’s variables by applying the inverse to the transformation (5.1). The Jacobians of the numerator and the denominator cancel with one another and we obtain
Consequently, the relation (5.7) takes the form
As the result
As before, the mean particle density at the disk border is finite for any and in the large- limit we have
In contrast to the case, the explicit value of is not known, but the important information is that this term is of order . Consequently,
Although the proof was made for the 2D OCP with a positive integer, it is reasonable to extend its validity to all real coupling values , within the fluid regime.
This work was motivated by the recent paper of Ferrero and Téllez Ferrero14 about the finite-size expansion of the grand potential for the 2D TCP at the coupling , in the presence of the impurity. The modification of the universal logarithmic term to the non-universal one, dependent on the charge of impurity, was there in agreement with the prediction (1.5) of non-neutral conformal field theories.
In this paper, we tested the prediction (1.5) of conformal field theory on the non-neutral 2D OCP. At the free-fermion coupling (Sect. 4), we reproduce the result (1.5). Using the mapping of the 2D OCP onto the anticommuting field theory formulated on the chain, we were able to extend the finite-size analysis to an arbitrary coupling constant . This is one of rare occasions when the exact results are obtained for a series of pair-integer values. The non-universal prediction of conformal field theory (1.6) is confirmed as well.
It would be interesting to generalize our results to the 2D TCP or even to Coulomb fluids with an arbitrary charge composition. For such systems, we miss techniques analogous to that for 2D OCP, but a phenomenological approach might reveal the general form of the finite-size logarithmic term.
This work might be a further motivation for specialists to use Coulomb fluids as practical realizations of conformal field theories to test their predictions.
Acknowledgements.The support received from Grant VEGA No. 2/0015/15 is acknowledged.
- (1) Affleck, I.: Universal term in the free energy at a critical point and the conformal anomaly. Phys. Rev. Lett. 56, 746–748 (1986)
- (2) Blöte, H.W.J., Cardy, J.L., Nightingale, M.P.: Conformal invariance, the central charge, and universal finite-size amplitudes at criticality. Phys. Rev. Lett. 56, 742–745 (1986)
- (3) Cardy, J.L., Peschel, I.: Finite-size dependence of the free energy in two-dimensional critical systems. Nucl. Phys. B 300, 377–392 (1988)
- (4) Cardy, J.L.: Conformal invariance and statistical mechanics. In: Brézin, E, Zinn-Justin, J. (eds.), Fields, Strings and Critical Phenomena, Les Houches 1988, Session XLIX, North-Holland, Amsterdam (1990)
- (5) Cornu, F., Jancovici, B.: On the two-dimensional Coulomb gas. J. Stat. Phys. 49, 33–56 (1987)
- (6) Di Francesco, P., Gaudin, M., Itzykson, C., Lesage, F.: Laughlin’s wave functions, Coulomb gases and expansions of the discriminant. Int. J. Mod. Phys. A 9, 4257–4351 (1994)
- (7) Deutsch, C., Lavaud, M.: Equilibrium properties of a two-dimensional Coulomb gas. Phys. Rev. A 9, 2598–2616 (1974)
- (8) Dotsenko, V.S.: Serie de cours sur la théorie conform. Université de Paris VI-VII (2004)
- (9) Erdélyi, A.: Higher Transcendental Functions. McGraw-Hill, New York (1953)
- (10) Ferrero, A., Téllez, G.: Screening of an electrically charged particle in a two-dimensional two-component plasma at . J. Stat. Mech., P11021 (2014)
- (11) Forrester, P.J.: Finite-size corrections to the free energy of Coulomb systems with a periodic boundary condition. J. Stat. Phys. 63, 491–504 (1991)
- (12) Forrester, P.J.: Exact results for two-dimensional Coulomb systems. Phys. Rep. 301, 235–270 (1998)
- (13) Friedman, H.L.: Ionic solution theory. Interscience, New York (1962)
- (14) Gaudin, M.: Critical isotherm of a lattice plasma. J. Phys. France 46 1027–1042 (1985)
- (15) Ginsparg, P.: Applied conformal field theory. In: Brézin, E, Zinn-Justin, J. (eds.), Fields, Strings and Critical Phenomena, Les Houches 1988, Session XLIX, North-Holland (1990)
- (16) Grimaldo, J.A.M., Telléz, G.: Relations among two methods for computing the partition function of the two-dimensional one-component plasma. J. Stat. Phys. 160, 4–28 (2015)
- (17) Jancovici, B.: Exact results for the two-dimensional one-component plasma. Phys. Rev. Lett. 46, 386–388 (1981)
- (18) Jancovici, B.: Classical Coulomb systems near a plane wall. I. J. Stat. Phys. 28, 43–65 (1982)
- (19) Jancovici, B.: Inhomogeneous two-dimensional plasmas. In: Henderson. D. (ed.) Inhomogeneous Fluids, pp. 201–237, Dekker, New York (1992)
- (20) Jancovici B., Manificat, G., Pisani, C.: Coulomb systems seen as critical systems: Finite-size effects in two dimensions. J. Stat. Phys. 76, 307–329 (1994)
- (21) Jancovici, B.: Classical Coulomb systems: Screening and correlations revisited. J. Stat. Phys. 80, 445–459 (1995)
- (22) Jancovici, B., Telléz, G.: Coulomb systems seen as critical systems: Ideal conductor boundaries. J. Stat. Phys. 82, 609–632 (1996)
- (23) Jancovici, B.: A sum rule for the two-dimensional two-component plasma. J. Stat. Phys. 100, 201–207 (2000)
- (24) Jancovici, B., Kalinay, P., Šamaj, L.: Another derivation of a sum rule for the two-dimensional two-component plasma. Physica A 279, 260–267 (2000)
- (25) Jancovici, B., Trizac, E.: Universal free energy correction for the two-dimensional one-component plasma. Physica A 284, 241–245 (2000)
- (26) Jancovici, B., Šamaj, L.: Coulomb systems with ideal dielectric boundaries: Free fermion point and universality. J. Stat. Phys. 104, 753–775 (2001)
- (27) Jancovici: Charge fluctuations in finite Coulomb systems. J. Stat. Phys. 110, 879–902 (2003)
- (28) Kalinay, P., Markoš, P., Šamaj, L., Travěnec, I.: The sixth-moment sum rule for the pair correlations of the two-dimensional one-component plasma: Exact result. J. Stat. Phys. 98, 639–666 (2000)
- (29) Lebowitz, J.L., Martin, Ph.A.: On potential and field fluctuations in classical charged systems. J. Stat. Phys. 34, 287–311 (1984).
- (30) Levesque, D., Weis, J.-J., Lebowitz, J.L.: Charge fluctuations in the two-dimensional one-component plasma. J. Stat. Phys. 100, 209–222 (2000).
- (31) Martin, Ph.A., Yalcin, T.: The charge fluctuations in classical Coulomb systems. J. Stat. Phys. 22, 435–463 (1980)
- (32) Salazar, R., Telléz, G.: Exact energy computation of the one component plasma on a sphere for even values of the coupling parameter, J. Stat. Phys. 164, 969–999 (2016)
- (33) Šamaj, L., Percus, J.K.: A functional relation among the pair correlations of the two-dimensional one-component plasma. J. Stat. Phys. 80, 811–824 (1995)
- (34) Šamaj, L.: Microscopic calculation of the dielectric susceptibility tensor for Coulomb fluids. J. Stat. Phys. 100, 949–967 (2000)
- (35) Šamaj, L., Jancovici, B.: Density correlations in the two-dimensional Coulomb gas. J. Stat. Phys. 106, 323–355 (2002)
- (36) Šamaj, L.: Is the two-dimensional one-component plasma exactly solvable? J. Stat. Phys. 117, 131–158 (2004)
- (37) Šamaj, L., Trizac, E.: Counter-ions at charged walls: Two-dimensional systems. Eur. Phys. J. E 34, 20 (2011)
- (38) Šamaj, L., Trizac, E.: Counter-ions between or at asymmetrically charged walls: 2D free-fermion point. J. Stat. Phys. 156, 932–947 (2014)
- (39) Šamaj, L.: Counter-ions near a charged wall: Exact results for disc and planar geometries. J. Stat. Phys. 161, 227–249 (2015)
- (40) Šamaj, L.: Amplitude function of asymptotic correlations along charged wall in Coulomb fluids. J. Stat. Phys. 164, 304–320 (2016)
- (41) Sari, R.R., Merlini, D., Calinon, R.: On the ground state of the one-component classical plasma. J. Phys. A: Gen. Phys. 9, 1539–1551 (1976)
- (42) Telléz, G., Forrester, P.J.: Exact finite-size study of the 2D OCP at and . J. Stat. Phys. 97, 489–521 (1999)
- (43) Telléz, G.: Two-dimensional Coulomb systems in a disk with ideal dielectric boundaries. J. Stat. Phys. 104, 945–970 (2001)
- (44) Telléz, G., Forrester, P.J.: Expanded Vandermonde powers and sum rules for the two-dimensional one-component plasma. J. Stat. Phys. 148, 824–855 (2012)