Dopant-controlled single-electron pumping through a metallic island

Dopant-controlled single-electron pumping through a metallic island


We investigate a hybrid metallic island / single dopant electron pump based on fully-depleted silicon on insulator technology. Electron transfer between the central metallic island and the leads is controlled by resonant tunneling through single phosphorus dopants in the barriers. Top gates above the barriers are used to control the resonance conditions. Applying radio frequency signals to the gates, non-adiabatic quantized electron pumping is achieved. A simple deterministic model is presented and confirmed by comparing measurements with simulations.


Accurate clocked control and transport of single electrons (SE) is a prerequisite for both the emerging field of electron quantum opticsGrenier et al. (2011) and the upcoming introduction of the quantum ampere by fixing the elementary charge .Mills et al. (2011) Transferring an integer number of electrons  in each cycle of a drive with frequency generates a current . Several approaches to clocked SE transport have been madePekola et al. (2013) of which tunable-barrier SE pumps are especially promising to reach sufficiently large currents at high accuracy for metrological needs.Kaestner and Kashcheyevs (2015) So far, GaAs based devices achieved the highest verified accuracy at large output currents,Stein et al. (2015); Giblin et al. (2012); Bae et al. (2015) however a high magnetic field and sub-kelvin temperature are required for accurate operation.

SE pumps made from silicon are promising candidates for less demanding operation requirements and even higher operation frequencyFujiwara et al. (2008) and have been realized using both adiabaticJehl et al. (2013) and non-adiabaticFujiwara et al. (2008); Rossi et al. (2014) pumping schemes. They benefit from mature CMOS technology with the potential of circuit integration, e.g. on chip drive electronicsClapera et al. (2015) or the charge detectorsYamahata et al. (2011); Tanttu et al. (2015) needed for a self-referenced quantized current source.Fricke et al. (2013, 2014) Additionally, CMOS technology offers excellent control of doping, opening the possibility to utilize the large addition energies of single dopants for single-electron control with a potential of increased accuracy. Up to now, the dopant or trap was used as the main quantum dot of a SE pump.Lansbergen et al. (2012); Roche et al. (2013); Tettamanzi et al. (2014); Yamahata et al. (2014) In this letter a new device scheme is examined: Single dopant states located in the barriers are used to control the tunnel coupling of a small metallic island to source and drain. While in other works similar states in barriers were an unwanted side-effect,Chan et al. (2011) here we actively utilize these states to perform SE pumping.

Figure 1: (a) Schematic cross section along the nanowire axis. Inset: SEM image of a similar device after gate and spacer patterning (scalebar 100 nm). (b) Schematic band structure of the investigated device. The central island is metallic. Under each of the gates, single dopants dominate the transport.

The investigated device is produced in a fully CMOS-compatible process using both deep UV and electron beam lithography. It is based on a silicon on insulator wafer with 145 nm of SiO, where the bulk silicon layer can be used as a back gate.Roche et al. (2012) The 20 nm thick Si-layer is n-doped with phosphorus ( cm). A 60 nm wide nanowire is etched out of the Si-layer. Two top gates (consisting of dielectrics and doped poly Si) are patterned around the nanowire with a length of 60 nm and a separation of 60 nm. The gates and a set of 15 nm wide self-aligned silicon nitride spacers act as a mask for the following doping and silicidation. Consequently, all regions of the nanowire not covered by gate or spacers will be metallic, including a central island formed between the two gates with a width of 30 nm. A schematic cross section of the device along the nanowire axis is shown in Fig. 1(a).

The simplified band structure of the device along the length of the nanowire is as follows [Fig. 1(b)]: Source and drain are metallic. By tuning the top gates to a negative voltage, barriers are formed below the top gates and isolate the metallic region in between the gates from the leads, causing discrete charge states to be formed. Under the gates, the average interdopant distance is approximately  nm, larger than the Bohr radius  nm for P dopants in Si, making these regions semiconducting.Tabe et al. (2010) By choosing the back gate voltage appropriately the device can be tuned such that for each gate only the energy level of a single dopant lies in the transport window and is sufficiently separated from any other state.Voisin et al. (2014); Verduijn et al. (2013) The states in the barrier under a single gate can be investigated individually by tuning the gate voltage of the other gate to  V, well above the threshold value. From Coulomb diamond measurements (not shown) the energy separations of the lowest state to the next state are found to be approximately  meV, comparable to the charging energy commonly found for single phosphorus dopants in Si nanostructures.Tan et al. (2010); Lansbergen et al. (2008) States at higher energies are more densely spaced.

