Purpose: This work proposes a novel RF pulse design for parallel transmit(pTx) systems to obtain uniform saturation of semisolid magnetization formagnetization transfer (MT) contrast in the presence of transmit field (B+1)inhomogeneities. The semisolid magnetization is usually modeled as beingpurely longitudinal, with the applied B+1 field saturating but not rotating itsmagnetization; thus, standard pTx pulse design methods do not apply.Theory and Methods: Pulse design for saturation homogeneity (PUSH) opti-mizes pTx RF pulses by considering uniformity of root-mean squared B+1 , Brms1 ,which relates to the rate of semisolid saturation. Here we considered designsconsisting of a small number of spatially non-selective sub-pulses optimized overeither a single 2D plane or 3D. Simulations and in vivo experiments on a 7T Terrasystem with an 8-TX Nova head coil in five subjects were carried out to study thehomogenization of Brms1 and of the MT contrast by acquiring MT ratio maps.Results: Simulations and in vivo experiments showed up to six and two timesmore uniform Brms1 compared to circular polarized (CP) mode for 2D and 3Doptimizations, respectively. This translated into 4 and 1.25 times more uniformMT contrast, consistently for all subjects, where two sub-pulses were enough forthe implementation and coil used.Conclusion: The proposed PUSH method obtains more uniform and higher MTcontrast than CP mode within the same specific absorption rate (SAR) budget.