Next Article in Journal
Biodegradation of Heterogeneous Industrial Multi-Walled Carbon Nanotubes by Pro-Inflammatory Macrophages
Previous Article in Journal
Unsupervised Learning for the Automatic Counting of Grains in Nanocrystals and Image Segmentation at the Atomic Resolution
Previous Article in Special Issue
A Quantum Mechanical MP2 Study of the Electronic Effect of Nonplanarity on the Carbon Pyramidalization of Fullerene C60
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Model and Energy Bounds for a Two-Dimensional System of Electrons Localized in Concentric Rings

1
Department of Physics, Prairie View A&M University, Prairie View, TX 77446, USA
2
Departament de Física and Institut d’Aplicacions Computacionals de Codi Comunitari (IAC3), University of the Balearic Islands, E-07122 Palma de Mallorca, Spain
3
CRISP—Centre de Recerca Independent de sa Pobla, E-07420 Mallorca, Spain
4
Mathematics Department, Faculty of Science, Sohag University, Sohag 82524, Egypt
5
Deanship of Graduate Studies and Research, Ahlia University, Manama P.O. Box 10878, Bahrain
6
Department of Civil Engineering, Faculty of Engineering, INTI International University, Nilai 71800, Malaysia
7
Department of Mathematics, Faculty of Science and Information Technology, Jadara University, Irbid P.O. Box 733, Jordan
*
Author to whom correspondence should be addressed.
Nanomaterials 2024, 14(20), 1615; https://doi.org/10.3390/nano14201615
Submission received: 15 September 2024 / Revised: 8 October 2024 / Accepted: 8 October 2024 / Published: 10 October 2024

Abstract

:
We study a two-dimensional system of interacting electrons confined in equidistant planar circular rings. The electrons are considered spinless and each of them is localized in one ring. While confined to such ring orbits, each electron interacts with the remaining ones by means of a standard Coulomb interaction potential. The classical version of this two-dimensional quantum model can be viewed as representing a system of electrons orbiting planar equidistant concentric rings where the kinetic energy may be discarded when one is searching for the lowest possible energy. Within this framework, the lowest possible energy of the system is the one that minimizes the total Coulomb interaction energy. This is the equilibrium energy that is numerically determined with high accuracy by using the simulated annealing method. This process allows us to obtain both the equilibrium energy and position configuration for different system sizes. The adopted semi-classical approach allows us to provide reliable approximations for the quantum ground state energy of the corresponding quantum system. The model considered in this work represents an interesting problem for studies of low-dimensional systems, with echoes that resonate with developments in nanoscience and nanomaterials.

1. Introduction

Nanoscience is the science of very small objects and materials that have sizes commensurate with those of atoms and molecules. It is one of the current topics of great interest in all of science, touching on physics, biology, chemistry as well as materials science and engineering. By advancing our understanding of materials at this scale, scientists are gaining a better understanding of how various properties emerge. This knowledge can be used to design better structures with great potential use. Research on nanostructures and nanomaterials has given rise to the development of novel technologies and has enhanced the progress of nanoscience towards unforeseen scientific frontiers. For instance, this progress is clearly seen in the world of microelectronics, where smaller has always meant more components per chip, a faster response, lower cost, lower power consumption, and improved device performance. This trend towards miniaturization and the production of microelectronic chips with dimensions in the nanometer range seems to continue unabated [1,2,3,4,5,6,7].
Within the framework of various nanostructures, the family of two-dimensional (2D) nanomaterials stands out because of the unique properties associated with their lower dimensionality. Typically, such materials are composed of thin layers that may have a thickness of one atomic layer. Differently from their bulk counterparts, these materials have a lot of atoms on their surface. These atoms behave differently from the atoms in bulk. As a result, a 2D material may behave quite differently from its bulk counterpart. One of the must famous examples of this category is graphene. As one of the most important 2D materials, graphene has unique properties that have led to its widespread use in various industries. In particular, the behavior of electrons in 2D graphene samples, as well as the properties of electrons in other low-dimensional structures such as 2D semiconductor systems, has been the object of intensive research for the past few decades. The competition between confinement, spin effects, delocalization, and strong Coulomb repulsion between electrons provides a fertile ground for observation of interesting physical phenomena pertaining to 2D systems of electrons under various conditions [8,9,10,11,12,13,14,15,16,17,18,19,20,21].
It is worth noting that the study of electron correlations in 2D systems has always been one of the main goals of quantum mechanics and condensed matter physics. Multi-pronged efforts in this direction run in parallel with new theories and models in this fast-evolving research area [22,23]. Recently, fabrication of GaAs/AlGaAs double concentric quantum rings [24,25] has enabled, for instance, the study of the Aharonov–Bohm effect [26] as well as the role of Coulomb interaction in magnetic systems [27]. The introduction of minimal models for these intricate quantum systems is a very desirable feature given the complexity of the phenomena that are observed in these structures. Quasi-exact quantum approaches have been successfully applied to different geometries of concentric [28] and parallel [29] quantum rings. However, quantum treatment becomes very difficult when one decides to carry out ab initio calculations and the number of particles is not small. For a situation like this, one may resort to a classical description that may provide very useful insights with less difficulty. For instance, classical studies of 2D systems of interacting electrons confined by a harmonic potential as well as electrons confined to a single circle [30,31,32,33] have shed light on their peculiar properties. For instance, the Thomson problem [34] has been extensively studied from very different perspectives in the literature [35,36,37,38].
In this work, we introduce a model of interacting electrons confined to a 2D planar system. The system consists of an arbitrary number of interacting electrons localized in equidistant circular rings. The electrons are considered spinless and each of them occupies one orbit. The novelty of this model is the fact that the orbits of the electrons are predefined. As a result, no external confinement is required. The ensuing model for the electrons, somehow, mimics that of a 2D classical atom, but without the presence of an attractive central force. The motivation for studying such a model from both a quantum and a classical perspective is two-fold: (i) we provide a brand new (to our knowledge) discrete problem to solve; and (ii) we shed light on the properties of an untractable quantum problem by using classical tools. As we show in this work, the corresponding equilibrium solutions of the 2D classical model under consideration are closely related to the counterpart quantum problem of searching for the ground state energy.
The present contribution is divided as follows: In Section 2, we introduce the classical model to study. In Section 3, we define lower and upper bounds to the ground state energy of the quantum model. The main results for the 2D classical model counterpart as well as related quantum energy bounds are shown in Section 4. An extension of the classical Coulomb model to include Hooke’s law interaction is discussed in Section 5. Finally, some brief conclusions are drawn in Section 6.

