Chapter 9 - Dynamical Decoupling Controls

Introduction

In the last chapter, it was shown that a symmetry in the system-bath Hamiltonian, if present, could be used to construct states immune to noise. Under certain conditions, it is possible to remove errors or to create a symmetry in the evolution of a quantum system. This is done by averaging away errors though repeated use of external controls which act on the system. These controls are often called "dynamical decoupling controls" due to their original objective of decoupling the system from the bath. They are quite generally useful controls to consider for the elimination and/or reduction of errors. In this chapter, a simple introduction to dynamical decoupling controls is given and some important concepts discussed.

General Conditions

As stated in Chapter 8 the Hamiltonian describing the evolution of a system and bath which are coupled together can always be written as

$H=H_{S}\otimes I_{B}+I_{S}\otimes H_{B}+H_{I},\,\!$

where $H_{S}\,\!$ acts only on the system, $H_{B}\,\!$ acts only on the bath, and

$H_{I}=\sum _{\alpha }S_{\alpha }\otimes B_{\alpha },\,\!$

is the interaction Hamiltonian with the $S_{\alpha }\,\!$ acting only on the system and the $B_{\alpha }\,\!$ acting only on the bath.

The idea is to modify the evolution of the system and bath such that the errors are reduced or eliminated using external control Hamiltonians. These controls are called dynamical decoupling controls since they are used to decouple (at least approximately decouple) the system from the bath. Since can be difficult to change states of a bath, indeed one often does not know details of the bath, the controls which are to be used for reducing errors should act on the system. As discussed previously, the errors arise from the system-bath interaction Hamiltonian $H_{I}\,\!$ and, in particular, the system operators $S_{\alpha }\,\!$ are the operators which describe the affect of the coupling on the system. In general the interaction Hamiltonian is time-dependent since the bath operators will change in time. However, for short times we may assume the interaction Hamiltonian is unchanged, or at least approximately constant. This is sometimes called the short-time assumption in dynamical decoupling.

The Magnus Expansion

A fairly good starting point to see how this is done is the so-called Magnus expansion. (See Blanes, et al. and references therein.) The general problem is that a time-dependent operation is to be applied to the Hamiltonian making the Hamiltonian itself time-dependent and one would like to solve the time-dependent Schrodinger equation:

	$i{\frac {\partial }{\partial t}}\left\vert \Psi (t)\right\rangle =H(t)\left\vert \Psi (t)\right\rangle ,\,\!$	(9.1)

which is sometimes written as

	$U^{\prime }(t)=-iH(t)U(t).\,\!$	(9.2)

The question is, what $U(t)\,\!$ will solve this equation? If $U(t)\,\!$ and $H(t)\,\!$ are just numbers, the solution would be

	$U^{\prime }(t)=\exp \left(-i\int _{0}^{t}H(t^{\prime })dt^{\prime }\right).\,\!$	(9.3)

However, when the Schrodinger equation is the equation to be solved, $U(t)\,\!$ and $H(t)\,\!$ are matrices. To be specific, $U(t)\,\!$ is a unitary matrix and $H(t)\,\!$ is a Hermitian matrix. The solution is often written in the form

	$U^{\prime }(t)={\mathcal {T}}\left[\exp \left(-i\int _{0}^{t}H(t^{\prime })dt^{\prime }\right)\right],\,\!$	(9.4)

where ${\mathcal {T}}\,\!$ denotes the time-ordered exponential. In this case, matrices do not commute so that the exponential must be handled with care. Operators must be ordered according to the time where they appear in the operation, and the solution Eq.(9.2) is not the solution to the problem unless $H(t)\,\!$ is a constant matrix.

The solution to this problem is the following,

	$U^{\prime }(T)=\exp \left(-i\Omega (T)T\right),\,\!$	(9.5)

where

	$\Omega (T)=\sum _{k=1}^{\infty }\Omega _{k},\,\!$	(9.6)

and

	${\begin{aligned}\Omega _{1}&={\frac {1}{T}}\int _{0}^{T}H(t_{1})dt_{1},\\\Omega _{2}&=-{\frac {1}{T^{2}}}{\frac {i}{2}}\int _{0}^{T}dt_{1}\int _{0}^{t_{1}}dt_{2}[H(t_{1}),H(t_{2})],\\\Omega _{3}&=-{\frac {1}{T^{3}}}{\frac {1}{6}}\int _{0}^{T}dt_{1}\int _{0}^{t_{1}}dt_{2}\int _{0}^{t_{2}}dt_{3}([H(t_{1}),[H(t_{2}),H(t_{3})]]+[H(t_{3}),[H(t_{2}),H(t_{1})]])\\&{\mbox{etc.}},\end{aligned}}\,\!$	(9.7)

where $T\,\!$ is some characteristic time scale.

$\,\!$

A First-Order Theory

Let us suppose that the

The Single-Qubit Case

The simplest case involves the elimination of an error on a single qubit state.