# Diagnostics of many-particle electronic states from non-stationary currents and residual charge

## Abstract

We propose the method for identifying many particle electronic states in the system of coupled quantum dots (impurities) with Coulomb correlations. We demonstrate that different electronic states can be distinguished by the complex analysis of localized charge dynamics and non-stationary characteristics. We show that localized charge time evolution strongly depends on the properties of initial state and analyze different time scales in charge kinetics for initially prepared singlet and triplet states. We reveal the conditions for existence of charge trapping effects governed by the selection rules for electron transitions between the states with different occupation numbers.

###### pacs:

## I Introduction

The control and diagnostics of electronic states in semiconductor nanostructures attracts a great deal of attention now a days. One of the key problems in this area is a development of efficient methods of detection of electronic states with different spin orientation as spin degrees of freedom are considered to play an important role in realizing new functions in modern nanoelectronic devices such as spin pumps (1),(2),(3),(4),(5) and turnstiles (6), (7), spin interference devices (8), quantum dot spin cellular automata (9),(10),(11) and devices for the qubit information (12),(13).

Double QDs are recently an attractive objects for spin-dependent transport analysis (14),(15),(16),(17). Electronic transport through the coupled QDs was considered both in the case of coupling to spin-polarized magnetic (15) and non-magnetic (17) leads. Coupled QDs can be applied for modern nanoelectronic devices creation due to the particular properties of charge and spin kinetics of individual localized states (18),(19),(20),(21). The possibility of QDs integration in a small size quantum circuits deals with careful analysis of relaxation processes and non-stationary effects influence on the electron transport through the dots system (22),(23),(24),(25),(26). Electronic transport in such systems is strongly governed by the presence of Coulomb correlations and by the ratio between the QDs coupling and interaction with the reservoir (27). Correct interpretation of quantum effects in nanoscale systems provides an opportunity to use them as a basis for high speed electronic and logic devices creation (28). Consequently, the problem of charge kinetics in correlated low-dimensional systems due to the coupling with reservoir is really vital. Moreover, non-stationary characteristics provide more information about the properties of nanoscale systems comparing to the stationary one.

In the present paper we propose the way of different many particle electronic states characterization in the system of two interacting quantum dots (impurity atoms) with Coulomb correlations by means of non-stationary current analysis and investigation of the charge trapping effects. Initial charge time evolution is analyzed in terms of pseudo particle technique with additional constraint on possible states (29),(30),(31),(32).

## Ii Theoretical model

We consider the problem of different many particle electronic states resolution in the system of two interacting correlated quantum dots (impurities) by analyzing its non-stationary characteristics. Coupled quantum dots are weakly connected to the substrate, so there is no charge transfer from the dots to the substrate. Charge transfer in the system is allowed only between the dots and due to the QDs coupling to the reservoir, switched to the system at (see Fig.1).

The Hamiltonian of the system

(1) |

is written as a sum of the QDs Hamiltonian

electronic reservoir Hamiltonian

(3) |

and the tunneling part, which describes transitions between the dots and reservoir

Here index labels continuous spectrum states in the reservoir. Localized charge time evolution depends on the way of coupling to reservoir as it was shown in (21). We’ll consider the symmetric coupling to reservoir and assume hopping amplitudes between the reservoir and QD with the energy to be independent on the momentum and spin, so further . Tunneling transfer amplitude between the dots is also considered to be independent on the momentum and spin. Operators are the creation/annihilation operators for the electrons in the continuous spectrum states . -localized state electron occupation numbers, where operator destroys electron with the spin on the single particle energy level . is the on-site Coulomb repulsion for the double occupation of the localized state. We consider low temperature regime when Fermi level is well defined and temperature is much lower than all the typical energy scales in the system. Let us further consider and elsewhere. As we consider the strong coupling between the QDs, the basis of exact eigenfunctions and eigenvalues of coupled QDs without interaction with the reservoir should be applied. Wave functions for the single- and multi-electronic states are well known:

Two single electron states with the wave function

(5) |

and energies

(6) |

exist in the system. Coefficients and are determined by the eigenvector of matrix:

(7) |

Six two electronic states exist in the system: two states with the same spin direction ; and four states with the opposite spins and wave function:

Two electron energies and coefficients , , and are determined by the eigenvalues and eigenvectors of matrix:

(9) |

These are low energy singlet and triplet states and excited singlet and triplet states and . Low energy triplet state with energy exists for any values of QDs energy levels and Coulomb interaction . Corresponding coefficients in Eq.(II) are and .

Two three electron states with the wave function

(10) |

exist in the system. Coefficients and and energies are determined by the eigenvectors and eigenvalues of matrix:

(11) |

Finally, single four-electronic state exists in the system with the wave function

(12) |

Coupling to reservoir leads to the changing of electrons number in the dots due to the tunneling processes. Kinetic properties are governed by the selection rules, which are determined by the matrix elements between the states with different number of electrons. Transitions between the states with different number of electrons can be analyzed in terms of pseudo-particle operators with constraint on the physical states (the number of pseudo-particles). Consequently, the electron operator can be written in terms of pseudo-particle operators:

