Current Generation in Double-Matrix Structure: A Theoretical Simulation

Peculiarities of kinetic characteristics in a supramolecular system, in particular, in a double-matrix structure observed at change of the guest content in a matrix-host are investigated. Results obtained within the framework of a time-independent one-dimensional Schrodinger equation with three barrier potential qualitatively explain experimental data. They indicate the importance of size quantization of a system, correlation between energy and geometric characteristics of both guest and host in this case.


Introduction
Recently, there are intensive both theoretical and experimental research in the field of new technologynanotechnology. In contrast to traditional technologies the nanotechnology operates with materials and systems whose size at least in one crystallographic direction is from one to a hundred nanometers. These materials possess a number of unique physical properties, which are quite perspective from the viewpoint of their practical applications in electronics.
One up-to-date technical problems is the creation of new efficient autonomous generate/conversion devices and energy storage. Today transition to nanoscale objects as energy or charge storages allows to solve this problem. Really, application of the nanodispersed FeS2 in an energy storage device with a lithium anode increases the specific capacity almost by 20% in comparison with a coarsegrained homologue [1], and the nanosized -Fe2O3 possesses the high recirculated capacity 200 mA hour/g and good cycling in a range 1,5-4,0 V regarding Li + /Li in comparison with the macrostructured α -Fe2O3, γ -Fe2O3 and α -Fe3O4 [2]. The concept of ʺelectrochemical graftingʺ [3,4] explains the effect by the features density of nanoobjects. From our point of view the nanostructural cathode material with a hierarchical architecture may be an efficient energy storage. The supramolecular and double-matrix structures belong to a broad class of nanoobjects. The supramolecular structure is structure in which two or more smaller components (guests) are introduced into larger one (host) without the formation of ordinary chemical compounds. It is possible in the case of a weak host -guest interaction. Clathrates [5], in particular relate to such structures, namely their ʺhost-guestʺ interaction is based on the principle of a molecular recognition ʺlock -keyʺ [6][7]. The nanostructural organization -the clathrate (or supramolecular) ensures not only vital new interfacial charge transfer, but also the efficient Faraday storage of energy.
The current-forming phenomenon there was investigated in the molecular lattice mesoporous regular structure on the SiO2 -МСМ-41-basis [8] in the hierarchical doublematrix structure МСМ-41<ТіО2> and in the supramolecular assembly of hierarchical architecture МСМ-41<hydroquinone>>.
The essence of a current-generation reaction in the nanoscale silica (similar in the structure TіO2) is the reaction The kinetics of an intercalation current generation was investigated by the impedance spectroscopy in the frequency range of 10 -3 ÷10 6 Hz. To establish the behavior of kinetic parameters in the current generation of synthesized double-matrix structures the Nyquist diagrams were applied, i.e., dependence of the imaginary part of total impedance on its real part was analyzed.
Obtained results are well simulated by the equivalent electrical circuit shown in Figure 1. Here Rs is the series equivalent resistance, which contains a resistance of electrolyte; the series electric circuits R1||C1|, R2||C2 and the Randles -Ershler circuit C3|| (R3-W1) are connected with the charge transfer through the matrix material, nematic liquid crystal, and matrix/nematic interface (molecular interlayer), respectively. Checking the adequacy of the model package of experimental data showed good results: the Kramers-Kronig coefficient did not exceed 3•10 -5 .
The computer-parametric identification of the diagram is shown in Figure 2. As can be seen, the parameters R2, R3, C2, C3 clearly show the unusual, oscillatory character which explicitly or implicitly (via the thermodynamic Wagner factor [9]) depend on the electron transport. There is reason to believe that such nonmonotonic dependence is caused, in particular, by the rearrangement of electronic states in the supramolecular assemblies due to the intercalant introduction. R3, which is connected with an ion transfer from the electrolyte to the guestʹs positions of cathode material, does not oscillate since it is less sensitive to the peculiarities of electron structure.
It was obtained that dependence of the lithium dissolution entropy on the guest content х in the МСМ-41<hydroquinone <Lix>> is nonmonotonic also.  The purpose of this paper is to show the quantum mechanical nature of the nonmonotonic character of current generation in the studied objects.

