Several recent studies have shown that very wide binary stars can potentially provide an interesting test for modified-gravity theories which attempt to emulate dark matter; these systems should be almost Newtonian according to standard dark-matter theories, while the predictions for MOND-like theories are distinctly different, if the various observational issues can be overcome. Here we explore an observational application of the test from the recent GAIA DR2 data release: we select a large sample of $\sim 24,000$ candidate wide binary stars with distance $< 200$ parsec and magnitudes $G < 16$ from GAIA DR2, and estimated component masses using a main-sequence mass-luminosity relation. We then compare the frequency distribution of pairwise relative projected velocity (relative to circular-orbit value) as a function of projected separation; these distributions show a clear peak at a value close to Newtonian expectations, along with a long `tail' which extends to much larger velocity ratios; the `tail' is considerably more numerous than in control samples constructed from DR2 with randomised positions, so its origin is unclear. Comparing the velocity histograms with simulated data, we conclude that MOND-like theories without an external field effect are strongly inconsistent with the observed data since they predict a peak-shift in clear disagreement with the data; testing MOND-like theories with an external field effect is not decisive at present, but has good prospects to become decisive in future with improved modelling or understanding of the high-velocity tail, and additional spectroscopic data.

The standard LambdaCDM model based on General Relativity (GR) including cold dark matter (CDM) is very successful at fitting cosmological observations, but recent non-detections of candidate dark matter (DM) particles mean that various modified-gravity theories remain of significant interest. The latter generally involve modifications to GR below a critical acceleration scale $\sim 10^{-10} \, m \, s^{-2}$. Wide-binary (WB) star systems with separations $> 5 \, kAU$ provide an interesting test for modified gravity, due to being in or near the low-acceleration regime and presumably containing negligible DM. Here, we explore the prospects for new observations pending from the GAIA spacecraft to provide tests of GR against MOND or TeVes-like theories in a regime only partially explored to date. In particular, we find that a histogram of (3D) binary relative velocities against circular velocity predicted from the (2D) projected separations predicts a rather sharp feature in this distribution for standard gravity, with an 80th (90th) percentile value close to 1.025 (1.14) with rather weak dependence on the eccentricity distribution. However, MOND/TeVeS theories produce a shifted distribution, with a significant increase in these upper percentiles. In MOND-like theories {\em without} an external field effect, there are large shifts of order unity. With the external field effect included, the shifts are considerably reduced to $\sim 0.04 - 0.08$, but are still potentially detectable statistically given reasonably large samples and good control of contaminants. In principle, followup of GAIA-selected wide binaries with ground-based radial velocities accurate to < 0.03 km/s should be able to produce an interesting new constraint on modified-gravity theories.