2. Classical Model

We consider a system of N 2 electrons with charge, e ( e > 0 ) confined to circular, concentric orbits that are equally spaced. In other words, the first electron is forced to move in a circular ring with radius R 1 = R ; the second electron moves in a circular orbit with radius R 2 = 2 R ; and, similarly, the N-th electron is confined to a ring with radius R N = N R . A schematic view of the model is shown in Figure 1.
The total Coulomb interaction energy of the system is written as
U N = 1 i < j N k e 2 R i 2 + R j 2 2 R i R j cos ( ϕ i ϕ j ) ,
where R i is the radial coordinate of the i-th electron, ϕ i is its azimuthal (polar) angle, and k is Coulomb’s electric constant. The equilibrium of this purely electrostatic system is reached when the total energy reaches a global minimum. Taking into account the fact that orbits have fixed radii ( R i = i R ; i = 1 , 2 , N ), we write a more compact dimensionless object to optimize:
U ˜ N U N k e 2 R = 1 i < j N 1 i 2 + j 2 2 i j cos ( ϕ i ϕ j ) .
For simplicity and for direct comparison of the results of this classical model to the quantum case counterpart, we assume that the radius of the first ring (namely, parameter R) is R = a B , where a B is the Bohr radius. This means that the chosen unit of energy is the atomic unit, k e 2 / a B (a Hartree), which is a commonly used unit in quantum atomic physics. As matter of fact, atomic units (with a B as unit of length and k e 2 / a B as unit of energy) are used for the quantum counterpart case. The complexity of having to deal with the expression in Equation (2) lies in the fact that the electrons are not positioned in the same ring. Nevertheless, there are certain bounds that we can obtain by looking at the functional form of Equation (2) even before embarking on the process of minimizing the total energy function. First, we notice that
i = 1 N 1 j = i + 1 N 1 j + i U ˜ N i = 1 N 1 j = i + 1 N 1 j i .
For the case of two electrons, we have 1 3 U ˜ 2 1 . The equilibrium solution in this case is trivial and coincides with the lower energy bound. The case of three particles reads as 1 3 + 1 4 + 1 5 U ˜ 3 1 + 1 2 + 1 . This result suggests that the optimal solution of Equation (2) should be found within the [ 0.783 , 2.5 ] range. The actual minimum energy value that we found numerically ( 0.8805 ) is quite close to the lower energy bound ( 0.783 ) . Thus, we reasonably expect the latter to be of major relevance as we proceed further with increasing the values of N. The lower energy bound does not seem to correspond to any easily identifiable physical configuration. On the other hand, the upper energy bound represents the maximum possible value of the quantity in Equation (2). This value corresponds to having all electrons aligned in the same direction, similar to the case of one-dimensional (1D) ionic crystals [39,40].
In addition to such observations for small values of N, we also point out that the energy bounds are related to the harmonic series, H M , truncated at some value M. Therefore, let us derive explicit bounds for the optimal energies of Equation (2) via Equation (3) and study the corresponding asymptotic behavior for a large N. It turns out that the lower energy bound can be written in terms of the digamma function, ψ ( N ) = H ( N 1 ) γ , with H ( N 1 ) being the summation of the harmonic series up to N 1 terms and γ the Euler constant. After some algebra, the lower energy bound in Equation (3) is obtained exactly. The upper energy bound is found to be exact as well. The final result reads
N H ( 2 N 1 ) ( N + 1 ) H ( N + 1 ) + 1 4 H ( N + 1 2 ) + 11 60 + ln 2 2 U ˜ N 1 + N [ H ( N 1 ) 1 ] .
At this juncture, let us divide by N all terms in the above inequality, taking into account that H ( N ) = ln N + γ + O ( 1 / N ) and that N 1 . Under these conditions, we have
ln 2 + ln N 4 N < U ˜ N N < ln N .
This means that the energy per particle will be no greater than a logarithmic function, ln N , whatever is the minimum of the classical energy. This is a common result for planar structures in electrostatics studies [41,42]. The significance of Equation (5) lies in the fact that it states the existence of a minimum (and a maximum) for the classical approach to the problem. However, the lower bound is not tight enough. If it was tight enough one might have argued that full optimization of the energy could be avoided, at least as a first crude step. However, as we found out, one must minimize the total energy numerically, in our case by means of a simulated annealing procedure, in order to obtain meaningful accurate results. Therefore, the minimum classical energy for systems with arbitrary values of N is obtained numerically and is denoted as E C m i n = m i n { U ˜ N } . The results obtained for both the classical and quantum model counterparts are discussed in more detail in Section 4.

