Cloud Hosting

Quantum computing implementation

Superconducting quantum computing is an implementation of a quantum computer in superconducting electronic circuits. Research in superconducting quantum computing is conducted by Google,[1] IBM,[2] BBN Technologies,[3] Rigetti,[4] and Intel.[5] as of May2016[update], up to nine fully controllable qubits are demonstrated in a 1D array,[6] up to sixteen in a 2D architecture.[2]

More than two thousand superconducting qubits are in a commercial product by D-Wave Systems, however these qubits implement quantum annealing instead of a universal model of quantum computation.

Classical computation models rely on physical implementations consistent with the laws of classical mechanics.[8] It is known, however, that the classical description is only accurate for specific cases, while the more general description of nature is given by quantum mechanics. Quantum computation studies the application of quantum phenomena, that are beyond the scope of classical approximation, for information processing and communication. Various models of quantum computation exist, however the most popular models incorporate the concepts of qubits and quantum gates. A qubit is a generalization of a bit - a system with two possible states, that can be in a quantum superposition of both. A quantum gate is a generalization of a logic gate: it describes the transformation that one or more qubits will experience after the gate is applied on them, given their initial state. The physical implementation of qubits and gates is difficult, for the same reasons that quantum phenomena are hard to observe in everyday life. One approach is to implement the quantum computers in superconductors, where the quantum effects become macroscopic, though at a price of extremely low operation temperatures.

In a superconductor, the basic charge carriers are pairs of electrons (known as Cooper pairs), rather than the single electrons in a normal conductor. The total spin of a Cooper pair is an integer number, thus the Cooper pairs are bosons (while the single electrons in the normal conductor are fermions). Cooled bosons, contrary to cooled fermions, are allowed to occupy a single quantum energy level, in an effect known as the Bose-Einstein condensate. In a classical interpretation it would correspond to multiple particles occupying the same position in space and having an equal momentum, effectively behaving as a single particle.

At every point of a superconducting electronic circuit (that is a network of electrical elements), the condensate wave function describing the charge flow is well-defined by a specific complex probability amplitude. In a normal conductor electrical circuit, the same quantum description is true for individual charge carriers, however the various wave functions are averaged in the macroscopic analysis, making it impossible to observe quantum effects. The condensate wave function allows designing and measuring macroscopic quantum effects. For example, only a discrete number of magnetic flux quanta penetrates a superconducting loop, similarly to the discrete atomic energy levels in the Bohr model. In both cases, the quantization is a result of the complex amplitude continuity. Differing from the microscopic quantum systems (such as atoms or photons) used for implementations of quantum computers, the parameters of the superconducting circuits may be designed by setting the (classical) values of the electrical elements that compose them, e.g. adjusting the capacitance or inductance.

In order to obtain a quantum mechanical description of an electrical circuit a few steps are required. First, all the electrical elements are described with the condensate wave function amplitude and phase, rather than with the closely related macroscopic current and voltage description used for classical circuits. For example, a square of the wave function amplitude at some point in space is the probability of finding a charge carrier there, hence the square of the amplitude corresponds to the classical charge distribution. Second, generalized Kirchhoff's circuit laws are applied at every node of the circuit network to obtain the equations of motion. Finally, the equations of motion are reformulated to Lagrangian mechanics and a quantum Hamiltonian is derived.

The devices are typically designed in the radio-frequency spectrum, cooled down in dilution refrigerators below 100mK and addressed with conventional electronic instruments, e.g. frequency synthesizers and spectrum analyzers. Typical dimensions on the scale of micrometers, with sub-micrometer resolution, allow a convenient design of a quantum Hamiltonian with the well-established integrated circuit technology.

A distinguishing feature of superconducting quantum circuits is the usage of a Josephson junction - an electrical element non existent in normal conductors. A junction is a weak connection between two leads of a superconducting wire, usually implemented as a thin layer of insulator with a shadow evaporation technique. The condensate wave functions on the two sides of the junction are weakly correlated - they are allowed to have different superconducting phases, contrary to the case of a continuous superconducting wire, where the superconducting wave function must be continuous. The current through the junction occurs by quantum tunneling. This is used to create a non-linear inductance which is essential for qubit design, as it allows a design of anharmonic oscillators. A quantum harmonic oscillator cannot be used as a qubit, as there is no way to address only two of its states.

