arxiv:cond-mat/ v1 27 Feb 1996

Similar documents
Single-Electron Parametron

Charge spectrometry with a strongly coupled superconducting single-electron transistor

Quantum Dot Structures Measuring Hamming Distance for Associative Memories

phys4.20 Page 1 - the ac Josephson effect relates the voltage V across a Junction to the temporal change of the phase difference

Noise effects and thermally induced switching errors in Josephson junctions

Analysis of flip flop design using nanoelectronic single electron transistor

Quasiadiabatic switching for metal-island quantum-dot cellular automata

arxiv:cond-mat/ v1 1 Mar 1996

"~l. >-T=SSsäÄs5i5ÄSfSSäsiwä AFRL-SR-BL-TR-98- REPORT DOCUMENTATION f'age "UNCLASSIFIED. ?inal 01 Jan 95 to 51 Mar

Information to energy conversion in an electronic Maxwell s demon and thermodynamics of measurements.

arxiv:chao-dyn/ v1 17 Feb 1999

arxiv:cond-mat/ v1 22 Nov 1993

arxiv:cond-mat/ v1 3 Mar 1994

A single-electron device and circuit simulator

arxiv:cond-mat/ v1 15 Jun 1994

Physics of switches. Luca Gammaitoni NiPS Laboratory

ROBUSTNESS AND POWER DISSIPATION IN QUANTUM-DOT CELLULAR AUTOMATA. A Dissertation. Submitted to the Graduate School. of the University of Notre Dame

EE650R: Reliability Physics of Nanoelectronic Devices Lecture 18: A Broad Introduction to Dielectric Breakdown Date:

The physics of information: from Maxwell s demon to Landauer. Eric Lutz University of Erlangen-Nürnberg