3. Quantum Model—Ground State Energy

Let us consider the quantum counterpart to the model introduced earlier. The problem setup assumes that the interacting electrons confined in equidistant 2D concentric rings are spinless (fully spin polarized). This is a simplification because the electrons possess two possible spin states. In the absence of interactions, it is evident that the wave function of each electron is confined in its respective ring, and thus, there is no overlap. In reality, the effect of the spin can be negligible in situations that are more realistic than our model as long as the inter-ring distance is large enough to lead to only a small overlap between the wave functions of respective electrons in neighboring rings. The exchange effect, and the effect of tunneling between the rings, can be neglected even in realistic systems when this condition is satisfied. The idea is to prove that there is an energy bound for the quantum case, which is also a bound for the minimum energy of the classical configuration. To this effect, let us look at the quantum mechanical problem for N = 2 electrons situated at circular rings with radii R 1 and R 2 , respectively. The quantum Hamiltonian is written as
H ^ ( R 1 , R 2 ) = 2 2 m e 1 2 + 2 2 + k e 2 | R 1 R 2 | ,
where R i = ( x i , y i ) is the 2D vector position of the i-th electron in the i-th ring, i 2 is the 2D Laplacian operator, m e is an electron’s mass, and is the reduced Planck’s constant.
As stated earlier, we use atomic units in this work. The stationary Schrödinger equation in atomic units reads
1 2 R 1 2 2 ϕ 1 2 Ψ ( ϕ 1 , ϕ 2 ) 1 2 R 2 2 2 ϕ 2 2 Ψ ( ϕ 1 , ϕ 2 ) + 1 d ( ϕ 1 , ϕ 2 ) Ψ ( ϕ 1 , ϕ 2 ) = E Ψ ( ϕ 1 , ϕ 2 ) ,
where ϕ 1 , 2 are the azimuthal (polar) angles and d ( ϕ 1 , ϕ 2 ) = | R 1 R 2 | = R 1 2 + R 2 2 2 R 1 R 2 cos ( ϕ 1 ϕ 2 ) is the separation distance between the pair of electrons localized in two different rings, R 1 and R 2 . One way to solve Equation (7) by preserving the periodicity is to span the unknown solution Ψ ( ϕ 1 , ϕ 2 ) in the basis functions of two non-interacting particles, one in each ring, and then, truncate the expansion to N + 1 terms (N even):
Ψ ( ϕ 1 , ϕ 2 ) = m = N 2 N 2 n = N 2 N 2 c m , n e i m ϕ 1 e i n ϕ 2 2 π R 1 R 2 .
By substituting Equation (8) into Equation (7), multiplying both sides by 1 2 π R 1 R 2 e i k ϕ 1 e i l ϕ 2 , and integrating over { ϕ 1 , ϕ 2 } one obtains
k = N 2 N 2 l = N 2 N 2 [ m 2 2 R 1 2 + n 2 2 R 2 2 δ k , m δ l , n + k l | 1 d | m n   E δ k , m δ l , n ] c k , l = 0 ,
for m , n = N 2 , . . , N 2 . Let us denote as H k l m n the term in the first line in Equation (9). Solving Equation (9) for c k , l is tantamount to providing an approximate solution to Equation (7) for the ground or excited states, with an accuracy that increases as the number of terms in the basis set increases. The matrix element for the Coulomb interaction term in Equation (9) reads explicitly as
k l | 1 d | m n = 1 4 π 2 0 2 π 0 2 π d ϕ 1 d ϕ 2 e i ( m k ) ϕ 1 e i ( n l ) ϕ 2 d ( ϕ 1 , ϕ 2 ) .
The set of Equation (9) for c k , l does not read yet as a standard eigenvalue–eigenstate problem. For this, we must transform H k l m n A i j and c k , l g j , with i , j = 1 , , ( N + 1 ) 2 using i = ( m + N 2 ) ( N + 1 ) + ( n + N 2 ) + 1 and j = ( k + N 2 ) ( N + 1 ) + ( l + N 2 ) + 1 ( k , l , m , n ) . With this transformation, we have the usual eigenvalue–eigenvector problem:
j = 1 ( N + 1 ) 2 A i j E δ i j g j = 0 ,
where i = 1 , 2 , , ( N + 1 ) 2 . The eigenvalues give the energy spectrum of the system. In order to find the eigenvectors, the inverse transformation g j c k , l can be proved to be unique. In other words, given j and N, we can find a sole couple ( k , l ).
Now, since we are looking for a bound to the quantum ground state energy, the kinetic energy terms in Equation (9) can be dropped (since they are always non-negative). The Coulomb energy term is bounded by virtue of the Cauchy–Schwarz inequality for integrals:
k l | 1 d | m n E Q M b o u n d 1 4 π 2 0 2 π 0 2 π d ϕ 1 d ϕ 2 1 d ( ϕ 1 , ϕ 2 ) .
Thus, the eigenvalue problem in Equation (9) gives rise to a total energy which is larger than the exact potential energy, which we denote as E Q M 0 . This is easily understood by recalling that the total energy includes non-negative kinetic energy terms. In other words, for the general N-particle problem, we can write that
E Q M 0 < E Q M b o u n d = 1 i < j N 1 4 π 2 0 2 π 0 2 π d ϕ i d ϕ j 1 d ( ϕ i , ϕ j ) .
The classical equilibrium energy for N electrons in concentric rings, E C m i n , is also bounded by E Q M b o u n d , because the latter is nothing more than the average over all inverses of the inter-particle distances. Since the minimum classical energy configuration is always less than the average, we have E C min < E Q M b o u n d . For the exactly tractable case of two electrons ( R 1 = a B and R 2 = 2 a B ), we have E C min = 1 3 = 0.333 and E Q M b o u n d = 0.536 , while E Q M 0 = 0.483 . The exact proof that the classical minimum energy is smaller than its exact quantum counterpart ( E C min < E Q M 0 ) is a little bit more involved, but it has been extensively checked numerically. The following lemma gives a sense of how the proof is obtained for the case of a one-particle quantum problem.
Lemma 1.
The classical (lowest) equilibrium energy is always a lower bound for the corresponding quantum energy of a system. To this effect, let us assume that we have a quantum system described by the quantum Hamiltonian (written in atomic units):
H ^ = 1 2 Δ + V ( x ) ,
where Δ is a d-dimensional Laplace operator on R d space and V ( x ) is a real-valued potential with x R d . Assume that V ( x ) is bounded below and that a minimum exists:
V m i n = min x R d V ( x ) .
Proof. 
For all normalized wave functions, ψ , one has
ψ , H ^ ψ V m i n .
To be more technical and consistent with the quantum rules, we must say that the wave function, ψ , should be in the domain of H ^ or the quadratic form domain of H ^ (where the energy of the system is finite), and the operator H ^ should be self-adjoint. Here, f , g = f ( x ) ¯ g ( x ) d x is the usual scalar product on the Hilbert space L 2 ( R d ) , where the overline symbol means complex conjugation.
Having reached this stage, we have
ψ , H ^ ψ = ψ , 1 2 Δ ψ     + ψ , V ψ = 1 2 | ψ ( x ) | 2 d x     + V ( x ) | ψ ( x ) | 2 d x     V ( x ) | ψ ( x ) | 2 d x .
To obtain the above results, we used the fact that the kinetic energy term is positive. At this juncture, we point out that | ψ ( x ) | 2 is non-negative. Therefore, for all x R d we have the following lower bound:
V ( x ) | ψ ( x ) | 2 V m i n | ψ ( x ) | 2 ,
which can be easily checked by recalling that V ( x ) V m i n . By using the fact that V m i n is a constant, we have
V ( x ) | ψ ( x ) | 2 d x V m i n | ψ ( x ) | 2 d x     = V m i n | ψ ( x ) | 2 d x = V m i n ,
since we are assuming that ψ is properly normalized. So, for any normalized quantum wave function we have
ψ , H ^ ψ V m i n .
The Rayleigh–Ritz principle now says that the ground state energy of the spectrum is given by
E 0 = min ψ L 2 ( R d ) , ψ 2 = 1 ψ , H ^ ψ V m i n ,
which concludes the proof. □
One can make sense of the main conclusion of this lemma by intuitively appealing to the role played by Heisenberg’s uncertainty principle for such a problem.