Figure 2: Conductance through the device as a function of top gate voltages. Anti-diagonal lines correspond to resonances of the central island. Horizontal and vertical lines are resonances of the states in the barriers. For  mV and  mV, transport through and from the central island is only possible via dopant states (extended with red dashed lines). Plot has been extended beyond measured values to show the line along which the measurement in Fig. 3(a) has been taken (red line) and RF amplitude used in that measurement (blue line).

Fig. 2 shows the measured conductance as a function of top gate voltages and . All measurements were taken at a base temperature of  mK. The top gates couple to the energy levels of the states in the barriers as well as to the energy levels of the central island. Multiple anti-diagonal lines indicate resonance of the central island with source and drain. They are visible for more positive gate voltages at which the barriers have sufficient transparency. Horizontal (vertical) lines of increased conductance appear, whenever a dopant’s state in the barrier under gate 1 (gate 2) is resonant.Golovach et al. (2011) Whenever the anti-diagonal lines intersect a horizontal or vertical resonance line anti-crossings appear, indicating that the central island is capacitively coupled to the states in both barriers. In the absence of these anti-crossings, resonance lines of the central island are evenly spaced as expected for a metallic dot.Jehl et al. (2013) The negligible slope of the horizontal and vertical resonances of the states under the gates demonstrate very small cross-coupling from gate 1 on states located under gate 2 and vice versa, caused by the strong screening provided by the gates and the central island. For more negative gate voltages, conductance is only appreciable in a small region at  mV and  mV, where the first relevant1 energy level of gate 1 and gate 2 ( and , respectively) are both aligned with the Fermi level , i.e. . These resonances are highlighted with dashed red lines in Fig. 2.

Figure 3: (a) Pumped current as a function of gate voltage along the path indicated in red in Fig. 2. (b) Stability diagram in the region where transport through the barriers is only possible via the lowest dopant states. Resonance lines of the central island (black) and the dopant states (red) are shown. The path followed by RF voltages is shown in blue (corresponds to blue path indicated in Fig. 2). This situation corresponds to the case marked with the red dot in (a), where electron is pumped. (c) Positions of energy levels during one cycle of RF modulation.

We will now demonstrate that using the dopant states in the barriers this new type of device allows to generate a quantized current . Sinusoidal radio frequency (RF) signals with 180 degree phase shift are applied to both gates to follow a certain path in the --plane as indicated by the blue line in Fig. 2 ( MHz for all shown measurements). When shifting the center of the blue path (red dot) along the solid red line, the produced current changes, as shown in Fig. 3(a). Starting at the most negative voltages, the current increases from to in steps of , clearly showing quantization of current. Subsequently, more steps between and follow.

To explain the structure of the current quantization steps we consider a simple deterministic model which treats the discrete charge states on the central island and under the barriers as energetically sharp at all times and allows transport only through the resonant states highlighted with red dashed lines in Fig. 2. The model assumes that the charge of the central island switches if and only if sequential tunneling through a dopant state under the barrier is energetically allowed. Further neglecting the dopant-island capacitive coupling we can reduce the model to the following simple rules. The th electron can only tunnel onto the metallic central island with electro-chemical potential if a gate’s energy level () is in the energy window of the Fermi energy and the dots potential: . Similarly, the th electron can only leave the island, if . The deterministic model is justified for large amplitude and intermediate frequency of modulation, such that on the one hand there is enough time for charge equilibration as the resonance energy passes through the bias window across a dot-lead barrier while on the other hand any off-resonance hopping events remain negligibly rare.

