Real arrays can be expected to have two main deviations from the ideal Hubbard model. Fitsr, in addition to the intradot Couloumb repulsion there will exist a certain degree of interdot repulsion. Second, individual dots will be invariably “detuned†from each other to some extent. Both these aspects are readily incorporated in our model, and we find that they have a noticeable effect on the two groups of peaks which we identified as the upper and the lower Hubbard bands. Deturning tends to localize the electron states, thus suppresses the conduction peaks. The effects of deturning have been presented in a separate publication. Consider a one-dimensional (1D) array of N coupled dots, indexed from left to right as 1─N, described by the Hamiltonian H,

In Eq (1.1) above, εka and ε in are energy levels in leads and the ith dot of the array, respectively, with an α being the spin index. Ui is the intradot repulsion of the ith dot, while Wi and ti are the interdot repulsions and the interdot coupling between the ith dot and its right neighbor [the (i+1)th dot]. The tunneling matrix element V L (V R ) connects dot 1 (dot N) to the left (right)
ka ka
lead. We assume two spin-degenerate levels on each dot.
The whole array were treated as a single quantum system and calculate its many-body eigenstates by exact diagonalization. The demand on computing power grows factorialy with the number of states. Arrays containing up ten dots (20 states) were studied. Once the eigenstates are known, the conductance is calculated from the relation.

In Eq (1.20), εn,i is the energy of many-body state (n,i), the ith if the n-particle states, Г L(R) are the transition rates between state (n,i) and (n-1,j) by losing or getting one electron through the left (right) lead, and P eq is the occupation probability of state (n,i) at equilibrium,