Since 2005 there has been a huge growth in the use of engineered control pulses to perform desired quantum operations in systems such as NMR quantum information processors. These approaches, which build on the original gradient ascent pulse engineering (GRAPE) algorithm, remain computationally intensive because of the need to calculate matrix exponentials for each time step in the control pulse. Here we discuss how the propagators for each time step can be approximated using the Trotter-Suzuki formula, and a further speed up achieved by avoiding unnecessary operations. The resulting procedure can give a substantial speed gain with negligible cost in propagator error, providing a more practical approach to pulse engineering.
PACS numbers:Quantum information processors encode information in two-level quantum systems (qubits) and manipulate this through a series of elementary unitary transformations (quantum logic gates) [1,2]. Quantum control seeks to implement some target unitary propagator U in a quantum system with background Hamiltonian H 0 by applying some time-dependent Hamiltonian H 1 (t). The resulting operator can be written aswhere T is the Dyson time-ordering operator. To make progress beyond this formal solution it is usually necessary to replace this continuously varying Hamiltonian by a piecewise constant form, so thatand to write the time-varying portion of the Hamiltonian as a weighted sum of a set of p distinct control fieldsAny particular control pulse can then be described by the corresponding set of amplitudes, a k j , and the time step δt, here taken as fixed. The quality of a control pulse can be measured by its fidelity with the desired operation Uwhere the Hilbert-Schmidt inner product is defined by U |V = tr(U † V ), possibly normalised by the dimension of the operators [2]. The optimal control problem is then to find the set of amplitudes which maximises this fidelity, usually in the presence of practical constraints on the magnitudes of the amplitudes and the total length * Electronic address: jonathan.jones@physics.ox.ac.uk of the sequence. This is computationally challenging as the dimension of the underlying Hilbert space rises exponentially with the number of qubits to be controlled, although this difficulty can in some cases be reduced by using subsystems to simplify the calculations [3]. One recent approach [4] is to use subsystem methods to find approximate control pulses and then optimise these directly using the quantum system itself. Whatever approach is adopted, it is important to perform any computations as efficiently as possible.