One of the tasks of quantum chemistry is determining the energy levels of an electron in a potential field, which is a superposition of centrally symmetric potentials with centers at the positions of atoms in a molecule. Solutions to quantum mechanical problems of this type are extremely difficult. Reliable results of ab initio calculations are only available for clusters with relatively small numbers of constituents. For example, the structure - and energy calculations of carbon cluster anions CN- are reported for only up to N = 18 in [1,2]. There exists, however, such a field, for which the problem of electron motion has a simple analytical solution. This is the field formed by zero-range potentials. Within the approach described in [3,4], the molecular constituents are replaced by a set of N zero-range potential wells located at the positions Rj (j =1, 2, 3…N) of the nuclei. The goal of the present communication is a demonstration that the model of zero-range potentials [3,4] can be successfully employed for the description of attached electrons in atomic and molecular anions, for example, negatively charged carbon clusters with an arbitrary number of constituents arranged in an arbitrary configuration. We briefly outline the model in the next paragraph.
The wave function of the electron
in a multicenter field, being a solution of the Schrödinger equation everywhere in space except the points Rj and decreasing at r → ∞, is a superposition of the wave functions of single carbon atoms
Where κ is the magnitude of the wave vector of the cluster. At the positions r = Rj, function (1) has to satisfy the boundary conditions
Where
,
is the binding energy of the electron in a single zero-range potential well. We use atomic units throughout the paper. Application of the boundary conditions (2) to the wave function (1) results in a system of homogeneous linear equations for the coefficients cj in equation (1)
The system of equations (3) has nontrivial solutions when its determinant, being a function of κ,
W(κ) = 0. (4)
Equation (4) is a transcendent equation connecting the wave vector κ with the wave vector α. It has several solutions that correspond to discrete energy levels (κ>0). The number of energy levels does not exceed the number of potential wells N. Other states of the negative ion belong to the continuum. The energies of the bound states of the molecular electron
.
To illustrate the capability of the above formulas, let us consider a family of carbon clusters consisting of two-, three- (equilateral triangle) and four (tetrahedron) carbon atoms equidistant from each other. These model carbon structures exhaust all possible configurations, in which all center-to-center distances are the same and equal to R. In terms of the notation
, the determinant W of the system of two equations reads
For the system of 3 equations, the determinant is
The system of 4 equations leads to the determinant
Evaluating the determinants (5) - (7), we use the inter-atomic distance R = 2.751 atomic units and
a.u. (here 1.27 eV is the electron affinity of С-). The graphs of the determinants versus the wave numbers κ are presented in Figure 1. The values of κ at which the curves Wi(κ) cross the X axis are collected in Table 1 along with the corresponding energy levels EN.
Figure 1: Determinants of the system of equations (3) vs the wave vector κ for carbon molecules C
2, C
3, C
4.
Table 1: Zeros of the determinant of the system of equations (3). |
Determinant |
Zeros of determinants, a.u. |
Energies
, eV |
W2(κ) |
Κ= 0.42 |
-2.40 eV |
vW3(κ) |
κ = 0.49 |
-3.26 eV |
W4(κ) |
Κ = 0.55 |
-4.11 eV |
W3(κ) chain |
Κ = 0.26 |
-0.919 eV |
W3(κ) chain |
Κ = 0.46 |
-2.878 eV |
Along with the equilateral triangular configuration, one can consider a linear chain of three carbon atoms with a distance of R/2 between the closest centers. The determinant corresponding to this configuration of atoms has the form [5]
Where A and B are the same as for the previous determinants, C = 2e-κR/2/R. Zeroes of the determinant (8) are shown in two lower lines in Table 1 and in Figure 1. New energy levels at κ = 0.26 a.u. and at κ =0.46 a.u. appear for the chain structure. As we can see, the number of energy levels and their positions are defined by the geometry of the molecule constituents.
Conclusion
In atomic physics, the simplest system that can be modeled by a particle in the field of the zero-range wells is a negative ion of an atom with an external weakly bounded electron in the s-state. Most of the time the outer electron spends outside the potential well, where it can be considered free. The zero-range potential approach opens perspectives for a study of the electron structure of more complex carbon structures, such as fullerenes and their derivatives. The wave functions and energies of the attached electron were calculated earlier in [6] for single-cage (CN)- and multi-cage (onion-like) (CM@CN@...)- singly charged fullerene anions using the model of the spherical potential well. This approach is only applicable for spherically symmetric structures. The model described here can be used for carbon structures of arbitrary geometries. The input information for such calculations includes the coordinates of carbon atoms (for example in C20 [7]) and the binding energy of a single carbon anion. The determinant of the system of equations (3) can be calculated numerically using standard software packages such as MATLAB [8].