  Recently a remarkable relation has been demonstrated between the observed
radial acceleration in disk galaxies and the acceleration predicted on the
basis of baryonic matter alone. Here we study this relation within the
framework of the modified gravity model MOND. The field equations of MOND
automatically imply the radial acceleration relation for spherically symmetric
galaxies, but for disk galaxies deviations from the relation are expected. Here
we investigate whether these deviations are of sufficient magnitude to bring
MOND into conflict with the observed relation. In the quasilinear formulation
of MOND, to calculate the gravitational field of a given distribution of
matter, an intermediate step is to calculate the "pristine field", which is a
simple nonlinear function of the Newtonian field corresponding to the same
distribution of matter. Hence, to the extent that the quasilinear gravitational
field is approximately equal to the pristine field, the radial acceleration
relation will be satisfied. We show that the difference between the quasilinear
and pristine fields obeys the equations of magnetostatics, the curl of the
pristine field serves as the source for the difference in the two fields, much
as currents serve as sources for the magnetic field. Using the magnetostatic
analogy we numerically study the difference between the pristine and
quasilinear fields for simple model galaxies with a Gaussian profile. Our
principal finding is that the difference between the fields is small compared
to the observational uncertainties and that quasilinear MOND is therefore
compatible with the observed radial acceleration relation.