The three superconducting qubit archetypes are the phase, charge and flux qubits, though many hybridizations exist (Fluxonium,[9] Transmon,[10] Xmon,[11] Quantronium[12]). For any qubit implementation, the logical quantum states { | 0 , | 1 } {displaystyle {|0rangle ,|1rangle }} are to be mapped to the different states of the physical system, typically to the discrete (quantized) energy levels or to their quantum superpositions. In the charge qubit, different energy levels correspond to an integer number of Cooper pairs on a superconducting island. In the flux qubit, the energy levels correspond to different integer numbers of magnetic flux quanta trapped in a superconducting ring. In the phase qubit, the energy levels correspond to different quantum charge oscillation amplitudes across a Josephson junction, where the charge and the phase are analogous to momentum and position correspondingly of a quantum harmonic oscillator. Note that the phase here is the complex argument of the superconducting wavefunction, also known as the superconducting order parameter, not the phase between the different states of the qubit.

In the table below, the three archetypes are reviewed. In the first row, the qubit electrical circuit diagram is presented. In the second, the quantum Hamiltonian derived from the circuit is shown. Generally, the Hamiltonian can be divided to a "kinetic" and "potential" parts, in analogy to a particle in a potential well. The particle mass corresponds to some inverse function of the circuit capacitance, while the shape of the potential is governed by the regular inductors and Josephson junctions. One of the first challenges in qubit design is to shape the potential well and to choose the particle mass in a way that the energy separation between specific two of the energy levels will differ from all other inter-level energy separations in the system. These two levels will be used as the logical states of the qubit. The schematic wave solutions in the third row of the table depict the complex amplitude of the phase variable. In other words, if a phase of the qubit is measured while the qubit is in a specific state, there is a non-zero probability to measure a specific value only where the depicted wave function oscillates. All three rows are essentially three different presentations of the same physical system.

Type

Aspect

A superconducting island (encircled with a dashed line) defined between the leads of a capacitor with capacitance C {displaystyle C} and a Josephson junction with energy E J {displaystyle E_{J}} is biased by voltage U {displaystyle U}

A superconducting loop with inductance L {displaystyle L} is interrupted by a junction with Josephson energy E J {displaystyle E_{J}} . Bias flux {displaystyle Phi } is induced by a flux line with a current I 0 {displaystyle I_{0}}

Josephson junction with energy parameter E J {displaystyle E_{J}} is biased by a current I 0 {displaystyle I_{0}}

H = E C ( N N g ) 2 E J cos {displaystyle H=E_{C}(N-N_{g})^{2}-E_{J}cos phi } ,where N {displaystyle N} is the number of Cooper pairs to tunnel the junction, N g = C V 0 / 2 e {displaystyle N_{g}=CV_{0}/2e} is the charge on the capacitor in units of Cooper pairs number, E C = ( 2 e ) 2 / 2 ( C J + C ) {displaystyle E_{C}=(2e)^{2}/2(C_{J}+C)} is the charging energy associated with both the capacitance C {displaystyle C} and the Josephson junction capacitance C J {displaystyle C_{J}} , and {displaystyle phi } is the superconducting wave function phase difference across the junction.

H = q 2 2 C J + ( 0 2 ) 2 2 2 L E J cos [ 2 0 ] {displaystyle H={frac {q^{2}}{2C_{J}}}+left({frac {Phi _{0}}{2pi }}right)^{2}{frac {phi ^{2}}{2L}}-E_{J}cos left[phi -Phi {frac {2pi }{Phi _{0}}}right]} ,where q {displaystyle q} is the charge on the junction capacitance C J {displaystyle C_{J}} and {displaystyle phi } is the superconducting wave function phase difference across the Josephson junction. {displaystyle phi } is allowed to take values greater than 2 {displaystyle 2pi } , and thus is alternatively defined as the time integral of voltage along the inductance L {displaystyle L} .

H = ( 2 e ) 2 2 C J q 2 I 0 0 2 E J cos {displaystyle H={frac {(2e)^{2}}{2C_{J}}}q^{2}-I_{0}{frac {Phi _{0}}{2pi }}phi -E_{J}cos phi } , where C J {displaystyle C_{J}} is the capacitance associated with the Josephson junction, 0 {displaystyle Phi _{0}} is the magnetic flux quantum, q {displaystyle q} is the charge on the junction capacitance C J {displaystyle C_{J}} and {displaystyle phi } is the phase across the junction.

The potential part of the Hamiltonian, E J cos {displaystyle -E_{J}cos phi } , is depicted with the thick red line. Schematic wave function solutions are depicted with thin lines, lifted to their appropriate energy level for clarity. Only the solid wave functions are used for computation. The bias voltage is set so that N g = 1 2 {displaystyle N_{g}={frac {1}{2}}} , minimizing the energy gap between | 0 {displaystyle |0rangle } and | 1 {displaystyle |1rangle } , thus making the gap different from other energy gaps (e.g. the gap between | 1 {displaystyle |1rangle } and | 2 {displaystyle |2rangle } ). The difference in gaps allows addressing transitions from | 0 {displaystyle |0rangle } to | 1 {displaystyle |1rangle } and vice versa only, without populating other states, thus effectively treating the circuit as a two-level system (qubit).

