Lecture Notes on Quantum Computing
Introduction
This lecture extends our discussion of quantum computing, focusing on the critical requirement of two-qubit gates for universal quantum computation. Having previously examined single-photon quantum gates, we now address the methods for coupling photons and other quantum systems to achieve this.
The topics covered in this lecture include:
Two-Photon Quantum Gates: Exploring cross-phase modulation as a mechanism for realizing two-qubit gates in optical quantum computing.
Interaction with External Fields: Introducing time-dependent perturbation theory and the concept of resonance to understand how external fields manipulate quantum systems.
Superconducting Qubits: Examining Josephson junctions and their application in charge and flux qubits for superconducting quantum computing.
Density Matrix and Superoperators: Introducing these tools for describing quantum system evolution, particularly for mixed states and open systems.
This lecture aims to provide a theoretical foundation and discuss practical challenges in the development of quantum computers.
Two-Photon Quantum Gates
For universal quantum computation, two-qubit gates capable of entangling qubits are indispensable. In optical quantum computing, this necessitates a mechanism for photons to interact. Utilizing non-linear media to mediate photon interactions is a promising strategy.
Cross-Phase Modulation (XPM)
Cross-phase modulation (XPM) is a non-linear optical phenomenon where the refractive index, and thus the phase of one optical beam, is altered by the intensity of another. In quantum computing, XPM can facilitate two-photon interactions.
The Hamiltonian for cross-phase modulation is given by: \[\mathcal{H}_{\text{XPM}} = \chi \hat{n}_A \hat{n}_B\] where \(\chi\) is the non-linear susceptibility, and \(\hat{n}_A = a^\dagger a\) and \(\hat{n}_B = b^\dagger b\) are the photon number operators for modes A and B, with \(a^\dagger, a, b^\dagger, b\) being the creation and annihilation operators.
The action of \(\mathcal{H}_{\text{XPM}}\) on the two-photon computational basis states is as follows:
\(\hat{n}_A \hat{n}_B \left|{00}\right\rangle = 0\)
\(\hat{n}_A \hat{n}_B \left|{01}\right\rangle = 0\)
\(\hat{n}_A \hat{n}_B \left|{10}\right\rangle = 0\)
\(\hat{n}_A \hat{n}_B \left|{11}\right\rangle = \left|{11}\right\rangle\)
The Hamiltonian only affects the \(\left|{11}\right\rangle\) state. The time evolution operator \(U(t) = e^{-i \frac{t}{\hbar} \mathcal{H}_{\text{XPM}}}\) introduces a phase shift only to the \(\left|{11}\right\rangle\) state: \[\begin{aligned} U(t) \left|{00}\right\rangle &= \left|{00}\right\rangle \\ U(t) \left|{01}\right\rangle &= \left|{01}\right\rangle \\ U(t) \left|{10}\right\rangle &= \left|{10}\right\rangle \\ U(t) \left|{11}\right\rangle &= e^{-i \frac{\chi t}{\hbar}} \left|{11}\right\rangle = e^{i\phi} \left|{11}\right\rangle\end{aligned}\] where \(\phi = -\frac{\chi t}{\hbar}\). By adjusting the interaction time \(t\) (or length of the non-linear medium), we can set \(\phi = \pi\). This results in a gate \(K\) where only the \(\left|{11}\right\rangle\) state acquires a phase of \(-1\): \[\begin{aligned} K \left|{00}\right\rangle &= \left|{00}\right\rangle \\ K \left|{01}\right\rangle &= \left|{01}\right\rangle \\ K \left|{10}\right\rangle &= \left|{10}\right\rangle \\ K \left|{11}\right\rangle &= -\left|{11}\right\rangle\end{aligned}\]
Definition: The \(K\) gate is a two-qubit gate that applies a phase shift of \(-1\) only to the \(\left|{11}\right\rangle\) state in the computational basis \(\{\left|{00}\right\rangle, \left|{01}\right\rangle, \left|{10}\right\rangle, \left|{11}\right\rangle\}\). Its action on the computational basis states is defined as: \[\begin{aligned} K \left|{00}\right\rangle &= \left|{00}\right\rangle \\ K \left|{01}\right\rangle &= \left|{01}\right\rangle \\ K \left|{10}\right\rangle &= \left|{10}\right\rangle \\ K \left|{11}\right\rangle &= -\left|{11}\right\rangle\end{aligned}\] In the computational basis \(\{\left|{00}\right\rangle, \left|{01}\right\rangle, \left|{10}\right\rangle, \left|{11}\right\rangle\}\), the \(K\) gate is represented by the matrix: \[K = \begin{pmatrix} 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & -1 \end{pmatrix}\]
CNOT Gate Construction
The \(K\) gate, combined with single-qubit gates, can implement a CNOT gate. The construction is given by \(CNOT = (I \otimes H) K (I \otimes H)\), where \(H\) is the Hadamard gate and \(I\) is the identity gate.
Definition: The Hadamard gate \(H\) is a single-qubit gate represented by the matrix: \[H = \frac{1}{\sqrt{2}} \begin{pmatrix} 1 & 1 \\ 1 & -1 \end{pmatrix}\]
and \(I \otimes H\) is:
Definition: The tensor product of the identity gate \(I\) and the Hadamard gate \(H\), denoted as \(I \otimes H\), is a two-qubit gate represented by the matrix: \[I \otimes H = \frac{1}{\sqrt{2}} \begin{pmatrix} 1 & 1 & 0 & 0 \\ 1 & -1 & 0 & 0 \\ 0 & 0 & 1 & 1 \\ 0 & 0 & 1 & -1 \end{pmatrix}\]
Theorem: The CNOT gate can be constructed using the \(K\) gate and Hadamard gates as follows: \[CNOT = (I \otimes H) K (I \otimes H)\] Verification: Let’s verify the product \((I \otimes H) K (I \otimes H)\): \[K (I \otimes H) = \frac{1}{\sqrt{2}} \begin{pmatrix} 1 & 1 & 0 & 0 \\ 1 & -1 & 0 & 0 \\ 0 & 0 & 1 & 1 \\ 0 & 0 & -1 & 1 \end{pmatrix}\] \[(I \otimes H) K (I \otimes H) = \frac{1}{2} \begin{pmatrix} 1 & 1 & 0 & 0 \\ 1 & -1 & 0 & 0 \\ 0 & 0 & 1 & 1 \\ 0 & 0 & 1 & -1 \end{pmatrix} \begin{pmatrix} 1 & 1 & 0 & 0 \\ 1 & -1 & 0 & 0 \\ 0 & 0 & 1 & 1 \\ 0 & 0 & -1 & 1 \end{pmatrix} = \begin{pmatrix} 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 0 & 1 \\ 0 & 0 & 1 & 0 \end{pmatrix}\] This resulting matrix is indeed the CNOT gate, with the first qubit as control and the second as target.
Practical Limitations of Non-linear Media
Despite the promise of cross-phase modulation, practical realization of two-photon gates using non-linear media is challenging due to:
Weak Non-linearities: Typical non-linear media exhibit weak non-linear effects. Achieving a \(\pi\) phase shift requires either high photon intensities or long interaction lengths.
Photon Absorption: Long interaction lengths increase photon absorption in the medium, leading to signal loss and reduced gate fidelity.
Higher-Order Non-linearities: High intensities can induce unwanted higher-order non-linear effects, complicating the desired XPM interaction.
These limitations necessitate exploring alternative methods to enhance photon interactions, such as cavity quantum electrodynamics (QED) or intermediary systems for indirect photon coupling.
Interaction of Quantum Systems with External Fields
External fields, particularly electromagnetic fields, are essential for manipulating quantum systems, enabling transitions between energy levels and control of quantum states. Understanding the response of quantum systems to time-dependent perturbations is crucial for quantum control and computation.
Time-Dependent Perturbation Theory
Time-dependent perturbation theory approximates solutions to the Schrödinger equation when the Hamiltonian is decomposed into a time-independent part \(H_0\) and a perturbation \(V(t)\), where \(V(t) \ll H_0\):
Definition: In time-dependent perturbation theory, the total Hamiltonian \(H(t)\) is decomposed into a time-independent Hamiltonian \(H_0\) and a time-dependent perturbation \(V(t)\), where the perturbation is assumed to be much smaller than the time-independent part: \[H(t) = H_0 + V(t)\] with \(V(t) \ll H_0\).
Let \(\{\left|{\phi_n}\right\rangle\}\) and \(E_n\) be the eigenstates and eigenvalues of \(H_0\), respectively, i.e.,
Definition: Let \(\{\left|{\phi_n}\right\rangle\}\) be the set of eigenstates and \(\{E_n\}\) be the corresponding eigenvalues of the time-independent Hamiltonian \(H_0\). They satisfy the eigenvalue equation: \[H_0 \left|{\phi_n}\right\rangle = E_n \left|{\phi_n}\right\rangle\]
We expand the time-dependent state \(\left|{\psi(t)}\right\rangle\) in the basis of \(\{\left|{\phi_n}\right\rangle\}\):
Definition: The time-dependent state \(\left|{\psi(t)}\right\rangle\) can be expanded in the basis of eigenstates \(\{\left|{\phi_n}\right\rangle\}\) of \(H_0\) as: \[\left|{\psi(t)}\right\rangle = \sum_n c_n(t) e^{-iE_n t/\hbar} \left|{\phi_n}\right\rangle\] where \(c_n(t)\) are time-dependent coefficients and \(e^{-iE_n t/\hbar}\) are phase factors associated with the eigenstates of \(H_0\).
Substituting this expansion into the time-dependent Schrödinger equation and projecting onto \(\left\langle{\phi_k}\right|\) yields the equation for the coefficients \(c_k(t)\):
Theorem: The time-dependent coefficients \(c_k(t)\) satisfy the following differential equation: \[i\hbar \dot{c}_k(t) = \sum_n \left\langle{\phi_k}\right| V(t) \left|{\phi_n}\right\rangle e^{-i(E_n - E_k) t/\hbar} c_n(t)\]
Defining
Definitions: Let \(\omega_{kn}\) be the transition frequency between eigenstates \(\left|{\phi_n}\right\rangle\) and \(\left|{\phi_k}\right\rangle\), and \(V_{kn}(t)\) be the matrix element of the perturbation \(V(t)\) between these eigenstates, defined as: \[\begin{aligned} \omega_{kn} &= (E_k - E_n)/\hbar \\ V_{kn}(t) &= \left\langle{\phi_k}\right| V(t) \left|{\phi_n}\right\rangle\end{aligned}\]
we have:
Theorem: Using the definitions of \(\omega_{kn}\) and \(V_{kn}(t)\), the equation for the coefficients \(c_k(t)\) can be rewritten as: \[i\hbar \dot{c}_k(t) = \sum_n V_{kn}(t) e^{i\omega_{kn} t} c_n(t)\]
To first order, assuming the system starts in state \(\left|{\phi_i}\right\rangle\) (i.e., \(c_i(0) = 1\) and \(c_n(0) = 0\) for \(n \neq i\)), and that \(c_i(t) \approx 1\) and \(c_n(t) \approx 0\) for \(n \neq i\) for short times, the equation for \(k \neq i\) becomes:
Theorem: Under the first-order approximation, the equation for the coefficients \(c_k^{(1)}(t)\) for \(k \neq i\) simplifies to: \[i\hbar \dot{c}_k^{(1)}(t) = V_{ki}(t) e^{i\omega_{ki} t}\]
Integrating from \(0\) to \(t\), we get the first-order coefficient:
Theorem: The first-order coefficient \(c_k^{(1)}(t)\) is given by the integral: \[c_k^{(1)}(t) = \frac{1}{i\hbar} \int_0^t V_{ki}(t') e^{i\omega_{ki} t'} dt'\]
This expression approximates the probability amplitude for transitioning from state \(\left|{\phi_i}\right\rangle\) to \(\left|{\phi_k}\right\rangle\) under the perturbation \(V(t)\).
Resonance with Oscillating Fields
Consider an oscillating perturbation
Definition: An oscillating perturbation \(V(t)\) with frequency \(\omega\) and amplitude \(V_0\) can be represented as: \[V(t) = V_0 \cos(\omega t) = \frac{V_0}{2} (e^{i\omega t} + e^{-i\omega t})\]
Then \(V_{ki}(t) e^{i\omega_{ki} t}\) becomes:
Theorem: For an oscillating perturbation \(V(t) = V_0 \cos(\omega t)\), the term \(V_{ki}(t) e^{i\omega_{ki} t}\) is given by: \[V_{ki}(t) e^{i\omega_{ki} t} = \frac{V_{0,ki}}{2} (e^{i(\omega_{ki} + \omega) t} + e^{i(\omega_{ki} - \omega) t})\] where \(V_{0,ki} = \left\langle{\phi_k}\right| V_0 \left|{\phi_i}\right\rangle\).
Substituting into the integral for \(c_k^{(1)}(t)\):
Theorem: Substituting the expression for \(V_{ki}(t) e^{i\omega_{ki} t}\) into the integral for \(c_k^{(1)}(t)\), we get: \[c_k^{(1)}(t) = \frac{V_{0,ki}}{2i\hbar} \int_0^t (e^{i(\omega_{ki} + \omega) t'} + e^{i(\omega_{ki} - \omega) t'}) dt'\]
Integrating yields:
Theorem: Integrating the expression above, the first-order coefficient \(c_k^{(1)}(t)\) for an oscillating perturbation is: \[c_k^{(1)}(t) = \frac{V_{0,ki}}{2\hbar} \left[ \frac{1 - e^{i(\omega_{ki} + \omega) t}}{\omega_{ki} + \omega} + \frac{1 - e^{i(\omega_{ki} - \omega) t}}{\omega_{ki} - \omega} \right]\]
The transition probability \(|c_k^{(1)}(t)|^2\) is maximized when the driving frequency \(\omega\) approaches the transition frequency \(\omega_{ki} = (E_k - E_i)/\hbar\).
Definition: The condition for resonance occurs when the driving frequency \(\omega\) of the oscillating perturbation is approximately equal to the transition frequency \(\omega_{ki} = (E_k - E_i)/\hbar\) between states \(\left|{\phi_i}\right\rangle\) and \(\left|{\phi_k}\right\rangle\): \[\omega \approx \omega_{ki}\]
This condition, \(\omega \approx \omega_{ki}\), is termed resonance. At resonance, the term \(\frac{1 - e^{i(\omega_{ki} - \omega) t}}{\omega_{ki} - \omega}\) dominates, leading to a significant transition probability.
This resonance phenomenon is analogous to efficiently pushing a swing; applying force at the swing’s natural frequency maximizes amplitude. In quantum systems, resonant oscillating fields efficiently drive transitions between energy levels, a principle fundamental to quantum control and spectroscopic techniques.
Superconducting Quantum Computing with Josephson Junctions
Superconducting quantum computing is a prominent approach for realizing quantum computers, leveraging superconducting circuits and Josephson junctions. These junctions exhibit unique quantum mechanical properties at extremely low temperatures, making them suitable for qubit implementation.
Josephson Junctions: The Core Element
A Josephson junction consists of two superconducting layers separated by a thin insulating barrier. At temperatures near absolute zero, Cooper pairs can tunnel across this barrier, giving rise to a dissipationless supercurrent—a phenomenon known as the Josephson effect.
Key properties of Josephson junctions are described by the Josephson relations:
Supercurrent:A supercurrent \(I\) flows through the junction up to a critical current \(I_c\) without any voltage drop.
Current-Phase Relation: The supercurrent is related to the phase difference \(\delta\) across the junction by \(I = I_c \sin(\delta)\).
Voltage-Phase Relation: The voltage \(V\) across the junction is related to the time derivative of the phase difference by \(V = \frac{\hbar}{2e} \frac{d\delta}{dt}\).
Definition: The Josephson relations describe the key properties of a Josephson junction:
Supercurrent: A supercurrent \(I\) flows through the junction up to a critical current \(I_c\) without any voltage drop.
Current-Phase Relation: The supercurrent \(I\) is related to the phase difference \(\delta\) across the junction by: \[I = I_c \sin(\delta)\]
Voltage-Phase Relation: The voltage \(V\) across the junction is related to the time derivative of the phase difference by: \[V = \frac{\hbar}{2e} \frac{d\delta}{dt}\]
where \(I_c\) is the critical current, \(\delta\) is the phase difference across the junction, \(\hbar\) is the reduced Planck constant, and \(e\) is the elementary charge.
These non-linear properties, combined with inherent inductance and capacitance, allow Josephson junctions to function as fundamental building blocks for quantum circuits and qubits.
Types of Superconducting Qubits
Josephson junctions can be configured to create various types of qubits. Two important examples are charge qubits and flux qubits.
Charge Qubits
Charge qubits encode quantum information in the charge states of a superconducting island. These qubits are formed by a superconducting island connected to a reservoir via a Josephson junction. The quantum states correspond to discrete numbers of Cooper pairs on the island. Control of charge qubits is achieved by applying a gate voltage, effectively tuning the electrostatic charging energy of the island. In essence, a capacitor is used to apply a static electric field, controlling the charge state.
Definition: Charge qubits are superconducting qubits that encode quantum information in the charge states of a superconducting island. They are formed by a superconducting island connected to a reservoir via a Josephson junction. The quantum states correspond to discrete numbers of Cooper pairs on the island. Control is achieved by applying a gate voltage to tune the electrostatic charging energy.
Flux Qubits
Flux qubits utilize superconducting loops containing Josephson junctions. The quantum states are defined by the clockwise or counterclockwise direction of persistent supercurrents circulating in the loop. These current directions correspond to quantized magnetic flux states through the loop. Control of flux qubits is realized by applying an external magnetic flux, typically using a coil to induce a current and generate a magnetic field that threads the loop, thereby manipulating the flux state.
Definition: Flux qubits are superconducting qubits that utilize superconducting loops containing Josephson junctions. Quantum information is encoded in the clockwise or counterclockwise direction of persistent supercurrents in the loop, corresponding to quantized magnetic flux states. Control is achieved by applying an external magnetic flux to manipulate the flux state.
Challenges and Practical Aspects
Superconducting quantum computers based on Josephson junctions are actively pursued due to their potential scalability and controllability. However, significant challenges remain:
Cryogenic Requirements: Operation necessitates extremely low temperatures, typically below 20 mK, requiring sophisticated and costly cryogenic infrastructure using liquid helium to cool down the system.
Coherence and Decoherence: Superconducting qubits are susceptible to decoherence from environmental noise, limiting their coherence times. Research focuses on mitigating decoherence to improve qubit performance.
Scalability and Complexity: Scaling up to large numbers of qubits while maintaining high fidelity and connectivity is a major engineering challenge. Fabricating and controlling complex multi-qubit circuits requires advanced nanofabrication and control techniques.
Control and Readout: Precise manipulation and measurement of qubit states are crucial. This involves developing sophisticated microwave pulse sequences for control and highly sensitive nanoelectronic devices based on quantum effects for state readout.
Remark: Superconducting quantum computing, while promising, faces several practical challenges including:
Cryogenic Requirements
Coherence and Decoherence
Scalability and Complexity
Control and Readout
Addressing these challenges is crucial for the advancement of this technology.
Despite these challenges, superconducting quantum computing is a leading technology. Real-world implementations involve complex, meter-sized apparatus to house and operate these quantum systems, as illustrated by existing quantum computing hardware. Current quantum computers are largely based on this technology, demonstrating its practical viability and promise for future advancements.
Density Matrix and Superoperators
While quantum states are typically represented by state vectors \(\left|{\psi}\right\rangle\) in Hilbert space, the density matrix formalism becomes essential when dealing with statistical mixtures of quantum states, known as mixed states, or when there is incomplete knowledge about the system’s state. Furthermore, the evolution of density matrices, particularly in open quantum systems, is described using superoperators.
Density Matrix: Describing Mixed States
Definition: The density matrix \(\rho\) is a positive semi-definite, Hermitian operator with trace equal to unity. It provides a comprehensive description of the quantum state, whether pure or mixed.
Definition: For a system in a pure state described by the state vector \(\left|{\psi}\right\rangle\), the density matrix is given by: \[\rho = \left|{\psi}\right\rangle\left\langle{\psi}\right|\]
Definition: A mixed state represents a statistical ensemble of pure states \(\{\left|{\psi_i}\right\rangle\}\) with corresponding probabilities \(\{p_i\}\). The density matrix for a mixed state is the ensemble average of the density matrices of the pure states: \[\rho = \sum_i p_i \left|{\psi_i}\right\rangle\left\langle{\psi_i}\right|\] where \(\sum_i p_i = 1\) and \(p_i \geq 0\).
The density matrix formalism is crucial for describing systems where we have statistical uncertainty about the quantum state, or when dealing with subsystems of larger entangled systems.
Superoperators: Evolution of Density Matrices
Definition: Superoperators are linear maps that act on operators, specifically on density matrices, transforming one density matrix into another. They describe transformations in the space of density matrices.
For a closed quantum system, the time evolution of the density matrix is governed by the von Neumann equation:
Theorem: The time evolution of the density matrix \(\rho\) for a closed quantum system is governed by the von Neumann equation: \[\frac{d\rho}{dt} = -\frac{i}{\hbar} [H, \rho]\] where \(H\) is the Hamiltonian of the system, and \([H, \rho] = H\rho - \rho H\) is the commutator.
where \(H\) is the Hamiltonian of the system, and \([H, \rho] = H\rho - \rho H\) is the commutator. This equation can be written using the Liouvillian superoperator \(\mathcal{L}\):
Definition: The Liouvillian superoperator \(\mathcal{L}\) is defined such that the von Neumann equation can be written as: \[\frac{d\rho}{dt} = \mathcal{L}(\rho), \quad \text{where} \quad \mathcal{L}(\rho) = -\frac{i}{\hbar} [H, \rho]\]
The formal solution to the von Neumann equation is given by:
Theorem: The formal solution to the von Neumann equation is: \[\rho(t) = e^{\mathcal{L}t} \rho(0) = U(t) \rho(0) U^\dagger(t)\] where \(U(t) = e^{-iHt/\hbar}\) is the unitary time evolution operator. The superoperator \(e^{\mathcal{L}t}\) is the Liouville propagator.
where \(U(t) = e^{-iHt/\hbar}\) is the unitary time evolution operator. In this case, the superoperator \(e^{\mathcal{L}t}\) is sometimes referred to as the Liouville propagator.
Superoperators are particularly important for describing open quantum systems, where the system interacts with an environment. In such scenarios, the evolution is generally non-unitary, and superoperators provide the necessary framework to describe processes such as:
Decoherence: Loss of quantum coherence due to interaction with the environment.
Dissipation: Energy loss from the system to the environment.
Quantum Channels: General transformations of quantum states, including noise and errors.
Remark: Superoperators are essential for describing open quantum systems and non-unitary processes such as:
Decoherence
Dissipation
Quantum Channels
These processes are crucial for understanding realistic quantum systems interacting with their environment.
The transition from using state vectors and operators in Hilbert space to density matrices and superoperators represents a shift to operator space (sometimes referred to as Liouville space or superoperator space). This generalization is essential for a complete description of quantum dynamics, especially when dealing with realistic quantum systems that are not perfectly isolated from their surroundings.
Conclusion
This lecture has covered advanced topics critical to quantum computing, starting with the challenge of realizing two-photon gates via cross-phase modulation and progressing to the practical domain of superconducting qubits using Josephson junctions. We emphasized the role of resonance and time-dependent perturbation theory in quantum manipulation and introduced the density matrix and superoperators as essential tools for describing quantum systems, particularly in complex and realistic scenarios.
Key insights from this lecture are:
Two-Qubit Gates: Essential for quantum universality, with cross-phase modulation offering a route for photonic implementations, albeit with significant practical hurdles.
Resonance: A fundamental principle for efficient quantum control using oscillating fields to drive transitions between quantum states.
Superconducting Qubits: Josephson junctions are the cornerstone of a promising quantum computing technology, balancing scalability and control, while facing challenges in cryogenic operation and decoherence mitigation.
Density Matrix and Superoperators: Provide a comprehensive framework for quantum mechanics beyond pure states, crucial for understanding mixed states and open quantum systems.
Remark: The key insights from this lecture are:
Two-Qubit Gates
Resonance
Superconducting Qubits
Density Matrix and Superoperators
These concepts are fundamental for understanding advanced topics in quantum computing.
Future studies should delve deeper into the mathematical intricacies of density matrices and superoperators, explore diverse physical implementations of quantum gates, and investigate advanced error correction and fault tolerance methods for robust quantum computation. The subsequent session will likely include exercises focused on applying the density matrix formalism to reinforce these concepts.