!"#$%&'()*"+,-./*($-"#+"0+'*"12%3+ (#3+4"#&'*"12%3+5'6+6)17-$%1$/)%8*

Quantum Cellular Automata: The Physics of Computing with Arrays of Quantum Dot Molecules

Quantum Size Effect of Two Couple Quantum Dots

ON BROWNIAN COMPUTATION

Maxwell s Demon. Kirk T. McDonald Joseph Henry Laboratories, Princeton University, Princeton, NJ (October 3, 2004; updated September 20, 2016)

DESIGN OF QCA FULL ADDER CIRCUIT USING CORNER APPROACH INVERTER

Environmental effects on Coulomb blockade in a small tunnel junction: A nonperturbative calculation

Bound states of two particles confined to parallel two-dimensional layers and interacting via dipole-dipole or dipole-charge laws

arxiv:cond-mat/ v2 [cond-mat.mes-hall] 11 Dec 1999

Observation of switching in a quantum-dot cellular automata cell

The physical basis of digital computing

arxiv:cond-mat/ v1 [cond-mat.mes-hall] 27 Nov 2001

2.0 Basic Elements of a Quantum Information Processor. 2.1 Classical information processing The carrier of information

Supporting Information

Final Report. Superconducting Qubits for Quantum Computation Contract MDA C-A821/0000

IMPLEMENTATION OF PROGRAMMABLE LOGIC DEVICES IN QUANTUM CELLULAR AUTOMATA TECHNOLOGY

ARTICLE IN PRESS. Physica E

THE INVERTER. Inverter

Physics of switches. Luca Gammaitoni NiPS Laboratory

Design Considerations for Integrated Semiconductor Control Electronics for a Large-scale Solid State Quantum Processor

ECE 438: Digital Integrated Circuits Assignment #4 Solution The Inverter

ANALYSIS OF AN INJECTION-LOCKED BISTABLE SEMICONDUCTOR LASER WITH THE FREQUENCY CHIRPING

Nanoelectronics. Topics

İlke Ercan, PhD Assistant Professor, Electrical & Electronics Eng. Dep. Boğaziçi University, İstanbul, Turkey

Demonstration of a functional quantum-dot cellular automata cell

Reversible Circuit Using Reversible Gate

Performance Enhancement of Reversible Binary to Gray Code Converter Circuit using Feynman gate

The mechanism of nuclear spontaneous time reversal symmetry breaking transfered to digital logic system

Maxwell s Demon. Kirk T. McDonald Joseph Henry Laboratories, Princeton University, Princeton, NJ (October 3, 2004)

Digital Integrated Circuits 2nd Inverter

2 Introduction The quantum measurement is an essential ingredient of investigations of quantum coherent eects. In particular, it is needed to probe ma

Conductance suppression due to correlated electron transport in coupled double quantum dots

MOS Transistors. Prof. Krishna Saraswat. Department of Electrical Engineering Stanford University Stanford, CA

Announcements. EE141- Fall 2002 Lecture 7. MOS Capacitances Inverter Delay Power

Thermal Bias on the Pumped Spin-Current in a Single Quantum Dot

Quantum efficiency of binary-outcome detectors of solid-state qubits

MOSFET: Introduction

Slanted coupling of one-dimensional arrays of small tunnel junctions

arxiv:cond-mat/ v1 29 Dec 1996

5.0 CMOS Inverter. W.Kucewicz VLSICirciuit Design 1

Power Dissipation. Where Does Power Go in CMOS?

Density of States in Superconductor -Normal. Metal-Superconductor Junctions arxiv:cond-mat/ v2 [cond-mat.mes-hall] 7 Nov 1998.

New Model for Drain and Gate Current of Single-Electron Transistor at High Temperature

A Theoretical Consideration of Disorder in a Finite 1D Metal Cluster Chain in a Nanoporous Solid

CMOS INVERTER. Last Lecture. Metrics for qualifying digital circuits. »Cost» Reliability» Speed (delay)»performance

Design for Manufacturability and Power Estimation. Physical issues verification (DSM)

Enhancing Otto-mobile Efficiency via Addition of a Quantum Carnot Cycle

Information entropy and thermal entropy: apples and oranges

Construction of a reconfigurable dynamic logic cell

Even if you're not burning books, destroying information generates heat.

NOVEL QCA CONTROL GATE AND NEW DESIGNING OF MEMORY ON THE BASIS OF QUANTUM DOT CELLULAR AUTOMATA WITH MINIMUM QCA BLOCKS

Supercondcting Qubits

Quantum annealing for problems with ground-state degeneracy

Low mass dileptons from Pb + Au collisions at 158 A GeV

A Tunable Kondo Effect in Quantum Dots

Superconducting Quantum Technologies ABSTRACTS. International Scientific Conference SUPERCONDUCTING QUANTUM TECHNOLOGIES 2018

Lecture 23. CMOS Logic Gates and Digital VLSI I

PHYS 414 Problem Set 4: Demonic refrigerators and eternal sunshine

Designing Information Devices and Systems II Fall 2017 Miki Lustig and Michel Maharbiz Homework 1. This homework is due September 5, 2017, at 11:59AM.

Chapter 5 MOSFET Theory for Submicron Technology

Design of an Online Testable Ternary Circuit from the Truth Table

SUPPLEMENTARY INFORMATION

(a) (b) Supplementary Figure 1. (a) (b) (a) Supplementary Figure 2. (a) (b) (c) (d) (e)

Basic Logic Gate Realization using Quantum Dot Cellular Automata based Reversible Universal Gate

Circuit for Revisable Quantum Multiplier Implementation of Adders with Reversible Logic 1 KONDADASULA VEDA NAGA SAI SRI, 2 M.

Topics. Dynamic CMOS Sequential Design Memory and Control. John A. Chandy Dept. of Electrical and Computer Engineering University of Connecticut

ESE 570: Digital Integrated Circuits and VLSI Fundamentals

Building blocks for nanodevices

Chapter 2: Capacitor And Dielectrics

arxiv:cond-mat/ v1 23 Feb 1993

! Crosstalk. ! Repeaters in Wiring. ! Transmission Lines. " Where transmission lines arise? " Lossless Transmission Line.

CSE140L: Components and Design Techniques for Digital Systems Lab. Power Consumption in Digital Circuits. Pietro Mercati

Control gates as building blocks for reversible computers

Quantum-dot cellular automata

PERFORMANCE ANALYSIS OF CLA CIRCUITS USING SAL AND REVERSIBLE LOGIC GATES FOR ULTRA LOW POWER APPLICATIONS

Realization of programmable logic array using compact reversible logic gates 1

Integrated Circuits & Systems

!"#"$#%&'(&)(*+,'#+-(."//+/,0( 1+#&-,#&'(1$#%&'(%'(2%/%$&'3&'3 %'4+/,#&0(."//4#

Switching in a reversible spin logic gate

Introduction to Theory of Mesoscopic Systems

QUANTUM CRYPTOGRAPHY QUANTUM COMPUTING. Philippe Grangier, Institut d'optique, Orsay. from basic principles to practical realizations.

Transcription:

Single-Electron Parametron: Reversible Computation in a Discrete State System Konstantin K. Likharev 1 and Alexander N. Korotkov 1,2 1 Department of Physics, State University of New York, arxiv:cond-mat/9602140v1 27 Feb 1996 Stony Brook, NY 11794-3800 and 2 Institute of Nuclear Physics, Moscow State University, Moscow 119899 GSP, Russia (September 23, 2018) Abstract We have analyzed energy dissipation in a digital device ( Single-Electron Parametron ) in which discrete degrees of freedom are used for presenting digital information. If the switching speed is not too high, the device may operate reversibly (adiabatically), and the energy dissipation E per bit may be much less than the thermal energy k B T. The energy-time product Eτ is, however, much larger than Planck s constant h, at least in the standard orthodox model of single-electron tunneling, which was used in our calculations. Typeset using REVTEX 1

Computation using any real physical system leads to the dissipation of energy, because of unavoidable coupling between the degrees of freedom which present information, and the environment. In most practical electronic devices (such as semiconductor transistor circuits), energy is dissipated at some rate even in static state, i.e. in the absence of information processing. Some prospective digital devices, such as various Single-Flux-Quantum 1,2 and Single-Electron 3,4 logics, 5 however, do not involve static power consumption, because they present conservative systems where digital information is coded by the choice of a local minimum of potential energy. In this case the energy dissipation is proportional to the number of logic operations. If such a conservative system is switched irreversibly (e.g., as in the RSFQ logic 2 ), the energy loss E per one logic operation is of the order of energy barrier W separating the states. The barrier should be sufficiently high to make the probability p exp( W/k B T) of thermally-induced errors low enough, so that at the physically irreversible computation E min k B T ln(1/p) k B T. Some conservative systems, e.g. the Parametric Quantron, 1,6 are flexible enough to allow independent adjustment ofthepotentialtiltandbarrier height. Ashasbeenshownearlier, 7,8 this flexibility allows physically reversible(adiabatic) switching of thesystem with E k B T, if the information content of the system remains intact. If, however, information is lost during the operation (informationally irreversible computation), the minimum energy loss is at least E min = k B T ln2 per each lost bit. Thus only at completely (physically and informationally) reversible computation, E may be made much less k B T. In this case, E scales as τ 1, where τ is the switching time, so that the product Eτ is fixed. 7 Eτ has the dimensionality of Planck s constant, so the natural question is whether quantum mechanics imposes any fundamental lower bound on this product. A quantitative analysis of a particular reversible system (the Parametric Quantron) has shown 9 that Eτ may be much less then h. The analysis has, however, used the assumption that the potential energy is a function of a continuous degree of freedom (in that particular case the Josephson phase φ). To our knowledge, a similar analysis has never before been carried out for any system presenting information with discrete states. 2

Recently we suggested 10 a system ( Single-Electron Parametron ) based on correlated single-electron tunneling (see, e.g. Ref. 3), which may be used for reversible computation using a discrete degree of freedom the electric charge Q. The goal of the present work was to find the minimum energy dissipation E for this system, and relate it to the switching time τ and to the error probability p. We have shown, in particular, that within the orthodox theory of single-electron tunneling, 11 Eτ is always considerably larger than h. Figure 1a shows the possible structure of a unit cell of the system. It consists of at least three small conducting islands (with capacitances C e 2 /k B T), with the middle island slightly shifted in y-direction. Tunnel barriers with small conductances (G e 2 / h) allow direct transfer of electrons only between the neighboring islands. The system is biased by a periodic clock electric field E(t) perpendicular to axis x. Let us consider the conceptually simplest case when the cell is charged as a whole by a single extra electron. (For practice, the operation using electron-hole pairs may be beneficial, 10,12 but for our present discussion both versions are identical.) When the vertical component E y (t) of the field is lower than a certain threshold value E t, the extra electron is kept inside the middle island. The energy diagram for this OFF state is shown at the top of Fig. 1b. As E y is increased, tunneling of the electron from the middle island into one of the edge islands becomes energy advantageous at E y > E t. If the system is completely x symmetric, this results in spontaneous symmetry breaking, so that the direction of the resulting electric dipole moment of the system is random: P x = Qd eff, d eff d, Q = ±e. However, if the symmetry is broken by a weak additional external field E x (say, induced by the dipole moment of a similar neighboring cell), the direction of the electron tunneling and hence the sign of P x is predetermined by this field. The middle frame in Fig. 1b shows this OFF ON switching stage. Finally, when E y is well above E t, the electron is trapped inside one of the edge islands, even if the signal field E x now favors its transfer in the opposite edge direction ( ON state, the bottom frame in Fig. 1b). In this state, electron transfer may only be achieved via a higher-order co-tunneling process; the probability of this process may be made negligibly small by either decreasing the tunnel 3

conductance or by inserting a few additional islands into the cell. 11 If this parasitic process is negligible, the cell has a fixed dipole moment and may serve as a robust source of signal field E x for similar neighboring cells. After this source has been used, the system is reset into the OFF state during the corresponding part of the clock cycle, when E(t) drops below E t again. It is evident that the operation of the cell is quite similar to that of the Parametric Quantron, 1 exceptnowtheinformationispresentedbyadiscretevariable,q = ±e. Similarly to the Parametric Quantron, the Single-Electron Parametron may be used for reversible transfer and processing of information. 10,12 For example, Figure 2 shows a possible structure of a shift register. In each neighboring cell, the extra charge sign alternates, while the direction of the middle island shift within plane yz is changed by Θ = π 2π/M, M > 2 (in Fig. 2, M = 3). The clock field E(t) has a fixed magnitude E > E t, but rotates within plane yz, providing periodic switching ON and OFF of the cells, with the phase shift 2π/M between the neighboring cells. At an appropriate choice of E and the distance between the cells, 12 the orientation of the dipole moment of the cells in ON state determines the direction of the field E x and hence the direction of electron tunneling in the neighboring cell which is being switched OFF ON. As a result, the information is being re-written from cell to cell, and thus transferred over M cells each clock period. Reversible logic operations may be implemented in a similar way, e.g. by using majority gates with additional output cells. 9 Within the orthodox theory 11 all properties of the system may be found from solution of the system of master equations for the probabilities p i (t) to find the extra electron in the middle (i = m), left (i = l), and right (i = r) islands: d dt p i = j=m,l,r(p j Γ ji p i Γ ij ), i p i (t) = 1, p m (0) = 1, (1) where in our case the tunneling rate matrix Γ has only four nonvanishing components: Γ ± mr = ±GW e 2 {1 exp[ W/k B T]}, Γ± ml = ±G(W ) e 2 {1 exp[ (W )/k B T]}. (2) Here Γ ij Γ ij,γ ij Γ ji, while W(t) E y (t)d const is the energy difference between the charge configurations with the extra electron in the middle and right islands. We will 4

accept that near the decision-making point of the ON OFF switching (t 0 in Fig. 1c) the difference is a linear function of time: W = αt. 2dE x is the energy difference between the left and right ON states. In order to operate with a low error probability p 1, at the decision-making moment this difference should be large enough. Without the loss of generality, we may assume > 0; then p can be found from the solution of the master equation as p l ( ), while the average energy E(t) dissipated by moment t can be calculated as E(t) = t {W(t) dp r dt [W(t) ] dp l }dt; (3) dt we will be mostly interested in the net dissipation E E( ). The solution of equations (1) (3) yields the following results. With an accuracy sufficient for our final result, the total error probability p may be calculated as a maximum of probabilities of the thermal and dynamic errors. The thermal error may occur due to thermally-activated tunneling to the wrong state (in our case, l), and its probability p therm may always be expressed as exp( /k B T). (Because of that, we will restrict our discussion to the limit k B T.) The dynamic error occurs when the switching speed α is too high, and the system remains in the initial (symmetric) state up to the moment when tunneling to the upper energy level becomes possible. If δ αe 2 /G k B T 1, the dynamic error dominates and its probability is given by expression p dyn = Kγexp( 1 2γ ), K = 1 π 2γ 1 exp( 4γ3/2 4γ )[1 Erf( 1 2 γ )] = 1 n=1 (2n1)! n! ( γ) n, (4) where γ αe 2 /G 2. In order to keep p dyn 1, γ should be much smaller than 1, so that one can use Eq. (4) with K = 1. Energy dissipation depends on another dimensionless parameter, β αe 2 /G(k B T) 2 = ( /k B T) 2 γ = ( /k B T)δ. Like γ and δ, parameter β is also proportional to the switching speed α, but is much larger than both of them (because k B T) and may be comparable to unity. In the low-speed limit β 1, the switching process is adiabatic. It consists of 5

numerous tunneling events (back and forth between m and r) taking place within the energy interval k B T around the point W(t) = 0. In this case E = κβk B T, where κ = e x (e x 1) dx 0.426, (5) x(1e x ) 3 so that for this (reversible) process E k B T. Notice that E = κβk B T = καe 2 /Gk B T decreases when temperature increases. Our model allows not only to calculate the net dissipation E, but also follow the time dynamics of energy transfer between the system and the environment ( heat bath ) during the switching process. During the first half of the process (when W(t) 0) the energy E 1 E(0) = T ln2 E is borrowed from the heat bath (which, hence, is cooled), while virtually the same amount E 2 E E(0) is returned back to the heat bath during the second half of the process (W(t) 0). This exchange is directly related as E(t) = T S(t) to the temporal increase and consequent decrease of the entropy corresponding to the degree of freedom used to code information (in this particular case, the polarization P x ). At the moment when W = 0, the system may be in either of two states (p m = p r = 1/2), i.e. S = k B ln2 has been acquired in the comparison with the definite initial state (p m = 1, p r = 0). By the end of the switching (W k B T) the entropy is restored to the initial value since the state is definite again (p m = 0, p r = 1). Finite switching speed decreases E 1 and increases E 2 (see the dotted lines in Fig. 3). In the limit β 1, the speed of energy change is so high that switching may take place only at W > 0, but within a much larger interval of energies: W β 1/2 k B T. The average energy dissipation for this (irreversible) process is of the same order, i.e. much larger than k B T and independent of temperature: E = (πβ/2) 1/2 k B T = (πe 2 α/2g) 1/2. The results of numerical calculation of E for intermediate values of β are presented by the solid line in Fig. 3. Taking into account that at p 1 the parameter τ = /α may be considered as the duration of the switching process (Fig. 1c), all our asymptotic results may be summarized as follows: 6

Eτ = h GR Q 0.67 ln 1, for δ,β 1, p 1.97β 1/2 ln 1, for δ 1 β, p 2.78(ln 1 2pln(1/p) )1/2, for 1 δ,β, (6) where R Q =π h/2e 2 6.45kΩ is the quantum unit of resistance. Since the orthodox theory is valid only at GR Q 1, within this theory Eτ h for any switching speed. To summarize, we have shown that reversible computation with the energy dissipation E per bit well below k B T may be implemented in a physical system with discrete states. The quantum bound for the product Eτ, obtained within our concrete model is, however, much higher than that obtained earlier for a system with continuous degrees of freedom. 9 Apparently the h-limit for Eτ may be overcome in the case of islands with discrete spectra of electron energies, 13,14 though this may require an exponentially high energy barrier W during the ON state of the cell. A quantitative analysis of this opportunity is in progress. Useful discussions with D. V. Averin and T. Usuki are gratefully acknowledged. The work was supported in part by ONR Grant #N00014-93-1-0880 and AFOSR Grant #F49620-95- 1-0044. 7

REFERENCES 1 K. K. Likharev, IEEE Trans. Magn. 13, 242 (1976). 2 K. K. Likharev and V. K. Semenov, IEEE Trans. Magn. 1, 3 (1991). 3 D. V. Averin and K. K. Likharev, in: Single Charge Tunneling, edited by H. Grabert and M. H. Devoret (Plenum, New York, 1992), p. 311. 4 A. N. Korotkov, Appl. Phys. Lett. 67, 2412 (1995). 5 In this work we consider only the case when the state of the system can be described as a classical variable. These classical systems (whose dynamics may nevertheless require quantum-mechanical description) may provide reliable computation even at substantial coupling to the environment, in contrast to the so-called quantum computing devices see, e.g., I. L. Chuang et al., Science 270, 1633 (1995). 6 This system was re-invented later as a Flux Quantum Parametron see K. F. Loe and E. Goto, IEEE Trans. Magn. 21, 884 (1985). 7 R. Landauer, IBM J. Res. Devel. 5, 183 (1961); R. W. Keyes and R. Landauer, ibid. 14, 152 (1970). 8 C. Bennett, IBM J. Res. Devel. 17, 525 (1973). 9 K. K. Likharev, Int. J. Teor. Phys. 21, 311 (1982). 10 K. K. Likharev and A. N. Korotkov, in: Abstr. of IWCE 95 (Tempe, AZ, 1995), p. 5. 11 D. V. Averin and K. K. Likharev, in: Mesoscopic Phenomena in Solids, edited by B. L. Altshuler, P. A. Lee, and R. A. Webb (Elsevier, Amsterdam, 1991), p. 173. 12 A. N. Korotkov and K. K. Likharev, to be published. 13 D. V. Averin and A. N. Korotkov, J. Low Temp. Phys. 80, 173 (1990). 14 M. Kastner, Rev. Mod. Phys. 64, 849 (1992). 8

FIGURES FIG. 1. The Single-Electron Parametron: (a) 3-island version of the system, (b) its energy diagram for three values of the clock field E y, and (c) energy of the extra electron in various islands as a function of time, close to the decision-making moment t 0. FIG. 2. Top (left) and side (right) views of a shift register using an array Single-Electron Parametron cells. Clock field E(t) rotates in yz plane. Digital bits are coded by positions of the extra charges in ON state of the cells, and are propagated from the top to the bottom, over M = 3 cells during one clock period. FIG. 3. Components of the energy exchange between the Parametron and the heat bath as functions of the process speed α = dw/dt. Dotted lines: average energy flow E 1 from the heat bath to the Parametron during the first half of the process (W 0) and the average flow E 2 from the device back into the heat bath during its second half (W 0), respectively. Solid line: net energy dissipation E = E 2 E 1. Dashed lines show the low-speed (adiabatic) and high-speed (diabatic) asymptotes of the function E(α) see formulas in the text. 9

island l -d y island m island r d' d x (a) tunnel barriers -e (b) State OFF (E <E, W<0) y t Switching OFF-->ON (E y ~E t, W~0) State ON (E y >E t, W>0) (c) island l 0 island m W(t) island r t Fig. 1

Binary 0 Binary 1 OFF Rotating Field E(t) State ON 1 - Switching OFF -> ON1 - State OFF - State ON 0 - - Switching OFF -> ON 0 State OFF - State ON 0 y y x z Fig. 2

10.0 ε/k B Tln2 1.0 ε 2 ε 1 ε=ε 2 -ε 1 0.1 0.1 1.0 10.0 100.0 β=αe 2 /G(k B T) 2 Fig. 3