The potential part of the Hamiltonian, ( 0 2 ) 2 2 2 L E J cos [ 2 0 ] {displaystyle left({frac {Phi _{0}}{2pi }}right)^{2}{frac {phi ^{2}}{2L}}-E_{J}cos left[phi -Phi {frac {2pi }{Phi _{0}}}right]} , plotted for the bias flux = 0 / 2 {displaystyle Phi =Phi _{0}/2} , is depicted with the thick red line. Schematic wave function solutions are depicted with thin lines, lifted to their appropriate energy level for clarity. Only the solid wave functions are used for computation. Different wells correspond to a different number of flux quanta trapped in the superconducting loops. The two lower states correspond to a symmetrical and an antisymmetrical superposition of zero or single trapped flux quanta, sometimes denoted as clockwise and counterclockwise loop current states: | 0 = [ | + | ] / 2 {displaystyle |0rangle =left[|circlearrowleft rangle +|circlearrowright rangle right]/{sqrt {2}}} and | 1 = [ | | ] / 2 {displaystyle |1rangle =left[|circlearrowleft rangle -|circlearrowright rangle right]/{sqrt {2}}} .

The so-called "washboard" potential part of the Hamiltonian, I 0 0 2 E J cos {displaystyle -I_{0}{frac {Phi _{0}}{2pi }}phi -E_{J}cos phi } , is depicted with the thick red line. Schematic wave function solutions are depicted with thin lines, lifted to their appropriate energy level for clarity. Only the solid wave functions are used for computation. The bias current is adjusted to make the wells shallow enough to contain exactly two localized wave functions. A slight increase in the bias current causes a selective "spill" of the higher energy state ( | 1 {displaystyle |1rangle } ), expressed with a measurable voltage spike - a mechanism commonly used for phase qubit measurement.

The GHz energy gap between the energy levels of a superconducting qubit is intentionally designed to be compatible with available electronic equipment, due to the terahertz gap - lack of equipment in the higher frequency band. In addition, the superconductor energy gap implies a top limit of operation below ~1THz (beyond it, the Cooper pairs break). On the other hand, the energy level separation cannot be too small due to cooling considerations: a temperature of 1K implies energy fluctuations of 20GHz. Temperatures of tens of mili-Kelvin achieved in dilution refrigerators allow qubit operation at a ~5GHz energy level separation. The qubit energy level separation may often be adjusted by means of controlling a dedicated bias current line, providing a "knob" to fine tune the qubit parameters.

An arbitrary single qubit gate is achieved by rotation in the Bloch sphere. The rotations between the different energy levels of a single qubit are induced by microwave pulses sent to an antenna or transmission line coupled to the qubit, with a frequency resonant with the energy separation between the levels. Individual qubits may be addressed by a dedicated transmission line, or by a shared one if the other qubits are off resonance. The axis of rotation is set by quadrature amplitude modulation of the microwave pulse, while the pulse length determines the angle of rotation.[14]

More formally, following the notation of,[14] for a driving signal

E ( t ) = E x ( t ) cos ( d t ) + E y ( t ) sin ( d t ) {displaystyle {mathcal {E}}(t)={mathcal {E}}^{x}(t)cos(omega _{d}t)+{mathcal {E}}^{y}(t)sin(omega _{d}t)}

of frequency d {displaystyle omega _{d}} , a driven qubit Hamiltonian in a rotating wave approximation is

H R / = ( d ) | 1 1 | + E x ( t ) 2 x + E y ( t ) 2 y {displaystyle H^{R}/hbar =(omega -omega _{d})|1rangle langle 1|+{frac {{mathcal {E}}^{x}(t)}{2}}sigma _{x}+{frac {{mathcal {E}}^{y}(t)}{2}}sigma _{y}} ,

where {displaystyle omega } is the qubit resonance and x , y {displaystyle sigma _{x},sigma _{y}} are Pauli matrices.

In order to implement a rotation about the X {displaystyle X} axis, one can set E y ( t ) = 0 {displaystyle {mathcal {E}}^{y}(t)=0} and apply the microwave pulse at frequency d = {displaystyle omega _{d}=omega } for time t g {displaystyle t_{g}} . The resulting transformation is

U x = exp { i 0 t g H R d t } = exp { i 0 t g E x ( t ) d t x / 2 } {displaystyle U_{x}=exp left{-{frac {i}{hbar }}int _{0}^{t_{g}}H^{R}dtright}=exp left{-iint _{0}^{t_{g}}{mathcal {E}}^{x}(t)dtcdot sigma _{x}/2right}} ,

