Abstract
The partition function for the random walk of an electrostatic field produced by several static parallel infinite charged planes in which the charge distribution could be either is obtained. We find the electrostatic energy of the system and show that it can be analyzed through generalized Dyck paths. The relation between the electrostatic field and generalized Dyck paths allows us to sum overall possible electrostatic field configurations and is used for obtaining the partition function of the system. We illustrate our results with one example.
1. Introduction
In 1961 Lenard considered the problem of a system of infinite charged planes free to move in one direction without any inhibition of free crossing over each other [1]. If all the charged planes carry a surface mass density of magnitude unity and the th one carries an electric surface charge density , the Hamiltonian of Lenard’s problem is given bywhere the system is overall neutral. With the Hamiltonian function at hand Lenard was able to obtain the exact statistical mechanics of the system with the technique of generating functions [1, 2].
From a physical point of view, Lenard’s model can be used in the study of a one-dimensional Coulomb gas [3]. The one-dimensional Coulomb gas is a statistical mechanical problem where particles of equal or opposite charges interact through the Coulomb potential [4]. The model has been extensively studied in the past and forms one of the classical exactly soluble problems in one dimension [5]. Also, the Coulomb gas model has also been used to describe the main features of ionic liquids [6, 7].
More recently, a very similar problem was proposed in [8] with the constraint that each of the charged planes has a fixed position in space and that the surface charge distribution in each plane could be either . This modified model gives rise to a random walk behavior of the electrostatic field that can be analyzed as a Markovian stochastic process.
The purpose of this paper is to give a statistical description of the electrostatic field generated by several static parallel infinite charged planes in which the surface charge distribution could be either . The key step in our formulation consists of the summation of all possible trajectories of the electrostatic field for every different charge configuration. We do this by showing that there is a one to one correspondence between every electrostatic field trajectory and a generalized Dyck path.
The article is organized as follows. In Section 2 we describe the model and derive a system of equations for obtaining the electrostatic energy of the system or Hamiltonian. In Section 3 we make a one to one correspondence between electrostatic field trajectories and generalized Dyck paths. In Section 4 the explicit expression for the partition function is given. In Section 5 a simple numerical analysis is given to show some features of the partition function. In the last section we summarize our conclusions.
2. Model of the System
Suppose that we are given a collection of infinite charged planes, half of them have a constant charge distribution , and the other half have a constant charge distribution ; that is, the system is overall neutral. If we randomly place the infinite charged planes parallel to each other along the -axis at position where , then the electrostatic field would evolve along the -axis making random jumps each time it crosses an infinite charged sheet.
For example, consider a system of four charged planes. We know that equals zero at the left and right side of the configuration due to the neutrality of the system [9]. After crossing the first charged plane, would increase or decrease in an amount of depending on whether the first charged plane had a positive or negative surface charge density. In Figure 1 we show all the possible electrostatic field configurations for the case when the charge density is not known.
Since we are concern with the statistical properties of the electrostatic field we have first to obtain the Hamiltonian of the field or the electrostatic energy of the system. We know from basic electrodynamics courses that the electrostatic energy is given by [9]where is the electrostatic energy density and is the volume inside the region between two consecutive parallel charged sheets. Denoting as the value of the electrostatic field in the th region between two charged sheets placed in and , then we have the following boundary condition:where is a random variable which takes one of the two possible values depending on whether the charged plane at is positively or negatively charged.
We can write down the boundary condition given in (3) in matrix form in the following way:where we have included the initial condition of the electrostatic field; that is, . Solving (4) for the electrostatic field we haveSubstituting (5) into (2) we obtain the electrostatic energy of the system which is given bywhere the Kronecker delta guarantees the neutrality of the system.
3. Electrostatic Field and Dyck Paths
In this section we present the relation between the electrostatic field and generalized Dyck paths. Let us consider for the sake of simplicity that all the charged planes carry a surface charge density of . Let be the total amount of charge to the left of the th region; then we can rewrite the electrostatic field given in (5) asNote that and . The electrostatic energy in terms of the new variable is given then bywhere .
It is convenient to investigate the nature of all the possible electrostatic field configurations by plotting . For the sake of simplicity we will graph versus and connect all the different values of the electrostatic field in each region by straight lines. An example of a graph is shown in Figure 2.
Note that every possible electrostatic field configuration consists of a closed path of length that starts at and ends at and is uniquely determined by a sequence of vertices , where is connected to by an edge [10, 11]. In the mathematical literature this is known as a closed simple walk of steps on the lattice of integers, where the discrete time that numbers the steps is measured in the horizontal axis, and the positions of the walker on the line are recorded in the vertical axis. Since our walks start from the origin, at time the walker is on site . A Dyck path is a simple closed walk with the following constraint for all [12]. The constraint on a Dyck path means that it never falls below the horizontal axis, but it may touch it at any number of integer points. A more general Dyck path or generalized Dyck paths are simple closed walks that are allowed to go below the horizontal axis [13]. We note then that all the electrostatic field configurations are the ensembles of generalized Dyck paths of length that start at and end at .
Let us denote by the ensemble of generalized Dyck paths of length and consider ; then the partition function of the electrostatic field will be given by [14]where , J/K is Boltzmann constant, and is the temperature. Equation (9) is the same as calculating a Feynman amplitude in the frame of generalized Dyck paths; that is, one sums overall generalized Dyck paths of length and to each whole path is then attributed a weight which is given by the Boltzmann factor.
4. Exact Partition Function
In order to calculate the sum of (9) we need to count overall possible generalized Dyck paths of length . A very powerful tool for counting walks is by means of the transfer matrix method [15, 16]. A transfer matrix method can be used to count all possible generalized Dyck paths of length . Let us introduce a set of matrices of the form which represent the number of simple walks in steps. Since the paths we are considering are made of steps , then the transfer matrix for this case is a bidiagonal one given bywhereIn this picture, the evaluation of consists in summing the weights of all walks of steps from to . Therefore, the partition function for the random walk of an electrostatic field generated by means of static parallel randomly charged planes where half of them have a charge surface density and the other half has a charge surface density is given bywhere represents the element in row and column of the transfer matrix raised to the power.
For example, consider the case of 4 parallel static randomly charged planes with surface charge distribution ; for this case and we can then write down the transfer matrix for this particular problem which istherefore, the partition function will be given byEquation (14) is given in the form , where represents the number of paths that share the same energy ; this means that for the case when all the possible electrostatic field configurations will have the same probability distribution; that is,and this was the case studied in [8].
We can calculate the expectation value for the energy of the system when ; that is,Equation (16) means that, for very low temperatures, only the paths with minimum electrostatic energy will contribute to the energy of the system. For this case, the system behaves like parallel plate capacitors connected in series.
5. Numerical Analysis
We can perform a direct evaluation of the partition function given in (12) for a given number of charged sheets using the numerical code given in the appendix. By using the numerical code for , we noticed that the minimum energy of the system is given by and the maximum energy of the system is given by . We have also found the entropic factor for the first two minimum energies and the maximum energy, which invites us to write the partition function asFrom (17) we see that the paths that contribute the most to the partition function are given by the first two terms, and we also see that the second term dominates over the first term ifEquation (18) tells us that when the paths that contribute the most to the partition energy belong to the first excited state.
6. Conclusions
We have shown that it is possible to obtain the exact partition function for the electrostatic field generated by means of several static parallel infinite charged planes in which the surface charge distribution is not explicitly known. We have worked out the special case where the charged planes have a constant surface charge distribution given by and the overall electrostatic system is neutral; that is, there is the same number of positive and negative charged planes. We use the connection between generalized Dyck paths and possible electrostatic field configurations to obtain the partition function of the system as a discrete Feynman path integral. A transfer matrix approach was used to count overall possible generalized Dyck paths. A simple numerical code is given to verify our results.
Appendix
In this appendix we present a simple Mathematica© code which evaluates the partition function given in (12). The only variable that needs to be specified is the number of positive or negative charged sheets .
Conflicts of Interest
The author declares that they have no conflicts of interest.
Acknowledgments
This work was supported by the program “Cátedras CONACYT.”