Heteronuclear decoupling by optimal tracking

The problem to design efficient heteronuclear decoupling sequences is studied using optimal control methods. A generalized version of the gradient ascent engineering (GRAPE) algorithm is presented that makes it possible to design complex non-periodic decoupling sequences which are characterized by tens of thousands of pulse sequence parameters. In contrast to conventional approaches based on average Hamiltonian theory, the concept of optimal tracking is used: a pulse sequence is designed that steers the evolution of an ensemble of spin systems such that at a series of time points, a specified trajectory of the density operator is tracked as closely as possible. The approach is demonstrated for the case of low-power heteronuclear decoupling in the liquid state for in vivo applications. Compared to conventional sequences, significant gains in decoupling efficiency and robustness with respect to offset and inhomogeneity of the radio-frequency field were found in simulations and experiments.