that is exactly the rotation operator R X ( ) {displaystyle R_{X}(theta )} by angle = 0 t g E x ( t ) d t {displaystyle theta =int _{0}^{t_{g}}{mathcal {E}}^{x}(t)dt} about the X {displaystyle X} axis in the Bloch sphere. An arbitrary rotation about the Y {displaystyle Y} axis can be implemented in a similar way. Showing the two rotation operators is sufficient for universality, as every single qubit unitary operator U {displaystyle U} may be presented as U = R X ( 1 ) R Y ( 2 ) R X ( 3 ) {displaystyle U=R_{X}(theta _{1})R_{Y}(theta _{2})R_{X}(theta _{3})} (up to a global phase, that is physically unimportant) by a procedure known as the X Y {displaystyle X-Y} decomposition.[15]

For example, setting 0 t g E x ( t ) d t = {displaystyle int _{0}^{t_{g}}{mathcal {E}}^{x}(t)dt=pi } results with a transformation

U x = exp { i 0 t g E x ( t ) d t x / 2 } = e i x / 2 = i x {displaystyle U_{x}=exp left{-iint _{0}^{t_{g}}{mathcal {E}}^{x}(t)dtcdot sigma _{x}/2right}=e^{-ipi sigma _{x}/2}=-isigma _{x}} ,

that is known as the NOT gate (up to the global phase i {displaystyle -i} ).

Coupling qubits is essential for implementing 2-qubit gates. Coupling two qubits may be achieved by connecting them to an intermediate electrical coupling circuit. The circuit might be a fixed element, such as a capacitor, or controllable, such as a DC-SQUID. In the first case, decoupling the qubits (during the time the gate is off) is achieved by tuning the qubits out of resonance one from another, i.e. making the energy gaps between their computational states different.[16] This approach is inherently limited to allow nearest-neighbor coupling only, as a physical electrical circuit is to be lay out in between the connected qubits. Notably, D-Wave Systems' nearest-neighbor coupling achieves a highly connected unit cell of 8 qubits in the Chimera graph configuration. Generally, quantum algorithms require coupling between arbitrary qubits, therefore the connectivity limitation is likely to require multiple swap operations, limiting the length of the possible quantum computation before the processor decoherence.

Another method of coupling two or more qubits is by coupling them to an intermediate quantum bus. The quantum bus is often implemented as a microwave cavity, modeled by a quantum harmonic oscillator. Coupled qubits may be brought in and out of resonance with the bus and one with the other, hence eliminating the nearest-neighbor limitation. The formalism used to describe this coupling is cavity quantum electrodynamics, where qubits are analogous to atoms interacting with optical photon cavity, with the difference of GHz rather than THz regime of the electromagnetic radiation.

One popular gating mechanism includes two qubits and a bus, all tuned to different energy level separations. Applying microwave excitation to the first qubit, with a frequency resonant with the second qubit, causes a x {displaystyle sigma _{x}} rotation of the second qubit. The rotation direction depends on the state of the first qubit, allowing a controlled phase gate construction.[17]

More formally, following the notation of,[17] the drive Hamiltonian describing the system excited through the first qubit driving line is

H D / = A ( t ) cos ( ~ 2 t ) ( x I J 12 z x + m 12 I x ) {displaystyle H_{D}/hbar =A(t)cos({tilde {omega }}_{2}t)left(sigma _{x}otimes I-{frac {J}{Delta _{12}}}sigma _{z}otimes sigma _{x}+m_{12}Iotimes sigma _{x}right)} ,

where A ( t ) {displaystyle A(t)} is the shape of the microwave pulse in time, ~ 2 {displaystyle {tilde {omega }}_{2}} is the resonance frequency of the second qubit, { I , x , y , z } {displaystyle {I,sigma _{x},sigma _{y},sigma _{z}}} are the Pauli matrices, J {displaystyle J} is the coupling coefficient between the two qubits via the resonator, 12 1 2 {displaystyle Delta _{12}equiv omega _{1}-omega _{2}} is the qubit detuning, m 12 {displaystyle m_{12}} is the stray (unwanted) coupling between qubits and {displaystyle hbar } is Planck constant divided by 2 {displaystyle 2pi } . The time integral over A ( t ) {displaystyle A(t)} determines the angle of rotation. Unwanted rotations due to the first and third terms of the Hamiltonian can be compensated with single qubit operations. The remaining part is exactly the controlled-X gate.

Architecture-specific readout (measurement) mechanisms exist. The readout of a phase qubit is explained in the qubit archetypes table above. A state of the flux qubit is often read by an adjust DC-SQUID magnetometer. A more general readout scheme includes a coupling to a microwave resonator, where the resonance frequency of the resonator is shifted by the qubit state.[18]

The list of DiVincenzo's criteria for a physical system to implement a logical qubit is satisfied by the superconducting implementation. The challenges currently faced by the superconducting approach are mostly in the field of microwave engineering.[18]

More here:
Superconducting quantum computing - Wikipedia