4. Results

4.1. Bounds to Classical Energy

Finding the minimum value, E C min , of the quantity in Equation (2) is a very difficult global minimization problem. However, this goal may be achieved numerically by using tools originating from molecular dynamics simulations. In our case, we employ the well-known simulated annealing approach [43]. This robust numerical approach contains a mechanism that allows for numerical searches to eventually escape local minima and thus, enables one to identify the global energy minimum in a reliable manner.
Results for the minimum total classical energies and configurations for systems with N = 3 20 electrons are shown in Table 1 and Figure 2, respectively. One notices that the electrons seem to spread homogeneously in the plane, but, while doing so, they preserve some parallel alignments to the diameters of the concentric circumferences.
Figure 3 depicts the evolution of the classical energy per particle, E C min / N , as a function of N for N = 2 20 electrons.
It is worth mentioning that we observed that the minimum energy per particle, E C min / N , obtained from the annealing simulations is closely followed from below by the energy bound, ln 2 + ln N 4 N , whereas the upper bound is not tight. Thus, it seems that the behavior of the total energy, E C min is to be found between N ln 2 + ln N 4 and N ln N . A stability analysis performed via Δ μ = E N + 1 + E N 1 2 E N does not reveal any special packing or magic number, at least for the system sizes currently considered. We also remark that finding the global energy minimum for this model turned out to be a very challenging simulation problem because of the peculiar constraints imposed in the model. This explains our computational limitations of not being able to go to very high values of N in our simulations.

4.2. Bounds to Quantum Energy

