Abstract

Designing optimal radiofrequency (RF) pulses for nuclear magnetic resonance (NMR) imaging can be formulated as an optimal control problem for the Bloch equation. Hardware constraints that restrict admissible pulse phase and amplitudes to a discrete set can be realized via a vector-valued "multi-bang penalty", which leads to a convex optimization problem that can be solved numerically using semismooth Newton methods.