(13) | |||||

where and - are pseudo-fermion creation (annihilation) operators for the electronic states with one and three electrons correspondingly. , and - are slave boson operators, which correspond to the states without any electrons, with two electrons or four electrons. Operators - describe system configuration with two spin up electrons and one spin down electron in the symmetric and asymmetric states.

Matrix elements , , , , and can be defined as:

(14) |

Finally one can easily express matrix elements through the matrixes (7), (9), (11) eigenvectors elements:

(15) |

For identical QDs and arbitrary values of Coulomb correlations charge trapping can occur due to the presence of one particle dark state with the energy . Matrix element, which corresponds to transitions between this state and empty states is equal to zero. This transition is forbidden by the symmetry of the tunneling Hamiltonian (II). If there is the way of allowed transitions from initial state to the dark state during relaxation processes, the residual charge is trapped in this one-particle state.

Conditions, which determine allowed and restricted transitions between the states with different number of electrons can be easily found:

I. Transitions between the state with zero electrons and single electron state:

(16) |

II. Transitions between single electron state and two electron state with:

(17) |

and

(18) |

III. Transitions between two electron state and three electron state:

(19) |

and

IV. Transitions between three electron state and four electron state:

For slightly different QDs () these matrix elements determine different time scales of the system dynamics. For rather large values of Coulomb interaction and low temperatures only single electron and low energy two-electron states can be considered, as all other states are separated by the Coulomb gap. So, the following non-stationary system of equations can be obtained for the pseudo particle filling numbers , , and :

(22) |

where

(23) |

and , - is the unperturbed density of states in the reservoir. Depending on the tunneling barrier width and height typical tunneling coupling strength can vary from (33) to and meV(34).

System of Eqs.(22) for the single and two-electron states can be solved both numerically and analytically with initial conditions , , and . Analytical expressions, which determine charge relaxation have the following form:

(24) |

where

(25) |

Index corresponds to the state with energy:

(26) |

and index - to the state with energy:

(27) |

Electron occupation numbers can be obtained from the pseudo particle occupation numbers considering spin degrees of freedom by the following expression:

We’ll consider charge time evolution from the singlet and triplet initial states. For singlet initial state coefficients , , and are determined as an eigenvector of matrix (9) corresponding to its minimal eigenvalue.

For the triplet initial state coefficients and .

Non-stationary behavior of the system occupation numbers depends on the initial conditions. Fig.2 demonstrates charge relaxation from initially occupied singlet (see solid lines in Fig.2) and triplet (see dashed lines in Fig.2) states for different system parameters.

When relaxation starts from the singlet state charge trapping effects are not present in the system even for identical QDs ( and ). Charge trapping is present for identical QDS ( and ) when relaxation starts from the triplet state (see black dashed line in the Fig.2). This is the direct manifestation of selections rules influence on the electrons transitions between QDs states and reservoir, which are determined by the matrix elements.

For different QDs, when conditions and are fulfilled, two different timescales for charge relaxation from the states with the energies and exist in the system. The presence of initial energy levels detuning leads to the appearance of two time scales and , related by the ration:

(29) |

All other relaxation rates, present in the system can be expressed through and :

(30) |

Charge relaxation from the triplet state to the single electron states does not depend on the value of Coulomb interaction . Moreover, relaxation rates from the single electron states also do not depend on the value of Coulomb correlations (see Fig.2). For identical QDs time evolution of initial triplet state leads to charge trapping. For slightly different QDs the second relaxation time scale appears in the system. It reveals in the slow charge relaxation instead of charge trapping (see red dashed line in the Fig.2) obtained in the absence of energy levels detuning. For initially singlet state the presence of Coulomb interaction slightly changes relaxation dynamics and do not influence the relaxation processes if initial state is a triplet one.

Obtained significant difference in non-stationary behavior of localized charge for singlet and triplet states gives us possibility to propose experimental scheme, which allows to distinguish different two-electronic states (see Fig.3).

Let us consider the situation, when unknown initial two-electronic state is prepared. There exists four different configurations , , or . To define the particular initial state one has to analyze non-stationary system characteristics and to control the value of residual charge . To measure non-stationary current (see and in the Fig.3) one has to use the system of two consecutive spin-polarized leads (see and in the Fig.3) with opposite directions of applied external magnetic field (see and in the Fig.3).

Proposed experimental scheme allows to distinguish different two-electronic states in the system of two interacting quantum dots with Coulomb correlations by means of the logic diagram shown in Fig.4.

States and can be resolved by means of residual charge control just after non-stationary current pulse has been registered after both leads. If measured charge is equal to zero , then initial state was . In the case of non-zero charge initial state was (see logic diagram in the Fig.4a). Initial states and can be resolved only by the measurements of non-stationary currents. If non-stationary current after both leads is absent but non-stationary current between the leads has non-zero value, then initial state was . Contrary, the absence of current registered between the leads means that initial state was (see logic diagram in the Fig.4b).