Theoretical model
Any solid body is a potential well or a system of potential wells for an electron. Therefore, its spectrum has a discrete character. In the wide potential well these discrete energy levels are so densely arranged that the spectrum can be regarded as continuous. As soon as a body size reduces the distance between the discrete energy levels so increases that the discreteness clearly manifests itself in different physical phenomena (e.g., in an optical spectra). This phenomenon of size quantization is especially important in the low-sized bodies, in particularly, in the nanoparticles or supramolecular systems.
is a Hamilton operator, U(x) is a potential energy of a particle with mass m (E is its total energy;  is the Planck constant).
The well-known quantum-mechanical problem ʺa particle in infinite wellʺ may serve as an illustration of nanoobject (d is width of the well). In this case the quantized levels of a particle are [10,12]: and the corresponding wave functions are: where n 1,2,3,...  From (3) it follows that the discreteness is greater for narrower wells.
Consider the electron tunneling in double-matrix or supramolecular structures. We use an one-dimensional time-independent Schrödinger equation (1) with the potential simulating such structures, namely: This potential is represented in Figure 3.
Here, the ranges 1-2, 5-6 describe the barrier created by the host shell and the range 3-4 with the size d2 is a quest introduced into the host (x2-x1 =c; x3-x2=d1; x4-x3 = d2; x5-x4 = d1; x6-x5 =c). Changing the dimensions 2 d and values of the potential c U it is possible to simulate the amount x of introduced guest. In the wells with finite thickness and height of barriers an electron is not fully localized in the wells, and there is a finite probability of its penetration outside of the wells. As a result, the energy levels have a finite width, and corresponding state is not purely stationary but quasistationary, n n n n1 n2 are the real and imaginary parts of n E , respectively).
Consider two problems: a) Electron energy states within the double-matrix structure -so-called virtual state. (By analogy with the problem of ʺa particle in infinite well ʺ one can predict the electron energy states as a set of levels), and b) an electron tunneling into the structure.
The solution of the 1-st problem is especially important for understanding of the 2-nd problem solution. We would remind you [12] that in the case of two-barrier structure it is probable so-called resonant tunneling. Its essence is increasing (up to 1) of the tunneling probability of an electron through both barriers if its energy coincides with any virtual levels. This situation is a result of the interference of incident and reflected waves in the interbarrier area. The phenomenon of interference should be expected in the double-matrix system. On its basis we can predict the effects of tunneling into the double-matrix structure, which is dominant in the current generation in the nanostructural electric power storages.
Let us analyze the consequences of difference between the potential (4) and the two-barrier potential for tunneling in case of the double-matrix structure.

Virtual states
The solutions of time-independent Schrödinger equation in any range of the rectangular potential are of the form: Let us consider the wave function (6) and its derivative Taking into account the continuity condition of wave functions in any point, matching such functions in the point l x (l = 1, 2, 3,..,6) we obtain: and from the continuity condition of their derivatives in this point we obtain: Recently, the alternative method for finding the energies and wave functions were widely used [13][14]. In [13] it was observed a formal coincidence of the analytical expressions of equations describing the quantummechanical problem with the rectangular potentials and the transmission line problem in classical electrodynamics. The equation can be represented as a recurrence relation in the matrix form:  Table 1 shows the values of ground and first three excited states with a U =0, b U =1 for 2 d = 1 nm and d2 = 8 nm on the guest height с U .   U c U c x, nm x, nm x, nm x, nm Let us compare effect of the guest presence in a well with the conclusions of the quantum-mechanical problem ʺa particle in infinite wellʺ. If in the problem ʺa particle in infinite well ʺ the square of the wave function of ground state has a single maximum in the center of well, in our case this characteristic has a hollow in the same point. Thus, the square of the ground state of wave function is a curve with two symmetrical (as a result of the symmetry of chosen potentials) located maxima. The square of the wave function of first excited state qualitatively agrees with an analogous one in the problem ʺa particle in infinite wellʺ .The only difference is a position of maximum in the first and second cases. Particularly evident changes at the consideration of the same problem take place at the higher values of the potential с U , namely at с U = 0.7 (see Figure 4e). Specific changes occur also with the 2 d width increasing of the barrier с U . At the same time the comparing of the results with the results of the problem ʺparticle in an infinitely wellʺ is less obvious (see Figure 4f).
The increasing of the barrier thickness 2 d simultaneously is accompanied by the narrowing of wells on both sides с U , reducing overlap of the wave functions between them. Therefore, these wells can be regarded as independent. The narrowing of wells itself causes a shift of the levels to the higher values.  Table 2 shows the values of levels depending on the height of the well a U at 2 d = 8 nm, and Figure 5 is a graphical representation of ground state level as a function of the height a U of the well. The absence of values is the fact that some electron states are above the range of our calculations, E 1  еV. The increasing a U leads to the increase of tunneling probability. In extreme case, when the position a U is minimal ( a U = 0) the tunneling grows. Obviously, the tunneling will be higher at lower c a U U  . In other words, a system of the wells can be regarded as the sum of the independent wells at a U = 0 and single well of the host value at a c U U  .
Such circumstances may explain it would seem the strange change   1 a E U . As the depth of wells decreases (with the increasing a U ) the levels in these wells are shifted. At some value of depth upper levels will go beyond observed range of calculation. However, further reduction of depth is accompanied by the increase of tunneling. Thus, the reducing of the depth of wells may be considered as the effective width of guest well that equals from about 1 2d at a U 0  to 1 Table 3 shows the imaginary part of ground state as a function of c U at a thickness of d = 1 nm.

Tunnelling
In the problem of the left to right motion of an electron its tunneling from the range 1 to the ranges 3 and 5 are described by the relations respectively.
The division of equations (8) and (9) The tendency of the dependences shown in Figure 7 is obvious: the tunneling increases when the electron energy tends to the value of barrier heights.
The nonmonotonic dependence of tunneling probability is obtained by the varying of incident electron energy relative to the fixed virtual states. A similar behavior of tunneling is possible with the shift of the virtual states relative the fixed incident electron energy. Such a shift is possible because the changes geometric and/or energy parameters of double-matrix structure. But such shift can be achieved, in particular, by change of the guest content in the structure.

Conclusions
The main result of the paper is the non-monotonic variation of electron tunneling into the host. Since the tunneling determines the current generation phenomenon of double-matrix structure, the result can be regarded as an explanation for the experimentally observed nonmonotonic behavior of current generation characteristics (see Figure 2). We especially note:  the importance of the size quantization in this effect;  the effect to a greater extent is determined by the ratio between the electron virtual states and the bombarding electrons energy than by only their individual parameters.