Using these assumptions, we will demonstrate how a quantized current of electron per cycle can be produced; equivalent to the red dot in Fig. 3(a). Figure 3(b) shows the resonances of energy levels with source and drain as a function of top gate voltages. Red dashed lines represent resonances of the states in the barriers (). Black lines represent resonances of the central island’s electrochemical potential (), which are at an angle of to the -axis. For clarity, the electron resonance is highlighted in cyan. By applying the RF signals to the top gates, the path shown in blue is produced (angle to -axis ). The RF signals on the gates are -shifted, therefore no area is enclosed in the --plane. Therefore, this single-parameter pumping process is non-adiabatic.Kaestner and Kashcheyevs (2015)

In Fig. 3(c) positions of the energy levels are sketched for certain points along the path. Starting at frame 1, , allowing the central island to be filled from the source to an occupancy of electrons. Next, frame 2 highlights the non-adiabaticity of the process:Kaestner and Kashcheyevs (2015) Even though , the occupancy of the central island remains unchanged, as no state in any barrier is in the appropriate energy window (). Finally, in frame 3, , allowing one electron to exit the central island to the drain. Following the path back to frame 1, no electrons enter or leave the dot, as none of the requirements for transport are met. During one cycle, electron has been transported from source to drain in this example.

Figure 4: (a) Measurement of pumped current as a function of DC gate voltages. Dashed red lines correspond to resonances of states in the barriers under the gates. Solid red line corresponds to the cross section shown in Fig. 3(a). (b) Simulation of the produced current using the same parameters as for the measurement in (a) [y-axis shared with (a)]. (c) Measurement of pumped current for smaller amplitude and closer to . Reversal of current can be observed. (d) Simulation corresponding to (c).

To change the number of electrons transported per cycle, the top gate DC voltages can be changed. Figure 4(a) shows a pump current measurement as a function of both top gate voltages, applying the same RF signals as before. The path along which the measurement in Fig. 3(a) has been taken is indicated by the solid red line. The situation described in Fig. 3(b) is highlighted by the red dot, where a current of is measured. Positions in gate voltage of the dopant resonances are marked with red dashed lines. The observed pattern features two main structures: Towards more negative , a set of diamonds with quantized current is visible. The first column has a value of , the second column a value of . Going towards more positive , the diamond shape is replaced with a set of anti-diagonal stripes on which the produced current stays constant. Their slope is found to be equal to the angle of the gate voltage path (blue line in Fig. 2).

To compare the measurement with the described model and explain the observed changes in the transferred number of electrons per cycle , a simple simulation has been made counting the number of electrons transferred during one pump cycle following the above assumptions for transport to and from the central island. The following experimentally gained parameters are used for the simulation (lever arms ): The resonance positions of the states in the barriers as function of gate voltages (from Fig. 2), the lever arms of gate voltage to the barrier state’s energy, and , and to the central island’s potential, and , and the charging energy of the central island  meV (all from Coulomb diamond measurements, not shown). As the tunneling rates are assumed to be infinite, electron transfer takes place whenever possible, and the simulation yields only integer values for .

The simulation is shown in Fig. 4(b). Red dashed lines indicate the resonances of the relevant states under each barrier. The produced pattern is almost point symmetric around the crossing point of those resonances, with the sign of the produced current being flipped for one side of the pattern. First, we will describe the region with positive current which fits well with the measurement. The other half will be discussed later. Point A is equivalent to the situation depicted in Figs. 3(b) and 3(c), where electron is transmitted. Moving towards more positive , point B is reached, where no electrons are transported. Here, the ejection of the electron in frame 3 can not take place, as the condition for ejection is not reached. Going from A into the opposite direction towards C the current increases to electrons. Here, an additional electron can tunnel out during frame 3, as . When is decreased further (point D), no more electrons are transported. Here, the time dependent gate voltage path does not cross the resonance any more, so no electrons are loaded into the central island during frame 1. When going from point C to E, an additional electron is loaded onto the central island, as during frame 1, leading to a total current of . Going further towards point F, the current switches between and several times while the diamond shape is lost. In this region, additional electrons tunnel onto the island during frame 1, but the number of electrons emitted does not increase further, as the resonance is crossed before frame 3. According to the model, the th electron can not be emitted if , which gives the anti-diagonal stripes of constant current.

