guide-tensor.rst

Using Tensor Products and Partial Traces

.. ipython::
   :suppress:

   In [1]: from qutip import *

Tensor products

To describe the states of multipartite quantum systems - such as two coupled qubits, a qubit coupled to an oscillator, etc. - we need to expand the Hilbert space by taking the tensor product of the state vectors for each of the system components. Similarly, the operators acting on the state vectors in the combined Hilbert space (describing the coupled system) are formed by taking the tensor product of the individual operators.

In QuTiP the function :func:`qutip.tensor.tensor` is used to accomplish this task. This function takes as argument a collection:

>>> tensor(op1, op2, op3)

or a list:

>>> tensor([op1, op2, op3])

of state vectors or operators and returns a composite quantum object for the combined Hilbert space. The function accepts an arbitray number of states or operators as argument. The type returned quantum object is the same as that of the input(s).

For example, the state vector describing two qubits in their ground states is formed by taking the tensor product of the two single-qubit ground state vectors:

.. ipython::

   In [1]: tensor(basis(2, 0), basis(2, 0))

or equivalently using the list format:

.. ipython::

   In [1]: tensor([basis(2, 0), basis(2, 0)])

This is straightforward to generalize to more qubits by adding more component state vectors in the argument list to the :func:`qutip.tensor.tensor` function, as illustrated in the following example:

.. ipython::

   In [1]: tensor((basis(2, 0) + basis(2, 1)).unit(), (basis(2, 0) + basis(2, 1)).unit(), basis(2, 0))

This state is slightly more complicated, describing two qubits in a superposition between the up and down states, while the third qubit is in its ground state.

To construct operators that act on an extended Hilbert space of a combined system, we similarly pass a list of operators for each component system to the :func:`qutip.tensor.tensor` function. For example, to form the operator that represents the simultaneous action of the \sigma_x operator on two qubits:

.. ipython::

   In [1]: tensor(sigmax(), sigmax())

To create operators in a combined Hilbert space that only act only on a single component, we take the tensor product of the operator acting on the subspace of interest, with the identity operators corresponding to the components that are to be unchanged. For example, the operator that represents \sigma_z on the first qubit in a two-qubit system, while leaving the second qubit unaffected:

.. ipython::

   In [1]: tensor(sigmaz(), qeye(2))

Example: Constructing composite Hamiltonians

The :func:`qutip.tensor.tensor` function is extensively used when constructing Hamiltonians for composite systems. Here we'll look at some simple examples.

Two coupled qubits

First, let's consider a system of two coupled qubits. Assume that both qubit has equal energy splitting, and that the qubits are coupled through a \sigma_x\otimes\sigma_x interaction with strength g = 0.05 (in units where the bare qubit energy splitting is unity). The Hamiltonian describing this system is:

.. ipython::

   In [1]: H = tensor(sigmaz(), qeye(2)) + tensor(qeye(2), sigmaz()) + 0.05 * tensor(sigmax(), sigmax())

   In [2]: H

Three coupled qubits

The two-qubit example is easily generalized to three coupled qubits:

.. ipython::

        In [1]: H = tensor(sigmaz(), qeye(2), qeye(2)) + tensor(qeye(2), sigmaz(), qeye(2)) + tensor(qeye(2), qeye(2), sigmaz()) + 0.5 * tensor(sigmax(), sigmax(), qeye(2)) +  0.25 * tensor(qeye(2), sigmax(), sigmax())

        In [2]: H

A two-level system coupled to a cavity: The Jaynes-Cummings model

The simplest possible quantum mechanical description for light-matter interaction is encapsulated in the Jaynes-Cummings model, which describes the coupling between a two-level atom and a single-mode electromagnetic field (a cavity mode). Denoting the energy splitting of the atom and cavity omega_a and omega_c, respectively, and the atom-cavity interaction strength g, the Jaynes-Cumming Hamiltonian can be constructed as:

>>> N = 10
>>> omega_a = 1.0
>>> omega_c = 1.25
>>> g = 0.05
>>> a = tensor(qeye(2), destroy(N))
>>> sm = tensor(destroy(2), qeye(N))
>>> sz = tensor(sigmaz(), qeye(N))
>>> H = 0.5 * omega_a * sz + omega_c * a.dag() * a + g * (a.dag() * sm + a * sm.dag())

Here N is the number of Fock states included in the cavity mode.

Partial trace

The partial trace is an operation that reduces the dimension of a Hilbert space by eliminating some degrees of freedom by averaging (tracing). In this sense it is therefore the converse of the tensor product. It is useful when one is interested in only a part of a coupled quantum system. For open quantum systems, this typically involves tracing over the environment leaving only the system of interest. In QuTiP the class method :func:`qutip.Qobj.ptrace` is used to take partial traces. :func:`qutip.Qobj.ptrace` acts on the :class:`qutip.Qobj` instance for which it is called, and it takes one argument sel, which is a list of integers that mark the component systems that should be kept. All other components are traced out.

For example, the density matrix describing a single qubit obtained from a coupled two-qubit system is obtained via:

.. ipython::

        In [1]: psi = tensor(basis(2, 0), basis(2, 1))

        In [2]: psi.ptrace(0)

        In [3]: psi.ptrace(1)

Note that the partial trace always results in a density matrix (mixed state), regardless of whether the composite system is a pure state (described by a state vector) or a mixed state (described by a density matrix):

.. ipython::

        In [1]: psi = tensor((basis(2, 0) + basis(2, 1)).unit(), basis(2, 0))

        In [2]: psi

        In [3]: psi.ptrace(0)

        In [4]: rho = tensor(ket2dm((basis(2, 0) + basis(2, 1)).unit()), fock_dm(2, 0))

        In [5]: rho

        In [6]: rho.ptrace(0)