The excitation k-space perspective on small-tip-angle selective excitation has facilitated RF pulse designs in a range of MR applications. In this paper, k-space-based design of multidimensional RF pulses is formulated as a quadratic optimization problem, and solved efficiently by the iterative conjugate-gradient (CG) algorithm. Compared to conventional design approaches, such as the conjugate-phase (CP) method, the new design approach is beneficial in several regards. It generally produces more accurate excitation patterns. The improvement is particularly significant when k-space is undersampled, and it can potentially shorten pulse lengths. A prominent improvement in accuracy is also observed when large off-resonance gradients are present. A further boost in excitation accuracy can be accomplished in regions of interest (ROIs) if they are specified together with "don't-care" regions. The design of RF pulses for multidimensional, small-tipangle selective excitation is facilitated by the excitation k-space perspective developed by Pauly et al. (1) under the small-tip-angle approximation to the Bloch equation. kSpace-based selective excitation has been used in a range of MR applications, such as functional MRI (fMRI) artifact correction (2), brain imaging with reduced field of view (FOV) (3), blood velocity measurement (4), parallel excitation using multiple transmit coils (5,6), and excitation inhomogeneity correction (7). The k-space perspective is popular because it provides a fairly accurate linear Fourier relationship between the time-varying gradient and RF waveforms, and the resulting transverse excitation pattern. The Fourier relationship can also be established for rotation angles (possibly large), provided that certain symmetry conditions are satisfied (8).A common approach to small-tip-angle RF pulse design is to predetermine the gradient waveforms and thus the k-space trajectory, and then obtain the complex-valued RF waveform by sampling the Fourier transform of the desired excitation pattern along the trajectory. Afterwards, the sample values are compensated for density variation in the trajectory. Some researchers have designed pulses by adopting the conjugate-phase (CP) approach from image reconstruction (9,10), which accounts for off-resonance effects and thus can correct for it to some extent (11,12).These conventional approaches are nonideal in several regards. First, in terms of minimizing excitation error, they generally produce pulses that are suboptimal even with respect to the linear design model. This design suboptimality is an eradicable source of excitation error, on top of the intractable amount of error due to the small-tip-angle approximation underlying excitation k-space. The excitation error due to design suboptimality is particularly large when the trajectory undersamples k-space, or when large spatial variations of off-resonance are present (11). To compensate for the effects of off-resonance gradients to some extent, Noll et al. (13) suggested the use of a sophisticated density c...