In view of the previous results that we found, it is reasonable to expect that the quantum ground state of N electrons, each in its respective fixed ring, may be well described by a many-electron wave function that has peaks around the equilibrium positions of the classical counterpart system. In order to find the quantum energy bound, E Q M b o u n d , defined in Equation (13), we calculate the mean value of the inverse distance, 1 d ( ϕ i , ϕ j ) , for all the N ( N 1 ) / 2 possible pairs of electrons. The resulting quantum energies per particle, E Q M b o u n d / N , are depicted in Figure 4 as a function of N (upper curve).
The classical equilibrium positions are obtained by the simulated annealing approach. While it is possible to find the global energy minimum for systems with N > 20 electrons, we found that this process comes at a huge computational cost. The constraints imposed in our model make the calculation of the global energy minimum difficult even for the case of N = 3 electrons (the case of N = 2 electrons is the only one where we can identify an analytical result). Since we wanted to display each of the individual equilibrium positions, a sort of Mendeleev table, we stopped at N = 20 , since the tendency is more or less clear. In this way we also avoided the pitfalls of additional costly simulations in which we would have been limited in accuracy by our computational power.
The bound energy curve for the quantum scenario does not correspond to a logarithm function. It is similar to the bound for classic equilibrium energies, namely, ln 2 + ln N 4 N , plus some additive terms. This fact would imply that, for large N the chemical potential for either the classical or the quantum system would be independent of N. In any case, an important result of the present work is the observation that the quantum energy of the system lies between two given precise bounds. This means that we have encountered a quantum system that allows the identification of distinct upper and lower bounds to the ground state energy.

5. Moshinsky-like Systems with Coulomb Interaction

The generalized Moshinsky system (or atom) [44] consists of N harmonically interacting particles confined by an overall external isotropic harmonic potential. Let us now consider a modification of the previous 2D planar ring model. The modification concerns the addition of a harmonic oscillator coupling between all electrons, m e 2 ω 2 ( R i R j ) 2 , where m e is the mass of an electron and ω is the angular frequency. This harmonic oscillator coupling is in addition to the already-present Coulomb interaction term. No external confinement is added since the electrons are kept confined in fixed rings in both the original or presently modified model. A quantum treatment of this modified problem is only partially possible for two electrons [44]. Therefore, let us deal here with lower energy bounds to the quantum case by employing classical tools. For such a situation, the total energy of the system would be written as
U N = 1 i < j N k e 2 d i j + m e 2 ω 2 d i j 2 ,
where d i j = R i 2 + j 2 2 i j cos ( ϕ i ϕ j ) is the usual pair separation distance between pairs of electrons. Let us write the total energy of the system in atomic units:
U ˜ N = 1 i < j N 1 d i j + 1 2 Ω 2 d i j 2 ,
where Ω is some rescaled dimensionless angular frequency. The rest is under the same conditions as before, namely, it is assumed that we have N electrons confined to circular, concentric rings that are equally spaced ( R 1 = R , R 2 = 2 R , . . , R N = N R ), with only one electron per ring, and we take R = a B , where a B is the Bohr radius. Equilibrium is reached when the quantity in Equation (23) has a global minimum.
There are two opposite regimes to consider depending on the strength of Ω . When Ω is small enough, a repulsive regime dominates. When Ω is big enough, all electrons seem to align along the same direction, having a total exact energy:
E N Ω = 1 2 Ω 2 R 2 N 2 ( N 2 1 ) 12 + 1 R 1 + N [ H ( N 1 ) 1 ] .
In order to gauge the behavior of the minimum equilibrium energy as a function of frequency, Ω , we performed numerical simulations using the simulated annealing approach to find the minimum energy corresponding to the expression in Equation (23). We express this quantity as a ratio to E N Ω and study its dependence as a function of Ω . For such a choice, a value close to one indicates that the full harmonic oscillator regime has become dominant over the repulsive Coulomb term. The results of the calculations are depicted in Figure 5 for N = 5 , 10 , and 15 electrons.
Note that as Ω 0 , the energies tend to the electrostatic equilibrium values. Also, notice how sudden the change is. Thus, for the case of this Moshinsky-like 2D classical system, the concomitant equilibrium values are very sensitive to the oscillator frequency. The question of what is any upper energy bound for the total or potential energy in the quantum setting of this case is quite intricate. The argument using the Cauchy–Schwarz inequality may not hold in this scenario. In the classical case, it looks like Equation (24) represents an upper bound (when all particles are aligned along the same direction). However, in the quantum setting, we do not think that there is such an upper limit that can be easily identified.

6. Conclusions

