MR acoustic radiation force imaging provides a promising method to monitor therapeutic ultrasound treatments. By measuring the displacement induced by the acoustic radiation force, MR acoustic radiation force imaging can locate the focal spot, without a significant temperature rise. In this work, the encoding gradient for MR acoustic radiation force imaging is optimized to achieve an enhanced accuracy and precision of the displacement measurement. By analyzing the sources of artifacts, bulk motion and eddy currents are shown to introduce errors to the measurement, and heavy diffusion-weighting is shown to result in noisy displacement maps. To eliminate these problems, a new encoding scheme is proposed, which utilizes a pair of bipolar gradients. Improved precision is achieved with robustness against bulk motion and background phase distortion, and improved accuracy is achieved with reduced diffusion-weighting and optimized encoding pulse width. The experiment result shows that the signal-to-noise ratio can be enhanced by more than 2-fold. These significant improvements are obtained at no cost of scan time or encoding sensitivity, enabling the detection of a displacement less than 0.1 gm in a gel phantom with MR acoustic radiation force imaging. Magn Reson Med 63:1050-1058, 2010. (C) 2010 Wiley-Liss, Inc.