Comparing the simulation in the region where positive current is expected with the measurement, one finds generally good agreement. While the maximum current in the measurement is , it is in the simulation. Further, the positive slope of the diamond shapes as well as the slope of the column of diamonds differ. The maximum current and both slopes depend on the lever arms used for the simulation. We have used the values of the latter measured in DC transport with either one of the barriers removed and the potential of the island fixed or both barriers sufficiently transparent even off dopant resonance, which is different from the conditions for pumping. A difference in screening conditions and the neglected capacitive coupling of the central island to the states in the barriers are possible sources of the discrepancy. Capacitive coupling is also likely to cause the zig-zag pattern at the bottom of the measured current map. These effects are subject to further investigation.

When the rules for electron transport are applied to the whole gate range, the simulation suggests that the pattern should be repeated with negative current in the region of more positive gate voltages. This reversal of current is not observed in the data in Fig. 4(a) due to the presence of additional transport channels that were neglected in the model. These transport channels become more effective at more positive gate voltages. Their effect can be minimized using RF signals with smaller amplitude and closer to , yielding a pump current that shows the expected reversal [Fig. 4(c)], although the produced current is now limited to . A simulation with a corresponding path shows good agreement [Fig. 4(d)]. The influence of additional transport channels is subject to further research.

In summary, we have investigated a hybrid quantum dot system consisting of a metallic island, whose tunnel coupling to source and drain is dominated by single dopant states located in the barriers under the top gates. A non-adiabatic pumping scheme using -shifted sinusoidal RF signals can be implemented to produce a quantized current which is strongly dependent on the DC gate voltages. With basic assumptions about the system the process can be explained and modeled. A predicted reversal of current, though obscured by additional transport channels at large RF amplitude, is observed for small RF amplitudes. This demonstrates the possibility of gate controlled current reversal for metrological applications. Future work will focus on the reduction of additional resonances in the barriers and on improving the model by including additional transport channels, finite rates and capacitive coupling between dopants and central island in order to investigate the dynamics of coupled single dopant-quantum dot systems.