To summarize, in this work we consider a 2D system of interacting electrons confined in equidistant planar circular rings. The electrons are considered spinless. A constraint is set pertaining to the localization of electrons in specific rings by mandating that one ring must be occupied by a single electron. While being confined to such rings, the electrons interact with each other via the standard Coulomb interaction potential. Both the quantum and the classical version of this model are considered. The intricacies of quantum mechanics and the presence of electron correlations does not permit analytic results. Therefore, a key idea of this work is to extract valuable information for the system by looking at energy bounds from both a quantum and a classical perspective. For example, in a quantum scenario, one may discard the quantum kinetic energy term when looking for lower energy bounds. At this juncture, it is worth discussing possible realistic settings where the present model, both the pure repulsive Coulomb interaction and the repulsive–attractive situation (with the added harmonic confinement that induces localization), can be useful. The discussed ring model is definitely useful for modeling a nanowire or a line dislocation since it can be made periodic solely in 1D. Furthermore, a collection of rings could be related to an array of 1D crystals. It is known that dealing with a classical counterpart to a 2D quantum model can be extremely useful [45]. In the present case, this approach leads to a global minimization problem for the energy of the system that can be handled numerically by using the simulated annealing method. By using numerical simulations we were able to accurately obtain the equilibrium position configuration and energy of the system. Then, we compared such classical energy bounds to the quantum case counterpart. The analogy between the classical discrete problem and its quantum counterpart is discussed from the perspective of energy bounds. Along these lines, we obtained analytic expressions for the energy bounds for both classical and quantum models. Furthermore, we also studied the concomitant behavior of the system as a function of the number of electrons. We believe that the correspondence found between equilibrium energies in classical systems and quantum energy bounds may open the door to studying complex quantum many-body systems by means of simpler tools in a way that leads to novel ideas.

Author Contributions

Conceptualization, O.C. and J.B.; Methodology, M.A.-A., M.A.H. and S.A.; Software, J.B.; Validation, O.C., M.A.-A., M.A.H. and S.A.; Formal analysis, O.C. and J.B.; Data curation, J.B.; Writing—review & editing, O.C. and J.B. All authors have read and agreed to the published version of the manuscript.

Funding

The research of O. Ciftja was supported in part by National Science Foundation (NSF) Grant No. DMR-2001980.

Data Availability Statement

The data presented in this study are available upon request from the authors.

Acknowledgments