## Iii Conclusion

We have proposed the method of how to distinguish different two-particle electronic states in the system of two interacting quantum dots (impurity atoms) with Coulomb correlations. This method is based both on the analysis of the system non-stationary current characteristics and on the control under the residual charge value. The resolution of two-particle electronic states becomes possible due to appearance of significantly different time scales in dynamics of initially singlet and triplet states. We also demonstrated the possibility of charge trapping in the proposed system due to the selection rules presence, which govern electron transitions between different electronic states. The method discussed in the paper can be used for the control under many-particle electronic states in modern nanoelectronic devices such as electronic pumps, turnstiles, logic and quantum information devices and new types of electronic devices based on non-equilibrium non-stationary currents.

This work was supported by RSF grant .

### References

- M.W. Keller, J.M. Martins, N.M. Zimmerman, A.H. Steinbach, Appl. Phys. Lett. 69 , 1804, (1996)
- H. Stoof, Yu. Nazarov, Phys. Rev. B. 53 , 1050, (1996)
- M. Covington, M.W. Keller, R.L. Kautz, J.M. Martins, Phys. Rev. Lett. 84 , 5192, (2000)
- P.I. Arseyev, N.S. Maslova, V.N. Mantsevich, JETP Lett. 95(10), 521, (2012)
- F. Rezoni, T. Brandes, Phys. Rev. Lett. 64 , 245301, (2001)
- J. P. Pekola, J. J. Vartiainen, M. Moðttoðnen, O.-P. Saira, M. Meschke, D. V. Averin, Nature Phys. 4 , 120, (2008)
- D. Averin, J. P. Pekola, Phys. Rev. Lett. 101 , 066801-1, (2008)
- T.-Z. Qian, Z.-B. Su, Phys. Rev. Lett. 72, 2311, (1994)
- A. Bayat, C.E. Creffield, J.H. Jefferson, M. Pepper, S. Bose, Semicond. Sci. Technol. 30, 105025, (2015)
- M.D. Schulman, O.E. Dial, S.P. Harvey, H. Bluhm, V. Umanski, A. Yacobi, Science 336, 202, (2012)
- E.P. Blair, C.S. Lent, J. Appl. Phys. 113, 124302, (2013)
- R. Hanson, J. Burkard, Phys. Rev. Lett. 98, 050502, (2007)
- I. van Weperen, B.D. Armstrong, E.A. Laird, J. Medford, C.M. Marcus, M.P. Hanson, A. Gossard, Phys. Rev. Lett. 107, 030506, (2011)
- V.N. Golovach, D. Loss, Phys. Rev. B, 69, (2004), 245327.
- R. Hornberger, S. Koller, G. Begemann et.al., Phys. Rev. B, 77, (2008), 245313.
- I. Weymann, Phys. Rev. B, 75, (2007), 195339.
- P.I. Arseyev, N.S. Maslova, V.N. Mantsevich, Eur. Phys. J. B 85, 410, (2012)
- C.A. Stafford, N. Wingreen, Phys. Rev. Lett., 76,(1996), 1916.
- B. Hazelzet, M. Wagewijs, H. Stoof et.al., Phys.Rev. B, 63, (2001), 165313.
- E. Cota, R. Aguadado, G. Platero, Phys.Rev Lett., 94, (2005), 107202.
- P.I. Arseyev, N.S. Maslova, V.N. Mantsevich, Eur. Phys. J. B 85, 249, (2012)
- S. J. Angus, A.J. Ferguson, A.S. Dzurak et. al., Nano Lett., 7, (2007), 2051.
- K. Grove-Rasmussen, H. Jorgensen, T. Hayashi et. al., Nano Lett., 8, (2008), 1055.
- S. Moriyama, D. Tsuya, E. Watanabe et. al., Nano Lett., 9, (2009), 2891.
- R. Landauer, Science, 272, (1996), 1914.
- D. Loss, D.P. DiVincenzo, Phys. Rev. A, 57, (1998), 120.
- V.N. Mantsevich, N.S. Maslova, P.I. Arseyev, JETP 118(1), 136, (2014)
- K.Y. Tan, K.W. Chan, M. Mottonen et. al., Nano Lett., 10, (2010), 11.
- P. Coleman, Phys. Rev. B, 29, (1984), 3035.
- P. Coleman, Phys. Rev. B, 35, (1987), 5072.
- N.S. Wingreen, Y. Meir, Phys. Rev. Lett., 49, (1994), 040.
- S.E. Barnes, J. Phys. F: Met. Phys., 6, (1976), 1375.
- S. Amaha, W. Izumida, T. Hatano, S. Teraoka, S. Tarucha, J. A. Gupta, and D. G. Austing, Phys. Rev. Lett., 110, (2013), 016803.
- J. Fransson Phys. Rev. B, 69, (2004), 201304.