We are grateful to R. Haug for fruitful discussions. We acknowledge financial support through the EC FP7-ICT initiative under Project SiAM No 610637 and from the Nanosciences foundation through the excellence chair program.


  1. The resonance of the lowest state under gate 1 (found at  mV) is very weakly coupled and cannot be resolved in this measurement.


  1. C. Grenier, R. Hervé, G. Fève,  and P. Degiovanni, Mod. Phys. Lett. B 25, 1053 (2011).
  2. I. M. Mills, P. J. Mohr, T. J. Quinn, B. N. Taylor,  and E. R. Williams, Phil. Trans. R. Soc. A 369, 3907 (2011).
  3. J. P. Pekola, O.-P. Saira, V. F. Maisi, A. Kemppinen, M. Möttönen, Y. A. Pashkin,  and D. V. Averin, Rev. Mod. Phys. 85, 1421 (2013).
  4. B. Kaestner and V. Kashcheyevs, Rep. Prog. Phys. 78, 103901 (2015).
  5. F. Stein, D. Drung, L. Fricke, H. Scherer, F. Hohls, C. Leicht, M. Götz, C. Krause, R. Behr, E. Pesel, K. Pierz, U. Siegner, F. J. Ahlers,  and H. W. Schumacher, Appl. Phys. Lett. 107, 103501 (2015).
  6. S. Giblin, M. Kataoka, J. Fletcher, P. See, T. Janssen, J. Griffiths, G. Jones, I. Farrer,  and D. Ritchie, Nat. Commun. 3, 930 (2012).
  7. M.-H. Bae, Y.-H. Ahn, M. Seo, Y. Chung, J. D. Fletcher, S. P. Giblin, M. Kataoka,  and N. Kim, Metrologia 52, 195 (2015).
  8. A. Fujiwara, K. Nishiguchi,  and Y. Ono, Appl. Phys. Lett. 92, 042102 (2008).
  9. X. Jehl, B. Voisin, T. Charron, P. Clapera, S. Ray, B. Roche, M. Sanquer, S. Djordjevic, L. Devoille, R. Wacquez,  and M. Vinet, Phys. Rev. X 3, 021012 (2013).
  10. A. Rossi, T. Tanttu, K. Y. Tan, I. Iisakka, R. Zhao, K. W. Chan, G. C. Tettamanzi, S. Rogge, A. S. Dzurak,  and M. Möttönen, Nano Lett. 14, 3405 (2014).
  11. P. Clapera, S. Ray, X. Jehl, M. Sanquer, A. Valentian,  and S. Barraud, Phys. Rev. Applied 4, 044009 (2015).
  12. G. Yamahata, K. Nishiguchi,  and A. Fujiwara, Appl. Phys. Lett. 98, 222104 (2011).
  13. T. Tanttu, A. Rossi, K. Y. Tan, K.-E. Huhtinen, K. W. Chan, M. Möttönen,  and A. S. Dzurak, New J. Phys. 17, 103030 (2015).
  14. L. Fricke, M. Wulf, B. Kaestner, V. Kashcheyevs, J. Timoshenko, P. Nazarov, F. Hohls, P. Mirovsky, B. Mackrodt, R. Dolata, T. Weimann, K. Pierz,  and H. W. Schumacher, Phys. Rev. Lett. 110, 126803 (2013).
  15. L. Fricke, M. Wulf, B. Kaestner, F. Hohls, P. Mirovsky, B. Mackrodt, R. Dolata, T. Weimann, K. Pierz, U. Siegner,  and H. W. Schumacher, Phys. Rev. Lett. 112, 226803 (2014).
  16. G. P. Lansbergen, Y. Ono,  and A. Fujiwara, Nano Lett. 12, 763 (2012).
  17. B. Roche, R.-P. Riwar, B. Voisin, E. Dupont-Ferrier, R. Wacquez, M. Vinet, M. Sanquer, J. Splettstoesser,  and X. Jehl, Nat. Commun. 4, 1581 (2013).
  18. G. C. Tettamanzi, R. Wacquez,  and S. Rogge, New J. Phys. 16, 063036 (2014).
  19. G. Yamahata, K. Nishiguchi,  and A. Fujiwara, Nat. Commun. 5, (2014).
  20. K. W. Chan, M. Möttönen, A. Kemppinen, N. S. Lai, K. Y. Tan, W. H. Lim,  and A. S. Dzurak, Appl. Phys. Lett. 98, 212103 (2011).
  21. B. Roche, B. Voisin, X. Jehl, R. Wacquez, M. Sanquer, M. Vinet, V. Deshpande,  and B. Previtali, Appl. Phys. Lett. 100, 032107 (2012).
  22. M. Tabe, D. Moraru, M. Ligowski, M. Anwar, R. Jablonski, Y. Ono,  and T. Mizuno, Phys. Rev. Lett. 105, 016803 (2010).
  23. B. Voisin, M. Cobian, X. Jehl, M. Vinet, Y.-M. Niquet, C. Delerue, S. de Franceschi,  and M. Sanquer, Phys. Rev. B 89, 161404 (2014).
  24. J. Verduijn, G. C. Tettamanzi,  and S. Rogge, Nano Lett. 13, 1476 (2013).
  25. K. Y. Tan, K. W. Chan, M. Möttönen, A. Morello, C. Yang, J. van Donkelaar, A. Alves, J.-M. Pirkkalainen, D. N. Jamieson, R. G. Clark,  and A. S. Dzurak, Nano Lett. 10, 11 (2010).
  26. G. P. Lansbergen, R. Rahman, C. J. Wellard, I. Woo, J. Caro, N. Collaert, S. Biesemans, G. Klimeck, L. C. L. Hollenberg,  and S. Rogge, Nat. Phys. 4, 656 (2008).
  27. V. N. Golovach, X. Jehl, M. Houzet, M. Pierre, B. Roche, M. Sanquer,  and L. I. Glazman, Phys. Rev. B 83, 075401 (2011).
  28. The resonance of the lowest state under gate 1 (found at  mV) is very weakly coupled and cannot be resolved in this measurement.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description