J. Batle acknowledges fruitful discussions with J. Rosselló, Maria del Mar Batle, and Regina Batle. J. Batle appreciates correspondence with Dirk Hundertmark.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Di Ventra, M.; Evoy, S.; Heflin, J.R., Jr. Introduction to Nanoscale Science and Technology, Series: Nanostructure Science and Technology; Springer: Berlin/Heidelberg, Germany, 2004. [Google Scholar]
  2. Theodore, L. Nanotechnology: Basic Calculations for Engineers and Scientists; Wiley: Hoboken, NJ, USA, 2005. [Google Scholar]
  3. Banyai, L.; Koch, S.W. Semiconductor quantum dots. In Series on Atomic, Molecular and Optical Physics; World Scientific: Singapore, 1993; Volume 2. [Google Scholar]
  4. Jacak, L.; Hawrylak, P.; Wojs, A. Quantum Dots, Nanoscience and Technology Series; Springer: Berlin/Heidelberg, Germany, 1998. [Google Scholar]
  5. Masumoto, Y.; Takagahara, T. Semiconductor Quantum Dots: Physics, Spectroscopy and Applications, Nanoscience and Technology Series; Springer: Berlin/Heidelberg, Germany, 2002. [Google Scholar]
  6. Michler, P. Single Quantum Dots: Fundamentals, Applications and New Concepts, Topics in Applied Physics; Springer: Berlin/Heidelberg, Germany, 2003; Volume 90. [Google Scholar]
  7. Bird, J.P. Electron Transport in Quantum Dots; Kluwer Academic Publishers: Dordrecht, The Netherlands, 2003. [Google Scholar]
  8. Ciftja, O. Quantum Dots: Applications, Synthesis, and Characterization; Nova Science Publishers: New York, NY, USA, 2012. [Google Scholar]
  9. Ciftja, O. Properties of confined small systems of electrons in a parabolic quantum dot. In Quantum Dots: Applications, Synthesis and Characterization; Nova Science Publishers: New York, NY, USA, 2012; pp. 1–12. [Google Scholar]
  10. Ciftja, O. Few-electron semiconductor quantum dots in magnetic field: Theory and methods. In Quantum Dots: Research, Technology and Applications; Nova Science Publishers: New York, NY, USA, 2008; pp. 1–46. [Google Scholar]
  11. Merkt, U.; Huser, J.; Wagner, M. Energy spectra of two electrons in a harmonic quantum dot. Phys. Rev. B 1991, 43, 7320. [Google Scholar] [CrossRef]
  12. Pfannkuche, D.; Gerhardts, R.R.; Maksym, P.A.; Gudmundsson, V. Theory of quantum dot helium. Phys. B 1993, 189, 6. [Google Scholar] [CrossRef]
  13. Ashoori, R.C.; Stormer, H.L.; Weiner, J.S.; Pfeiffer, L.N.; Baldwin, K.W.; West, K.W. N-electron ground state energies of a quantum dot in magnetic field. Phys. Rev. Lett. 1993, 71, 613. [Google Scholar] [CrossRef] [PubMed]
  14. Kumar, A.; Laux, S.E.; Stern, F. Electron states in a GaAs quantum dot in a magnetic field. Phys. Rev. B 1990, 42, 5166. [Google Scholar] [CrossRef]
  15. Pfannkuche, D.; Gudmundsson, V.; Maksym, P.A. Comparison of a Hartree, a Hartree-Fock, and an exact treatment of quantum-dot helium. Phys. Rev. B 1993, 47, 2244. [Google Scholar] [CrossRef] [PubMed]
  16. Fujito, M.; Natori, A.; Yasunaga, H. Many-electron ground states in anisotropic parabolic quantum dots. Phys. Rev. B 1996, 53, 9952. [Google Scholar] [CrossRef] [PubMed]
  17. Müller, H.M.; Koonin, S.E. Phase transitions in quantum dots. Phys. Rev. B 1996, 54, 14532. [Google Scholar] [CrossRef]
  18. Reusch, B.; Grabert, H. Unrestricted Hartree-Fock for quantum dots. Phys. Rev. B 2003, 68, 045309. [Google Scholar] [CrossRef]
  19. Koskinen, M.; Manninen, M.; Reimann, S.M. Hund’s rules and spin density waves in quantum dots. Phys. Rev. Lett. 1997, 79, 1389. [Google Scholar] [CrossRef]
  20. Steffens, O.; Rössler, U.; Suhrke, M. Generalized Hund’s rule in the addition spectrum of a quantum dot. Europhys. Lett. 1998, 42, 529. [Google Scholar] [CrossRef]
  21. Ferconi, M.; Vignale, G. Current-density-functional theory of quantum dots in a magnetic field. Phys. Rev. B 1994, 50, 14722. [Google Scholar] [CrossRef] [PubMed]
  22. Parr, R.G.; Yang, W. Density Functional Theory for Atoms and Molecules; Oxford University Press: Oxford, UK, 1989. [Google Scholar]
  23. Helgaker, T.; Jorgensen, P.; Olsen, J. Molecular Electronic-Structure Theory; John Wiley & Sons: Hoboken, NJ, USA, 2000. [Google Scholar]
  24. Szafran, B.; Peeters, F.M. Few-electron eigenstates of concentric double quantum rings. Phys. Rev. B 2005, 72, 155316. [Google Scholar] [CrossRef]
  25. Mano, T.; Kuroda, T.; Sanguinetti, S.; Ochiai, T.; Tateno, T.; Kim, J.; Noda, T.; Kawabe, M.; Sakoda, K.; Kido, G.; et al. Self-assembly of concentric quantum double rings. Nano Lett. 2005, 5, 425–428. [Google Scholar] [CrossRef] [PubMed]
  26. Aharonov, Y.; Bohm, D. Significance of electromagnetic potentials in the quantum theory. Phys. Rev. 1959, 115, 485. [Google Scholar] [CrossRef]
  27. Mühle, A.; Wegscheider, W.; Haug, R.J. Coupling in concentric double quantum rings. Appl. Phys. Lett. 2007, 91, 133116. [Google Scholar] [CrossRef]
  28. Loos, P.-F.; Gill, P.M.W. Exact wave functions for concentric two-electron systems. Phys. Lett. A 2014, 378, 329. [Google Scholar] [CrossRef]
  29. Guo, G.-J.; Ren, Z.-Z.; Zhou, B.; Guo, X.-Y. Quasi-exactly analytical solutions of two electrons in vertically coupled rings. Int. J. Mod. Phys. B 2012, 26, 1250201. [Google Scholar] [CrossRef]
  30. Yu, E.L.; Mandelshtam, V.A. Classical and quantum melting of a Coulomb cluster in a trap. Phys. Lett. A 1992, 165, 469. [Google Scholar]
  31. Bolton, F.; Rössler, U. Classical model of a Wigner crystal in a quantum dot. Superlatt. Microstruct. 1993, 13, 139. [Google Scholar] [CrossRef]
  32. Bedanov, V.M.; Peeters, F.M. Ordering and phase transitions of charged particles in a classical finite two-dimensional system. Phys. Rev. B 1994, 49, 2667. [Google Scholar] [CrossRef]
  33. Kong, M.; Partoens, B.; Matulis, A.; Peeters, F.M. Structure and spectrum of two-dimensional clusters confined in a hard wall potential. Phys. Rev. E 2004, 69, 036412. [Google Scholar] [CrossRef] [PubMed]
  34. Thomson, J.J. On the structure of the atom: An investigation of the stability and periods of oscillation of a number of corpuscles arranged at equal intervals around the circumference of a circle; with application of the results to the theory of atomic structure. Phil. Mag. 1904, 7, 237. [Google Scholar] [CrossRef]
  35. Mughal, A.; Moore, M.A. Topological defects in the crystalline state of one-component plasmas of nonuniform density. Phys. Rev. E 2007, 76, 011606. [Google Scholar] [CrossRef] [PubMed]
  36. Hardin, D.P.; Saff, E.B. Discretizing manifolds via minimum energy points. Notices Am. Math. Soc. 2004, 51, 1186. [Google Scholar]
  37. Levin, Y.; Arenzon, J.J. Why charges go to the surface: A generalized Thomson problem. Europhys. Lett. 2003, 63, 415. [Google Scholar] [CrossRef]
  38. Mughal, A. Packing of softly repulsive particles in a spherical box - A generalised Thomson problem. Forma 2014, 29, 13. [Google Scholar] [CrossRef]
  39. Ciftja, O.; Rossel, S.; Smith, S.; Thomas, P. Results for the energy of a finite one-dimensional ionic crystal. Results Phys. 2017, 7, 3696. [Google Scholar] [CrossRef]
  40. Ciftja, O. Equivalence of an infinite one-dimensional ionic crystal to a simple electrostatic model. Results Phys. 2019, 13, 102325. [Google Scholar] [CrossRef]
  41. Batle, J.; Vlasiuk, O.; Ciftja, O. Correspondence between electrostatics and contact mechanics with further results in equilibrium charge distributions. Ann. Phys. 2023, 535, 2300269. [Google Scholar] [CrossRef]
  42. Vučić, F. Cylindrical multipole expansion of electrostatic potential in form of triple sum of associated Legendre polynomials. J. Electrostat. 2024, 129, 103932. [Google Scholar] [CrossRef]
  43. Kirkpatrick, S.; Gelatt, C.D., Jr.; Vecchi, M.P. Optimization by simulated annealing. Science 1983, 220, 671. [Google Scholar] [CrossRef] [PubMed]
  44. Moshinsky, M. How good is the Hartree-Fock approximation. Am. J. Phys. 1968, 36, 52. [Google Scholar] [CrossRef]
  45. Batle, J.; Ciftja, O. Zigzag phase transition of electrons confined within a thin annulus region. Phys. Lett. A 2024, 517, 129665. [Google Scholar] [CrossRef]
Figure 1. (Color online) Relative position of the rings, each one containing one electron. See text for details.
Figure 1. (Color online) Relative position of the rings, each one containing one electron. See text for details.
Nanomaterials 14 01615 g001
Figure 2. (Color online) Equilibrium configuration for systems with N = 3 20 electrons. Distances are given in atomic units, where it is assumed that R = a B . See text for details.
Figure 2. (Color online) Equilibrium configuration for systems with N = 3 20 electrons. Distances are given in atomic units, where it is assumed that R = a B . See text for details.
Nanomaterials 14 01615 g002
Figure 3. (Color online) Minimum classical energy per particle E C min / N (upper curve) as a function of N for systems ranging from N = 2 to N = 20 electrons. Its lower bound from Equation (5) is plotted as the lower curve for comparison. Note that the lower curve steadily tends to ln 2 as N increases. The energy is given in atomic units of k e 2 / a B , where R = a B is assumed. See text for details.
Figure 3. (Color online) Minimum classical energy per particle E C min / N (upper curve) as a function of N for systems ranging from N = 2 to N = 20 electrons. Its lower bound from Equation (5) is plotted as the lower curve for comparison. Note that the lower curve steadily tends to ln 2 as N increases. The energy is given in atomic units of k e 2 / a B , where R = a B is assumed. See text for details.
Nanomaterials 14 01615 g003
Figure 4. (Color online) Energy per particle, E Q M b o u n d / N , as a function of N (upper curve). The classical counterpart E C min / N is also shown for comparison (lower curve). The energy is given in atomic units of k e 2 / a B , where R = a B is assumed. See text for details.
Figure 4. (Color online) Energy per particle, E Q M b o u n d / N , as a function of N (upper curve). The classical counterpart E C min / N is also shown for comparison (lower curve). The energy is given in atomic units of k e 2 / a B , where R = a B is assumed. See text for details.
Nanomaterials 14 01615 g004
Figure 5. (Color online) Equilibrium energies as U ˜ N / E N Ω versus Ω for N = 5 (upper curve), N = 10 (middle curve), and N = 15 (lower curve). The change in regime occurs quite abruptly. See text for details.
Figure 5. (Color online) Equilibrium energies as U ˜ N / E N Ω versus Ω for N = 5 (upper curve), N = 10 (middle curve), and N = 15 (lower curve). The change in regime occurs quite abruptly. See text for details.
Nanomaterials 14 01615 g005
Table 1. Minimum total classical energy, E C m i n , for a given number of electrons. The energy is given in atomic units of k e 2 / a B , where R = a B is assumed. See text for details.
Table 1. Minimum total classical energy, E C m i n , for a given number of electrons. The energy is given in atomic units of k e 2 / a B , where R = a B is assumed. See text for details.
N E C min N E C min
30.8805128.16254
41.5266139.06952
52.24421410.0028
63.01761510.9343
73.81001611.8688
84.64701712.8314
95.50581813.7806
106.37441914.7164
117.25432015.7123
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ciftja, O.; Batle, J.; Abdel-Aty, M.; Hafez, M.A.; Alkhazaleh, S. Model and Energy Bounds for a Two-Dimensional System of Electrons Localized in Concentric Rings. Nanomaterials 2024, 14, 1615. https://doi.org/10.3390/nano14201615

AMA Style

Ciftja O, Batle J, Abdel-Aty M, Hafez MA, Alkhazaleh S. Model and Energy Bounds for a Two-Dimensional System of Electrons Localized in Concentric Rings. Nanomaterials. 2024; 14(20):1615. https://doi.org/10.3390/nano14201615

Chicago/Turabian Style

Ciftja, Orion, Josep Batle, Mahmoud Abdel-Aty, Mohammad Ahmed Hafez, and Shawkat Alkhazaleh. 2024. "Model and Energy Bounds for a Two-Dimensional System of Electrons Localized in Concentric Rings" Nanomaterials 14, no. 20: 1615. https://doi.org/10.3390/nano14201615

APA Style

Ciftja, O., Batle, J., Abdel-Aty, M., Hafez, M. A., & Alkhazaleh, S. (2024). Model and Energy Bounds for a Two-Dimensional System of Electrons Localized in Concentric Rings. Nanomaterials, 14(20), 1615. https://doi.org/10.3390/nano14201615

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop